51nod1126【矩阵快速幂】

思路:
自己的一点心得:中间矩阵为最终矩阵。
搞出来很简单的;

#include <bits/stdc++.h>
using namespace std;

const int N=1e2+10;
const int mod=7;

struct asd{
    int num[3][3];
};


asd mul(asd a,asd b)
{
    asd ans;
    memset(ans.num,0,sizeof(ans.num));
    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++)
            for(int k=0;k<2;k++)
                ans.num[i][j]=(ans.num[i][j]+a.num[i][k]*b.num[k][j])%mod;
    return ans;
}
asd quickmul(int g,asd x)
{
    asd ans;
    for(int i=0;i<2;i++)
        for(int j=0;j<2;j++)
        {
            if(i==j)
                ans.num[i][j]=1;
            else
                ans.num[i][j]=0;
        }
    while(g)
    {
        if(g%2) ans=mul(ans,x);
        x=mul(x,x);
        g>>=1;
    }
    return ans;
}
int main()
{
    int c,a,b,n;
    asd x,ans;
    scanf("%d%d%d",&a,&b,&n);
    if(n==1||n==2)
    {
        puts("1");
        return 0;
    }
    x.num[0][0]=a;x.num[0][1]=b;
    x.num[1][0]=1;x.num[1][1]=0;
    ans=quickmul(n-2,x);
    c=ans.num[0][0]+ans.num[0][1];
    while(c<0)
        c+=mod;
    printf("%d\n",c%mod);
    return 0;
}
posted @ 2016-09-05 23:24  see_you_later  阅读(94)  评论(0编辑  收藏  举报