BZOJ 3195: [Jxoi2012]奇怪的道路 状压dp

3195: [Jxoi2012]奇怪的道路

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 479  Solved: 307
[Submit][Status][Discuss]

Description

小宇从历史书上了解到一个古老的文明。这个文明在各个方面高度发达，交通方面也不例外。考古学家已经知道，这个文明在全盛时期有n座城市，编号为1..n。m条道路连接在这些城市之间，每条道路将两个城市连接起来，使得两地的居民可以方便地来往。一对城市之间可能存在多条道路。
据史料记载，这个文明的交通网络满足两个奇怪的特征。首先，这个文明崇拜数字K，所以对于任何一条道路，设它连接的两个城市分别为u和v，则必定满足1 <=|u - v| <= K。此外，任何一个城市都与恰好偶数条道路相连（0也被认为是偶数）。不过，由于时间过于久远，具体的交通网络我们已经无法得知了。小宇很好奇这n个城市之间究竟有多少种可能的连接方法，于是她向你求助。
方法数可能很大，你只需要输出方法数模1000000007后的结果。

Input

输入共一行，为3个整数n，m，K。

Output

输出1个整数，表示方案数模1000000007后的结果。

Sample Input

输入样例1
3 4 1
输入样例2
4 3 3

Sample Output

输出样例1
3
输出样例2
4
数据规模

HINT

100%的数据满足

1<= n <= 30, 0 <= m <= 30, 1 <= K <= 8.

【题目说明】
两种可能的连接方法不同当且仅当存在一对城市，它们间的道路数在两种方法中不同。
在交通网络中，有可能存在两个城市无法互相到达。

转载

考试的时候我真的没有想到……看到了数据范围比较小，但是我想的是深搜而不是状压……

我们发现k的数据范围很小，所以我们考虑状压，首先考虑状态数组定义。

首先，f数组肯定有2维表示枚举到的点数和边数。我们发现，题目中对每个点的奇偶性有限制，而奇偶是两种相反的概念，

因此我们尝试再用一位表示可以和这个点连边的点的出度奇偶性状态，1表示奇；但是我们发现，如果这个设置为所有它能转移到的点的奇偶性（一共16位），

不仅时间复杂度变大，转移时候的讨论也会很复杂，而我们考虑一下，其实只弄一边的状态就可以了，从现在的点往后转移相当于从后面的点往前转移

因此我们只保存i点以及前面k个点的状态即可，我们再设置一维l，表示我们已经考虑到了前面k+1个点中的第l个点(我们设高位离着i点最近，第0位是点i-k，第k位是i点，一共k+1位）

于是我们得到了一个四维的状态数组，f[35][35][(1<<9)[10]，接下来我们考虑转移

考虑对第i-k+l个点与第i个点，如果我们不再加边，就直接转移到下一个f[i][j][state][l+1]

如果我们再还合法的情况下加边，会同时改变i点和i-k+l这两个的出度奇偶性，也就是转移到f[i][j+1][state^(1<<k)^(1<<l)][l]

然后我们考虑不同i之间的转移。如果当前状态要转移的话，我们现在第i-k个点永远不会再被考虑，因此转移的前提之一是第i-k位奇偶性是0，即state&1==0

并且，我们还可以发现一个很有用的性质：对于点i+1，它初始state就是第i个点的state>>1，这一点很显然，从我们第3维的定义中就能看出这一操作可以实现

这样的话，就有跨第一维的转移，即向f[i+1][j][u>>1][0]转移

再最后的时候，输出f[n+1][m][0][0]即可，代码见下：

#include<cmath>#include<ctime>#include<cstdio>#include<cstring>#include<cstdlib>#include<iostream>#include<algorithm>#include<iomanip>#include<vector>#include<bitset>#include<string>#include<queue>#include<set>#include<map>using namespace std;inline int read(){int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}return x*f;}void print(int x){if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}const int mod=1000000007,N=35;int f[N][N][1<<9][10],bin[16];int main(){int n=read(),m=read(),K=read();register int i,j,k,p;bin[0]=1;for(i=1;i<=15;++i)bin[i]=bin[i-1]<<1;f[1][0][0][0]=1;for(i=1;i<n;++i)for(j=0;j<=m;++j)for(k=0;k<(bin[K+1]);++k){for(p=0;p<K;++p)if(f[i][j][k][p]){(f[i][j][k][p+1]+=f[i][j][k][p])%=mod;if(j<m&&i-K+p>=0)(f[i][j+1][k^bin[p]^bin[K]][p]+=f[i][j][k][p])%=mod;}if((k&1)==0&&f[i][j][k][p])f[i+1][j][k>>1][0]=f[i][j][k][K];}print(f[n][m][0][0]);puts("");return 0;}/*3 4 134 3 34*/