POJ 3233 矩阵快速幂&二分

题意：
给你一个n*n的矩阵 让你求S:
思路：
只知道矩阵快速幂 然后nlogn递推是会TLE的。

所以呢 要把那个n换成log

那这个怎么搞呢

二分！

当k为偶数时：

当k为奇数时：

就按照这么搞就能搞出来了
（我是看的题解才A的，，， 中间乱搞的时候犯了一些脑残的错误）

// by SiriusRen
#include <cstdio>
#include <cstring>
using namespace std;
int n,mod,k;
struct matrix{int a[33][33];void init(){memset(a,0,sizeof(a));}}first;
matrix mul(matrix a,matrix b){
matrix temp;temp.init();
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++)
temp.a[i][j]=(temp.a[i][j]+a.a[i][k]*b.a[k][j])%mod;
return temp;
}
matrix add(matrix a,matrix b){
matrix temp;temp.init();
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++)
temp.a[i][j]=(a.a[i][j]+b.a[i][j])%mod;
return temp;
}
matrix pow(matrix a,int x){
matrix temp;
x--;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
temp.a[i][j]=a.a[i][j];
while(x){
if(x&1)temp=mul(temp,a);
a=mul(a,a),x>>=1;
}
return temp;
}
matrix recursive(int x){
if(x==1)return first;
matrix temp=recursive(x/2);
if(x&1){
matrix jy=pow(first,x/2+1);
temp=add(temp,mul(temp,jy));
return add(jy,temp);
}
else{
matrix jy=pow(first,x/2);
return add(mul(jy,temp),temp);
}
}
int main(){
scanf("%d%d%d",&n,&k,&mod);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++){
scanf("%d",&first.a[i][j]);
first.a[i][j]=first.a[i][j]%mod;
}
matrix jy=recursive(k);
for(int i=1;i<=n;i++){
for(int j=1;j<n;j++)
printf("%d ",jy.a[i][j]);
printf("%d\n",jy.a[i]
);
}
}