【BZOJ1879】【SDOI2009】Bill的挑战 [状压DP]

Bill的挑战

Time Limit: 4 Sec  Memory Limit: 64 MB
[Submit][Status][Discuss]

Description

　　

Input

　　第一行：一个整数T，表示数据的个数。 

　　对于每组数据： 

　　　　第一行：两个整数，N和K（含义如题目表述）。 

　　　　接下来N行：每行一个字符串。

Output

　　T行，每行一个整数表示答案

Sample Input

　　1
　　2 1
　　a?
　　?b

Sample Output

　　50

HINT

　　T ≤ 5，M ≤ 15，字符串长度≤ 50。

Solution

　　我们运用状压DP，令 g[i][c] 表示第 i 位，用 字符c 来匹配可行的串的集合。

　　然后显然就可以DP啦！f[i][opt] 表示做到了第 i 位，匹配集合为opt的方案数。

Code

 #include<iostream>
 #include<string>
 #include<algorithm>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<cmath>
 #include<queue>
 using namespace std;
 typedef long long s64;

 const int ONE = 4e5 + ;
 const int MOD = ;

 int n, m, T;
 int g[][], f[][];
 char s[][];
 int Ans;

 int get()
 {
         int res=,Q=;char c;
         while( (c=getchar())< || c> )
         if(c=='-')Q=-;
         res=c-;
         while( (c=getchar())>= && c<= )
         res=res*+c-;
         return res*Q;
 }

 void Deal()
 {
         memset(g, , sizeof(g));
         memset(f, , sizeof(f));
         n = get();  m = get();
         for(int i = ; i <= n; i++)
             scanf("%s", s[i] + );

         int len = strlen(s[] + );
         for(int i = ; i <= len; i++)
             for(int c = ; c <= ; c++)
                 for(int j = ; j <= n; j++)
                     if(s[j][i] == '?' || s[j][i] == c + 'a' - )
                         g[i][c] |=  << j - ;

         int total = ( << n) - ;
         f[][total] = ;
         for(int i = ; i <= len; i++)
             for(int opt = ; opt <= total; opt++)
                 if(f[i][opt])
                     for(int c = ; c <= ; c++)
                         (f[i + ][opt & g[i][c]] += f[i][opt]) %= MOD;

         Ans = ;
         for(int opt = ; opt <= total; opt++)
         {
             int num = ;
             for(int j = ; j <= n; j++)
                 if(opt & ( << j - )) num++;
             if(num == m) Ans = (Ans + f[len + ][opt]) % MOD;
         }

         printf("%d\n", Ans);
 }

 int main()
 {
         T = get();
         while(T--)
             Deal();
 }