【BZOJ3143】【Hnoi2013】 游走 期望DP 高斯消元

卧槽这题BZOJ上的数据用了dcmp就会WA不讲道理。。。裸的才能A。。。

除了第一个点以外,每个点到达的期望次数是与它相连的点的期望次数除以度数之和,一条边的期望次数是他非n端点(因为到n就不再走出来了)的期望次数除以度数之和,

排序后赋值即可。

#include<cstdlib>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
#define maxn 500005
#define oo 999999999
using namespace std;
int n,m;
int g[505][505];
double G[505][505];
struct data 
{
	int x,y;
	double n;
}E[250005];
bool cmp(data x,data y)
{
	return x.n>y.n;
}
void Debug()
{
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n+1;j++)
		{
			printf("%.3lf ",G[i][j]);
		}
		putchar('\n');
	}
	putchar('\n');
	return ;
}
void Init()
{
	scanf("%d%d",&n,&m);
	memset(g,0,sizeof(g));
	int x,y;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&y);
		E[i].x=x; E[i].y=y;
		g[x][y]=g[y][x]=1;
		g[x][0]++; g[y][0]++;
	}

	for(int i=1;i<=n;i++)
	{
		G[i][i]=1.0;if(i==1)G[i][n+1]=1.0;

		for(int j=1;j<n;j++)if(g[i][j])
		{
			G[i][j]=-(1.0)/(double)g[j][0];
		}
	}
	//Debug();
	return ;
}
double Abs(double x)
{
	return x<0 ? -x : x;
}
double f,eps=0.00000000001;
void Gauss()
{
	int ni=1,nj=1,r;
	
	while(ni<=n && nj<= n)
	{
		r=ni;
		for(int i=ni+1;i<=n;i++)
		{
			if(Abs(G[i][nj])>Abs(G[r][nj]))
			{
				r=i;
			}
		}
		if(Abs(G[r][nj])>=eps)
		{
			for(int i=1;i<=n+1;i++)
			{
				swap(G[r][i],G[ni][i]);
			}
			for(int i=1;i<=n;i++)if(i!=ni)
			{
				f=G[i][nj]/G[ni][nj];
				for(int j=1;j<=n+1;j++)
				{
					G[i][j]-=G[ni][j]*f;
				}
			}
			ni++;
		}
		nj++;
	}
	
	for(int i=1;i<=n;i++)
	{
		G[i][n+1]=G[i][n+1]/G[i][i];
		G[i][i]=1.0;
	}
	//Debug();
	return ;
}
int main()
{
	//freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	Init();
	Gauss();
	int x,y;
	for(int i=1;i<=m;i++)
	{
		x=E[i].x; y=E[i].y;
		E[i].n=0.0;
		if(x!=n && g[x][0]!=0)
		{
			E[i].n+=G[x][n+1]/(double)g[x][0];
		}
		if(y!=n && g[y][0]!=0)
		{
			E[i].n+=G[y][n+1]/(double)g[y][0];
		}
	}
	sort(E+1,E+1+m,cmp);
	double ans=0;
	for(int i=1;i<=m;i++)
	{
		ans+=E[i].n*(double)(i);
	}
	printf("%.3lf\n",ans+eps);
	return 0;
}


 

 

你可能感兴趣的:(【BZOJ3143】【Hnoi2013】 游走 期望DP 高斯消元)