content/math/linearRecurrenceOld.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33

constexpr ll mod = 1'000'000'007;
vector<ll> modMul(const vector<ll>& a, const vector<ll>& b,
                  const vector<ll>& c) {
	ll n = ssize(c);
	vector<ll> res(n * 2 + 1);
	for (int i = 0; i <= n; i++) { //a*b
		for (int j = 0; j <= n; j++) {
			res[i + j] += a[i] * b[j];
			res[i + j] %= mod;
	}}
	for (int i = 2 * n; i > n; i--) { //res%c
		for (int j = 0; j < n; j++) {
			res[i - 1 - j] += res[i] * c[j];
			res[i - 1 - j] %= mod;
	}}
	res.resize(n + 1);
	return res;
}

ll kthTerm(const vector<ll>& f, const vector<ll>& c, ll k) {
	assert(ssize(f) == ssize(c));
	vector<ll> tmp(ssize(c) + 1), a(ssize(c) + 1);
	tmp[0] = a[1] = 1; //tmp = (x^k) % c

	for (k++; k > 0; k /= 2) {
		if (k & 1) tmp = modMul(tmp, a, c);
		a = modMul(a, a, c);
	}

	ll res = 0;
	for (int i = 0; i < ssize(c); i++) res += (tmp[i+1] * f[i]) % mod;
	return res % mod;
}