标签：xk USACO10HOL 期望 int 题解 303 洛谷 x1 高斯消

貌似没有卡我精度？

这道题跟这道的思路和做法都挺像的，也是期望+高斯消元

设fu为一个点期望的经过次数，那么我们可以发现，炸弹在每个点爆炸的概率其实就是fu∗p/q，求出每个点的fi即可得到最终的答案，显然,每个点的期望是由相连的点的期望决定的，dux为点x的度数,点 x 1 x1 x1 , x 2 x2 x2 , x 3 x3 x3 … x k xk xk 与点x相邻，则

∑ i − 1 k \sum_{i-1}^k ∑i−1k​ f x i d u x i \frac{fxi}{duxi} duxifxi​

最后用高斯消元解一下每个点的期望即可，上代码

#include<bits/stdc++.h>
using namespace std;
int n,m,t,x,y;
double p,q,gai,chu,du[303],a[303][303],b[303],ans[303];
vector<int>l[303];
int main()
{
	scanf("%d%d%lf%lf",&n,&m,&p,&q);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&x,&y);
		l[x].push_back(y),l[y].push_back(x),du[x]++,du[y]++;
	}
	gai=p/q,b[1]=-1;
	for(int i=1;i<=n;i++)
	{
		a[i][i]=-1;
		for(int j=0;j<l[i].size();j++)
			a[i][l[i][j]]=(1-gai)*(1/du[l[i][j]]);			
	}
	for(int i=1;i<=n;i++)
	{
		t=i;
		for(int j=i+1;j<=n;j++)
			if(a[j][i]>a[j][t])
				t=j;
		if(t!=i)
		{
			for(int j=i;j<=n;j++)
				swap(a[i][j],a[t][j]);
			swap(b[i],b[t]);
		}
		if(a[i][i]!=0)
			for(int j=i+1;j<=n;j++)
			{
				chu=a[j][i]/a[i][i];
				for(int k=i;k<=n;k++)
					a[j][k]-=a[i][k]*chu;
				b[j]-=b[i]*chu;
			}
	}
	for(int i=n;i>=1;i--)
	{
		ans[i]=b[i]/a[i][i];
		for(int j=1;j<i;j++)
			b[j]-=a[j][i]*ans[i];
	}
	for(int i=1;i<=n;i++)
		printf("%.9lf\n",fabs(ans[i]*gai));
	return 0;
}

标签：xk,USACO10HOL,期望,int,题解,303,洛谷,x1,高斯消
来源： https://blog.csdn.net/guajiProMax/article/details/110713074

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。