bzoj 1879: [Sdoi2009]Bill的挑战【状压dp】

石乐志写容斥……其实状压dp就行
设f[i][s]表示前i个字母，匹配状态为s，预处理g[i][j]为第i个字母是j的1~n的集合，转移的时候枚举26个字母转移，最后答案加上正好有k个的方案即可

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int mod=1000003;
int T,n,m,len,t,f[55][50005],g[55][27],ans;
char c[20][55];
void jia(int &x,int y)
{
x+=y;
if(x>mod)
x-=mod;
}
int main()
{
scanf("%d",&T);
while(T--)
{
memset(f,0,sizeof(f));
memset(g,0,sizeof(g));
ans=0;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
scanf("%s",c[i]+1);
len=strlen(c[1]+1);
f[0][(1<<n)-1]=1;
for(int i=1;i<=len;i++)
for(int j=0;j<26;j++)
for(int k=1;k<=n;k++)
if(c[k][i]=='a'+j||c[k][i]=='?')
g[i][j]|=1<<(k-1);
for(int i=1;i<=len;i++)
for(int j=0,l=(1<<n);j<l;j++)
if(f[i-1][j])
for(int k=0;k<26;k++)
jia(f[i][g[i][k]&j],f[i-1][j]);
for(int i=0;i<(1<<n);i++)
{
int sm=0;
for(int j=0;j<n;j++)
if(i&(1<<j))
sm++;
if(sm==m)
jia(ans,f[len][i]);
}
printf("%d\n",ans);
}
return 0;
}