【Ctsc2011】幸福路径

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2306

给定一张有向图，每个点有权值，蚂蚁从某个节点出发，初始体力值为$1$，每走一条边$体力值*=p$，每经过一个点会获得幸福值为$点权*体力值$，求最大幸福值。即求

${\sum _{i=0}^{\infty }w[i]*p^{i}}$且${U(i-1,i)=1}$其中${U(A,B)}$表示是否存在A到B这样一条路径。

正解是一个叫做倍增Floyed的东西。

其实就是说考虑到步数是可以走无限步的，我们只需要知道一个近似值，而精度要求也还是比较高的，考虑用倍增的方法快速确定精度。
令${f[i][j][k]}$表示从第$i$个点走到第$j$个点，期间走了${2^{k}}$步。

那么可以利用Floyed进行如下转移：
$${f[i][j][k]=max(\forall (f[i][t][k-1]+f[t][j][k-1]*p^{2^{k}}))}$$

注意$k$这一维是可以滚动的。

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<vector>
 5 #include<cstdlib>
 6 #include<cmath>
 7 #include<cstring>
 8 using namespace std;
 9 #define maxn 110
10 #define inf 0x7fffffff
11 #define llg long long
12 #define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
13 llg n,m,S;
14 double p,ans,a[maxn],f[maxn][maxn],g[maxn][maxn],tmp;
15 int main()
16 {
17     yyj("a");
18     cin>>n>>m;
19     for (llg i=1;i<=n;i++) scanf("%lf",&a[i]);
20     cin>>S>>p;
21     for (llg i=1;i<=n;i++)
22         for (llg j=1;j<=n;j++)
23             if (i!=j) f[i][j]=inf*-1;
24     for (llg i=1;i<=m;i++)
25     {
26         llg x,y;
27         scanf("%lld%lld",&x,&y);
28         f[x][y]=a[y];
29     }
30     llg T=50;
31     double tmp=p;
32     while (T--)
33     {
34         //tmp*=tmp;
35         for (llg i=1;i<=n;i++)
36             for (llg j=1;j<=n;j++)
37                 if (i!=j) g[i][j]=inf*-1;
38         for (llg k=1;k<=n;k++)
39             for (llg i=1;i<=n;i++)
40                 for (llg j=1;j<=n;j++)
41                     g[i][j]=max(g[i][j],f[i][k]+f[k][j]*tmp);
42         memcpy(f,g,sizeof f);
43         tmp*=tmp;
44     }
45     for (llg i=1;i<=n;i++) ans=max(ans,f[S][i]);
46     printf("%.1lf\n",ans*p+a[S]);
47     return 0;
48 }

转载于:https://www.cnblogs.com/Dragon-Light/p/6391345.html

【Ctsc2011】幸福路径相关推荐

【bzoj2360】【cstc2011】【幸福路径】【倍增+floyd】
Description 有向图 G有n个顶点 1, 2, -, n,点i 的权值为 w(i).现在有一只蚂蚁,从给定的起点 v0出发,沿着图 G 的边爬行.开始时,它的体力为 1.每爬过一条边 ...
【传递闭包】【倍增】幸福路径（P4308）
正题 P4308 题目大意给一个图,每个点有一个权值 axa_xax,当你经过一条边时体力×p\times p×p(体力初始为1),每到一个点造成的贡献为体力×ax\times a_x×ax,问 ...
Noip前的大抱佛脚----赛前任务
赛前任务 tags:任务清单前言现在xzy太弱了,而且他最近越来越弱了,天天被爆踩,天天被爆踩题单不会在作业部落发布,所以可(yi)能(ding)会不及时更新省选前的练习莫名其妙地成为了Noi ...
2019.4.summary
2019.4.1 BZOJ1061: [Noi2008]志愿者招募真心有点难QAQ https://www.byvoid.com/zhs/blog/noi-2008-employee 看void爷的 ...
rest spring_带有Spring的REST的ETag
rest spring 1.概述本文将重点介绍ETag – Spring支持,RESTful API的集成测试以及带有curl的使用场景. 这是关于使用Spring 3.1和Spring Secur ...
使用Spring Security对RESTful服务进行身份验证
1.概述本文重点介绍如何针对提供安全服务的安全REST API进行身份验证 -主要是RESTful用户帐户和身份验证服务. 2.目标首先,让我们看一下参与者-典型的启用了Spring Securi ...
带有Spring的REST的ETag
1.概述本文将重点介绍ETags-Spring支持,RESTful API的集成测试以及带有curl的使用场景. 这是关于使用Spring 3.1和Spring Security 3.1和基于Jav ...
ios中位置权限_iOS 13中的位置权限
ios中位置权限 In iOS 13, Apple made a big changes in location permission's behaviour, especially for &quo ...
《幸福就在你身边》第八课、幸福比成功更重要【哈佛大学幸福课精华】
一.财富与幸福 "财富不是幸福之源".什么会让你感到幸福,一所大房子,一辆好车,还是一位更性感或者更善解人意的伴侣?积极心理学权威米哈伊·西卡森特米哈伊曾问过一个问题:" ...

【Ctsc2011】幸福路径

【Ctsc2011】幸福路径相关推荐

最新文章

热门文章