poj 3155 Hard Life 分数规划+最大权闭合图

分类： 网络流2013-10-25
23:54 274人阅读 评论(0) 收藏 举报

这个就是完全按照论文中的做法来搞了，输出割边的时候就dfs判断一下在残余网络中该点是否可达汇点，不可达则是割边。

输出答案的证明可以百度到。

[cpp] view
plaincopy

#include <iostream>

#include <cstdio>

#include <cstring>

using namespace std;

const int maxn=1e2+9,inf=1e9;

const double epx=1e-6;

int n,m;

int d[maxn][maxn],level[maxn+1300],que[maxn+1300];

int visit[maxn+1300];

int S=0,T=1200;

int head[maxn+1300],lon;

struct

{

int next,to;

double c;

}e[111111];

void edgeini()

{

memset(head,-1,sizeof(head));

lon=-1;

}

void edgemake(int from,int to,double c)

{

e[++lon].c=c;

e[lon].to=to;

e[lon].next=head[from];

head[from]=lon;

}

void makegraph(double ret)

{

edgeini();

int tot=n;

for(int i=1;i<=n;i++)

{

edgemake(i,T,ret);

edgemake(T,i,0);

}

for(int i=1;i<=n;i++)

for(int j=1;j<=n;j++)

if(d[i][j])

{

++tot;

edgemake(S,tot,1);

edgemake(tot,S,0);

edgemake(tot,i,inf);

edgemake(i,tot,0);

edgemake(tot,j,inf);

edgemake(j,tot,0);

}

}

bool makelevel(int s,int t)

{

memset(level,0,sizeof(level));

int front=1,end=0;

que[++end]=s;

level[s]=1;

while(front<=end)

{

int u=que[front++];

// cout<<u<<endl;

if(u==t) return true;

for(int k=head[u];k!=-1;k=e[k].next)

{

int v=e[k].to;

if(e[k].c>0&&!level[v])

{

level[v]=level[u]+1;

que[++end]=v;

}

}

}

return false;

}

double dfs(int now,int t,double maxf)

{

double ret=0;

if(now==t) return maxf;

for(int k=head[now];k!=-1;k=e[k].next)

{

int u=e[k].to;

if(level[u]==level[now]+1&&e[k].c>0)

{

double f=dfs(u,t,min(maxf-ret,e[k].c));

e[k].c-=f;

e[k^1].c+=f;

ret+=f;

}

}

return ret;

}

double maxflow(int s,int t)

{

double ret=0;

while(makelevel(s,t))

{

ret+=dfs(s,t,inf);

}

return ret;

}

bool chk()

{

double ans=maxflow(S,T);

// cout<<m-ans<<endl;

return (m-ans)>epx;

}

bool dfs(int now,int t)

{

// printf("%d %d\n",now,t);

visit[now]=1;

if(now==t) return true;

for(int k=head[now];k!=-1;k=e[k].next)

{

int u=e[k].to;

if(e[k].c>0&&!visit[u])

if(dfs(u,t)) return true;

}

return false;

}

void prin()

{

// for(int i=0;i<=T;i++)

// for(int k=head[i];k!=-1;k=e[k].next)

//// if(e[k].c>0)

// printf("%d %d %f\n",i,e[k].to,e[k].c);

//

// dfs(4,T);

int sum=0;

for(int i=1;i<=n;i++)

{

memset(visit,0,sizeof(visit));

if(!dfs(i,T))

{

sum++;

}

}

cout<<sum<<endl;

for(int i=1;i<=n;i++)

{

memset(visit,0,sizeof(visit));

if(!dfs(i,T))

{

printf("%d\n",i);

}

}

}

int main()

{

// freopen("in.txt","r",stdin);

while(scanf("%d%d",&n,&m)!=EOF)

{

memset(d,0,sizeof(d));

for(int i=1,from,to;i<=m;i++)

{

scanf("%d%d",&from,&to);

d[from][to]=1;

}

double l=0,r=n,mid;

while(r-l>epx)

{

mid=(l+r)/2;

makegraph(mid);

// cout<<mid<<endl;

if(chk()) l=mid;

else r=mid;

//prin();

}

// printf("%f\n",epx);

// cout<<l<<endl;

makegraph(l);

chk();

prin();

}

return 0;

}