【Vijos 1966】 夜夜的旅游计划|期望DP|高斯消元|信心题QAQ

做高斯消元忽然想起以前的一个比赛 便拿来水水啦

同时也对期望DP又有啦一个新的认识!!!!抓狂

提醒下最后一个点一直wa的


#include <cmath>
#include <cstdio>
#include <iostream>
#include <algorithm> 
using namespace std;
const double EPS=1e-7;
const int MAXN=200+10;
int n,m,d[MAXN],v[MAXN][MAXN];
double ans[MAXN],f[MAXN][MAXN];
void gauss()
{
	for(int i=1;i<n;i++)
	{
		int t;
		for(int j=i;j<n;j++) 
			if(abs(f[i][j])>EPS) 
			{
				t=i;
				break;	
			}
		if(t!=i) 
			for(int j=i;j<=n;j++) 
				swap(f[i][j],f[t][j]);
		for(int j=i+1;j<n;j++)
		{
			double tmp=f[j][i]/f[i][i];
			for(int k=i;k<=n;k++) 
				f[j][k]-=tmp*f[i][k];
		}
	}
	for(int i=n-1;i>=1;i--)
	{
		ans[i]=f[i][n]/f[i][i];
		for(int j=i-1;j>=1;j--)
			f[j][n]-=f[j][i]*ans[i];
	}
}
int main()
{
	cin >>n >>m;
	for (int i=1;i<=m;i++)
	{
		int x,y,z;
		scanf("%d %d %d",&x,&y,&z);
		d[x]++;d[y]++;
		v[x][y]=v[y][x]=z;
	}
	for(int i=1;i<n;i++)
		for(int j=1;j<n;j++)
		{
			if(i==j) f[i][j]=1*d[i];
			else if(v[i][j]!=0)
			{
				f[i][j]=-1.0;
				f[i][n]+=(double)v[i][j];
			}
		}
	for(int i=1;i<n;i++) 
		if(v[i][n]!=0)
			f[i][n]+=(double)v[i][n];
	gauss();
	printf("%.1lf",ans[1]);
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值