BZOJ3270 博物馆

建一张新图,有n^2个点,新图中一个点代表发生第一个男孩在点x,第二个男孩在点y这种情况的期望次数,两个人都在一个点的情况就是终点

然后就是高斯消元了

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<map>
#include<set>
#include<bitset>
#include<queue>
#include<stack>
using namespace std;
#define MAXN 30
#define MAXM 410
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
#define ll long long
bool mp[MAXN][MAXN];
double MP[MAXM][MAXM];
int d[MAXN];
int n,m,s1,s2,N;
double p[MAXN];
double a[MAXM][MAXM];
int pt(int x,int y){
	return (x-1)*n+y;
}
void gs(){
	int i,j,k;
	for(i=1;i<=N;i++){
		if(fabs(a[i][i])<eps){
			for(j=i+1;j<=N;j++){
				if(fabs(a[j][i])>eps){
					for(k=1;k<=N+1;k++){
						swap(a[i][k],a[j][k]);
					}
					break;
				}
			}
		}
		if(fabs(a[i][i])>eps){
			for(j=1;j<=N;j++){
				if(i!=j){
					double t=a[j][i]/a[i][i];
					for(k=1;k<=N+1;k++){
						a[j][k]-=a[i][k]*t;
					}
				}
			}
		}
	}
}
int main(){
	int i,j,k,l,x,y;
	scanf("%d%d%d%d",&n,&m,&s1,&s2);
	N=n*n;
	for(i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		mp[x][y]=mp[y][x]=1;
		d[x]++;
		d[y]++;
	}
	for(i=1;i<=n;i++){
		scanf("%lf",&p[i]);
	}
	for(i=1;i<=n;i++){
		for(j=1;j<=n;j++){
			if(i==j){
				continue ;
			}
			for(k=1;k<=n;k++){
				if(!mp[i][k]&&i!=k){
					continue ;
				}
				for(l=1;l<=n;l++){
					if(!mp[j][l]&&j!=l){
						continue ;
					}
					double t=1;
					if(k==i){
						t*=p[i];
					}else{
						t*=(1.0-p[i])/d[i];
					}
					if(l==j){
						t*=p[j];
					}else{
						t*=(1.0-p[j])/d[j];
					}
					MP[pt(i,j)][pt(k,l)]=t;
				}
			}
		}
	}
	for(i=1;i<=N;i++){
		a[i][i]=1;
		for(j=1;j<=N;j++){
			a[i][j]-=MP[j][i];
		}
	}
	a[pt(s1,s2)][N+1]=1;
	gs();
	for(i=1;i<=n;i++){
		printf("%.6lf ",a[pt(i,i)][N+1]/a[pt(i,i)][pt(i,i)]);
	}
	printf("\n");
	return 0;
}

/*

*/


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值