【最短路】UVA11090 Going in Cycle!!解题报告

link

题目大意：

给出一个$N$个点$M$条边的有向图，求边权平均值最小的连通块。

题目分析：

首先，每个连通块的边权平均值为：$x=\frac{\sum _{k=i}^j{a}_k}{j-i+1}$ ($i$到$j$联通,$a_i$表示边权)。题目要求平均值最小，即，$ans>=\frac{\sum _{k=i}^j{a}_k}{j-i+1}$($i$到$j$联通) 。同时，$ans$是具有单调性的，所以我们考虑二分答案。

我们二分出一个$ans$。考虑对原式子进行化简：

$ans>=\frac{\sum _{k=i}^j{a}_k}{j-i+1}$($i$到$j$联通)

$ans\ast (j-i+1)>=\sum _{k=i}^j{a}_k$

$ans\ast (j-i+1)>=a_i+a_{i+1}+...+a_{j-1}+aj$

$(a_i-ans)+(a_{i+1}-ans)+...+(a_{j-1}-ans)+(a_j-ans)<=0$ ①

即，当我们二分出的一个答案$ans$在①式成立的情况下是合法的。

现在我们考虑如何维护①式。

观察可知，我么可以将原来的边权$a_i$转化为$a_i-ans$，然后判断在这个新图里面是否有负环，如果有，那么①式肯定成立。即我们二分出来的这个答案$ans$合法。

注意事项：

1.用$SPFA$判断是否存在负环是通过记录每个点入队的次数来实现的。如果一个点被超过$n$个点更新，则一定存在负环。

2.这个图有可能不连通。（比如存在$1->2,3->2$这种情况）所以在$SPFA$之前，需要将所有的点加入队列。

3.$SPFA$可以用双端队列或者优先队列优化。

4.由于我们二分的是一个实数，所以需要$eps$来确定精度。

5.由于是多组数据，所以要注意初始化。

6.如果我们$ans=Maxedge+1$仍没有负环，则说明该图没有环。即输出"No cycle found."

代码如下：

#include <bits/stdc++.h>
#define inf 1e8+7.0
#define eps 1e-8 
using namespace std;
const int maxn=1e4+5;


double dis[maxn];
int used[maxn],head[maxn],num[maxn];
double l,r;
int size=0,n,m,cas,t;

struct EDGE
{
	int u,v,nxt;
	double w;
}edge[maxn];

void add(int u,int v,double w)
{
	edge[size].u=u;
	edge[size].v=v;
	edge[size].w=w;
	edge[size].nxt=head[u];
	head[u]=size++;
}


bool SPFA()
{
	deque<int>q;
	for(int i=1;i<=n;i++)
	{
		used[i]=0;
		num[i]=0;
	}
	for(int i=1;i<=n;i++)
	q.push_back(i),dis[i]=0;
	
	while(!q.empty())
	{
		int u=q.front();
		q.pop_front();
		used[u]=0;
		for(int i=head[u];~i;i=edge[i].nxt)
		{
			int v=edge[i].v;
			if(dis[v]-(dis[u]+edge[i].w)>eps)
			{
				dis[v]=dis[u]+edge[i].w;
				if(used[v])continue;
				used[v]=0;
				if(!q.empty())
				{
					if(dis[q.front()]-dis[v]>eps)q.push_front(v);
					else q.push_back(v);
				}
				else q.push_back(v);
				num[v]++;
				if(num[v]>n)return true;
			}
		}
	}
	return false;
}

bool check(double mid)
{
	for(int i=0;i<size;i++)
	{
		edge[i].w-=mid;
	}
	bool flag=SPFA();
	for(int i=0;i<size;i++)
	{
		edge[i].w+=mid;
	}
	return flag;
}

void init()
{
	freopen("in.in","r",stdin);
}

void readdata()
{
	memset(head,-1,sizeof(head));
	memset(edge,0,sizeof(edge));
	size=0;
	l=inf,r=-inf;
	scanf("%d%d",&n,&m);
	int u,v;
	double w;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%lf",&u,&v,&w);
		add(u,v,w);
		r=max(r,w);
		l=min(l,w);
	}
	printf("Case #%d: ",++cas);
}

void work()
{
	if(!check(r+1))
	{
		printf("No cycle found.\n");  
        return;
	}	
	while(r-l>=eps)
	{
		double mid=l+(r-l)/2;
		if(check(mid))
		{
			r=mid;
		}
		else 
		{
			l=mid;
		}
	}
	printf("%.2lf\n",r);
}


int main()
{
	init();
	scanf("%d",&t);
	while(t--)
	{
		readdata();
		work();
	}
	return 0;
}