[BZOJ] 幸福路径

bzoj 2306

最短路水题，多次迭代就可以求得答案。
我在这里迭代了100下，应该很安全了orz

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int N = 110;

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

int main(){
	int i,j,k,T,x,y;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;i++) scanf("%lf",&a[i]);
	scanf("%d %lf",&s,&p);
	memset(f,0xc2,sizeof f);
	for(i=1;i<=n;i++) f[i][i]=0;
	for(i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		f[x][y]=a[y];
	}
	double temp = p;
	for(T = 0 ; T <= 100 ; T++,temp*=temp){
		memset(g,0xc2,sizeof g);
		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]*temp);
		memcpy(f,g,sizeof f);
	}
	for(i=1;i<=n;i++)
		ans = max(ans , f[s][i]);
	printf("%.1lf\n",ans * p +a[s]);
	return 0;
}