UVA 11090 判负圈问题

题目链接http://vjudge.net/problem/viewProblem.action?id=34650

题目大意：

给定n个点m条边的加权有向图，求平均权值最小的回路。
平均权值=路径权值之和/路径边数

我们可以通过找到他其中的最小和最大值，然后通过二分不断查找满足的点，然后尽可能的取到它的最大值，因为这里保留两位有效小数，所以设立

while(la-st>0.001)即可

找到一个满足的值是要建立在图上的每一条线段减去这个值后便能得到一个负圈，我们通常用spfa判负圈。

这个spfa只是用来作判断并不是算最短路径，为了防止出现多个连通分量，你只从一个点开始可能遍历整个图，所以最开始就把节点全放入队列中，dp[i]的

值全设定为0，如果出现负圈，则会为了找到最小值一直循环更新，我们用cnt[]数组，当某个点被访问了n次以上时，说明出现了负圈。

具体spfa函数如下所示：

bool spfa()
{
    for(int i=1;i<=n;i++) cnt[i]=0,visit[i]=0;
    visit[1]=1;
    queue<int> q;
    for(int i=1;i<=n;i++) q.push(i),dp[i]=0;
    while(!q.empty()){
        int u=q.front();
        visit[u]=0;
        q.pop();
        for(int i=first[u];i!=-1;i=path[i].next){
            if(dp[path[i].y]>dp[u]+path[i].d){
                dp[path[i].y]=dp[u]+path[i].d;
                if(!visit[path[i].y]) {
                    q.push(path[i].y);
                    if(++cnt[path[i].y]>n) return true;
                }
            }
        }
    }
    return false;
}

我们每次减去一个要找到的mid值，那我们一定要在判断结束后给它加回来以免下次判断出意外

总代码如下：

#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
using namespace std;
#define N 52
int first[N],visit[N],k,cnt[N],n;
double dp[N],maxn,minn;
struct Path{
    int y,next;
    double d;
}path[100010];

void add(int a,int b,double c)
{
    path[k].y=b,path[k].d=c,path[k].next=first[a];
    first[a]=k++;
}

bool spfa()
{
    for(int i=1;i<=n;i++) cnt[i]=0,visit[i]=0;
    visit[1]=1;
    queue<int> q;
    for(int i=1;i<=n;i++) q.push(i),dp[i]=0;
    while(!q.empty()){
        int u=q.front();
        visit[u]=0;
        q.pop();
        for(int i=first[u];i!=-1;i=path[i].next){
            if(dp[path[i].y]>dp[u]+path[i].d){
                dp[path[i].y]=dp[u]+path[i].d;
                if(!visit[path[i].y]) {
                    q.push(path[i].y);
                    if(++cnt[path[i].y]>n) return true;
                }
            }
        }
    }
    return false;
}

bool findMid(double mid)
{
    for(int i=0;i<k;i++) path[i].d-=mid;
    bool temp=spfa();
    for(int i=0;i<k;i++) path[i].d+=mid;
    return temp;
}

int main()
{
    int T,m,a,b;
    double c,st,la,mid,ans=-1;
    scanf("%d",&T);
    for(int j=1;j<=T;j++){
        memset(first,-1,sizeof(first));
        scanf("%d%d",&n,&m);
        k=0;
        maxn=0,minn=10000005;
        for(int i=0;i<m;i++)
        {
            scanf("%d%d%lf",&a,&b,&c);
            add(a,b,c);
            maxn=max(maxn,c);
            minn=min(minn,c);
        }
        st=minn-1,la=maxn;
        printf("Case #%d: ",j);
        if(!findMid(la+1)){printf("No cycle found.
");continue;}
        while(la-st>0.001){
            mid=st+(la-st)/2;
            if(findMid(mid)) la=mid;
            else ans=mid,st=mid;
        }
        printf("%.2f
",ans);
    }
    return 0;
}