文章目錄

• 題意：
• 思路：

傳送門

題意：

你需要從$1$走到$n$，初始速度是$t$，某些地方有自行車，每個位置自行車有$p_i$的概率是壞掉的，如果自行車沒壞可以騎上自行車，速度是$r$，可以一直騎著到終點。

$1\le t\le r\le 1e4,1\le n,m\le 1e5,0\le k\le 18,1\le a_i\le n,0\le p_i\le 100$

思路：

注意到$k$很小，可以選擇狀壓一下到哪些有自行車的位置，設$f[state][j]$表示當前選擇的自行車位置集合為$state$，最后一次停在$j$的時候，到終點的期望，顯然我們需要倒著推，轉移方程：
$$f[state][i]=min(f[state][i],(f[state|(1<<j)][j]+dis[i][a[j]]/t)?p[i]+dis[i][n]?(1?p[i])/r)$$

讓后選擇記憶化或者循環都可以，這個題由于有邊界問題，顯然選$dfs$更好寫。

記憶化：

#include<bits/stdc++.h> #define X first #define Y second #define L (u<<1) #define R (u<<1|1) #define Mid (tr[u].l+tr[u].r>>1) #define pb push_back using namespace std;const int N=100010,INF=0x3f3f3f3f,mod=1e9+7; typedef long long LL; typedef pair<int,int> PII;int t,r; int n,m,k; vector<PII>v[N]; int dis[21][N]; int a[N],p[N]; bool st[N];void dijkstra(int s) {priority_queue<PII,vector<PII>,greater<PII>>q;memset(st,0,sizeof(st));memset(dis[s],0x3f,sizeof(dis[s]));dis[s][a[s]]=0;q.push({0,a[s]});while(q.size()) {auto u=q.top(); q.pop();int id=u.Y;if(st[id]) continue;st[id]=1;for(auto x:v[id]) {if(dis[s][x.X]>dis[s][id]+x.Y) {dis[s][x.X]=dis[s][id]+x.Y;q.push({dis[s][x.X],x.X});}}} }double f[1<<20][20],P[1010]; LL d[1<<20][20];double dfs(int state,int pos) {if(f[state][pos]!=-1) return f[state][pos];double tmp=P[pos]*dis[pos][n]/t+(1-P[pos])*dis[pos][n]/r;for(int i=0;i<k;i++) {if(i==pos) continue;if(state&(1<<i)) continue;tmp=min(tmp,P[pos]*(dfs(state|(1<<i),i)+1.0*dis[pos][a[i]]/t)+(1-P[pos])*dis[pos][n]/r);}return f[state][pos]=tmp; }void solve() {scanf("%d%d%d%d",&t,&r,&n,&m);while(m--) {int a,b,c; scanf("%d%d%d",&a,&b,&c);v[a].pb({b,c});v[b].pb({a,c});}scanf("%d",&k); a[k]=1; p[k]=100;for(int i=0;i<k;i++) scanf("%d%d",&a[i],&p[i]);for(int i=0;i<=k;i++) dijkstra(i);if(dis[k][n]==INF) {puts("-1");return;}for(int i=0;i<=k;i++) P[i]=1.0*p[i]/100;for(int i=0;i<1<<20;i++) for(int j=0;j<20;j++) f[i][j]=-1;printf("%.8f\n",dfs(0,k)); }int main() {int _=1;while(_--) {solve();}}

循環二進制

#include<bits/stdc++.h> #define X first #define Y second #define L (u<<1) #define R (u<<1|1) #define Mid (tr[u].l+tr[u].r>>1) #define pb push_back using namespace std;const int N=100010,INF=0x3f3f3f3f,mod=1e9+7; typedef long long LL; typedef pair<int,int> PII;int t,r; int n,m,k; vector<PII>v[N]; int dis[21][N]; int a[N],p[N]; bool st[N];void dijkstra(int s) {priority_queue<PII,vector<PII>,greater<PII>>q;memset(st,0,sizeof(st));memset(dis[s],0x3f,sizeof(dis[s]));dis[s][a[s]]=0;q.push({0,a[s]});while(q.size()) {auto u=q.top(); q.pop();int id=u.Y;if(st[id]) continue;st[id]=1;for(auto x:v[id]) {if(dis[s][x.X]>dis[s][id]+x.Y) {dis[s][x.X]=dis[s][id]+x.Y;q.push({dis[s][x.X],x.X});}}} }double f[1<<20][20],P[1010];void solve() {scanf("%d%d%d%d",&t,&r,&n,&m);while(m--) {int a,b,c; scanf("%d%d%d",&a,&b,&c);v[a].pb({b,c});v[b].pb({a,c});}scanf("%d",&k); a[k]=1; p[k]=100;for(int i=0;i<k;i++) scanf("%d%d",&a[i],&p[i]);for(int i=0;i<=k;i++) dijkstra(i);if(dis[k][n]==INF) {puts("-1");return;}for(int i=0;i<=k;i++) P[i]=1.0*p[i]/100;for(int i=0;i<1<<20;i++) for(int j=0;j<20;j++) f[i][j]=1e18;for(int i=(1<<k)-1;i>=1;i--) {for(int x=0;x<k;x++) {//當前走到的點if(i&(1<<x)) {f[i][x]=(1-P[x])*dis[x][n]/r+P[x]*dis[x][n]/t;for(int y=0;y<k;y++) {//下一個要到的點if(x==y) continue;if((i&1<<y)) continue;f[i][x]=min(f[i][x],(1-P[x])*dis[x][n]/r+P[x]*(1.0*dis[x][a[y]]/t+f[i|(1<<y)][y]));}}}}double ans=1.0*dis[k][n]/t;for(int y=0;y<k;y++) {int x=1<<y,i=0;ans=min(ans,1.0*dis[k][a[y]]/t+f[i|(1<<y)][y]);}printf("%.8f\n",ans);}int main() {int _=1;while(_--) {solve();}}

總結

以上是生活随笔為你收集整理的2021 ccpc 哈尔滨 G. Damaged Bicycle 状压 + 期望dp的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。