bzoj 2306: [Ctsc2011]幸福路径 倍增+floyed

       这道题果然是倍增卡精度。。。。

       首先可以得到走2^k步时，f[i][j]表示i~j的最大收益；那么第2^(k+1)步就有f[i][j]=min{f[i][k]+f[k][j]*ρ^(2^(k+1))}，然后当ρ<eps的时候退出即可。

AC代码如下：

#include<iostream>#include<cstdio>#include<cstring>#define N 105using namespace std;int n,m; double a[N],f[N][N],g[N][N];int main(){scanf("%d%d",&n,&m); int i,j,k;for (i=1; i<=n; i++) scanf("%lf",&a[i]);int sta; scanf("%d",&sta);for (i=1; i<=n; i++)for (j=1; j<=n; j++) if (i!=j) f[i][j]=-1e100;int x,y; double p,ans=0; scanf("%lf",&p);while (m--){scanf("%d%d",&x,&y); f[x][y]=a[y]*p;}for (; p>1e-10; p*=p){for (i=1; i<=n; i++) for (j=1; j<=n; j++) g[i][j]=-1e100;for (k=1; k<=n; k++)for (i=1; i<=n; i++)for (j=1; j<=n; j++)g[i][j]=max(g[i][j],f[i][k]+f[k][j]*p);memcpy(f,g,sizeof(g));}for (i=1; i<=n; i++) ans=max(ans,f[sta][i]);printf("%.1f\n",ans+a[sta]);return 0;}

by lych

2016.4.16