hihocoder #1162 矩阵加速dp

#1162 : 骨牌覆盖问题·三

时间限制:10000ms

单点时限:1000ms

内存限制:256MB

描述

前两周里，我们讲解了2xN,3xN骨牌覆盖的问题，并且引入了两种不同的递推方法。
这一次我们再加强一次题目，对于给定的K和N，我们需要去求KxN棋盘的覆盖方案数。

提示：KxN骨牌覆盖

输入

第1行：2个整数N。表示棋盘宽度为k，长度为N。2≤K≤7，1≤N≤100,000,000

输出

第1行：1个整数，表示覆盖方案数 MOD 12357

样例输入

2 62247088

样例输出

1399

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int MOD=;
int m,n;
struct Lu { int mx[][]; }L1,L2;
void dfs(int now,int next,int col)
{
    if(col==m){
        L1.mx[now][next]=;//表示now状态能够转移到next状态
        return;
    }
    dfs((now<<)+,next<<,col+);//不放
    dfs(now<<,(next<<)+,col+);//竖放
    if(col+<=m) dfs((now<<)+,(next<<)+,col+);//横放
}
Lu multi(Lu a,Lu b)
{
    Lu c;
    for(int i=;i<(<<m);i++){
        for(int j=;j<(<<m);j++){
            c.mx[i][j]=;
            for(int k=;k<(<<m);k++)
                c.mx[i][j]=(c.mx[i][j]+a.mx[i][k]*b.mx[k][j])%MOD;
        }
    }
    return c;
}
int solve(int x)
{
    memcpy(L2.mx,L1.mx,sizeof(L1.mx));
    while(x){
        if(x&) L1=multi(L1,L2);
        L2=multi(L2,L2);
        x>>=;
    }
    return L1.mx[(<<m)-][(<<m)-];//最后放满，转移到末状态
}
int main()
{
    while(scanf("%d%d",&m,&n)==){
        memset(L1.mx,,sizeof(L1.mx));
        dfs(,,);
        int ans=solve(n-);
        printf("%d\n",ans);
    }
    return ;
}