BZOJ 3143 [HNOI2013] 游走期望dp+高斯消元

最新推荐文章于 2019-09-03 19:46:00 发布

SigmaQuadrant

最新推荐文章于 2019-09-03 19:46:00 发布

阅读量124

点赞数

分类专栏：比赛题解

本文链接：https://blog.csdn.net/qq_40859782/article/details/100116784

版权

比赛题解专栏收录该内容

218 篇文章 0 订阅

订阅专栏

题意给出一个点为 $1 - n$ ，边为 $1 - m$ 的无向连通图，从 $1$ 出发走到 $n$ 停下来，走过一条边的时候要加上这条边的标号，然后要求你给所有的边重新标号，求最小的期望和。
根据排序不等式，我们按照每条边的期望从大到小编从小到大的编号。假设这条边的两个端点分别是 $u, v$ ，并且记从u出发的期望次数为 $E (u)$ ，这个点的度数为 $d e g (u) ，$ 那么这条边被走过的期望为 $\frac{E(u)}{deg(u)}+\frac{E(v)}{deg(v)}$ 。
现在考虑计算每个点出发的期望和，由于是无向图无法直接递推。考虑列式高斯消元。注意 $E (1)$ 要多 $+ 1$ ， $E (n)$ 为 $0$ 。
$G a u s s$ 重新抄了 $r q y$ 的板子，非常好用。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const ll INF=LONG_LONG_MAX;
const int N=5e2+7;

int n,m; 
bool mp[N][N]; 
int u[N*N],v[N*N],deg[N];
double a[N][N];
int Gauss() {
	for(int i=1;i<=n;i++) {
		int j=i;
		for(int k=i+1;k<=n;k++) if(abs(a[k][i])>abs(a[j][i])) j=k;
//		if(abs(a[j][i])<eps) return -1;
		for(int k=i;k<=n+1;k++) swap(a[i][k],a[j][k]);
		for(int j=i+1;j<=n;j++) 
			for(int k=n+1;k>=i;k--)
				a[j][k]-=a[i][k]*a[j][i]/a[i][i];	
	}
	for(int i=n;i>0;i--) {
		a[i][n+1]/=a[i][i];
		for(int j=i-1;j>0;j--)
			a[j][n+1]-=a[i][n+1]*a[j][i];
	}
}
int main() {
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++) {
		scanf("%d%d",&u[i],&v[i]);
		mp[u[i]][v[i]]=mp[v[i]][u[i]]=1;
		++deg[u[i]];
		++deg[v[i]];
	}
	for(int i=1;i<=n;i++) {
		a[i][i]=1.0;
		if(i==n) continue;
		for(int j=1;j<=n;j++) {
			if(j!=i&&mp[i][j]==1) {
				a[i][j]=-1.0/deg[j];
			}
		}
		if(i==1) {
			a[1][n+1]=1.0;
		}
	}
	Gauss();
	vector<double> ans;
	for(int i=1;i<=m;i++) {
		ans.push_back(a[u[i]][n+1]/deg[u[i]]+a[v[i]][n+1]/deg[v[i]]);
	}
	sort(ans.begin(),ans.end());
	double E=0.0;
	for(int i=0;i<ans.size();i++) 
		E+=ans[i]*(m-i); 
	printf("%.3lf\n",E);
	return 0;
}