BZOJ1879bill的挑战——状压DP

最新推荐文章于 2020-08-09 00:23:58 发布

stevensonson

最新推荐文章于 2020-08-09 00:23:58 发布

阅读量212

点赞数 1

分类专栏： BZOJ

本文链接：https://blog.csdn.net/stevensonson/article/details/80227043

版权

BZOJ 专栏收录该内容

94 篇文章 1 订阅

订阅专栏

Description
这里写图片描述
Input

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

如题

Sample Input

3 3

???r???

???????

3 4

???????

?????a?

???????

3 3

???????

?a??j??

????aa?

3 2

a??????

???????

3 2

???????

???a???

????a??
Sample Output

914852

871234

67018

这道题一看就是状压DP，然而我们很难直接推出。我们可以换一个思路，先预处理出数组g。

数组g[i][j]表示字符串到第i位，并且这位上选择小写字母j所符合的个数。（用二进制表示）那么我们可以推出g[i][j]的表示方式：

for(int i=0;i<len;i++)
         for(int j=0;j<26;j++)
          for(int k=1;k<=n;k++) if(s[k][i]=='?'||s[k][i]==(j+'a')) g[i][j]|=1<<k-1;

我们预处理出g数组之后就可以发现一个DP了：

f[(1<<n)-1][0]=1;
        for(int i=0;i<len;i++)
         for(int j=0;j<(1<<n);j++)
          if(f[j][i]) for(int k=0;k<26;k++) f[j&g[i][k]][i+1]=(f[j&g[i][k]][i+1]+f[j][i])%md;
        ans=0;

这个DP的f[i][j]表示二进制状态为i，dp到第j位时的方案数，最后输出所有的f[i][len]，其中i满足二进制下有k个1即可。

#include<bits/stdc++.h>
#define md 1000003
using namespace std;
int read(){
    char c;int x;while(c=getchar(),c<'0'||c>'9');x=c-'0';
    while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
} 
int T,n,m,len,ans,g[105][105],f[1<<15][105];
string s[105];
int main()
{
    T=read();
    while(T--){
        n=read();m=read();
        memset(g,0,sizeof(g));memset(f,0,sizeof(f));
        for(int i=1;i<=n;i++) cin>>s[i];
        len=s[1].length();
        for(int i=0;i<len;i++)
         for(int j=0;j<26;j++)
          for(int k=1;k<=n;k++) if(s[k][i]=='?'||s[k][i]==(j+'a')) g[i][j]|=1<<k-1;
        f[(1<<n)-1][0]=1;
        for(int i=0;i<len;i++)
         for(int j=0;j<(1<<n);j++)
          if(f[j][i]) for(int k=0;k<26;k++) f[j&g[i][k]][i+1]=(f[j&g[i][k]][i+1]+f[j][i])%md;
        ans=0;
        for(int i=0;i<(1<<n);i++){
            int p=i,sum=0;
            while(p){sum+=p&1;p>>=1;}
            if(sum==m) ans=(ans+f[i][len])%md;
        }
        printf("%d\n",ans);
    }
    return 0;
}