这题给我坑坏了
构造矩阵乘法,时刻要注意%mod 还有n==0时输出0
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const long long maxn=10;
long long n,m,k;
const long long mod=1000000007;
typedef struct
{
long long m[maxn][maxn];
} Matrax;
Matrax a,per;
Matrax multi(Matrax a,Matrax b,long long m) {
Matrax c;
for(long long i=0;i<n;i++)
{
for(long long j=0;j<n;j++)
{
c.m[i][j]=0;
for(long long k=0;k<n;k++)
{
c.m[i][j]+=(a.m[i][k]*b.m[k][j])%mod;
// c.m[i][j]%=m;
if(c.m[i][j]>mod) c.m[i][j]-=mod;
}
}
}
return c;
}
Matrax power(Matrax a,long long k,long long m)
{
Matrax p,ans;
p=a;
for(long long i=0;i<n;i++)
for(long long j=0;j<n;j++) ans.m[i][j]=(i==j);
while(k)
{
if(k&1)
{
ans=multi(ans,p,m);
k--;
}
k/=2;
p=multi(p,p,m);
}
return ans;
}
Matrax zero;
int main()
{
std::ios::sync_with_stdio(false);
memset(zero.m,0,sizeof(zero.m));
long long nn;
n=8;
while(cin>>nn)
{
long long ao,ax,ay,bo,bx,by;
cin>>ao>>ax>>ay>>bo>>bx>>by;
//scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&ao,&ax,&ay,&bo,&bx,&by);
ao%=mod;
ax%=mod;
ay%=mod;
bo%=mod;
bx%=mod;
by%=mod;
Matrax a;
memset(a.m,0,sizeof(a));
int fuck=0;
if(nn==0)
{
//cout<<(ao*bo+(ao-ay)/ax*(bo-by)/bx)<<endl;
cout<<fuck<<endl;
continue;
}
a.m[0][0]=(ax*bx)%mod; a.m[0][1]=0;a.m[0][2]=0;a.m[0][3]=0;
a.m[1][0]=(ay*bx)%mod;a.m[1][1]=bx%mod;a.m[1][2]=0;a.m[1][3]=0;
a.m[2][0]=(ax*by)%mod;a.m[2][1]=0;a.m[2][2]=ax%mod;a.m[2][3]=0;
a.m[3][0]=(ay*by)%mod;a.m[3][1]=by%mod;a.m[3][2]=ay%mod;a.m[3][3]=1;
for(long long i=0;i<8;i++)
{
for(long long j=4;j<8;j++)
if(i==j-4||i==j)
a.m[i][j]=1;
}
for(long long i=4;i<8;i++)
{
for(long long j=0;j<4;j++)
a.m[i][j]=0;
}
a=power(a,nn,mod);
long long ans=0;
ans=ao*bo%mod*a.m[0][4];
ans=ans%mod;
ans+=bo*a.m[1][4];
ans=ans%mod;
ans=ans+ao*a.m[2][4]%mod;
ans%=mod;
ans+=a.m[3][4]%mod;
ans%=mod;
cout<<(ans+mod)%mod<<endl;
// printf("%I64d\n",&ans);
}
}