poj 3621 二分+spfa判負環

http://poj.org/problem?id=3621

求一個環的{點權和}除以{邊權和}，使得那個環在所有環中{點權和}除以{邊權和}最大。

0/1整數劃分問題

令在一個環裏，點權爲v[i]，對應的邊權爲e[i]， 
即要求：∑（i=1，n）v[i]/∑（i=1，n)e[i]最大的環（n爲環的點數）， 
設題目答案爲ans， 
即對於所有的環都有 ∑（i=1，n）（v[i]）/∑（i=1，n)（e[i]）<=ans 
變形得ans* ∑（i=1，n）（e[i]）>=∑（i=1，n）（v[i]） 
再得 ∑（i=1，n）（ans*e[i]-v[i]） >=0 
稍分析一下，就有： 
當k<ans時，就存在至少一個環∑（i=1，n）（k*e[i]-v[i]）<0，即有負權迴路(邊權爲k*e[i]-v[i])； 
當k>=ans時，就對於所有的環∑（i=1，n）（k*e[i]-v[i]）>=0，即沒有負權迴路。 
然後我們就可以使新的邊權爲k*e[i]-v[i]，用spfa來判斷負權迴路,二分ans。

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <cstring>
#include <algorithm>
#define RD(x) scanf("%d",&x)
#define RD2(x,y) scanf("%d%d",&x,&y)
#define RD3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define clr0(x) memset(x,0,sizeof(x))
#define clr1(x) memset(x,-1,sizeof(x))
using namespace std;
typedef long long LL;
const int maxn = 1005,maxm = 105;
const double eps = 1e-3;
const int inf = 0x7fffffff;
int v[maxn];
struct edge{
    int v,w,next;
    edge(){};
    edge(int vv,int ww,int nnext):v(vv),w(ww),next(nnext){};
}e[maxn*5];
int head[maxn],vis[maxn],_v[maxn],cnt[maxn],ecnt,n,m;
double dist[maxn];
void add(int u,int v,int w)
{
    e[ecnt] = edge(v,w,head[u]);
    head[u] = ecnt++;
}
void init()
{
    for(int i = 1;i <= n;++i)
        RD(_v[i]);
    ecnt = 0;
    clr1(head);//判負環的初始化
    int u,v,w;
    while(m--){
        RD3(u,v,w);
        add(u,v,w);
    }
    return ;
}
bool spfa(double mid)
{
    clr0(vis),clr0(cnt);
    fill(dist,dist + n + 1,inf);
    dist[1] = 0;
    queue<int> q;
    q.push(1);
    while(!q.empty()){
        int u = q.front();
        q.pop();
        vis[u] = false;
        cnt[u]++;
        if(cnt[u] > n)
            return true;
        for(int i = head[u];i != -1;i = e[i].next){
            int v = e[i].v;
            double tmp = mid*e[i].w - _v[v];//"邊權"
            if(dist[u] + tmp < dist[v]){
                dist[v] = dist[u] + tmp;
                if(!vis[v]){
                    vis[v] = true;
                    q.push(v);
                }
            }
        }
    }
    return false;
}
void work()
{
    double l = 0,r = 10000,mid,ans;
    while(r - l > eps){
        mid = (l + r)/2;
        if(spfa(mid)){
            ans = mid;
            l = mid - 0.000001;
        }else
            r = mid + 0.000001;
    }
    printf("%.2f\n",ans);
}
int main()
{
    while(~RD2(n,m)){
        init();
        work();
    }
    return 0;
}