hdu 4686 Arc of Dream

矩阵快速幂：

1,A0%N,B0%N,A0*B0%N,A0*B0%N;

1,AY%N,BY%N,AY*BY%N,AY*BY%N;

AX%N, 0,AX*BY%N,AX*BY%N;

BX%N,BX*AY%N,BX*AY%N;

AX*BX%N,AX*BX%N;

1;

#include <cstdio>
#include <cstring>
#define N 1000000007
#define LL long long
struct matrix
{
    LL num[5][5];
};
matrix p,m,q;
matrix cal(matrix x,matrix y)
{
    for(int i=0; i<5; ++i)
        for(int j=0; j<5; ++j)
        {
            q.num[i][j]=0;
            for(int k=0; k<5; ++k)
            q.num[i][j]=(q.num[i][j]+x.num[i][k]*y.num[k][j]%N)%N;
        }
    return q;
}
void solve(LL n)
{
    while(n)
    {
        if(n&1) m=cal(m,p);
        n>>=1;
        p=cal(p,p);
    }
}
int main()
{
//    freopen("in.txt","r",stdin);
    LL AX,AY,A0,BX,BY,B0,n;
    while(scanf("%I64d%I64d%I64d%I64d%I64d%I64d%I64d",&n,&A0,&AX,&AY,&B0,&BX,&BY)!=EOF)
    {
        if(n==0)
        {
            printf("0\n");
            continue;
        }
        memset(m.num,0,sizeof(m.num));
        memset(p.num,0,sizeof(p.num));
        m.num[0][0]=1,m.num[0][1]=A0%N,m.num[0][2]=B0%N,m.num[0][3]=A0*B0%N,m.num[0][4]=A0*B0%N;

        p.num[0][0]=1,p.num[0][1]=AY%N,p.num[0][2]=BY%N,p.num[0][3]=AY*BY%N,p.num[0][4]=AY*BY%N;
        p.num[1][1]=AX%N,p.num[1][3]=AX*BY%N,p.num[1][4]=AX*BY%N;
        p.num[2][2]=BX%N,p.num[2][3]=BX*AY%N,p.num[2][4]=BX*AY%N;
        p.num[3][3]=AX*BX%N,p.num[3][4]=AX*BX%N;
        p.num[4][4]=1;
        solve(n-1);
        printf("%I64d\n",m.num[0][4]);
    }
    return 0;
}