bzoj 2306: [Ctsc2011]幸福路径

题意:在一个有向图中,一只蚂蚁在上面爬,每爬过一条边体力就*ρ,经过一个点获得点权*体力的收益,求最大收益。
题解:倍增floyd。因为蚂蚁可以无限走下去,直接求貌似不现实,但是体力下降得非常快,所以走足够多次就可以进入答案的误差允许范围内了。
f[i][j][l]表示从i点走到j点走2^l步的收益。
f[i][j][l]=max{f[i][k][l1]+f[k][j][l1]ρl}
滚一下就可以了。
代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

int n,m,s;
double f[110][110],g[110][110],p,a[110];

int main()
{
    memset(f,0xc2,sizeof(f));
//  printf("%lf",f[0][0][0]);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    scanf("%lf",&a[i]);
    scanf("%d%lf",&s,&p);
    for(int i=1;i<=n;i++)
    f[i][i]=0;
    for(int i=0;i<m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        f[x][y]=a[y];
    }
    double hh=p;
    for(int l=0;l<=70;l++)
    {
        memset(g,0xc2,sizeof(g));
        for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        g[i][j]=max(g[i][j],f[i][k]+f[k][j]*hh);
        hh*=hh;
        memcpy(f,g,sizeof(f));
    }
    double ans=0;
    for(int i=1;i<=n;i++)
    ans=max(ans,f[s][i]);
    printf("%.1lf",ans*p+a[s]);
}
阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/herobrine_tkj/article/details/77922930
个人分类: bzoj 倍增 Floyd
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭
关闭