首先一个结论是任意一个长度为 x x x的串的概率为 1 2 x \frac{1}{2^x} 2x1
考虑设所有不是终止答案的状态为
S
S
S
f
i
f_i
fi表示
i
i
i赢的概率
那么可以直接填特定
n
n
n个得到
i
i
i
即
f
i
←
1
2
l
e
n
S
f_i\leftarrow\frac{1}{2^{len}}S
fi←2len1S
当然有可能
S
S
S结尾可以直接拼出一段
或者拼的过程中先拼出其他的一个串
j
j
j
可以发现这种情况
只有当
j
j
j的一个后缀等于
i
i
i的一个前缀
那么
f
i
=
1
2
m
−
∑
j
=
1
n
∑
k
=
1
m
[
p
r
e
[
i
]
[
k
]
=
s
u
f
[
j
]
[
k
]
]
1
2
k
f
j
f_i=\frac{1}{2^m}-\sum_{j=1}^{n}\sum_{k=1}^{m}[pre[i][k]=suf[j][k]]\frac{1}{2^k}f_j
fi=2m1−∑j=1n∑k=1m[pre[i][k]=suf[j][k]]2k1fj
O ( n 3 ) O(n^3) O(n3)高斯消元即可
不过精度十分爆炸
高斯消元的时候把每列最大的拿来消
#include<bits/stdc++.h>
using namespace std;
const int RLEN=1<<20|1;
inline char gc(){
static char ibuf[RLEN],*ib,*ob;
(ob==ib)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
return (ob==ib)?EOF:*ib++;
}
#define gc getchar
inline int read(){
char ch=gc();
int res=0,f=1;
while(!isdigit(ch))f^=ch=='-',ch=gc();
while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
return f?res:-res;
}
#define ll long long
#define re register
#define pii pair<int,int>
#define fi first
#define se second
#define pb push_back
#define cs const
#define bg begin
#define poly vector<int>
#define chemx(a,b) ((a)<(b)?(a)=(b):0)
#define chemn(a,b) ((a)>(b)?(a)=(b):0)
cs int N=305;
cs double eps=1e-10;
cs int bas1=1331,mod1=1291245795,bas2=233,mod2=719831223;
double f[N][N],pw[N];
int n,m;
char s[N];
int pre1[N][N],pt1[N],pre2[N][N],pt2[N];
inline bool comp(double x){
return x>eps||x<-eps;
}
inline void Gauss(int n){
for(int i=1;i<=n;i++){
int pos=i;
for(int j=i;j<=n;j++)if(comp(f[pos][i]-f[j][i]))pos=j;
if(pos!=i)swap(f[pos],f[i]);
for(int j=i+1;j<=n;j++){
double p=f[j][i]/f[i][i];
for(int k=i;k<=n+1;k++)f[j][k]-=f[i][k]*p;
}
}
for(int i=n;i;i--){
for(int j=i+1;j<=n;j++)
f[i][n+1]-=f[i][j]*f[j][n+1];
f[i][n+1]/=f[i][i];
}
}
int main(){
n=read(),m=read();
pt1[0]=pt2[0]=1;
for(int i=1;i<=m;i++)pt1[i]=1ll*pt1[i-1]*bas1%mod1,pt2[i]=1ll*pt2[i-1]*bas2%mod2;
for(int i=1;i<=n;i++){
scanf("%s",s+1);
for(int j=1;j<=m;j++)
pre1[i][j]=(1ll*pre1[i][j-1]*bas1+(s[j]-'A'))%mod1,pre2[i][j]=(1ll*pre2[i][j-1]*bas2+(s[j]-'A'))%mod2;
}
pw[0]=1;
for(int i=1;i<=m;i++)pw[i]=pw[i-1]*0.5;
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++)
for(int k=1;k<=m;k++)
if((pre1[i][k]==((pre1[j][m]-1ll*pre1[j][m-k]*pt1[k]%mod1)+mod1)%mod1)&&(pre2[i][k]==((pre2[j][m]-1ll*pre2[j][m-k]*pt2[k]%mod2)+mod2)%mod2))
f[i][j]+=pw[m-k];
f[i][n+1]=-pw[m];
}
for(int i=1;i<=n;i++)f[n+1][i]=1;
f[n+1][n+2]=1;
Gauss(n+1);
for(int i=1;i<=n;i++)printf("%.8lf\n",f[i][n+2]);
}