BZOJ-1879 Bill的挑战 状态压缩DP

MD....怎么又是状压.......

1879: [Sdoi2009]Bill的挑战
Time Limit: 4 Sec Memory Limit: 64 MB
Submit: 537 Solved: 280
[Submit][Status][Discuss]

Description

Input
本题包含多组数据。 第一行：一个整数T，表示数据的个数。 对于每组数据： 第一行：两个整数，N和K（含义如题目表述）。 接下来N行：每行一个字符串。

Output
1
2 1
a?
?b

Sample Input
50

Sample Output
对于30%的数据，T ≤ 5，M ≤ 5，字符串长度≤ 20；
对于70%的数据，T ≤ 5，M ≤ 13，字符串长度≤ 30；
对于100%的数据，T ≤ 5，M ≤ 15，字符串长度≤ 50。

HINT

Source
Day2

一看数据范围，50？感觉不能状压啊...哦，M<=15.....

思路比较简单：
f【i】【j】表示 匹配到第i位，时状态为j的方案数；
具体的转移：
f[i+1][j&(g[i][l])]+=f[i][j],f[i+1][j&(g[i][l])]%=p;
用g来存储状态，枚举状态即可；

PS：当天晚上调了20多分钟没调完....被YveH神犇叫去颓废了，Carry了他一盘后睡觉去了...第二天调完1A辣

code：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
int read()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
#define p 1000003
int t,n,k;
char s[20][60];

int f[60][1<<15],g[60][1<<5];
void DP()
{
    memset(f,0,sizeof(f)); memset(g,0,sizeof(g));
    int len=strlen(s[1]);
    for (int i=0; i<len; i++)
        for (int j=1; j<=n; j++)
            for (int l=0; l<26; l++)
                if (s[j][i]=='?' || s[j][i]=='a'+l)
                    g[i][l]|=1<<(j-1);
    f[0][(1<<n)-1]=1;
    for (int i=0; i<len; i++)
        for (int j=0; j<(1<<n); j++)
            if (f[i][j]!=0)
                for (int l=0; l<26; l++)
                    f[i+1][j&(g[i][l])]+=f[i][j],f[i+1][j&(g[i][l])]%=p;
    int ans=0;
    for (int i=0; i<(1<<n); i++)
        {
            int now=i,tmp=0;
            while (now) tmp+=now&1,now>>=1;
            if (tmp==k) ans=(ans+f[len][i])%p;
        }
    printf("%d\n",ans);
}

int main()
{
    t=read();
    while (t--)
        {
            n=read(),k=read();
            for (int i=1; i<=n; i++)
                scanf("%s",s[i]);
            DP();
        }
    return 0;
}