魔力手环，黑科技的矩阵快速幂优化（n3变n2）规律吗，只要愿找，总还是有的！

题目链接

优化：矩阵很好推

因为就是两两和

n=3的时候就是

1 0 1

1 1 0

0 1 1

还有就是这个矩阵进行求幂次的话，你会发现不可描述的规律，（自己算下，不要偷懒）；又因为初始矩阵除了第一行其他全是0，也是n2复杂度，轻松过题

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<algorithm>
#include<stack>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<string>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f;
struct Matrix
{
    int a[205][205];
    Matrix(){memset(a,0,sizeof(a));}
};
Matrix tmp,pre;
int num[205];
int n,k;
void multi1(Matrix p1,Matrix p2)
{
    Matrix ans;
        for(int j=1;j<=n;j++)
        {
            for(int m=1;m<=n;m++)
            {
                ans.a[1][j]=(ans.a[1][j]+p1.a[1][m]*p2.a[m][j])%100;
            }
        }
    pre=ans;
}
void multi2(Matrix p1,Matrix p2)
{
    Matrix ans;
    for(int j=1;j<=n;j++)
    {
        for(int i=1;i<=n;i++)
        {
            ans.a[1][j]=(ans.a[1][j]+p1.a[i][j]*p2.a[1][i])%100;
        }
    }
    for(int i=2;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(j==1)ans.a[i][j]=ans.a[i-1][n];else ans.a[i][j]=ans.a[i-1][j-1];
        }
    }
    //return ans;
    tmp=ans;
}
void qpow()
{
    while(k)
    {
        if(k&1)
            multi1(pre,tmp);
        k>>=1;
        multi2(tmp,tmp);
    }
}

int main()
{
    scanf("%d%d",&n,&k);
    for(int i=1;i<=n;i++){if(i==1)tmp.a[1][n]=tmp.a[1][1]=1;else tmp.a[i][i]=tmp.a[i][i-1]=1;}
    for(int i=1;i<=n;i++){scanf("%d",&num[i]);pre.a[1][i]=num[i];}
    qpow();
    for(int i=1;i<=n;i++){printf("%d%c",pre.a[1][i],i==n?'\n':' ');}
    return 0;
}