voj 1067 经典矩阵7 递推+矩阵快速幂

题目：https://vijos.org/p/1067

感觉网上其他题解有点奇怪

这个图写错了，矩阵写到前面去。

#include<bits/stdc++.h>using namespace std;typedef long long int ll;const ll mod = 7777777;const int maxn=11;struct matrix{    ll arr[maxn][maxn];    matrix operator*(matrix b){        matrix ans;        ll tmp;        for(int i=0; i<maxn; i++)        for(int j=0; j<maxn; j++){            ans.arr[i][j] = 0;            for(int k=0; k<maxn; k++){                tmp = (arr[i][k]*b.arr[k][j])%mod;                ans.arr[i][j] = (ans.arr[i][j] + tmp)%mod;            }        }        return ans;    }};matrix quick_pow(matrix a,ll N){    matrix ans;    memset(ans.arr,0,sizeof(ans.arr));    for(int i=0; i<maxn; i++)        ans.arr[i][i] = 1;    while(N){        if(N&1)            ans = ans*a;        a = a*a;        N /= 2;;    }    return ans;}int main(){    ll k,n,tmp[11];    matrix a;    scanf("%lld%lld",&k,&n);    memset(a.arr,0,sizeof(a.arr));    for(int i=0;i<k;i++){        if(i!=k-1)            a.arr[i][i+1]=1;        else            for(int j=0;j<k;j++)                a.arr[i][j]=1;    }    for(int i=0;i<k;i++)        tmp[i]=pow(2,i);    a=quick_pow(a,n-k);    ll ans=0;    for(int i=0;i<k;i++){        //cout<<tmp[i]<<" ";        ans+=a.arr[k-1][i]*tmp[i];        ans%=mod;    }    printf("%lld\n",ans);    /*for(int i=0;i<k;i++){        for(int j=0;j<k;j++){            printf("%lld ",a.arr[i][j]);        }        cout<<endl;    }*/    return 0;}