UVA 11090 Going in Cycle!! SPFA判断负环+二分

题目链接:

https://vjudge.net/problem/UVA-11090

题意:

给你一个有向图，问你定义一个环的平均值为这个环上所有边的平均值，问你最小的环的平均值是多少

题解:

二分答案：若当前的二分值是mid，让所有的边都减去这个值，如果此时图中出现负环，则说明有环的平均值比这个更小

假设一个包含k条边的回路，回路上各条边的权值为w1,w2……wk，那么平均值小于mid意味着 w1+w2+……wk< k* mid即：
（w1-mid)+(w2-mid)+……(wk-mid)<0
有负环，当前平均值就大了

注意图可能不连通，每个没访问过得点都spfa判一遍

代码:

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define MS(a) memset(a,0,sizeof(a))
#define MP make_pair
#define PB push_back
const int INF = 0x3f3f3f3f;
const ll INFLL = 0x3f3f3f3f3f3f3f3fLL;
inline ll read(){
    ll x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
//
const int maxn = 50+10;
double eps = 1e-6;

int n,m;

vector<pair<int,double> > G[maxn];
int vis[maxn],inq[maxn],cnt[maxn];
double d[maxn];

bool spfa(int u){
    MS(inq);MS(cnt);
    for(int i=0; i<=n; i++) d[i] = INF;
    queue<int> q;
    q.push(u),inq[u]=1,d[u]=0;
    while(!q.empty()){
        int now = q.front(); q.pop(); inq[now] = 0;
        vis[now] = 1;
        for(int i=0; i<(int)G[now].size(); i++){
            int v = G[now][i].first; double w = G[now][i].second;
            if(d[v] > d[now]+w){
                d[v] = d[now]+w;
                if(inq[v]) continue;
                inq[v] = 1;
                q.push(v);
                cnt[v]++;
                if(cnt[v] > n) return true;
            }
        }
    }
    return false;
}

bool check(double x){
    for(int i=1; i<=n; i++)
        for(int j=0; j<(int)G[i].size(); j++)
            G[i][j].second -= x; // 让所有的边都减去mid，如果此时图中出现负环，则说明有环的平均值比这个更小
    MS(vis);
    bool flag;
    for(int i=1; i<=n; i++){
        if(!vis[i]) // 图可能不连通，每个没访问过得点都判一遍
            flag = spfa(i);
        if(flag) break;
    }

    for(int i=1; i<=n; i++)
        for(int j=0; j<(int)G[i].size(); j++)
            G[i][j].second += x;

    return flag; // flag为true 就是有负环，当前平均值就大了
}


int main(){
    int T = read();
    for(int cas=1; cas<=T; cas++){
        for (int i = 0; i <= n; i++)G[i].clear();
        scanf("%d%d",&n,&m);
        for(int i=1; i<=m; i++){
            int u,v; double w;
            scanf("%d%d%lf",&u,&v,&w);
            G[u].push_back(MP(v,w));
        }

        double L=0,R=1e7,ans=1e7;
        while((R-L)>eps){
            double mid = (L+R)/2;
            if(check(mid)) ans=mid,R=mid-eps;
            else L=mid+eps;
        }
        if(abs(ans-1e7)<eps) printf("Case #%d: No cycle found.\n",cas);
        else printf("Case #%d: %.2lf\n",cas,ans);
    }

    return 0;
}