【SDOI2017】硬币游戏（高斯消元）（AC自动机）

传送门

---

因为不取模，所以高消的时候可能会在一些奇奇怪怪的地方炸精度。

都2017的省选了，在这种题目上为什么不取模啊

对计数和概率取模这种套路不是在2014年的时候就在各种OJ的比赛里面普及过了吗。。。

这道题的严格推导要用到无限集合的理论，我自己尝试证明了一遍，有点复杂，而且对于不懂高等概率的人来说可能看不懂，这里只给一个口胡式的解答。

题解：

首先我们明白一点，在接下来的$t$次随机中强行指定一个长度为$t$的串出现的概率为$2^t$，这里的出现是指接下来$t$次立即出现。

设答案为$p_i,i=1,2\cdots n$

假设在进行了足够多次随机之后（中间就算出现了某一个串也不停下）我们强行停止这个过程（就算一个串都没有出现也停止），设$x_i(i=1,2,\cdots n)$表示在该串中第一个出现的串为$i$的概率，一个串都没有出现的概率为$x_0$。

注意，在进行了足够多次乃至无限次随机之后，有$x_1:x_2:\cdots :x_n=p_1:p_2:\cdots :p_n$。没好好学过微积分请感性理解，学过微积分请用极限理解。

所以我们求出$x_i=k{x}_0$之后利用${\sum }_{i=1}^n{x}_i=1$来解方程就行了。

现在请忘掉上面那个关于足够多次的构造。

设$N$表示还没有结束的局面集合，我们考虑在$N$后面加一个$s_i$来让游戏强行结束。

但是注意，这时候获胜的不一定是$s_i$。。

如果$s_i$的前缀和另一个串$s_j$的后缀相同，那么就可能会出事。。。

出事的概率是所有$i$的前缀等于$j$的后缀的长度的$\frac{1}{2}$的幂次之和。

对这个玩意列方程求解即可。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

using std::cerr;
using std::cout;

cs int N=3e2+7,M=1e5+7;
cs double eps=1e-10,P=.5;

int n,m;
int son[M][2],cnt,fail[M],d[M],ps[M];
std::vector<int> nd[M];
double p[M],ans[M],A[N][N];

inline void ins(int id){
	static char s[N];
	scanf("%s",s+1);
	int cur=0;
	for(int re i=1;i<=m;++i){
		int c=s[i]=='H';
		if(!son[cur][c]){son[cur][c]=++cnt;d[cnt]=d[cur]+1;}
		nd[cur=son[cur][c]].push_back(id);
	}ps[id]=cur;
}

inline void get_fail(){
	std::queue<int> q;
	if(son[0][0])q.push(son[0][0]);
	if(son[0][1])q.push(son[0][1]);
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int re i=0;i<2;++i)
		if(son[u][i])fail[son[u][i]]=son[fail[u]][i],q.push(son[u][i]);
		else son[u][i]=son[fail[u]][i];
	}
}

inline void get(int x){
	for(int re u=ps[x];u;u=fail[u])
	for(int re v:nd[u])A[v][x]+=p[m-d[u]];
}

inline void Gauss(int n){
	for(int re i=1;i<=n;++i){
		int p;for(p=i;p<=n;++p)if(fabs(A[p][i])>eps)break;
		if(p>n)p=i;
		if(p!=i)for(int re j=i;j<=n+1;++j)std::swap(A[p][j],A[i][j]);
		double t=A[i][i];
		for(int re j=i;j<=n+1;++j)A[i][j]/=t;
		for(int re j=i+1;j<=n;++j)if(fabs(A[j][i])>eps){
			double t=A[j][i];
			for(int re k=i;k<=n+1;++k)A[j][k]-=t*A[i][k];
		}
	}
	for(int re i=n;i;--i){
		ans[i]=A[i][n+1];
		for(int re j=i+1;j<=n;++j)ans[i]-=A[i][j]*ans[j];
	}
}

signed main(){
#ifdef zxyoi
	freopen("coin.in","r",stdin);
#endif
	scanf("%d%d",&n,&m);
	p[0]=1;for(int re i=1;i<=m;++i)p[i]=p[i-1]*P;
	for(int re i=1;i<=n;++i)ins(i);get_fail();
	for(int re i=1;i<=n;++i)get(i);
	for(int re i=1;i<=n;++i)A[i][n+1]=-p[m],A[n+1][i]=1;
	A[n+1][n+2]=1;Gauss(n+1);
	for(int re i=1;i<=n;++i)printf("%.10f\n",ans[i]);
	return 0;
}