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 105
using 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