[Luogu P2973&BZOJ 1778][USACO10HOL]赶小猪DOtP(高斯消元+期望)
Posted Zars19
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[Luogu P2973&BZOJ 1778][USACO10HOL]赶小猪DOtP(高斯消元+期望)相关的知识,希望对你有一定的参考价值。
Description
一个无向图,节点1有一个炸弹,在每个单位时间内,有可能在这个节点炸掉,也有p/q的概率随机选择一条出去的路到其他的节点上。问最终炸弹在每个节点上爆炸的概率。
Solution
炸弹最终爆炸的概率是1(可能是在无穷无尽的时间之后…),于是到达某一点的概率除以到达所有点的概率就是在这个节点上爆炸的概率
f[u]=∑f[v]*(1-p/q)/d[v]
#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> using namespace std; int n,m,d[303],head[303],cnt=0; double a[303][303],f[303],p,q; struct Node { int next,to; }Edges[100005]; void addedge(int u,int v) { Edges[++cnt].next=head[u]; head[u]=cnt; Edges[cnt].to=v; } void Gauss() { for(int i=1;i<=n;i++) { int maxline=i; for(int j=i;j<=n;j++) if(a[j][i]>a[maxline][i])maxline=j; if(maxline!=i) for(int j=i;j<=n+1;j++) swap(a[maxline][j],a[i][j]); for(int j=i+1;j<=n;j++) { double t=a[j][i]/a[i][i]; for(int k=i;k<=n+1;k++) a[j][k]-=t*a[i][k]; } } for(int i=n;i>0;i--) { for(int j=n;j>i;j--) a[i][n+1]-=f[j]*a[i][j]; f[i]=a[i][n+1]/a[i][i]; } } int main() { memset(head,-1,sizeof(head)); scanf("%d%d%lf%lf",&n,&m,&p,&q); for(int i=1;i<=m;i++) { int u,v; scanf("%d%d",&u,&v); addedge(u,v); addedge(v,u); d[u]++,d[v]++; } for(int i=1;i<=n;i++) { a[i][i]+=1; for(int j=head[i];~j;j=Edges[j].next) { int v=Edges[j].to; a[v][i]-=(1-p/q)/d[i]; } } a[1][n+1]+=1-p/q; Gauss(); double sum=0; for(int i=1;i<=n;i++)sum+=f[i]; for(int i=1;i<=n;i++) printf("%.9lf\\n",f[i]/sum); return 0; }
以上是关于[Luogu P2973&BZOJ 1778][USACO10HOL]赶小猪DOtP(高斯消元+期望)的主要内容,如果未能解决你的问题,请参考以下文章
Luogu2973 [USACO10HOL]Driving Out the Piggies G
BZOJ1229 & 洛谷2917:[USACO2008 NOV]toy 玩具 & 洛谷4480:[BJWC2018]餐巾计划问题——题解