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

Posted fcwww

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了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);
}

 

以上是关于BZOJ_1486_[HNOI2009]最小圈_01分数规划的主要内容,如果未能解决你的问题,请参考以下文章

bzoj 1486: [HNOI2009]最小圈

bzoj1486HNOI2009最小圈

BZOJ1486 [HNOI2009]最小圈

bzoj1486[HNOI2009]最小圈

[BZOJ1486][HNOI2009]最小圈

[BZOJ1486][HNOI2009]最小圈