题意:宽度从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;
}