[BZOJ1486][HNOI2009]最小圈(01分数规划图+负环)

题目:

我是超链接

题解:

相当于是 wij1<L ∑ w i j ∑ 1 < L 要求更小的啊,然后 wijL<0 ∑ w i j − L < 0
那就是一个图的01分数规划了,如果有负环说明可以更小
注意求负环的时候dis数组要设为0,这样快很多

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int N=3005;
const int M=10005;
int tot,nxt[M],point[N],v[M],n;
double c[M],dis[N];bool vis[N],fff;
void addline(int x,int y,double z){++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y; c[tot]=z;}
void spfa(int now,double mid)
{
    vis[now]=1;
    for (int i=point[now];i;i=nxt[i])
      if (dis[v[i]]>dis[now]+c[i]-mid && !fff)
      {
        dis[v[i]]=dis[now]+c[i]-mid;
        if (vis[v[i]]) {fff=1;return;}
        spfa(v[i],mid);
      }
    vis[now]=0;
}
bool check(double mid)
{
    fff=0;
    memset(dis,0,sizeof(dis));
    memset(vis,0,sizeof(vis));
    for (int i=1;i<=n;i++) 
    {
        spfa(i,mid);
        if (fff) return 1;
    }
    return 0;
}
int main()
{
    int m;scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++)
    {
        int x,y;double w;scanf("%d%d%lf",&x,&y,&w);
        addline(x,y,w);
    }
    double l=-10000000,r=10000000;
    while (r-l>1e-10)
    {
        double mid=(l+r)/2;
        if (check(mid)) r=mid;
        else l=mid; 
    }
    printf("%.8lf",l);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值