1. 程式人生 > >【POJ3155】Hard Life 分數規劃+最小割

【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; }