UESTC - 1652 递推方程

最新推荐文章于 2018-11-20 19:17:40 发布

转载最新推荐文章于 2018-11-20 19:17:40 发布 · 55 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/caturra/p/8292113.html

本文介绍了一种使用概率动态规划解决特定路径选择问题的方法。通过定义状态转移方程，解决了一个涉及多条路径选择的概率计算问题。文章通过具体代码实现了这一算法，并详细解释了各个部分的功能。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

方程很简单，每一公里往上推就行
WA了2发，忘了单通道时的特判，还有n m傻傻分不清，忘了fixed什么的我好弱啊QAQ..

#include<bits/stdc++.h>
#define rep(i,j,k) for(int i=j;i<=k;i++)
using namespace std;
double dp[2][30010];
vector<int> vec[1003];
int n,m,p,k,a,b;
int main(){
    ios::sync_with_stdio(0);
    while(cin>>m>>k>>n>>p){
        memset(dp,0,sizeof dp);dp[0][p]=1;
        memset(vec,0,sizeof vec);
        rep(i,1,k){
            cin>>a>>b;
            vec[b].push_back(a);
        }
        rep(i,1,n) sort(vec[i].begin(),vec[i].end());
        vector<int>::iterator it;
        if(m>=2) rep(i,0,n-1){
            memset(dp[i+1&1],0,sizeof dp[i+1&1]);
            rep(j,1,m) if( ( it = find(vec[i+1].begin(),vec[i+1].end(),j) ) !=vec[i+1].end() ) dp[i+1&1][j]=-6666666;
            rep(j,1,m){
                if(dp[i&1][j]<=0) continue;
                if(j==1){
                    if(dp[i+1&1][j+1]>=0) dp[i+1&1][j+1]+=1.0/2.0*dp[i&1][j];
                    if(dp[i+1&1][j]>=0) dp[i+1&1][j]+=1.0/2.0*dp[i&1][j];
                }
                else if(j==m){
                    if(dp[i+1&1][j-1]>=0) dp[i+1&1][j-1]+=1.0/2.0*dp[i&1][j];
                    if(dp[i+1&1][j]>=0) dp[i+1&1][j]+=1.0/2.0*dp[i&1][j];
                }
                else{
                    if(dp[i+1&1][j]>=0) dp[i+1&1][j]+=1.0/3.0*dp[i&1][j];
                    if(dp[i+1&1][j-1]>=0) dp[i+1&1][j-1]+=1.0/3.0*dp[i&1][j];
                    if(dp[i+1&1][j+1]>=0) dp[i+1&1][j+1]+=1.0/3.0*dp[i&1][j];
                }
            }
        }
        else{
            bool flag=0;
            rep(i,1,n)  if((it=find(vec[i].begin(),vec[i].end(),1))!=vec[i].end()) flag=1;
            if(flag==1) cout<<"0.000000"<<endl;
            else cout<<"1.000000"<<endl;
            continue;
        }
        double ans=0;
        rep(j,1,m) if( ( it = find(vec[n].begin(),vec[n].end(),j) ) !=vec[n].end() ) dp[n&1][j]=-1;
        rep(i,1,m) ans+=(dp[n&1][i]<=0?0:dp[n&1][i]);
        cout<<fixed<<setprecision(6)<<ans<<endl;
    }
    return 0;
}

转载于:https://www.cnblogs.com/caturra/p/8292113.html