JZOJ4702&codeforces534F 状压DP

题意就不说了。。
不说个毛啊！！！wori我就是题目没看懂所以搞了半天结果搞错了好吗QAQ。。
这里的块如果空间允许一定是两个，如果剩下的空间不足两个才能为1.
。。。。
怎么暴力怎么来。
设f[i][j][a][b][c][d][e]表示第i列，状态为j，5行的块数分别为a,b,c,d,e。
那么转移显然，把之前的全部加上就好了。
对于一个长度为m的行，我们最多能取的块数是m/2+m%2.
那么假设最终要取n个，那么这行最少取n-(n-m)/2-(n-m)%2.
然后就可以按照上面两个公式来枚举a,b,c,d,e。
由于每一列的限制是固定的，那么我们可以预处理出可行的状态，优化时间复杂度，我是设g[i][j]表示块数为i时，第j种状态。。。
细节比较多。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int mo=1e9+7;
int f[21][23][11][11][11][11][11],mi[6],a[6],c[6],b[21],g[11][33];
int n,m,ans,tot;
inline bool check(int y,int x)
{
    int last=0,ans=0;
    while (x)
    {
        if (last==0&&x%2==1)ans++;
        last=x%2;
        x/=2;
    }
    if (ans==y)return 1;
    return 0;
}
inline int calc(int x)
{
    return x/2+x%2;
}
int main()
{
    scanf("%d%d",&n,&m);
    mi[0]=1;
    fo(i,1,n)mi[i]=mi[i-1]<<1;
    fo(i,1,n)scanf("%d",&a[i]);
    fo(i,1,m)scanf("%d",&b[i]);
    fo(i,0,10)
    fo(j,0,mi[n]-1)
    if (check(i,j))g[i][++g[i][0]]=j;
    f[0][0][0][0][0][0][0]=1;
    fo(i,0,m-1)
    fo(j,0,mi[n]-1)
    fo(_1,max(a[1]-calc(m-i),0),min(a[1],calc(i)))
    fo(_2,max(a[2]-calc(m-i),0),min(a[2],calc(i)))
    fo(_3,max(a[3]-calc(m-i),0),min(a[3],calc(i)))
    fo(_4,max(a[4]-calc(m-i),0),min(a[4],calc(i)))
    fo(_5,max(a[5]-calc(m-i),0),min(a[5],calc(i)))
    if (f[i][j][_1][_2][_3][_4][_5])
    fo(k,1,g[b[i+1]][0])
    {
        int x=g[b[i+1]][k];//status k
        bool bz=0;
        c[1]=_1;
        c[2]=_2;
        c[3]=_3;
        c[4]=_4;
        c[5]=_5;
        fo(l,1,n)
        if (!(j&mi[l-1])&&(x&mi[l-1]))
        {
            c[l]++;
            if (c[l]>a[l])
            {
                bz=1;
                break;
            }
        }
        if (bz) continue;
        f[i+1][x][c[1]][c[2]][c[3]][c[4]][c[5]]+=f[i][j][_1][_2][_3][_4][_5];
        f[i+1][x][c[1]][c[2]][c[3]][c[4]][c[5]]%=mo;
    }
    fo(i,1,g[b[m]][0])
    ans=(ans+f[m][g[b[m]][i]][a[1]][a[2]][a[3]][a[4]][a[5]])%mo;
    printf("%d\n",ans);
}