bzoj 4820: [Sdoi2017]硬币游戏【kmp+高斯消元】

有点神，按照1444的做法肯定会挂
注意到它的概率是相同的，所以可以简化状态
详见http://www.cnblogs.com/candy99/p/6701221.html
https://www.geek-share.com/detail/2707064841.html
总之就是靠在kmp中的ne数组上跳来找到对当前串概率有影响的串然后把概率加起来……

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int N=305;
int n,m,ne

;
double a

,p
;
char c

;
void gtne(char s[],int ne[])
{
for(int i=2,j=0;i<=m;i++)
{
while(j&&s[j+1]!=s[i])
j=ne[j];
if(s[j+1]==s[i])
j++;
ne[i]=j;
}
}
void gaosi(int n)
{
for(int i=1;i<=n;i++)
{
int nw=i;
for(int j=i+1;j<=n;j++)
if(fabs(a[j][i])>fabs(a[nw][i]))
nw=j;
for(int j=i;j<=n+1;j++)
swap(a[i][j],a[nw][j]);
for(int j=i+1;j<=n+1;j++)
a[i][j]/=a[i][i];
a[i][i]=1;
for(int j=i+1;j<=n;j++)
{
for(int k=i+1;k<=n+1;k++)
a[j][k]-=a[j][i]*a[i][k];
a[j][i]=0;
}
}
for(int i=n-1;i>=1;i--)
{
for(int j=i+1;j<=n;j++)
a[i][n+1]-=a[j][n+1]*a[i][j];
a[i][n+1]/=a[i][i];
}
}
int main()
{
scanf("%d%d",&n,&m);
p[0]=1;
for(int i=1;i<=m;i++)
p[i]=p[i-1]/2;
for(int i=1;i<=n;i++)
scanf("%s",c[i]+1),a[i][n+1]=-p[m],a[n+1][i]=1;
a[n+1][n+2]=1;
for(int i=1;i<=n;i++)
gtne(c[i],ne[i]);
for(int y=1;y<=n;y++)
for(int x=1;x<=n;x++)
{
int j=0;
for(int i=1;i<=m;++i)
{
while(j&&c[y][j+1]!=c[x][i])
j=ne[y][j];
if(c[y][j+1]==c[x][i])
j++;
}
while(j)
a[y][x]+=p[m-j],j=ne[y][j];
}
gaosi(n+1);
for(int i=1;i<=n;i++)
printf("%.10f\n",a[i][n+2]);
return 0;
}