CF482C Game with Strings (状压DP+期望DP)

Posted 2021-01-08 guapisolo

题目大意：甲和乙玩游戏，甲给出n(n<=50)个等长的字符串(len<=20)，然后甲选出其中一个字符串，乙随机询问该字符串某一位的字符(不会重复询问一个位置)，求乙能确定该串是哪个字符串的询问次数的期望值

这题不看题解好难想......(感谢zhx和zhx两位大佬的题解)

len很小，考虑状压DP，显然我们要状压询问，要定义两个状态，f[]和num[]

1表示询问，0表示未询问

那么，我们定义f[s]表示询问状态s距离确定一个字符串所需要的期望值。

 定义tot是s状态剩余的询问的次数，那么显然因为询问剩余的每一位的概率是相等的

而一个询问并不一定能确定唯一一个字符串，可能是好几个，所以我们要处理这个询问状态能确定几个字符串，即num[s]

而num[s]并不能通过暴力枚举得到，所以我们只能通过枚举它的父集来获取它的状态

我们把字符串两两匹配，定义为s状态不能确定的串(用二进制表示哪些串不能确定)，接着我们从大到小枚举状态，每个子集都从它的父集更新，取所有的并集，其中1的数量即为num

tip1：而如果1的数量为1，该状态仍然可以确定所有的串，所以此时num是0

最后我们得到一个很玄学的方程

最后的f[0]即为答案

tip2：当n==1时，出题人貌似想告诉我们，因为不论问不问都只有一个串，所以不用问也知道是哪个......

 1 #include <cstdio>
 2 #include <algorithm>
 3 #include <cstring>
 4 #define N 55
 5 #define maxn (1<<20)+100
 6 #define ll long long 
 7 #define dd double
 8 #define seed 13131
 9 using namespace std;
10 
11 int n,m;
12 char str[55][22];
13 dd f[maxn];
14 ll unidf[maxn];
15 int num[maxn];
16 
17 int main()
18 {
19     //freopen("aa.in","r",stdin);
20     scanf("%d",&n);
21     if(n==1) {printf("0.0000000000
");return 0;}
22     for(int i=0;i<n;i++)
23         scanf("%s",str[i]);
24     m=strlen(str[1]);
25     for(int i=0;i<n;i++)
26         for(int j=i+1;j<n;j++)
27         {
28             int pos=0;
29             for(int k=0;k<m;k++){
30                 if(str[i][k]==str[j][k])
31                    pos|=(1<<k);
32             }
33             unidf[pos]|=(1ll<<j)|(1ll<<i);
34         }
35     for(int s=(1<<m)-1;s>=0;s--)
36     {
37         for(int i=0;i<m;i++)
38             if(!(s&(1<<i))) unidf[s]|=unidf[s|(1<<i)];
39         for(int j=0;j<n;j++) 
40             if(unidf[s]&(1ll<<j)) num[s]+=1;
41         if(num[s]==1) num[s]=0;
42     }
43     f[(1<<m)-1]=0;
44     for(int s=(1<<m)-2;s>=0;s--)
45     {
46         if(!num[s]) {f[s]=0;continue;}
47         int tot=0;
48         for(int i=0;i<m;i++) if(!(s&(1<<i)))tot++;
49         for(int i=0;i<m;i++)
50         {
51             if(((1<<i)&s)) continue;
52             f[s]+=(f[s|(1<<i)]/(dd)tot*(dd)num[s|(1<<i)]/(dd)num[s]); 
53         }
54         f[s]+=1.0;
55     } 
56     printf("%.10lf
",max(f[0],1.00000000000));
57     return 0;
58 }