usaco

Clean implementations of solutions to USACO problems

  1. 1
  2. 2
  3. 3
  4. 4
  5. 5
  6. 6
  7. 7
  8. 8
  9. 9
  10. 10
  11. 11
  12. 12
  13. 13
  14. 14
  15. 15
  16. 16
  17. 17
  18. 18
  19. 19
  20. 20
  21. 21
  22. 22
  23. 23
  24. 24
  25. 25
  26. 26
  27. 27
  28. 28
  29. 29
  30. 30
  31. 31
  32. 32
  33. 33
  34. 34
  35. 35
  36. 36
  37. 37
  38. 38
  39. 39
  40. 40
  41. 41
  42. 42
  43. 43
  44. 44
  45. 45
  46. 46
  47. 47
  48. 48
  49. 49
  50. 50
  51. 51
  52. 52
  53. 53
  54. 54
  55. 55
  56. 56
  57. 57
  58. 58
  59. 59
  60. 60
  61. 61
  62. 62
  63. 63
  64. 64
  65. 65
  66. 66
  67. 67
  68. 68
  69. 69
  70. 70
  71. 71
  72. 72
  73. 73
  74. 74
  75. 75
  76. 76
  77. 77
  78. 78
  79. 79
  80. 80
  81. 81
  82. 82
  83. 83
  84. 84
  85. 85
  86. 86
  87. 87
  88. 88
  89. 89
  90. 90
  91. 91
  92. 92
  93. 93
  94. 94
  95. 95
  96. 96
  97. 97
  98. 98
  99. 99
  100. 100
  101. 101
  102. 102
  103. 103
#include <bits/stdc++.h>
#define init_io(pname) ifstream cin((string)pname+".in"); ofstream cout((string)pname+".out"); ios_base::sync_with_stdio(false); cin.tie(NULL)
#define FOR(i, a, b, in) for (int i = (a); i < (b); i += in)
#define REP(i, a, b) for (int i = (a); i < (b); i++)
#define RFOR(i, a, b, in) for (int i = (a) - 1; i >= (b); i -= in)
#define RREP(i, a, b) for (int i = (a) - 1; i >= (b); i--)
#define trav(a, x) for (auto& a : x)
#define mp make_pair
#define pb push_back
#define f first
#define s second
#define lb lower_bound
#define ub upper_bound
#define sz(x) (int)x.size()
#define all(x) begin(x), end(x)
#define rsz resize
#define mem(a, b) memset(a, (b), sizeof(a))
using namespace std;
typedef string str;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> ii; typedef pair<ll, ll> pl; typedef pair<ld, ld> pd;
typedef vector<int> vi; typedef vector<ll> vl; typedef vector<ld> vd;
typedef vector<ii> vii; typedef vector<pl> vpl; typedef vector<pd> vpd;
constexpr auto INF = (int)1e9;
constexpr auto LINF = (ll)1e18;
constexpr auto MOD = 1000000007;

typedef vector<vi> matrix;

int d = 0;
map<str, int> var; // Variables
vector<str> program; // Program
stack<int> loop_st; // Loop stack
vector<matrix> mat_st; // Matrix stack

void mult(matrix& A, matrix& B) { // B = A * B;
	matrix res(d + 1, vi(d + 1, 0));
	for (int i = 0; i <= d; i++) {
		for (int j = 0; j <= d; j++) {
			for (int k = 0; k <= d; k++) res[i][j] = ((ll)A[i][k] * (ll)B[k][j] + res[i][j]) % MOD;
		}
	}
	B = res;
}

void pow(matrix& A, int n) { // A = A ^ n
	matrix res(d + 1, vi(d + 1, 0));
	for (int i = 0; i <= d; i++) res[i][i] = 1;
	while (n) {
		if (n % 2 == 1) mult(A, res);
		mult(A, A);
		n /= 2;
	}
	A = res;
}

int main() {
	init_io("cowbasic");

	str input;
	while (getline(cin, input)) program.push_back(input);
	for (auto& line : program) { // Record variables
		stringstream ss(line);
		str n; ss >> n;
		if (islower(n[0]) && var.find(n) == var.end()) var[n] = d++;
	}
	mat_st.push_back({});
	mat_st.back().resize(d + 1, vi(d + 1, 0));
	for (int i = 0; i <= d; i++) mat_st.back()[i][i] = 1;
	for (auto& line : program) { // Process each line
		stringstream ss(line);
		str n; ss >> n;
		if (islower(n[0])) { // Variable assignment
			matrix M(d + 1, vi(d + 1, 0));
			for (int i = 0; i <= d; i++) {
				if (i != var[n]) M[i][i] = 1;
			}
			str tmp;
			while (ss >> tmp) {
				if (islower(tmp[0])) M[var[n]][var[tmp]]++;
				else if (isdigit(tmp[0])) M[var[n]][d] += stoi(tmp);
			}
			mult(M, mat_st.back());
		}
		else if (isdigit(n[0])) { // Begin loop
			mat_st.push_back({});
			mat_st.back().resize(d + 1, vi(d + 1, 0));
			for (int i = 0; i <= d; i++) mat_st.back()[i][i] = 1;
			loop_st.push(stoi(n));
		}
		else if (n[0] == '}') { // End loop
			pow(mat_st.back(), loop_st.top());
			mult(mat_st.back(), mat_st[mat_st.size() - 2]);
			mat_st.pop_back();
			loop_st.pop();
		}
		else { // Return
			ss >> n;
			cout << mat_st.back()[var[n]][d] << endl;
		}
	}
}