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 longstruct 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;}