hdu 4686 构造矩阵乘法 待更新

这题给我坑坏了 

构造矩阵乘法，时刻要注意%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);
    }
}