1. 程式人生 > >bzoj 1879: [Sdoi2009]Bill的挑戰【狀壓dp】

bzoj 1879: [Sdoi2009]Bill的挑戰【狀壓dp】

處理 %s while ++ const ring sin amp namespace

石樂誌寫容斥……其實狀壓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;
}

bzoj 1879: [Sdoi2009]Bill的挑戰【狀壓dp】