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

歡迎訪問 生活随笔!

生活随笔

當(dāng)前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

CF605E-Intergalaxy Trips【期望dp】

發(fā)布時間:2023/12/3 编程问答 26 豆豆
生活随笔 收集整理的這篇文章主要介紹了 CF605E-Intergalaxy Trips【期望dp】 小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.

正題

題目鏈接:https://www.luogu.com.cn/problem/CF605E


題目大意

給出nnn個點(diǎn)的一張完全有向圖,每一天iiijjj的路徑有pi,jp_{i,j}pi,j?的概率出現(xiàn)。
詢問從111出發(fā)走到nnn在最優(yōu)策略下的期望天數(shù)。

1≤n≤103,pi,i=11\leq n\leq 10^3,p_{i,i}=11n103,pi,i?=1


解題思路

設(shè)gig_igi?表示從iiinnn的答案,那么ggg滿足式子
gi=∑j=1ngj∏gk<gj(1?pj,k)×pi,j×∏gk<gj(1?pi,k)g_i=\sum_{j=1}^{n}\frac{g_{j}}{\prod_{g_k<g_j}(1-p_{j,k})}\times p_{i,j}\times \prod_{g_k<g_j}(1-p_{i,k})gi?=j=1n?gk?<gj??(1?pj,k?)gj??×pi,j?×gk?<gj??(1?pi,k?)
注意到里面有條件gk<gjg_k<g_jgk?<gj?,但是我們并不知道ggg,好像就死循環(huán)了。

但是我們可以知道gng_ngn?肯定是最小的而且是000,然后剩下第二小的只會受到gng_ngn?的影響,也就是這個里面有一個的ggg是完整的。

可以猜測這個第二小的肯定是剩下的ggg里面最小的那個,因?yàn)槿绻皇?#xff0c;那么顯然不滿足我們最小化的條件,從一個更大的走到了本應(yīng)該可以更小的。

所以我們每次剩下的點(diǎn)中找到一個最小的gj∏gk<gj(1?pj,k)\frac{g_{j}}{\prod_{g_k<g_j}(1-p_{j,k})}gk?<gj??(1?pj,k?)gj??來更新它對其他點(diǎn)的影響即可。

時間復(fù)雜度O(n2)O(n^2)O(n2)


code

#include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int N=1100; int n;bool v[N]; double p[N][N],prod[N],f[N]; int main() {scanf("%d",&n);for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){scanf("%lf",&p[i][j]);p[i][j]/=100.0;}}if(n==1)return puts("0")&0;for(int i=1;i<n;i++)prod[i]=1-p[i][n],f[i]=1;v[n]=1;while(1){double low=1e9;int x;for(int i=1;i<=n;i++)if(f[i]/(1-prod[i])<low&&!v[i])low=f[i]/(1-prod[i]),x=i;if(x==1)return printf("%.10lf",low)&0;v[x]=1;for(int i=1;i<=n;i++){if(v[i])continue;f[i]+=low*p[i][x]*prod[i];prod[i]*=(1-p[i][x]);}}return 0; }

總結(jié)

以上是生活随笔為你收集整理的CF605E-Intergalaxy Trips【期望dp】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。

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