[SDOI2017]硬币游戏

最新推荐文章于 2023-08-06 20:26:45 发布

Tomorrow......

最新推荐文章于 2023-08-06 20:26:45 发布

阅读量363

点赞数

文章标签：游戏算法 c++

本文链接：https://blog.csdn.net/m0_74821245/article/details/128274119

版权

#include<bits/stdc++.h>

#define N 310

using namespace std;

const unsigned int base=19260817;

int n,m;
char s[N][N];
unsigned int poww[N],sum[N][N];
double div2[N],x[N],a[N][N];

unsigned int get(int i,int l,int r)
{
if(l>r) return 0;
return sum[i][r]-sum[i][l-1]*poww[r-l+1];
}

void Gauss()
{
   for(int i=1;i<=n+1;i++)
   {
for(int j=i;j<=n+1;j++)
       {
if(!a[j][i]) continue;
for(int k=i;k<=n+2;k++)
               swap(a[i][k],a[j][k]);
break;
}
for(int j=i+1;j<=n+2;j++)
       {
if(!a[j][i]) continue;
double tmp=a[j][i]/a[i][i];
for(int k=i;k<=n+2;k++)
               a[j][k]-=a[i][k]*tmp;
}
}
x[n+1]=a[n+1][n+2]/a[n+1][n+1];
for(int i=n;i>=1;i--){
for(int j=i+1;j<=n+1;j++)
           x[i]-=a[i][j]*x[j];
x[i]/=a[i][i];
}
}

int main()
{
   scanf("%d%d",&n,&m);
   poww[0]=div2[0]=1;
   for(int i=1;i<=m;i++)
       poww[i]=poww[i-1]*base,div2[i]=div2[i-1]/2.0;
   for(int i=1;i<=n;i++)
   {
       scanf("%s",s[i]+1);
       for(int j=1;j<=m;j++)
           sum[i][j]=sum[i][j-1]*base+s[i][j];
   }
   for(int i=1;i<=n;i++)
   {
       for(int j=1;j<=n;j++)
       {
           double ans=0;
           for(int k=1;k<=m;k++)
               if(get(i,1,k)==get(j,m-k+1,m))
                   ans+=div2[m-k];
           a[i][j]=ans;
       }
       a[i][n+1]=-div2[m];
   }
   for(int i=1;i<=n;i++)
       a[n+1][i]=1;
   a[n+1][n+2]=1;
   Gauss();
   for(int i=1;i<=n;i++)
       printf("%.10lf\n",x[i]);
   return 0;
}
/*
2 2
TTH
HTT
*/