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[5][5];
};
ma mul(ma a,ma b)
{
    ma ans;
    for(int i=0;i<5;i++)
    for(int j=0;j<5;j++){
        ans.a[i][j]=0;
        for(int k=0;k<5;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&1) ans=mul(ans,a);
        n>>=1;
        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=0;i<5;i++)
        for(int j=0;j<5;j++){
                p.a[i][j]=0;
                ans.a[i][j]=0;
        }
    p.a[0][0]=p.a[1][0]=p.a[4][4]=1;
    while(scanf("%I64d",&n)!=EOF){
        scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&a0,&ax,&ay,&b0,&bx,&by);
        p.a[1][1]=(ax*bx)%mod;p.a[3][1]=(ax*by)%mod;
        p.a[2][1]=(ay*bx)%mod;p.a[4][1]=(ay*by)%mod;
        p.a[2][2]=bx%mod;p.a[4][2]=by%mod;
        p.a[3][3]=ax%mod;p.a[4][3]=ay%mod;
        ans.a[0][1]=a0*b0%mod;ans.a[0][3]=a0%mod;
        ans.a[0][2]=b0%mod;ans.a[0][4]=1;
        re=pow(p,n,ans);
        printf("%I64d
",re.a[0][0]%mod);
    }
    return 0;
}