hdu 4686 Arc of Dream

思路：构造矩阵

a[i]*b[i]=ax*bx*a[i-1]*b[i-1]+ax*by*a[i-1]+ay*bx*b[i-1]+ay*by

代码如下：

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cmath>
 #include<algorithm>
 #include<map>
 #define ll __int64
 #define mod 1000000007
 #define phi 1000000006
 using namespace std;
 struct ma
 {
     ll a[][];
 };
 ma mul(ma a,ma b)
 {
     ma ans;
     for(int i=;i<;i++)
     for(int j=;j<;j++){
         ans.a[i][j]=;
         for(int k=;k<;k++)
             if(a.a[i][k]&&b.a[k][j])
                 ans.a[i][j]=(ans.a[i][j]+a.a[i][k]*b.a[k][j])%mod;
         ans.a[i][j]%=mod;
     }
     return ans;
 }
 ma pow(ma a,ll n,ma ans)
 {
     while(n){
         if(n&) ans=mul(ans,a);
         n>>=;
         a=mul(a,a);
     }
     return ans;
 }
 int main()
 {
     ll n,an,s,ab,a,b;
     ll a0,b0,ax,ay,bx,by;
     ma p,ans,re;
     for(int i=;i<;i++)
         for(int j=;j<;j++){
                 p.a[i][j]=;
                 ans.a[i][j]=;
         }
     p.a[][]=p.a[][]=p.a[][]=;
     while(scanf("%I64d",&n)!=EOF){
         scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&a0,&ax,&ay,&b0,&bx,&by);
         p.a[][]=(ax*bx)%mod;p.a[][]=(ax*by)%mod;
         p.a[][]=(ay*bx)%mod;p.a[][]=(ay*by)%mod;
         p.a[][]=bx%mod;p.a[][]=by%mod;
         p.a[][]=ax%mod;p.a[][]=ay%mod;
         ans.a[][]=a0*b0%mod;ans.a[][]=a0%mod;
         ans.a[][]=b0%mod;ans.a[][]=;
         re=pow(p,n,ans);
         printf("%I64d\n",re.a[][]%mod);
     }
     return ;
 }