hdu 6050 Funny Function

 /*  打表找到公式，矩阵快速幂加速- -  马虎：ll 类型  写成int , 一直 wrong T.T  1-7^n / 1-7 f(n) =2^(n-1)+f(n-2) */ #include<iostream>#include<cstdio>#include<algorithm>#include<cstring>#include<cmath>using namespace std;typedef long long ll;const int mod =1e9+7;struct mex{    int a[4][4];}m1[6];mex Plus(mex a,mex b){    mex c;ll d;    memset(c.a,0,sizeof(c.a));    for(int i=1;i<=3;++ i)     for(int j=1;j<=3;++ j)      for(int k=1;k<=3; ++ k)        {                        d=a.a[i][k];            d *=b.a[k][j];            d %=mod;            c.a[i][j] =(c.a[i][j]+d)%mod;        }      return c;}ll quickm(ll n,ll m)  {      ll a=1,x=n;      while(m)      {          if(m&1) a =(a*x)%mod;          x=(x*x)%mod;          m /=2;      }      return a%mod;  }  void gcd(ll a,ll b,ll &d,ll &x,ll &y)  {      if(!b) {d=1;x=1;y=0;}      else      {          gcd(b,a%b,d,y,x); y-= x*(a/b);      }  }  int inc(int a,int mod)  {      ll x,y,d;      gcd(a,mod,d,x,y);      return d==1?(x+mod)%mod:-1;  }void init(){    m1[1].a[1][1]=1;    m1[1].a[2][1]=4;    m1[1].a[3][1]=2;    memset(m1[2].a,0,sizeof(m1[2].a));    m1[2].a[1][3]=1;m1[2].a[3][1]=1;m1[2].a[3][2]=1;    m1[2].a[2][2]=2;}int main(){    int t;    ll n,m;    scanf("%d",&t);    while(t--)    {        scanf("%lld%lld",&n,&m);        ll n0,nn=n;        if(m==1||n==1) { printf("1\n"); continue;}        init();        n0=n-1;n -=2;        while(n)        {                            if(n&1) m1[1]=Plus(m1[2],m1[1]);            m1[2]=Plus(m1[2],m1[2]);            n /=2;        }        int a1=m1[1].a[3][1],a0=0;        if(n0==1) a0=1;        if(!a0)        {            n=n0;n -=2;            init();            while(n)           {              if(n&1) m1[1]=Plus(m1[2],m1[1]);                            m1[2]=Plus(m1[2],m1[2]);              n /=2;           }           a0=m1[1].a[3][1];        }        if(m==2)         {            printf("%d\n",a1);            continue;         }        int dos=(a0+a1)%mod;        ll d=quickm(dos,m-2);        ll ans=(ll)d*a1%mod;        if(nn%2)        {            int ink=inc(dos-1,mod);                     ll dd =ink;dd *=d-1; dd%=mod; dd *=a0;dd %=mod;            ans -=dd;            ans =(ans+mod)%mod;        }        printf("%lld\n",ans);    }    return 0;}