speroc的博客

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int mod = 10000;int n;struct mat{ long long a[5][5];};mat operator *(mat x,mat y){

快速幂int ksm(int a,int b){int i=1,x=a;while(b!=0){if(b&1)i*=x;x*=x;b>>=1;}return i;}矩阵乘法struct mat{ //定义矩阵结构体 long long a[maxn][maxn];};mat operator*(mat x,mat y){ //重载乘号(定义