BZOJ3195: [Jxoi2012]奇怪的道路（状压DP）

传送门

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

题解：
一眼以为是状压裸题，仔细想想好像复杂度不对，不过还是厚着脸皮交了，居然过了，好像比正解慢二三十倍。。又水了一道题。

考虑$f[i][j][m]$表示第$i$个点，前$K$个点度数奇偶性为$j$，现在一共连接了$m$条边。转移$f[i-1][j_2][m_2]\rightarrow f[i][j][m]$要满足$j_2$的第一位不能是奇数（因为后面的点无法改变前面的奇偶性）。同时$count(j\,\verb '^'\, j_2)\verb '&'1 == (m-m_2)\verb'&'1$（$count$表示二进制下为1的位数，这里要处理第 K <script type="math/tex" id="MathJax-Element-187">K</script>位和最后一位，注意细节）。显然的是先给每一个奇偶性改变的点连一条边，之后可以给任意一个点连两条边，这是一个组合数问题。直接暴力DP就好了。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
using namespace std;
const int Maxn=(1<<9)+20,Maxm=35;
const int Mod=1e9+7;
int n,m,k,f[2][Maxn][Maxm],cnt[Maxn],now,mxpos[Maxn],C[150][150];

int main(){
    scanf("%d%d%d",&n,&m,&k);
    f[now][0][0]=1;
    static int lim=(1<<k)-1;
    for(int i=1;i<=lim;i++){
        int res=0;
        for(int j=0;(1ll<<j)<=i;++j){
            if(i&(1ll<<j))++res,mxpos[i]=j;
        }
        cnt[i]=res;
    }
    for(int i=0;i<=m;i++)C[i][0]=1;
    for(int i=1;i<=100;i++){
        for(int j=1;j<=100;j++){
            (C[i][j]=C[i-1][j-1]+C[i-1][j])%=Mod;
        }  
    }
    mxpos[0]=-1;f[now][0][0]=1;
    for(int i=2;i<=n;i++){
        now^=1;
        for(int j=0;j<=lim;j++){
            memset(f[now][j],0,sizeof(int)*(m+1));
            if(mxpos[j]>=i)continue;
            for(int l=0;l<=lim;l++){
                if(mxpos[l]>=i-1)continue;
                int c=cnt[(l^(l&(1<<k)))^(j>>1)];
                if((c&1)!=(j&1))continue;
                for(int t=c;t<=m;t++){
                    for(int o=(t-c)&1;o<=t-c;o+=2){
                        (f[now][j][t]+=1ll*f[now^1][l][o]*C[(t-o-c)/2+min(i-2,k-1)][min(i-2,k-1)]%Mod)%=Mod;
                    }
                }
            }
        }
    }
    printf("%d\n",f[now][0][m]);
}