BZOJ 3195: [Jxoi2012]奇怪的道路 【状压DP】

题目描述：

n个点编号1~n，两个点u,v之间可以连边当且仅当1<=|u-v|<=K，两点之间可以连任意条边，要求每个点都恰好与偶数条道路相连。
给出n，K，以及总共的边数m，求连边的方案数
1<=n<=30，0<=m<=30，1<=K<=8

题目分析：

K很小，可以状压当前K个点已经连好边的奇偶性。
考虑两种方法:

• 逐行转移，设$f[i][j][s]$表示现在考虑第i个点向后连边，已经用了j条边，从i到i+K-1已经连好边的奇偶性为s
  我们需要枚举当前需要用的边x，下一个状态y，预处理用x条边使得s转移到y的方案数
  但是这样的做法时间复杂度是$O(2^K{2}^K{m}^{2}n)$的，跑不过

• 逐格转移，设$f[i][j][s][l]$表示现在考虑第i个点向前连边，已经用了j条边，从i到i-K的奇偶性为s，现在考虑i与i-K+l是否连边

  • 当前不连边，转移到f[i][j][s][l+1]
  • 连边，转移到f[i][j+1][s^(1<<K)(1<<l)][l]
  • 当第i-K位为偶时，f[i][j][s][K]转移到f[i+1][j][s>>1][0]

  惊奇地发现时间复杂度为$O(nmK{2}^K)$
  这种上一个状态和下一个状态没有限制关系的DP逐格转移的复杂度就是状态数

#include<cstdio>
const int mod = 1e9+7;
int n,m,K,f[33][33][1<<9][9];
int main()
{
	scanf("%d%d%d",&n,&m,&K);
	f[2][0][0][0]=1;
	for(int i=2;i<=n;i++)
		for(int j=0;j<=m;j++)
			for(int k=0;k<(1<<(K+1));k++)
			{
				for(int l=0;l<K;l++) if(f[i][j][k][l])
				{
					(f[i][j][k][l+1]+=f[i][j][k][l])%=mod;
					if(j<m&&i-K+l>0) (f[i][j+1][k^(1<<K)^(1<<l)][l]+=f[i][j][k][l])%=mod;
				}
				if((k&1)==0) (f[i+1][j][k>>1][0]+=f[i][j][k][K])%=mod;
			}
	printf("%d\n",f[n+1][m][0][0]);
}