1、本题考查状压DP+矩阵乘法优化。先用dfs求出状态转移的可达矩阵,然后用矩阵快速幂求出矩阵的n次方,an【0】【0】即为答案。
2、更加深刻地理解了状态转移。所谓状态转移就是某个状态通过一系列操作到达另一个状态(或者说“局面”,“情况”)。比如1011,有没有办法让它变为0100呢?有的,答案就是在那个唯一的0的位置摆上一个骨牌,那么新的一行也就只有那个位置是填满的,其余位置都是空的,即0100。其余情况依次类推,用dfs计算。
3、首先用dfs以0——15为终点求出的一个状态转移可达矩阵,状态转移,就相当于在有向图上连一条边,矩阵的n次幂的作用就相当于求从i到j的长度为n的道路有多少条(离散数学课本有证明,这里不赘述)。对应到本题,就是求从0000开始,经过n列的排放后,回到初始状态(0000)的方法数,即本题的答案。
4、注意,&的优先级低于==,要适当加括号。
5、这道题做得真辛苦,所幸一次AC~
#include<cstdio>
#include<cstring>
using namespace std;
int n,m;
const int MAXN=16;
const int MAXM=16;
struct Matrix{
int n,m;
int a[MAXN][MAXM];
void clear(){
n=m=0;
memset(a,0,sizeof(a));
}
Matrix operator * (const Matrix &b) const{
Matrix temp;
temp.clear();
temp.n=n;temp.m=b.m;
for(int i=0;i<n;i++)
for(int j=0;j<b.m;j++)
for(int k=0;k<m;k++)
temp.a[i][j]+=a[i][k]*b.a[k][j];
return temp;
}
Matrix operator % (const int mod) const{
Matrix temp;
temp.n=n;temp.m=m;
for(int i=0;i<n;i++)
for(int j=0;j<m;j++)
temp.a[i][j]=a[i][j]%mod;
return temp;
}
};
Matrix a,an;
Matrix pow_mod(Matrix A,int i,int mod){
if(i==0){
Matrix u;
u.clear();
u.n=MAXN;u.m=MAXM;
for(int i=0;i<u.n;i++)
u.a[i][i]=1%mod;
return u;
}
Matrix temp=pow_mod(A,i>>1,mod);
temp=temp*temp%mod;
if(i&1) temp=temp*A%mod;
return temp;
}
void dfs(int now,int final,int bit){
if(bit>4) return;
if(bit==4) {a.a[now][final]=1;return;}
int temp=(1<<bit);
if(final&temp) dfs(now,final,bit+1);
else{
dfs(now|temp,final,bit+1);
if((final&(temp<<1))==0) dfs(now,final,bit+2);
}
return;
}
int main(){
a.clear();
a.n=MAXN;a.m=MAXM;
for(int i=0;i<16;i++)
dfs(0,i,0);
while(scanf("%d%d",&n,&m)==2){
if(n==0&&m==0) break;
an=pow_mod(a,n,m);
printf("%d\n",an.a[0][0]);
}
return 0;
}