BZOJ2306/Ctsc2011 幸福路径

最新推荐文章于 2020-06-06 11:43:00 发布

sherlock_zhuang

最新推荐文章于 2020-06-06 11:43:00 发布

阅读量775

点赞数

分类专栏： my-code 文章标签：倍增 Floyd

本文链接：https://blog.csdn.net/sherlock_zhuang/article/details/50811286

版权

my-code 专栏收录该内容

28 篇文章 0 订阅

订阅专栏

思路：
还是比较巧妙的一种求期望的方法，以前没遇到过。其实只要Floyd倍增求即可逼近近似值。具体看其他人的题解。

//采用floyd多次可以得到近似值 
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
using namespace std;
const int imax=229;
const double inf=1000000229.229;
int n,m,S; double val[imax],p;
double f[62][imax][imax];


void iread()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%lf",&val[i]);
    scanf("%d",&S); scanf("%lf",&p);


    for(int k=0;k<=60;k++)
        for(int i=0;i<=n;i++) 
            for(int j=0;j<=n;j++) f[k][i][j]=-inf;
    for(int i=0;i<=n;i++) f[0][i][i]=0;//有可能一开始不动就比较大    

    int u,v;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&u,&v);
        f[0][u][v]=val[v]*p; //不包括起点 不然会算重复; 
    }
}

void iwork()
{
    for(int s=1;s<=60;s++)
    {
        for(int k=1;k<=n;k++)
            for(int i=1;i<=n;i++)
                for(int j=1;j<=n;j++)
                    f[s][i][j]=max(f[s][i][j],f[s-1][i][k]+f[s-1][k][j]*p);
        p*=p;   
    }
    double ans=-inf;
    for(int i=1;i<=n;i++) ans=max(ans,f[60][S][i]);
    printf("%.1lf\n",ans+val[S]);   
}

int main()
{
//  freopen("p3.in","r",stdin);
    iread();
    iwork();
    return 0;
}