BZOJ2306-[Ctsc2011]幸福路径

最新推荐文章于 2019-02-11 20:39:00 发布

JackflyDC

最新推荐文章于 2019-02-11 20:39:00 发布

阅读量272

点赞数 2

分类专栏： OI

本文链接：https://blog.csdn.net/qq_34531807/article/details/81211693

版权

OI 专栏收录该内容

222 篇文章 0 订阅

订阅专栏

题解：
倍增 $Floyd$ 。
令 $f[t][i][j]$ 表示走 $2^t$ 步，从 $i$ 走到 $j$ 获得的最大幸福度。
$f[t][i][j]=max(f[t-1][i][k]+f[t-1][k][j]*p^{2^{t-1}})$
当 $t$ 足够大时， $f$ 得到的就近似为最大的幸福值。
$Code：$

#include<bits/stdc++.h>
#define N 229
#define inf 1000000000
using namespace std;
int n,m,s,u,v;
double a[N],p,f[62][N][N];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%lf",&a[i]);
    scanf("%d%lf",&s,&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;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&u,&v);
        f[0][u][v]=a[v]*p;
    }
    for(int t=1;t<=60;t++)
    {
        for(int k=1;k<=n;k++)
            for(int i=1;i<=n;i++)
                for(int j=1;j<=n;j++)
                    f[t][i][j]=max(f[t][i][j],f[t-1][i][k]+f[t-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+a[s]);
    return 0;
}