BZOJ1879: [Sdoi2009]Bill的挑战

字符串不是很多，可以考虑状态压缩
f[i][j]表示到匹配到第i位，匹配状态为j，j在二进制下就表示和这n个串的匹配情况
然后加个滚动
然而貌似有很多废的状态，直接DP复杂度好像也过不去，于是我就记录了当前这一层有哪些状态

---

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 16;
const int maxl = 60;
const int Mod = 1000003;

int f[2][1<<maxn];
int num[2][1<<maxn];
char st[maxn][maxl];
int n,m;

int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
            scanf("%s",st[i]);
        int len = strlen(st[1]);

        int now=1;
        memset(f[now],0,sizeof f[now]);
        memset(num[now],0,sizeof num[now]);
        f[now][(1<<n)-1]=1; num[now][0]=1; num[now][1]=(1<<n)-1;
        for(int i=0;i<len;i++)
        {
            now^=1;
            memset(f[now],0,sizeof f[now]);
            memset(num[now],0,sizeof num[now]);
            for(int j=1;j<=num[1^now][0];j++)
            {
                int ki=num[1^now][j];
                for(int l=0;l<26;l++)
                {
                    int ti=0,ni=0;
                    for(int k=1;k<=n;k++)
                    if(((1<<k-1)&ki)&&(st[k][i]=='?'||st[k][i]-'a'==l))
                        ni++,ti|=(1<<k-1);
                    if(ni>=m)
                    {
                        if(!f[now][ti])
                        {
                            num[now][0]++;
                            num[now][num[now][0]]=ti;
                            f[now][ti]=f[1^now][ki];
                        }
                        else
                        {
                            f[now][ti]+=f[1^now][ki];
                            f[now][ti]=(f[now][ti]-1)%Mod+1;
                        }
                    }
                }
            }
        }
        int r=0;
        for(int i=0;i<1<<n;i++)
        {
            int ni=0;
            for(int j=1;j<=n;j++) if(i&(1<<j-1)) ni++;
            if(ni==m) (r+=f[now][i])%=Mod;
        }
        printf("%d\n",r);
    }

    return 0;
}