输入的数据会大于mod，当n=0时输出0.
#include#include #include #include #include #include #include #include using namespace std;typedef long long ll;const ll mod = 1000000007;struct Matrix{ ll da[5][5]; Matrix (){memset(da,0,sizeof(da));} Matrix operator *(const Matrix &a) { Matrix res; for(int i=0;i<5;i++) { for(int j=0;j<5;j++) { res.da[i][j]=0; for(int k=0;k<5;k++) res.da[i][j]=(res.da[i][j]+da[i][k]*a.da[k][j])%mod; } } return res; }}base,I;Matrix f(ll n){ if(n==0)return I; if(n%2==0) { Matrix res=f(n/2); return res*res; }else { Matrix res=f(n/2); return res*res*base; }}ll n,AX,AY,BX,BY,A0,B0;int main(){ while(scanf("%I64d",&n)!=EOF) { scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&A0,&AX,&AY,&B0,&BX,&BY); if(n==0){printf("0\n");continue;} base.da[0][0]=AX*BX%mod; base.da[0][1]=AX*BY%mod; base.da[0][2]=BX*AY%mod; base.da[0][3]=AY*BY%mod; base.da[0][4]=A0*B0%mod; base.da[1][1]=AX%mod; base.da[1][3]=AY%mod; base.da[1][4]=A0%mod; base.da[2][2]=BX%mod; base.da[2][3]=BY%mod; base.da[2][4]=B0%mod; base.da[3][3]=base.da[4][4]=base.da[3][4]=1; for(int i=0;i<5;i++) I.da[i][i]=1; Matrix res=f(n-1); printf("%I64d\n",((((res.da[0][0]*A0)%mod)*B0)%mod+res.da[0][1]*A0+res.da[0][2]*B0+res.da[0][3]+res.da[0][4])%mod); } return 0;}