リングの確率はガウス関数によって排除できます
マッチング状況はある文字列から別の文字列に転送される可能性があるため、転送関係のグラフを作成する必要があります
ACオートマトンを構築できますが、ノードの数はnmです。未知数を設定してから、いくつかの方程式をまとめることができます。
Nを誰も勝てない確率とします(プレフィックスは許可されます)
その後、N + A文字列/ B文字列/...n文字列は確実に停止します。
確実に停止する文字列はN + n種類のみです。
若A = HTH B = THT
N + HTH =(... A + TH)+(... A)+(... B + H)
最後にA + B = 1;
次に、方程式を直接解きます。ガウスの消去法では、各項の最大値をその上に配置する必要があることに注意してください。最大値は大きな三角形を形成します。そうしないと、精度が失われます。。
コード:
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
#define eps 0.0000000000001
int n,i,j,l,m,k,sp[505];
long double f[505][505],t,ans[505],er[505];
char ch[505][505];
void gauss()
{
int i,j,k;
for(i=1;i<=n;i++)//枚举方程
{
for(j=i;j<=n;j++)if(fabs(f[j][i])>fabs(f[i][i]))for(k=1;k<=n+1;k++) swap(f[i][k],f[j][k]);
t=f[i][i];
for(j=1;j<=n+1;j++)
f[i][j]/=t;
for(j=1;j<=n;j++)
{
if(j==i)continue;
t=f[j][i];
for(l=1;l<=n+1;l++)
f[j][l]-=f[i][l]*t;
}
}
}
int main()
{
scanf("%d%d",&n,&m);
er[0]=1;
for(i=1;i<=m;i++)
{
er[i]=er[i-1]/2;
}
for(i=1;i<=n;i++)
{
for(j=1;j<=m;j++)
{do{
scanf("%c",&ch[i][j]);
} while(ch[i][j]!='H'&&ch[i][j]!='T');
}
}
for(i=1;i<=n;i++)
{
k=0;
for(j=2;j<=m;j++)//适配
{
while(ch[i][k+1]!=ch[i][j]&&k!=0)k=sp[k];
if(ch[i][k+1]==ch[i][j])k++;
sp[j]=k;
}
for(j=1;j<=n;j++)
{
l=0;
for(k=1;k<=m;k++)
{
while(ch[i][l+1]!=ch[j][k]&&l!=0)l=sp[l];
if(ch[i][l+1]==ch[j][k])l++;
}
while(l)
{
f[i][j]+=er[m-l];
l=sp[l];
}
f[i][n+1]=-er[m];
}
}
for(i=1;i<=n;i++)
f[n+1][i]=1;
f[n+1][n+2]=1;
n++;
gauss();
for(i=1;i<=n; i++)
{
for(j=1;j<=n;j++)
{
if(f[i][j]-eps>0||f[i][j]+eps<0)ans[j]=f[i][n+1];
}
}
for(i=1;i<=n-1;i++)
{
printf("%.10lf\n",double(ans[i]));
}
}