uva11090 Bellman-Ford 运用

给定一一个n个点m条边的加权有向图， 平均值最小的回路。

二分答案，对于每个二分的mid 做一次Bellman-Fprd , 假设有k条边组成的回路。 回路上各条边的权值为  w1 , w2 ..wk ，

那么平均值小于mid意味着w1+w2+w3..+wk< k*mid 即：

（w1 - min）+(w2-mid)+...+(w2-mid)<0；

也就是说 这k条边能组成 一个负环，用 Bellman_Ford 来检查

 #include <iostream>
 #include <cstdio>
 #include <algorithm>
 #include <string.h>
 #include <vector>
 #include <queue>
 #include <cmath>
 using namespace std;
 const int maxn = ;
 int cmp(double a ,double b){
      if(fabs(a-b)<=0.00000001) return ;
      return a-b>?:-;
 }
 struct Edge{
    int from,to;
    double dist;
 };
 struct BellmanFord{
     int n,m;
     vector<Edge> edges;
      vector<int> G[maxn];
      bool inq[maxn];
      double d[maxn];
      int p[maxn];
      int cnt[maxn];
      void inti(int n){
          m=;
         this->n = n;
         for(int i=; i<n; ++i ) G[i].clear();
         edges.clear();
      }
      void AddEdge(int form, int to, double dist){
         edges.push_back((Edge){form,to,dist});
         m =edges.size();
         G[form].push_back(m-);
      }
      bool negativeCycle(){
          queue<int> Q;
          memset(inq, , sizeof(inq));
          memset(cnt, , sizeof(cnt));
          for(int i=; i < n; ++i)  {  d[i] =; inq[i] = true; Q.push(i);}
          while(!Q.empty()){
             int u = Q.front(); Q.pop();
             inq[u] = false;
             for(int i=; i<G[u].size(); i++){
                 Edge &e = edges[G[u][i]];
                 if(cmp(d[e.to] , d[u] + e.dist)>){
                      d[e.to] = d[u] +e.dist;
                      p[e.to] = G[u][i];
                      if(!inq[e.to]){
                          Q.push(e.to); inq[e.to] = true;
                          if(++cnt[e.to]>n) return true;
                      }
                 }
             }
          }
          return false;
      }
 }solver;
 bool test(double x){
      for(int i=; i<solver.m; i++){
         solver.edges[i].dist-=x;
      }
      bool ret = solver.negativeCycle();
      for(int i= ; i<solver.m; i++)
          solver.edges[i].dist+=x;
      return ret;
 }
 int main()
 {
    int T;
    scanf("%d",&T);
    for(int kase =; kase<=T; kase++){
        int n,m;
        scanf("%d%d",&n,&m);
        solver.inti(n);
        int ub =;
        while(m--){
           int u,v,w;
           scanf("%d%d%d",&u,&v,&w); u--,v--; ub= max(ub,w);
           solver.AddEdge(u,v,w);
        }
        printf("Case #%d: ",kase);
        double ans = ub;
        if(!test(ub+)){
          printf("No cycle found.\n");
        }else{
             double L =,R= ub;
             while(R-L>1e-){
                 double M = L+(R-L)/;
                 if(test(M)){
                         R=M;
                 }else {
                     L=M;
                 }
             }
             printf("%.2lf\n",L);
        }

    }
     return ;
 }