日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

bzoj3143,P3232-[Hnoi2013]游走【数学期望,高斯消元,贪心】

發(fā)布時間:2023/12/3 编程问答 28 豆豆
生活随笔 收集整理的這篇文章主要介紹了 bzoj3143,P3232-[Hnoi2013]游走【数学期望,高斯消元,贪心】 小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.

正題

題目鏈接:
https://www.lydsy.com/JudgeOnline/problem.php?id=3143
https://www.luogu.org/problem/P3232


題目大意

一張無向圖nnn個點mmm條邊,然后給每條邊附上1~m1\sim m1m的權值(不能重復),求1走到nnn的最小期望值。


解題思路

我們可以計算出每條邊經(jīng)過的期望次數(shù),在此之前我們要先計算出每個點經(jīng)過的期望次數(shù)。
我們發(fā)現(xiàn)很容易列出方程fx=∑y?>xfyinxf_x=\sum_{y->x} \frac{f_y}{in_x}fx?=y?>x?inx?fy??
但是我們發(fā)現(xiàn)該方程有后效性所以我們用高斯消元求解。

然后對于一條邊(x,y)(x,y)(x,y)經(jīng)過的期望次數(shù)就是
fxinx+fyiny\frac{f_x}{in_x}+\frac{f_y}{in_y}inx?fx??+iny?fy??
然后貪心選取邊權即可。


codecodecode

#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> using namespace std; const int N=510; int n,m,x[N*N],y[N*N],in[N]; double a[N][N],q[N*N],ans; int main() {scanf("%d%d",&n,&m);for(int i=1;i<=m;i++){scanf("%d%d",&x[i],&y[i]);a[x[i]][y[i]]=a[y[i]][x[i]]=1;in[y[i]]++;in[x[i]]++;}memset(a[n],0,sizeof(a[n]));for(int i=1;i<n;i++){for(int j=1;j<n;j++)a[i][j]=-a[i][j]/in[j];a[i][n]=0;a[i][i]=1;}a[1][n]=1;for(int i=1;i<n;i++){int z=i;for(int j=i+1;j<n;j++)if(fabs(a[j][i])>fabs(a[z][i])) z=j;for(int j=1;j<=n;j++)swap(a[i][j],a[z][j]);for(int j=n;j>=i;j--)a[i][j]/=a[i][i];for(int j=1;j<n;j++){if(i==j) continue;double rate=a[j][i]/a[i][i];for(int k=i;k<=n;k++)a[j][k]-=rate*a[i][k];}}for(int i=1;i<=m;i++){if(x[i]!=n)q[i]+=a[x[i]][n]*(1.0/in[x[i]]);if(y[i]!=n)q[i]+=a[y[i]][n]*(1.0/in[y[i]]);}sort(q+1,q+1+m);for(int i=1;i<=m;i++)ans+=q[i]*((m-i+1)*1.0);printf("%.3lf",ans); }

總結

以上是生活随笔為你收集整理的bzoj3143,P3232-[Hnoi2013]游走【数学期望,高斯消元,贪心】的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內容還不錯,歡迎將生活随笔推薦給好友。