本质上就是倍增Floyd

由于体力越来越小，所以在走一定次数之后剩下的部分就可以忽略不计了，我们可以搞出邻接矩阵，然后自乘 50 50 50 次即可。
代码如下：

#include <cstdio>
#include <cstring>
#include <iostream>
 
using namespace std;
double ans;
double v[110],bas;
int n,m,S;
double f[110][110],g[110][110];
int main()
{
	scanf("%d%d",&n,&m);
	int i,j,k,l,a,b;
	for(i=1;i<=n;i++)	for(j=1;j<=n;j++)	f[i][j]=g[i][j]=-99999999;
	for(i=1;i<=n;i++)	scanf("%lf",&v[i]),f[i][i]=0;
	scanf("%d%lf",&S,&bas);
	for(i=1;i<=m;i++)	scanf("%d%d",&a,&b),f[a][b]=v[b]*bas;
	for(l=0;l<=50;l++)
	{
		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]+bas*f[k][j]);
		for(i=1;i<=n;i++)	for(j=1;j<=n;j++)	f[i][j]=g[i][j],g[i][j]=-99999999;
		bas=bas*bas;
	}
	for(i=1;i<=n;i++)	ans=max(ans,v[S]+f[S][i]);
	printf("%.1lf",ans);
	return 0;
}