记录编号	469454	评测结果	AAAAAAAAAAAAAAAA
题目名称	[USACO Hol10] 臭气弹	最终得分	100
用户昵称	Shirry	是否通过	通过
代码语言	C++	运行时间	0.133 s
提交时间	2017-11-03 10:12:04	内存使用	1.59 MiB

#include<cstdio>
#include<vector>
#include<cstring>
#include<cmath>
#include<algorithm>
#define eps 1e-13
using namespace std;
const int maxn=410;
int n,m,row,in[maxn];
double p,q,pro,untipro,A[maxn][maxn];
vector<int>s[maxn];
void Gauss(){
	int i,j,k;
	for(i=1;i<=n;i++){
		if(fabs(A[i][i])<=eps){
			for(j=i+1;j<=n;j++){if(A[j][i]>eps)break;}
			for(k=1;k<=m;k++)swap(A[i][k],A[j][k]);
		}
		for(j=1;j<=n;j++){
			if(j!=i&&fabs(A[j][i])>eps){
				double t=A[j][i]/A[i][i];
				for(k=1;k<=m;k++)A[j][k]-=t*A[i][k];
			}
		}
	}
	for(i=1;i<=n;i++)A[i][m]/=A[i][i];
}
int main(){
	freopen("dotp.in","r",stdin);
	freopen("dotp.out","w",stdout);
	scanf("%d%d%lf%lf",&n,&m,&p,&q);
	pro=p/q,untipro=1.0-pro;
	int u,v;
	for(int i=1;i<=m;i++)scanf("%d%d",&u,&v),s[u].push_back(v),s[v].push_back(u),in[u]++,in[v]++;
	m=n+1;
	for(int i=1;i<=n;i++){
		A[i][i]=1.0;
		for(int j=0;j<(int)s[i].size();j++){
			u=s[i][j];
			A[i][u]=-untipro/(double)in[u];
		}
	}
	A[1][m]=1.0;
	Gauss();
	for(int i=1;i<=n;i++)printf("%.9lf\n",A[i][m]*pro);
	return 0;
}