hit2255矩阵连乘求Fibonacci的和

   1 1 0 s(n-1) s(n)

(  0 p q  )   *   (  f(n)   )     =         (   f(n+1)     )

   0 1 0   f(n-1)   f(n)

注意不能因为m[0][2]=0在求结果的时候就不乘它。

另外这个题应该用long long。。。

还有，因为在题中有quickpow（s-1） 所以应该分为3中情况。。。。

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
#define MAX 3
#define Mod 10000000
typedef struct
{
    long long m[MAX][MAX];
}Matrix;
Matrix I={1,0,0,0,1,0,0,0,1};
Matrix P={1,1,0,0,0,0,0,1,0};
long long a,b,p,q,e,s;
Matrix matrixmul(Matrix a,Matrix b)
{
    Matrix c;
    int i,j,k;
    for(i=0;i<MAX;i++)
        for(j=0;j<MAX;j++)
        {
            c.m[i][j]=0;
            for(k=0;k<MAX;k++)
            {
                c.m[i][j] = ((c.m[i][j] + a.m[i][k] * b.m[k][j] % Mod) + Mod) % Mod;
            }
        }
    return c;
}
Matrix quickpow(long long n)
{
    Matrix b=I,m=P;
    while(n>=1)
    {
        if(n&1)
            b=matrixmul(b,m);
        n=n>>1;
        m=matrixmul(m,m);
    }
    return b;
}
int main()
{
    int t;
    long long ans1,ans2,ans,tmp1,tmp2,tmp3,tmp4,tmp5,tmp6;
    while(scanf("%d",&t)!=EOF){
    while(t--)
    {
        scanf("%lld%lld%lld%lld%lld%lld",&a,&b,&p,&q,&s,&e);
        P.m[1][1]=p;P.m[1][2]=q;
        Matrix x,y;
        x=quickpow(e);
        if(s>1)
        y=quickpow(s-1);

        tmp1=(x.m[0][0]*a)%Mod;
        tmp2=(x.m[0][1]*b)%Mod;
        tmp5=(x.m[0][2]*a)%Mod;
        ans1=(tmp1+tmp2+tmp5)%Mod;
        if(s>1){
        tmp3=(y.m[0][0]*a)%Mod;
        tmp4=(y.m[0][1]*b)%Mod;
        tmp6=(y.m[0][2]*a)%Mod;
        ans2=(tmp3+tmp4+tmp6)%Mod;
        }
        else if(s==0)   ans2=0;
        else ans2=a;
       // ans2=(ans2+Mod)%Mod;
        //cout<<"ans1="<<ans1<<" ans2="<<ans2<<endl;
        ans=(ans1-ans2)%Mod;
        ans=(ans+Mod)%Mod;
        printf("%lld\n",ans);
    }
    }
    return 0;
}