POJ 3420 Quad Tiling

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;
}