codevs 5966 [SDOI2017]硬币游戏

输入描述 Input Description

输入输出数据精度为1e-10

[题解]

#include<cstdio>
using namespace std;
const int N=1e3+;
char s[N][N];
long double ec[N],f[N][N];
int va[N],vb[N];
int n,m;
void MakeMatrix(){
    ec[m]=;
    for(int i=m-;i;i--) ec[i]=ec[i+]*0.5;
    for(int i=;i<=n;i++){
        va[]=;
        for(int k=,p=;k<=m;k++){
            while(p&&s[i][k]!=s[i][p+]) p=va[p];
            if(s[i][k]==s[i][p+]) p++;
            va[k]=p;
        }
        for(int j=;j<=n;j++){
            vb[]=;
            for(int k=,p=;k<=m;k++){
                while(p&&s[j][k]!=s[i][p+]) p=va[p];
                if(s[j][k]==s[i][p+]) p++;
                vb[k]=p;
            }
            for(int now=vb[m];now;now=va[now]) f[i][j]+=ec[now];
        }
        f[i][n+]=;
    }
}
void Guass(){
    for(int i=;i<=n;i++){
        for(int j=i+;j<=n+;j++) f[i][j]/=f[i][i];
        f[i][i]=;
        for(int j=i+;j<=n;j++){
            for(int k=i+;k<=n+;k++){
                f[j][k]-=f[i][k]*f[j][i];
            }
        }
    }
    for(int i=n;i;i--){
        for(int j=i+;j<=n;j++){
            f[i][n+]-=f[i][j]*f[j][n+];
        }
    }
    long double sum=,ans=;
    for(int i=;i<=n;i++) sum+=f[i][n+];
    for(int i=;i<=n;i++){
        ans=f[i][n+];ans/=sum;
        printf("%.10f\n",(double)ans);
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++) scanf("%s",s[i]+);
    MakeMatrix();
    Guass();
    return ;
}