题目地址:http://poj.org/problem?id=3621
思路:
1.求点权值除以边权值最大的环。
2.二分答案ans,sigma(v)/sigma(e)<=ans,即sigma(ans*e-v)>=0。对原图中的边,设其边权值为ans*e-v,判断是否存在负环即可。
3.由于图有可能不连通,设一点0,向所有点连权值为0的边,从0点开始跑spfa判断是否存在负环即可。
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define debug
using namespace std;
const int maxn=1000+50;
const double eps=1e-6;
const int INF=0x3f3f3f3f;
struct Edge
{
int nt,to;
double w;
};
struct Node
{
int u,v;
double w;
};
int n,m,tot;
queue<int> q;
Node e[maxn*10];
Edge edge[maxn*10];
double dist[maxn],a[maxn];
int head[maxn],vis[maxn],cnt[maxn];
void addEdge(int u,int v,double w)
{
edge[tot].to=v,edge[tot].w=w;
edge[tot].nt=head[u],head[u]=tot++;
}
int check(int s,double x)
{
tot=0;
memset(head,-1,sizeof(head));
for(int i=1;i<=n;i++)
{
addEdge(0,i,0);
}
for(int i=1; i<=m; i++)
{
addEdge(e[i].u,e[i].v,x*e[i].w-a[e[i].u]);
}
while(!q.empty()) q.pop();
for(int i=0; i<=n; i++)
{
dist[i]=INF;
vis[i]=cnt[i]=0;
}
q.push(s),vis[s]=1,dist[s]=0,cnt[s]=1;
while(!q.empty())
{
int now=q.front();
q.pop(),vis[now]=0;
for(int i=head[now]; ~i; i=edge[i].nt)
{
int nt=edge[i].to;
double w=edge[i].w;
if(dist[nt]>dist[now]+w)
{
dist[nt]=dist[now]+w;
if(!vis[nt])
{
vis[nt]=1;
q.push(nt);
if(++cnt[nt]>=n)
{
return 1;
}
}
}
}
}
return 0;
}
int main()
{
#ifdef debu
freopen("in.txt","r",stdin);
#endif // debug
scanf("%d%d",&n,&m);
for(int i=1; i<=n; i++)
{
scanf("%lf",&a[i]);
}
for(int i=1; i<=m; i++)
{
int u,v;
double w;
scanf("%d%d%lf",&e[i].u,&e[i].v,&e[i].w);
}
double l=0.0,r=100000.0,ans=0.0;
while(fabs(r-l)>eps)
{
double mid=(l+r)/2;
if(!check(0,mid))
{
ans=mid;
r=mid;
}
else
{
l=mid;
}
}
if(fabs(ans)<eps) printf("0\n");
else printf("%.2f\n",ans);
return 0;
}