UVA11090 Going in Cycle!!

analysis

二分+spfa判负环，绝了！

为什么会想到二分？因为问题可以转化为一个单调函数：

平均值越大，存在这种环的可能性就越小

也就是说可以二分平均值然后判断是否存在就可以了

判断如何存在的话，可以考虑这样的回路的性质：

如果对于一个平均值x，存在这样的一个回路A使得A的平均值小于x，那么

$\frac{{\sum }_{i\in A}len[i]}{|A|}<x$

${\sum }_{i\in A}len[i]<|A|x$

${\sum }_{i\in A}(len[i]-x)<0$

也就是说，如果我们把每条边的权值都减去x，那么如果图中有负环，就说明这个平均值偏大

于是就可以套spfa判负环的模板了

但是注意，如果没有特别要求，spfa判负环时要将所有点加入队列然后判负环，否则在不连通的图中会出现漏掉负环的情况

code

#include<bits/stdc++.h>
using namespace std;
#define loop(i,start,end) for(register int i=start;i<=end;++i)
#define anti_loop(i,start,end) for(register int i=start;i>=end;+--i)
#define clean(arry,num) memset(arry,num,sizeof(arry))
#define ll long long
template<typename T>void read(T &x){
	x=0;char r=getchar();T neg=1;
	while(r>'9'||r<'0'){if(r=='-')neg=-1;r=getchar();}
	while(r>='0'&&r<='9'){x=(x<<1)+(x<<3)+r-'0';r=getchar();}
	x*=neg;
}
int n,m;
int NFG=0;
const int maxn=50+10; 
const int maxm=maxn*maxn;
const double eps=1e-5;
struct node{
	int e;
	double w;
	int nxt;
}edge[maxm<<1];
int head[maxn];
int cnt=0;
inline void addl(int u,int v,double w){
	edge[cnt].e=v;
	edge[cnt].w=w;
	edge[cnt].nxt=head[u];
	head[u]=cnt++;
}
double dis[maxn];
int in[maxn];
queue<int>q;
inline void change(double x){
	loop(i,0,cnt-1)
		edge[i].w+=x;
}
inline bool spfa_circle(){
	clean(dis,0x3f);
	clean(in,0);
	while(q.empty()==false)q.pop();
	dis[1]=0;
	loop(i,1,n){
		q.push(i);
		dis[i]=0;
	}
	while(q.empty()==false){
		int f=q.front();
		q.pop();
		for(int i=head[f];i!=-1;i=edge[i].nxt){
			int v=edge[i].e;
			if(dis[v]>dis[f]+edge[i].w){
				dis[v]=dis[f]+edge[i].w;
				++in[v];
				if(in[v]>n)
					return true;
				q.push(v);
			}
		}
	}
	return false;
}
inline void bin(double L,double R){
	while(R-L>eps){
		double mid=L+(R-L)/2;
		change(-mid);
		if(!spfa_circle())
			L=mid;
		else R=mid;
		change(mid);
	}
	printf("Case #%d: %.2lf\n",NFG,L);
}
int main(){
	int T;
	read(T);
	while(T--){
		++NFG;
		clean(head,-1);
		cnt=0;
		read(n);
		read(m);
		double L=1e7,R=0;
		loop(i,1,m){
			int ui,vi;
			ll wi;
			read(ui);
			read(vi);
			read(wi);
			L=min(L,(double)wi);
			R=max(R,(double)wi);
			addl(ui,vi,(double)wi);
		}
		change(-(R+1));
		if(!spfa_circle()){
			printf("Case #%d: No cycle found.\n",NFG);
			continue;
		}
		change(R+1);
		bin(L,R);
	}
	return 0;
}