bzoj2306 幸福路径 倍增 Floyd

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

题意：一张有向图，每一个点有一个权值$w(x)$，给出路径起点求出最大$f(x)=sigma(w(x)*p)$，其中，$p$初始值为$1$，每走一步这个值都会乘上另外一个给出的常量。

因为这个题精度要求极低（只有$1e-1$），因此咱们直接迭代求值便可。

可是若是咱们这么一步一步搞确定会$T$……这时候咱们就须要用一些黑科技：倍增。咱们每次倍增前进，前进的时候在每一层作一次$Floyd$，作完以后合并再到下一层。

这样就能够愉快的解决啦~最后不要忘记加上起点时候的权值~

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 using namespace std;
 6 const int maxn=105,maxm=1005;
 7 int n,m;double a[maxn],f[maxn][maxn],g[maxn][maxn];
 8 int haha()
 9 {
10     scanf("%d%d",&n,&m);
11     for(int i=1;i<=n;i++)scanf("%lf",&a[i]);
12     int s;scanf("%d",&s);
13     for(int i=1;i<=n;i++)
14         for(int j=1;j<=n;j++)f[i][j]=i==j?0:-1e100;
15     double p;scanf("%lf",&p);
16     for(int i=1;i<=m;i++)
17     {
18         int x,y;scanf("%d%d",&x,&y);
19         f[x][y]=a[y]*p;
20     }
21     for(;p>1e-10;p*=p)
22     {
23         for(int i=1;i<=n;i++)
24             for(int j=1;j<=n;j++)g[i][j]=-1e100;
25         for(int k=1;k<=n;k++)
26             for(int i=1;i<=n;i++)
27                 for(int j=1;j<=n;j++)g[i][j]=max(g[i][j],f[i][k]+f[k][j]*p);
28         memcpy(f,g,sizeof(g));
29     }
30     double ans=0;
31     for(int i=1;i<=n;i++)ans=max(ans,f[s][i]);
32     printf("%0.1lf\n",ans+a[s]);
33 }
34 int sb=haha();
35 int main(){;}

bzoj2306