BZOJ_1486_[HNOI2009]最小圈_01分数规划

BZOJ_1486_[HNOI2009]最小圈_01分数规划

Description

Input

Output

Sample Input

4 5
1 2 5
2 3 5
3 1 5
2 4 3
4 1 3

Sample Output

3.66666667

 

---

二分答案，边权减去答案，判负环即可。

然而spfa判负环会T掉，于是我使用了dfs判负环的方法。

dfs判负环代码：

void dfs(int x)
{
    vis[x]=1;
    int i;
    for(i=head[x];i&&!ok;i=nxt[i]){
        if(dis[to[i]]>dis[x]+cost[i]){
            dis[to[i]]=dis[x]+cost[i];
            if(vis[to[i]]){
                ok=1;return ;
            }
            dfs(to[i]);
        }
    }
    vis[x]=0;
}

不断更新最小值，直到更新了一圈回来，则说明有负环存在。

总之是比spfa快到不知哪里去。

dis数组清零即可。

 

代码：

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 3050
#define M 10050
#define du double
#define eps 1e-9
int head[N],to[M],nxt[M],val[M],cnt,n,m;
int dep[N],inq[N],Q[N],l,r,ok,vis[N];
du dis[N],cost[M];
inline void add(int u,int v,int w)
{
    to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;val[cnt]=w;
}
void dfs(int x)
{
    vis[x]=1;
    int i;
    for(i=head[x];i&&!ok;i=nxt[i]){
        if(dis[to[i]]>dis[x]+cost[i]){
            dis[to[i]]=dis[x]+cost[i];
            if(vis[to[i]]){
                ok=1;return ;
            }
            dfs(to[i]);
        }
    }
    vis[x]=0;
}
bool check(du x)
{
    int i;
    for(i=1;i<=m;i++) cost[i]=val[i]-x;
    for(i=1;i<=n;i++) dis[i]=vis[i]=0;
    for(i=1,ok=0;i<=n;i++)
    {
        dfs(i);
        if(ok) return 1;
    }
    return 0;
}
int main()
{
    scanf("%d%d",&n,&m);
    int i,x,y,z;
    for(i=1;i<=m;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z);
    }
    du L=-100000,R=100000;
    while(R-L>eps)
    {
        du mid=(L+R)*0.5;
        if(check(mid)) R=mid;
        else L=mid;
    }
    printf("%.8lf\n",L);
}