BZOJ 4820 [Sdoi2017]硬币游戏 ——期望DP 高斯消元

做法太神了,理解不了。

自己想到的是建出AC自动机然后建出矩阵然后求逆计算,感觉可以过$40%$

用一个状态$N$表示任意一个位置没有匹配成功的概率和。

每种匹配不成功的情况都是等价的。

然后我们强制在后面加上长度为m的01串,那么这个串的概率是一定的。

然后考虑加上的这些字符还能拼成什么串,因为状态$N$的末尾是不确定的。

如果另外一个串的后缀等于这个串的前缀的话,是可能带来影响的。

所以计算出影响的概率,然后高斯消元即可。

然而有一个问题,N的概率最后消出来代表什么意思啊,是指期望的长度吗?

希望各位dalao不吝赐教。

#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
#define ll long long
#define mp make_pair int n,m,s[305][305],str[605],fail[605];
char ss[305];
double a[305][305],pw[305],ans[305]; void kmp()
{
// F(i,1,2*m) printf("%d ",str[i]); printf("\n");
str[0]=-1;
memset(fail,0,sizeof fail);
for (int i=2,j=0;i<=2*m;++i)
{
while (j&&str[i]!=str[j+1]) j=fail[j];
if (str[j+1]==str[i]) j++;
fail[i]=j;
}
// F(i,1,2*m) printf("%d ",fail[i]);printf("\n");
} void solve(int x)
{
a[x][x]=1;F(i,1,m) str[i]=s[x][i];
F(y,1,n)// if (y!=x)
{
F(i,1,m) str[i+m]=s[y][i]; kmp();
int now=fail[2*m];
// printf("now is %d\n",now);
while (now>=m) now=fail[now];
while (now)
{
// printf("Can %d\n",now);
a[x][y]+=pw[m-now];
now=fail[now];
}
}
a[x][n+1]=-pw[m]; a[x][n+2]=0;
} void Gauss()
{
F(i,1,n+1)
{
int tmp=i;
F(j,i+1,n+1) if (fabs(a[j][i])>fabs(a[i][i])) tmp=j;
if (tmp!=i) F(j,1,n+2) swap(a[i][j],a[tmp][j]);
F(j,1,n+1) if (j!=i)
{
double t=a[j][i]/a[i][i];
F(k,1,n+2) a[j][k]-=t*a[i][k];
}
// F(i,1,n+1){F(j,1,n+2)printf("%.3f ",a[i][j]);printf("\n");}
}
F(i,1,n+1) ans[i]=a[i][n+2]/a[i][i];
F(i,1,n) printf("%.10lf\n",ans[i]);
} int main()
{
// freopen("in.txt","r",stdin);
scanf("%d%d",&n,&m);
pw[0]=1;F(i,1,m)pw[i]=pw[i-1]*0.5;
F(i,1,n){scanf("%s",ss+1);F(j,1,m) s[i][j]=(ss[j]=='H');}
F(i,1,n) solve(i);
F(i,1,n) a[n+1][i]=1; a[n+1][n+2]=1;
// F(i,1,n+1){F(j,1,n+2)printf("%.3f ",a[i][j]);printf("\n");}
Gauss();
}

  

上一篇:[Android][Framework] 添加系统服务


下一篇:初识Avro