数学专项matrix：UVa 10689

快速幂模板题。

#include <iostream>#include <cstdio>#include <cstring>using namespace std;typedef int matrix[3][3];int mod;int a,b,n,m;void mat_mul(matrix A,matrix B,matrix res){    matrix C;    memset(C,0,sizeof(C));    for(int i=0;i<2;i++)        for(int j=0;j<2;j++)            for(int k=0;k<2;k++)                C[i][j]=(C[i][j]+A[i][k]*B[k][j])%mod;    memcpy(res,C,sizeof(C));}void mat_pow(matrix A,int n,matrix res){    matrix a,r;    memset(r,0,sizeof(r));    memcpy(a,A,sizeof(a));    for(int i=0;i<2;i++) r[i][i]=1;    while(n)    {        if(n&1) mat_mul(r,a,r);        n>>=1;        mat_mul(a,a,a);    }    memcpy(res,r,sizeof(r));}int main(){    int T;    scanf("%d",&T);    while(T--)    {        scanf("%d%d%d%d",&a,&b,&n,&m);        mod=1;for(int i=0;i<m;i++) mod*=10;        a%=mod;b%=mod;        if(n==0)        {            printf("%d\n",a);            continue;        }        matrix A;        A[0][0]=0;A[0][1]=1;A[1][1]=1;A[1][0]=1;        mat_pow(A,n-1,A);        int ans=(A[1][0]*a%mod+A[1][1]*b%mod)%mod;        printf("%d\n",ans);    }    return 0;}