[USACO Hol10] 臭气弹图上期望概率dp 高斯

记住一开始和后来的经过是两个事件因此概率可以大于一

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<vector>
#include<algorithm>
#define MAXN 333
using namespace std;
typedef double D;
D a[MAXN][MAXN],ans[MAXN];
int p,q,n,m;
vector<int> Link[MAXN];
int in[MAXN];
inline D abs(D x)
{return x<0?0.0-x:x;
}
inline void swap(D &x,D &y)
{D temp=x;x=y;y=temp;
}
void gauss()
{ for(int i=1,k=1;i<=n;i++,k++){int t=i;D h=abs(a[i][k]);for(int j=i+1;j<=n;j++)if(abs(a[j][k])>abs(a[t][k])){t=j;h=abs(a[j][k]);}if(t!=i){for(int j=k;j<=n+1;j++)swap(a[i][j],a[t][k]);}for(int j=i+1;j<=n;j++){h=a[j][k]/a[i][k];for(int l=k;l<=n+1;l++)a[j][l]-=a[i][l]*h;}}for(int i=n;i>0;i--){for(int j=i+1;j<=n;j++)a[i][n+1]-=ans[j]*a[i][j];ans[i]=a[i][n+1]/a[i][i];}
}
int main()
{freopen("dotp.in","r",stdin);freopen("dotp.out","w",stdout);scanf("%d%d%d%d",&n,&m,&p,&q);for(int i=1;i<=m;i++){int x,y;scanf("%d%d",&x,&y);Link[x].push_back(y);in[x]++;Link[y].push_back(x);in[y]++;}for(int i=1;i<=n;i++){a[i][i]=-1.0;for(int j=0;j<Link[i].size();j++)a[i][Link[i][j]]=(D)(1.0-(D)p/q)*(1.0/in[Link[i][j]]);a[i][1+n]=0.0;}a[1][1+n]=-1.0;gauss();for(int i=1;i<=n;i++){ans[i]=ans[i]*p/q;if(ans[i]==-0.0)ans[i]=0.0;printf("%.9lf\n",ans[i]);}return 0;
}

转载于:https://www.cnblogs.com/TSHugh/p/7072765.html

[USACO Hol10] 臭气弹图上期望概率dp 高斯相关推荐

期望/概率dp 学习报告
期望dp 期望的基本性质离散型随机变量和连续型随机变量例题1:红包发红包期望的基本运算性质例题2:Little Pony and Expected Maximum 期望的一些经典题型图上期望 ...
BZOJ 3270: 博物馆 1778: 驱逐猪猡【概率DP+高斯消元】
题目描述: 中文题面,不多解释.1778传送门 3270 传送门 (博物馆)题目分析: 也许很多人做概率题的时候都有种虚幻感..感觉莫名其妙就得出一个期望.概率,一知半解... 所以我在这里仔细地剖析 ...
BZOJ_1778_[Usaco2010 Hol]Dotp 驱逐猪猡_概率DP+高斯消元
BZOJ_1778_[Usaco2010 Hol]Dotp 驱逐猪猡_概率DP+高斯消元题意: 奶牛们建立了一个随机化的臭气炸弹来驱逐猪猡.猪猡的文明包含1到N (2 <= N <= 3 ...
【bzoj1778】[Usaco2010 Hol]Dotp 驱逐猪猡矩阵乘法+概率dp+高斯消元
题目描述奶牛们建立了一个随机化的臭气炸弹来驱逐猪猡.猪猡的文明包含1到N (2 <= N <= 300)一共N个猪城.这些城市由M (1 <= M <= 44,850)条由两 ...
BZOJ3270 博物館概率DP 高斯消元
BZOJ3270 博物館概率DP 高斯消元 @(XSY)[概率DP, 高斯消元] Description 有一天Petya和他的朋友Vasya在进行他们众多旅行中的一次旅行,他们决定去参观一座城堡博 ...
LightOJ - 1038 Race to 1 Again 基础期望概率 dp
传送门刚刚学习期望&概率我们设数X的期望改变次数为P[X] 如果要求X的期望,很容易想到找x的因子; 可以得到下式 ,cnt为X因子个数,ai为X的因子可以这么理解,当因子ai为1时, ...
【BZOJ3640】JC的小苹果概率DP+高斯消元
[BZOJ3640]JC的小苹果 Description 让我们继续JC和DZY的故事. "你是我的小丫小苹果,怎么爱你都不嫌多!" "点亮我生命的火,火火火火火!&qu ...
bzoj 4820: [Sdoi2017]硬币游戏概率dp+高斯消元+KMP
题意周末同学们非常无聊,有人提议,咱们扔硬币玩吧,谁扔的硬币正面次数多谁胜利. 大家纷纷觉得这个游戏非常符合同学们的特色,但只是扔硬币实在是太单调了. 同学们觉得要加强趣味性,所以要找一个同学扔很多 ...
hdu4418 概率dp+高斯消元（近日吐槽）
因为过年,好久没更博客了,最近要狂补题解了== 1.做了几天的概率dp,稍后还要写一篇概率dp和期望的总结,当然,下面几天的专题还是概率dp 2.接着还要整理一个漂亮的高斯消元模板= 3.最近打了几场 ...

[USACO Hol10] 臭气弹图上期望概率dp 高斯

[USACO Hol10] 臭气弹图上期望概率dp 高斯相关推荐

最新文章

热门文章

[USACO Hol10] 臭气弹 图上期望概率dp 高斯

[USACO Hol10] 臭气弹 图上期望概率dp 高斯相关推荐

最新文章

热门文章

[USACO Hol10] 臭气弹图上期望概率dp 高斯

[USACO Hol10] 臭气弹图上期望概率dp 高斯相关推荐