Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 459  Solved: 114

Description

在一家公司中,人事部经理与业务部经理不和。一次,总经理要求人事部从公司的职员中挑选出一些来帮助业务部经理完成一项任务。人事部经理发现,在公司的所有职员中,有一些人相处得很不好。如果把他们同时放在一个工作小组中,他们将会给主管带来不小的麻烦。人事部经理还知道,在一个工作小组中,主管人员的麻烦程度可以认为是(带来麻烦的人的对数/总人数) 。于是,人事部经理决定选取这样的一些人,使得业务部经理的麻烦度最大。你的任务是帮助人事部经理达到他的目的。 在这样的一个公司中,保证公司的人数小于100 人,带来麻烦的员工不会多于1000 对。

Input

第一行给出N,M代表共有N个人,M对人之间有怨念 下面M行,每行两个数,代表这两个人有怨念

Output

输出最多能找出多少个人,人数>=1

Sample Input

5 6
1 5
5 4
4 2
2 5
1 2
3 1

Sample Output

4

HINT

Source

图论 最大密度子图 01分数规划 网络流

将人抽象成点,关系抽象成边,则“带来麻烦的人的对数/总人数”就是一个子图中边数/点数,这玩意儿叫做图的密度。

原题POJ3155 要求输出密度最大时的任一方案

这里要求输出密度最大时最多能选出多少个人。

设密度r=边数/点数,显然是一个01分数规划问题。

解法1:

  将每条边看做一个点,对于原图中的一条无向边<u,v>,从代表<u,v>的点向点u和点v各连一条边,容量为INF;

  从S向<u,v>连边,容量为1

  从u和v向T连边,容量为r

  ↑原问题转化成了最大权闭合子图问题。

解法2:

  证明见胡伯涛《最小割模型在信息学竞赛中的应用》

  U=m

  从S向u连边,容量为U

  从v向T连边,容量为U+2*r-deg[v]

  对于边<u,v>从u向v连边,容量为1

  若$(m*n-maxflow)/2>0$说明r可以扩大

这样就求出了最大的r

再用这个r建一遍图,跑最大流,在残量网络上DFS就可以找出所有可选的点。

就可以过POJ3155

把输出方案去掉就可以过Bzoj1312

但是博主傻傻不能理解为什么这样贪心一定能选出最多的人

 /*by SilverN*/
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<queue>
#define LL long long
using namespace std;
const int INF=0x3f3f3f3f;
const double eps=1e-;
const int mxn=;
int read(){
int x=,f=;char ch=getchar();
while(ch<'' || ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>='' && ch<=''){x=x*+ch-'';ch=getchar();}
return x*f;
}
struct EG{
int x,y;
}eg[mxn<<];
int deg[mxn];
//
struct edge{
int v,nxt;
double f;
}e[mxn<<];
int hd[mxn],mct=;
inline void add_edge(int u,int v,double f){
e[++mct].v=v;e[mct].nxt=hd[u];e[mct].f=f;hd[u]=mct;return;
}
void insert(int u,int v,double f){
add_edge(u,v,f); add_edge(v,u,);
return;
}
//
int n,m,S,T,U;
int d[mxn];
bool BFS(){
memset(d,,sizeof d);
queue<int>q;
d[S]=;
q.push(S);
while(!q.empty()){
int u=q.front();q.pop();
for(int i=hd[u];i;i=e[i].nxt){
int v=e[i].v;
if(!d[v] && e[i].f>){
d[v]=d[u]+;
q.push(v);
}
}
}
return d[T];
}
double DFS(int u,double lim){
if(u==T)return lim;
double f=,tmp;
for(int i=hd[u];i;i=e[i].nxt){
int v=e[i].v;
if(d[v]==d[u]+ && e[i].f>eps && (tmp=DFS(v,min(lim,e[i].f)))){
e[i].f-=tmp;
e[i^].f+=tmp;
lim-=tmp;
f+=tmp;
if(fabs(lim)<eps)return f;
}
}
d[u]=;
return f;
}
double Dinic(){
double res=;
while(BFS())res+=DFS(S,INF);
return res;
}
void Build(double r){
memset(hd,,sizeof hd);mct=;
S=;T=n+;int i;
for(i=;i<=n;i++){
insert(S,i,U);
insert(i,T,U+*r-deg[i]);
}
for(i=;i<=m;i++){
add_edge(eg[i].x,eg[i].y,);
add_edge(eg[i].y,eg[i].x,);
}
return;
}
bool use[mxn];int cnt=;
void DFS(int u){
use[u]=;++cnt;
for(int i=hd[u];i;i=e[i].nxt){
if(!use[e[i].v] && e[i].f>eps){DFS(e[i].v);}
}
return;
}
void solve(){
double l=,r=m; U=m;
double X=1.0/n/n;
while(r-l>X){
double mid=(l+r)/;
Build(mid);
// printf("mid:%.3f\n",mid);
if((m*n-Dinic())/>=eps){
l=mid;
}else r=mid;
}
Build(l);Dinic();
return;
}
void init(){
memset(hd,,sizeof hd);mct=;
memset(use,,sizeof use);cnt=;
return;
}
int main(){
while(scanf("%d%d",&n,&m)!=EOF){
if(!m){printf("1\n1\n");continue;}
init();
for(int i=;i<=m;i++){
eg[i].x=read();eg[i].y=read();
++deg[eg[i].x];
++deg[eg[i].y];
}
solve();
DFS(S);
printf("%d\n",cnt-);
for(int i=;i<=n;i++)if(use[i])printf("%d\n",i);
break;
}
return ;
}

Bzoj1312 / POJ3155 Neerc2006 Hard Life的更多相关文章

  1. POJ3155 Hard Life

    Time Limit: 8000MS   Memory Limit: 65536K Total Submissions: 8482   Accepted: 2461 Case Time Limit:  ...

  2. bzoj2503&poj3387[NEERC2006]IdealFrame

    其实只是把别人的题解强行扩写了 写这篇题解之前我不会的预备知识: 欧拉通路:从图中一个点出发不重复地遍历所有边的路径(可以停在另一个点) 欧拉回路:从图中一个点出发不重复地遍历所有边的回路(必须回到出 ...

  3. 最大密集子图(01分数规划+二分+最小割)POJ3155

    题意:给出一副连通图,求出一个子图令g=sigma(E)/sigma(V); h[g]=sigma(E)-g*sigma(V):设G是最优值 则当h[g]>0:g<G h[g]<0, ...

  4. bzoj1312

    忘写题解了,经典的最大密度子图 可以类似分数规划的做,二分密度,然后转化为最大权闭合子图做,判断是否大于0 注意方案的输出 const eps=1e-6; lim=1e-12; inf=; type ...

  5. POJ3155 Hard Life [最大密度子图]

      题意:最大密度子图 #include<iostream> #include<cstdio> #include<cstring> #include<algo ...

  6. BZOJ.1312.[Neerc2006]Hard Life(分数规划 最大权闭合子图)

    BZOJ 最大密度子图. 二分答案\(x\),转为求是否存在方案满足:\(边数-x*点数\geq 0\). 选一条边就必须选两个点,所以可以转成最大权闭合子图.边有\(1\)的正权,点有\(x\)的负 ...

  7. poj3155 最大密度子图

    求最大密度子图 记得在最后一次寻找的时候记得将进入的边放大那么一点点,这样有利于当每条边都满流的情况下会选择点 #include <iostream> #include <algor ...

  8. 【POJ3155】生活的艰辛Hard Life

    题面 Description ADN公司内部共 n个员工,员工之间可能曾经因为小事有了过节,总是闹矛盾.若员工u和员工 v有矛盾,用边(u, v)表示,共 m个矛盾.最近,ADN公司内部越来越不团结, ...

  9. poj分类 很好很有层次感。

    初期: 一.基本算法:      (1)枚举. (poj1753,poj2965)      (2)贪心(poj1328,poj2109,poj2586)      (3)递归和分治法.      ( ...

随机推荐

  1. java实现几种简单的排序算法

    public class SimpleAri { public static void main(String[] args) { int[] t = {11, 21, 22, 1, 6, 10, 3 ...

  2. ACM 第六天

    图论 网络流 最大流 INF(初始值) 路径上权值最小的边,决定流量大小. 流量网络的三个特性: ①流量控制 ②反对称性 ③流量守恒 残余网络:保留了c(e)容量<f(e)流量[可以继续流,因为 ...

  3. <Effective C++>读书摘要--Implementations<一>

    1.For the most part, coming up with appropriate definitions for your classes (and class templates) a ...

  4. PAT 甲级 1101 Quick Sort

    https://pintia.cn/problem-sets/994805342720868352/problems/994805366343188480 There is a classical p ...

  5. Cacti自定义脚本,监测Docker信息(Script/Command方式)

    一 环境背景 监控主机A:192.168.24.231:被监控主机B:192.168.24.233 A/B主机,通过公私钥建立ssh连接 [操作B主机时不需要输入密码,详见笔记:http://app. ...

  6. c# 调用c++ 使用指针传递的时候

    http://www.cnblogs.com/warensoft/archive/2011/12/09/warenosoft3d.html 上面这篇文章很好解释了. 简单记录一下: 1. 声明  注意 ...

  7. 从一个简单的main方法执行谈谈JVM工作机制

    本来JVM的工作原理浅到可以泛泛而谈,但如果真的想把JVM工作机制弄清楚,实在是很难,涉及到的知识领域太多.所以,本文通过简单的mian方法执行,浅谈JVM工作原理,看看JVM里面都发生了什么. 先上 ...

  8. 移除 ios 上 input 的默认样式

    input{ -webkit-appearance:none; }

  9. 【Json】C#格式化JSON字符串

    很多时候我们需要将json字符串以 {     "status": 1,     "sum": 9 }这种方式显示,而从服务端取回来的时候往往是这样 {&quo ...

  10. byte字节的输入流 建议使用字节数组形式接受内容 因为当输入流是汉字时候 会超出字节的范围 出现无法读取的现象

    byte字节的输入流 建议使用字节数组形式接受内容 因为当输入流是汉字时候 会超出字节的范围 出现无法读取的现象