POJ3420状压dp和矩阵快速幂的结合

相当好的一道题。
首先这道题排除轮廓线的方法，很明显不可能，数据这么大。
不过我们可以从这里开始思考。
只要上一排的顺序排列好，那么下一排的方法是固定的，下一排的方法固定，下下一排就是固定的，便是机械的固定操作。
我们这里知道了，一共就只有4列，其实可以dp出4列，两行的所有情况。
那么往后dp，便是一系列一摸一样的操作，这一行与上一行匹配，那么下一行就出来了.假设一个dp数组，代表能不能转移。1表示能，0表不能，res数组表示方法数目。
于是就有res[pre][new]+=dp[pre][new1]*dp[new1][new]
可以看出来，等号左边就类似于矩阵乘法，一个行不变，一个列不变，相乘做相加运算。
所以就可表达成为res=dpⁿ。dp是一个16乘以16的矩阵，因为一共有0~1111种方法。
最后取res[15][15]因为相当于考虑在边框的外面一行是全1行，这样才能保证我们的第一行是正确的，再加上最后一行是满的，就可以保证我们的做法是合理的做法，因为dfs里面每一个都是合法操作。
精髓在于dp的值，只有1和1相乘才能为1.
相当于找，从1111到1111的有向图中，路径长为n的路径数目。（离散数学）

#include<iostream>
#include<string>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<climits>
#include<stack>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<cstdio>
#define up(i,a,b)  for(int i=a;i<b;i++)
#define dw(i,a,b)  for(int i=a;i>b;i--)
#define upd(i,a,b) for(int i=a;i<=b;i++)
#define dwd(i,a,b) for(int i=a;i>=b;i--)
//#define local
typedef long long ll;
const double esp = 1e-6;
const double pi = acos(-1.0);
const long long INF = 0x3f3f3f3f;
using namespace std;
typedef pair<int, int> pir;
int n, m;
struct mat {
int  a[16][16];
};
mat mod;
mat mat_mul(mat x, mat y)
{
	mat res;
	memset(res.a, 0, sizeof(res.a));
	up(i, 0, 16)
	{
		up(j, 0, 16)
		{
			up(k, 0, 16)
			{
				res.a[i][j] = (res.a[i][j] + (x.a[i][k] * y.a[k][j])%m) % m;
			}
		}
	}
	return res;
}
mat mat_pow()
{
	mat res;
	memset(res.a, 0, sizeof(res.a));
	up(i, 0, 16)
		res.a[i][i] = 1;
	while (n)
	{
		if (n & 1)res = mat_mul(res, mod);
		mod = mat_mul(mod, mod);
		n >>= 1;
	}
	return res;
}
void dfs(int i,int news,int old)
{
	if (i > 4)return;
	if (i == 4)
	{
		mod.a[old][news] = 1;
	}
	dfs(i + 1, news << 1, old << 1 | 1);
	dfs(i + 1, news << 1 | 1, old << 1);
	dfs(i + 2, news << 2 | 3, old << 2|3);
}
int main()
{
	while (cin >> n >> m&&n&&m)
	{
		memset(mod.a, 0, sizeof(mod.a));
		dfs(0, 0, 0);
		mat res = mat_pow();
		cout << res.a[15][15] << endl;
	}
	return 0;
}

这道题目a完，认真思考一下，其实很多dp都可以用矩阵来优化。
因为我们可以发现，这种具有滚动数组性质的，和前后状态不变的，都可以用
res+=dp[new][pre]*dp[pre][new1]大致来刻画。
比如类似dp[i]=dp[i-a]+dp[i-b]的，都可以用矩阵来刻画，只是系数不一样。