码迷,mamicode.com
首页 > 其他好文 > 详细

【POJ3155】Hard Life 分数规划+最小割

时间:2015-06-10 10:28:04      阅读:161      评论:0      收藏:0      [点我收藏+]

标签:poj3155   hard-life   分数规划   最小割   

链接:

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

题解:

如题。先算出那个分数值,然后看有哪些人还与源点相连。
最小割建图:原图每个点对应一个点,原图每条边对应一个点。每条边对应点向两端点对应点连边,注意要单向边。

这道题卡精度:

所以一些细节问题扒代码吧Qwq
eps:1e-5
因为是double网络流,所以二分上界别太大,边数就好。
另外这份代码在BZOJ并不能AC、
而在BZOJ能AC的代码也并不一定能在POJ上AC、
什么鬼。

代码:

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 12000
#define M 70000
#define eps 1e-6
#define inf 0x3f3f3f3f
using namespace std;
struct ELi
{
    int v,n;
    double l;
}e[M];
int head[N],cnt;
inline void add(int u,int v,double l)
{
    e[++cnt].v=v;
    e[cnt].l=l;
    e[cnt].n=head[u];
    head[u]=cnt;
}
inline void Add(int u,int v,double l)
{add(u,v,l),add(v,u,0);}
int s,t,dep[N];
queue<int>q;
bool bfs()
{
    while(!q.empty())q.pop();
    memset(dep,0,sizeof dep);
    q.push(s),dep[s]=1;
    int i,u,v;
    while(!q.empty())
    {
        u=q.front(),q.pop();
        for(i=head[u];i;i=e[i].n)
        {
            if(!dep[v=e[i].v]&&e[i].l>eps)
            {
                dep[v]=dep[u]+1;
                if(v==t)return 1;
                q.push(v);
            }
        }
    }
    return 0;
}
double dinic(int x,double flow)
{
    if(x==t)return flow;
    int i,v;
    double remain=flow,k;
    for(i=head[x];i&&remain>eps;i=e[i].n)
    {
        if(dep[v=e[i].v]==dep[x]+1&&e[i].l>eps)
        {
            k=dinic(v,min(remain,e[i].l));
            if(k<eps)dep[v]=0;
            e[i].l-=k,e[i^1].l+=k;
            remain-=k;
        }
    }
    return flow-remain;
}
int U[N],V[N],n,m;
bool check(double mid)
{
    int i;
    memset(head,0,sizeof head);

    for(cnt=i=1;i<=m;i++)
    {
        Add(n+i,U[i],inf);
        Add(n+i,V[i],inf);
    }
    for(i=1;i<=n;i++)Add(i,t,mid);
    for(i=1;i<=m;i++)Add(s,n+i,1);

    double maxflow=m;
    while(bfs())
        maxflow-=dinic(s,inf);
    return maxflow>eps;
}
int main()
{
//  freopen("test.in","r",stdin);
    int i;

    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(!m)
        {
            puts("1\n1");
            continue;
        }
        s=0,t=n+m+1;
        for(i=1;i<=m;i++)scanf("%d%d",&U[i],&V[i]);

        double l=0,r=m,mid;
        while(l+1e-5<r)
        {
            mid=(l+r)/2.0;
            if(check(mid))l=mid;
            else r=mid;
        }
        check(l),bfs();

        int ans=0;
        for(i=1;i<=n;i++)if(dep[i])ans++;
        printf("%d\n",ans);

        for(i=1;i<=n;i++)if(dep[i])
            printf("%d\n",i);
    }
    return 0;
}

【POJ3155】Hard Life 分数规划+最小割

标签:poj3155   hard-life   分数规划   最小割   

原文地址:http://blog.csdn.net/vmurder/article/details/46437961

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!