Codeforce 1182 E. Product Oriented Recurrence（递推 + 矩阵快速幂 + 欧拉降幂）

---

$$f_x=c^{2x-6}\ast f_{x-1}\ast f_{x-2}\ast f_{x-3}$$$$c^x{f}_x=c^{x-1}{f}_{x-1}\ast c^{x-2}{f}_{x-2}\ast c^{x-3}{f}_{x-3}$$令$F_x=c^x\ast f_x$$$F_4=F_1\ast F_2\ast F_3$$$$F_5=F_4\ast F_3\ast F_2=F_1\ast F_2^2\ast F_3^2$$$$F_6=F_5\ast F_4\ast F_3=F_1^2\ast F_2^3\ast F_3^4$$$$F_7=F_6\ast F_5\ast F_4$$即任意一项$F_x$都可以转化为 $F_1^{x_1}\ast F_2^{x_2}\ast F_3^{x_3}$，第7项及以后的项,$x_1=(x-1{)}_1+(x-2{)}_1+(x-3{)}_1$，对于 $x_2$和$x_3$也同样满足。
也就是说它的幂次是一个递推式，可以用矩阵快速幂进行加速，根据欧拉降幂来取模。根据递推式构造3阶递推矩阵，以及幂次的初始矩阵：$F_4,F_5,F_6$的各项幂次构成初始矩阵。
这题就做完了

---

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 1e9 + 7;
ll a[5][5],tmp[5][5],c[5][5],ans[5][5];
ll b1[5][5],b2[5][5],b3[5][5];
ll n,f1,f2,f3,C;
void mul(ll a[5][5],ll b[5][5],ll p,int n,int m,int c) {
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= c; j++)
			tmp[i][j] = 0;
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= m; j++)
			for(int k = 1; k <= c; k++)
				tmp[i][k] += 1ll * a[i][j] * b[j][k] % p;
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= c; j++)
			a[i][j] = tmp[i][j];				
}
void mfpow(ll a[5][5],ll b,ll p,int n) {
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= n; j++)
			ans[i][j] = i == j;	
	while(b) {
		if(b & 1)
			mul(ans,a,p,n,n,n);
		b >>= 1;
		mul(a,a,p,n,n,n);
	}	
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= n; j++)
			a[i][j] = ans[i][j];
}
ll fpow(ll a,ll b) {
	ll r = 1;
	while(b) {
		if(b & 1) r = r * a % mod;
		a = a * a % mod;
		b >>= 1;
 	}
 	return r;
}
void solve(ll a[5][5],ll b[5][5],ll p,ll n,ll m,ll c) {
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= m; j++)
			ans[i][j] = a[i][j];
	mul(ans,b,p,n,m,c);
	for(int i = 1; i <= n; i++)
		for(int j = 1; j <= c; j++)
			b[i][j] = ans[i][j];
}
int main() {
	c[1][1] = c[1][2] = c[1][3] = 1;
	c[2][1] = c[3][2] = 1;
	b1[1][1] = 2,b1[2][1] = b1[3][1] = 1;
	b2[1][1] = 3,b2[2][1] = 2,b2[3][1] = 1;
	b3[1][1] = 4,b3[2][1] = 2,b3[3][1] = 1;
	ll a1,a2,a3;
	scanf("%lld%lld%lld%lld%lld",&n,&f1,&f2,&f3,&C);
	f1 = f1 * C % mod;
	f2 = f2 * fpow(C,2) % mod;
	f3 = f3 * fpow(C,3) % mod;
	ll p = fpow(C,n);
	ll inv = fpow(p,mod - 2);
	n -= 3;
	if(n > 3) {
		mfpow(c,n - 3,mod - 1,3);
		solve(c,b1,mod - 1,3,3,1);
		solve(c,b2,mod - 1,3,3,1);
		solve(c,b3,mod - 1,3,3,1);
		//printf("%lld %lld %lld\n",tmp[1][1],b2[1][1],b3[1][1]);
		a1 = fpow(f1,b1[1][1]);
		a2 = fpow(f2,b2[1][1]);
		a3 = fpow(f3,b3[1][1]);
	}
	else {
		a1 = fpow(f1,b1[3 - n + 1][1]);
		a2 = fpow(f2,b2[3 - n + 1][1]);
		a3 = fpow(f3,b3[3 - n + 1][1]);		
	}
	printf("%lld\n",a1 * a2 % mod * a3 % mod * inv % mod);	
	return 0;
}