acdream 1214 矩阵快速幂

题意：宽度从1到5的矩阵 然后进行染色 要求染色完以后不存在2X2的格子 完全同色的所有方案数

矩阵的长度小于10的100次

解法：看到长度的范围就应该很清楚了 矩阵快速幂啊 然后比赛中也是往这个方向去想的 然后讨论1～5各自宽度的矩阵

好吧 其实只要按照行考虑就可以了 每一行有的方案数就是2的次方啊 那么宽度5所需要的矩阵的就是32*32了

好吧 手推这个是不现实的 打表或者直接暴力处理就可以了

突然发现 这题需要大数 然后vector实现真是不能再赞了  想到自己一直没有掌握的快速幂和大数 这是什么一种心情呢？

感觉一下子就回到原点  然后又变成那个被吊打的白痴 这趟远行没有目的：这个世上没有什么偶然，有的只有必然

#include<cstdio>
#include<cstring>
#include<iostream>
#include<vector>
using namespace std;
int zt,p,n,m;
char s[111];
vector<int>g;
struct matrix{
	int mat[33][33];
	void clc(){memset(mat,0,sizeof(mat));}
}a,b;
// 00 11
// 00 11
int judge(int x,int y){
    for(int i=1;i<m;++i){
        if((x%4==0&&y%4==0)||(x%4==3&&y%4==3))return 0;
        x>>=1;y>>=1;
    }
    return 1;
}
matrix mul(matrix a,matrix b){
    matrix c;c.clc();
    for(int i=0;i<zt;++i)
        for(int j=0;j<zt;++j)
            for(int k=0;k<zt;++k)
                c.mat[i][j]=(c.mat[i][j]+a.mat[i][k]*b.mat[k][j])%p;
    return c;
}
void gao(){
    for(int i=(int)g.size()-1;i;--i){
        if(g[i]&1)g[i-1]+=10;
        g[i]/=2;
    }g[0]/=2;
    if(!g[g.size()-1])g.pop_back();
}
int main(){
    
    scanf("%s%d%d",s,&m,&p);
    
    g.clear();zt=(1<<m);
    int sz=(int)strlen(s);
    
    for(int i=sz-1;i>=0;--i)g.push_back(s[i]-'0');
    
//    ---to pow(m-1)
    --g[0];
    for(int i=0;i<g.size();++i)if(g[i]<0)g[i]+=10,g[i+1]--;
    if(!g[g.size()-1])g.pop_back();
    
    a.clc(),b.clc();
    for(int i=0;i<zt;++i)a.mat[i][i]=1;
    for(int i=0;i<zt;++i)
        for(int j=0;j<zt;++j)
            b.mat[i][j]=judge(i,j);
    
    while(!g.empty()){
        if(g[0]&1)a=mul(a,b);
        gao();b=mul(b,b);
    }

    int ans=0;
    for(int i=0;i<zt;++i)
        for(int j=0;j<zt;++j)
            ans=(ans+a.mat[i][j])%p;
    
    printf("%d\n",ans);
    
	return 0;
}