POJ 3155 Hard Life

最新推荐文章于 2020-10-29 12:27:41 发布

nlj1999

最新推荐文章于 2020-10-29 12:27:41 发布

阅读量368

点赞数

分类专栏：最小割分数规划 !@#$%^&*()

本文链接：https://blog.csdn.net/nlj1999/article/details/50723336

版权

!@#$%^&*() 同时被 3 个专栏收录

20 篇文章 0 订阅

订阅专栏

最小割

4 篇文章 0 订阅

订阅专栏

分数规划

4 篇文章 0 订阅

订阅专栏

首先很明显是最大密度子图（感觉没学过啊怎么一下子就脑补出来这个名词了）

然后……

不会

于是去Orz 07 年胡伯(#)涛的论文，然后就可做了。

结果……精度这东西明显是玄学啊，WA了快10次了

不过看到了Discuss里面wty大爷的小tip，就过了。。。。玄学！玄学！

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<algorithm>
using namespace std;
const double eps=1e-9;
const double inf=1e9;
const double U=1e4;
struct Edge{int to,next;double v,cap;}e[200005];
int head[105],cnt,d[105];
void ins(int u,int v,double w){
	e[++cnt]=(Edge){v,head[u],w,w};head[u]=cnt;
}
void insert1(int u,int v,double w){
	ins(u,v,w);ins(v,u,0);
}
void insert2(int u,int v,double w){
	ins(u,v,w);ins(v,u,w);
}
int dcmp(double x){
	if(fabs(x)<eps)return 0;
	return x<0?-1:1;
}
bool bfs(int s,int t){
	memset(d,-1,sizeof(d));
	d[s]=0;queue<int>q;q.push(s);
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int i=head[u];i;i=e[i].next)
		if(dcmp(e[i].v)&&d[e[i].to]==-1){
			d[e[i].to]=d[u]+1;
			q.push(e[i].to);
		}
	}
	return d[t]!=-1;
}
double dfs(int x,double f,int t){
	if(x==t)return f;
	double flow=0,w;
	for(int i=head[x];i;i=e[i].next)
	if(dcmp(e[i].v)&&d[e[i].to]==d[x]+1){
		w=dfs(e[i].to,min(f-flow,e[i].v),t);
		e[i].v-=w;e[i^1].v+=w;flow+=w;
		if(!dcmp(f-flow))return flow;
	}
	if(!dcmp(flow))d[x]=-1;
	return flow;
}
double dinic(int s,int t){double ans=0;while(bfs(s,t))ans+=dfs(s,inf,t);return ans;}
void init(){
	memset(head,0,sizeof(head));
	cnt=1;
}
int S,T,du[105];
int u[1005],v[1005];
int n,m;
bool check(double g){
	init();
	for(int i=1;i<=n;i++){
		insert1(S,i,U);
		insert1(i,T,U+2*g-du[i]);
	}
	for(int i=1;i<=m;i++)
	insert2(u[i],v[i],1);
	double ans=dinic(S,T);
	return dcmp(U*n-ans)>0;
}
int ans[105],tot;
bool vis[105];
void dfs(int u){
	vis[u]=true;
	if(u>=1&&u<=n)ans[++tot]=u;
	for(int i=head[u];i;i=e[i].next)
	if(e[i].v>0&&!vis[e[i].to])
	dfs(e[i].to);
}
int main(){
	//freopen("a.in","r",stdin);
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)scanf("%d%d",&u[i],&v[i]),du[u[i]]++,du[v[i]]++;
	S=0;T=n+1;
	double l=0,r=m,lim=1.0/n/n;
	while(r-l>lim){
		double mid=(l+r)/2;
		if(check(mid))l=mid;
		else r=mid;
	}
	check(l-(1e-4));
	dfs(S);
	sort(ans+1,ans+1+tot);
	printf("%d\n",tot);
	for(int i=1;i<=tot;i++)
	printf("%d\n",ans[i]);
	return 0;
}