[HDU] A Simple Math Problem

矩阵快速幂裸题，可以当做模板题练手。

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cmath>

using namespace std;
const int N = 15;

int m;

struct matrix{
    int M[N][N];
    int l,r;
    
    void operator * (matrix x){
        int i,j,k; matrix t;
        for(i=1;i<=l;i++)
            for(j=1;j<=r;j++) t.M[i][j]=0;
        for(i=1;i<=r;i++)
            for(j=1;j<=x.l;j++)
                for(k=1;k<=l;k++)
                    t.M[j][i]=(t.M[j][i]+M[k][i]*x.M[j][k]) % m;
        for(i=1;i<=l;i++)
            for(j=1;j<=r;j++) M[i][j]=t.M[i][j];
    }
        
}c,res;

int main(){
    int i,j,n,k;
    while(scanf("%d%d",&k,&m)!=EOF){
        for(i=10;i>=1;i--)
            scanf("%d",&c.M[10][i]);
        for(i=1;i<=9;i++)
            for(j=1;j<=10;j++)
                c.M[i][j]=((j-1)==i);
        c.l=10; c.r=10;
        /*for(i=1;i<=10;i++){
            for(j=1;j<=10;j++) printf("%d",c.M[i][j]);
            cout<<endl;
        }*/
        for(i=0;i<=9;i++) res.M[i+1][1]=i;
        res.l=10; res.r=1;
        for(; k ; k >>= 1){
            if(k & 1){
                res * c;
            }
            c * c;
        }
        printf("%d\n",res.M[1][1]);
    }
    return 0;
}

转载于:https://www.cnblogs.com/notseefire/p/8010609.html