【BZOJ1486】【HNOI2009】最小圈 分数规划 dfs判负环。

7 篇文章 0 订阅
2 篇文章 0 订阅

链接:

#include <stdio.h>
int main()
{
    puts("转载请注明出处[辗转山河弋流歌 by 空灰冰魂]谢谢");
    puts("网址:blog.csdn.net/vmurder/article/details/46348771");
}

题解:

分数规划Qwq。
然而它卡判点入n次的那种spfa判断负环。
于是有了一种黑科技:
我们从枚举点 i 开始 dfs ,然后扫到点 j 时,保持 i~j 这一条链上的点被标记,然后强行判断再扫一个点 k 时,是否会到这个链上,然后是不是能重新更新此点 k 与 i 的距离。。。
这个东西是指数级别时间复杂度的,然而却可以过这道题。

代码:

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 3010
#define M 10100
#define eps 1e-8
using namespace std;
double mid;
struct Eli
{
    int v,n;
    double f,l;
    void re(){l=f-mid;}
}e[M];
int head[N],cnt;
inline void add(int u,int v,double l)
{
    e[++cnt].v=v;
    e[cnt].f=l;
    e[cnt].n=head[u];
    head[u]=cnt;
}

int n,m;
bool vis[N];
double dist[N];
bool dfs(int x)
{
    vis[x]=1;
    int i,v;
    for(i=head[x];i;i=e[i].n)
    {
        v=e[i].v;
        if(dist[v]>dist[x]+e[i].l)
        {
            if(vis[v])
            {
                vis[x]=0;
                return 1;
            }
            else {
                dist[v]=dist[x]+e[i].l;
                if(dfs(v))
                {
                    vis[x]=0;
                    return 1;
                }
            }
        }
    }
    vis[x]=0;
    return 0;
}

bool check()
{
    memset(dist,0,sizeof dist);
    for(int i=1;i<=m;i++)e[i].re();
    for(int i=1;i<=n;i++)if(dfs(i))return 1;
    return 0;
}
int main()
{
    int i,a,b;
    double c;

    scanf("%d%d",&n,&m);
    for(i=1;i<=m;i++)
    {
        scanf("%d%d%lf",&a,&b,&c);
        add(a,b,c);
    }
    double l=-1e7,r=1e7;
    for(i=60;i--;)
    {
        mid=(l+r)/2.0;
        if(check())r=mid;
        else l=mid;
    }
    printf("%.8lf\n",l);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值