poj3155 最大密度子图

最新推荐文章于 2020-04-12 20:44:49 发布
最新推荐文章于 2020-04-12 20:44:49 发布
阅读量498 收藏
点赞数
分类专栏： ACM/ICPC
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u012915516/article/details/46884053
版权

给定一个无向图，要求它的一个子图，使得子图中边数|E|与点数|V|的比值最大。
解法一：
假设答案为k ，则要求解的问题是：选出一个合适的点集V和边集E，令(|E| - k * |V|) 取得最大值。所谓“合适”是指满足如下限制：若选择某条边，则必选择其两端点。
建图:以原图的边作为左侧顶点，权值为1；原图的点作为右侧顶点，权值为-k(相当于 支出k)。
若原图中存在边(u,v)，则新图中添加两条边([uv]–>u), ([uv]–>v)，转换为最大权闭合子图。
解法二：
把原图中的无向边转换成两条有向边，容量为1。
设一源点，连接所有点，容量为U（取m）。
设一汇点，所有点连接汇点，容量为 U+2g-dv 。
二分枚举最大密度g，其中dv为v的度。
判断(U*n-MaxFlow)/2.>=0。
最后跳出的L就是最大密度。
拿这个L再重新建图，求最大流。
然后从s出发bfs或者dfs，走残留容量大于0的边，所有能到达的点就是答案。

方法1code:

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <math.h>
#include <vector>
#include <cstdio>
#include <string>
#include<string.h>
#include <fstream>
#include <iostream>
#include <algorithm>
using namespace std;
#define exp 1e-8
#define INF 0x3f3f3f3f
#define ll long long
#define mm(a,b) memset(a,b,sizeof(a));
#define for1(a,b) for(int a=1;a<=b;a++)//1---(b)
#define rep(a,b,c) for(int a=b;a<=c;a++)//b---c
void bug(string st="bug")
{cout<<st<<endl;}
template<typename __ll>
inline void READ(__ll &m){
    __ll x=0,f=1;char ch=getchar();
    while(!(ch>='0'&&ch<='9')){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    m=x*f;
}
template<typename __ll>
inline void read(__ll &m){READ(m);}
template<typename __ll>
inline void read(__ll &m,__ll &a){READ(m);READ(a);}
template<typename __ll>
inline void read(__ll &m,__ll &a,__ll &b){READ(m);READ(a);READ(b);}
const int maxn=20000*5; //边的数量，应该比正常边的数量多两倍
const int MAXN=10000;//节点数

int s,t,tot;
int cur[MAXN];//当前弧优化
int pre[MAXN];//先祖
int level[MAXN];//高度
int gap[MAXN];//gap优化
int dat[maxn][3];
int d[MAXN];
int n,m;
struct EDGE
{
    int to,next;
    double cost;
}edge[maxn];
int head[maxn],cnt;
inline void addedge(int a,int b,double c)
{
    edge[cnt].to=b;
    edge[cnt].cost=c;
    edge[cnt].next=head[a];
    head[a]=cnt++;

    edge[cnt].to=a;
    edge[cnt].cost=0;
    edge[cnt].next=head[b];
    head[b]=cnt++;
}
void build_graph(int n,int m,double mid,double &sum) //节点数(1->n)，边数
{
    sum=0;
    s=0,t=n+m+1,tot=t+1,cnt=0;//cnt=0为了更好的异或
    memset(head,-1,sizeof(head));
    for1(i,m)addedge(s,i,1),addedge(i,dat[i][1]+m,INF),addedge(i,dat[i][2]+m,INF),sum++;
    for1(i,n)addedge(i+m,t,mid);
}
double sap()  //无需改动
{
    memset(pre,-1,sizeof(pre));
    memset(gap,0,sizeof(gap));
    memset(level,0,sizeof(level));
        gap[0]=tot;
    int u=s;
    double ans=0.,inf=INF;
    for(int i=0; i<tot; i++)//初始化当前弧为第一条弧
        cur[i]=head[i];
    while(level[s]<tot)
    {
        bool flag=0;
        for(int &i=cur[u];i!=-1;i=edge[i].next) //读取边，并更新当前弧
        {
            int v=edge[i].to;
            if(edge[i].cost>0&&level[u]==level[v]+1)
            {
                flag=1;
                if(edge[i].cost<inf)
                    inf=edge[i].cost;
                pre[v]=u;
                u=v;
                if(u==t)
                {
                    while(u!=s) //结束后u==s
                    {
                        u=pre[u];
                        edge[cur[u]].cost-=inf;
                        edge[cur[u]^1].cost+=inf;//异或是找与其配对的边
                    }
                    ans+=inf;
                    inf=INF;
                }
                break;
            }
        }
        if(!flag)
        {
            if(--gap[level[u]]==0)//gap优化
                break;
            int minn=tot;
            for(int i=head[u];i!=-1;i=edge[i].next)
            {
                int v=edge[i].to;
                if(edge[i].cost>0&&level[v]<minn)
                {
                    minn=level[v];
                    cur[u]=i;//更新当前弧，让其从度最小的线段开始
                }
            }
            level[u]=minn+1;
            gap[minn+1]++;
            if(u!=s)
                u=pre[u];
        }
    }
    return ans;
}
bool ok(double mid)
{
    double sum=0;
    build_graph(n,m,mid,sum);
    return sum-sap()>=exp;
}
int vis[MAXN];
void dfs(int idx)
{
    vis[idx]=1;
    for(int i=head[idx];i!=-1;i=edge[i].next)
        if(!vis[edge[i].to]&&edge[i].cost)dfs(edge[i].to);
}
int main()
{
    while(scanf("%d %d",&n,&m)!=EOF)
    {
        mm(d,0);
        for1(i,m)read(dat[i][1],dat[i][2]);
        for1(i,m)d[dat[i][1]]++,d[dat[i][2]]++;
        if(m==0){puts("1\n1\n");continue;}
        double l=0,r=m;
        while(r-l>1./n/n)
        {
            double mid=(r+l)/2.;
            if(ok(mid))l=mid;
            else r=mid;
        }
        ok(l);
        mm(vis,0);
        dfs(s);
        int ans=0;
        rep(i,m+1,n+m)ans+=vis[i];
        cout<<ans<<endl;
        rep(i,m+1,n+m)if(vis[i])
            cout<<i-m<<endl;
    }
    return 0;
}

方法2code:

#include <map>
#include <set>
#include <queue>
#include <stack>
#include <math.h>
#include <vector>
#include <cstdio>
#include <string>
#include<string.h>
#include <fstream>
#include <iostream>
#include <algorithm>
using namespace std;
#define exp 1e-8
#define INF 0x3f3f3f3f
#define ll long long
#define mm(a,b) memset(a,b,sizeof(a));
#define for1(a,b) for(int a=1;a<=b;a++)//1---(b)
#define rep(a,b,c) for(int a=b;a<=c;a++)//b---c
void bug(string st="bug")
{cout<<st<<endl;}
template<typename __ll>
inline void READ(__ll &m){
    __ll x=0,f=1;char ch=getchar();
    while(!(ch>='0'&&ch<='9')){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    m=x*f;
}
template<typename __ll>
inline void read(__ll &m){READ(m);}
template<typename __ll>
inline void read(__ll &m,__ll &a){READ(m);READ(a);}
template<typename __ll>
inline void read(__ll &m,__ll &a,__ll &b){READ(m);READ(a);READ(b);}
#define INF 999999999
const int maxn=20000*5; //边的数量，应该比正常边的数量多两倍
const int MAXN=1000;//节点数

int s,t,tot;
int cur[MAXN];//当前弧优化
int pre[MAXN];//先祖
int level[MAXN];//高度
int gap[MAXN];//gap优化
int dat[maxn][3];
int d[MAXN];
int n,m;
struct EDGE
{
    int to,next;
    double cost;
}edge[maxn];
int head[maxn],cnt;
inline void addedge(int a,int b,double c)
{
    edge[cnt].to=b;
    edge[cnt].cost=c;
    edge[cnt].next=head[a];
    head[a]=cnt++;

    edge[cnt].to=a;
    edge[cnt].cost=0;
    edge[cnt].next=head[b];
    head[b]=cnt++;
}
void build_graph(int n,int m,double mid) //节点数(1->n)，边数
{
    s=0,t=n+1,tot=n+2,cnt=0;//cnt=0为了更好的异或
    memset(head,-1,sizeof(head));
    for1(i,n) addedge(s,i,m);
    for1(i,n) addedge(i,t,m+2*mid-d[i]);
    for1(i,m) addedge(dat[i][1],dat[i][2],1),addedge(dat[i][2],dat[i][1],1);
}
double sap()  //无需改动
{
    memset(pre,-1,sizeof(pre));
    memset(gap,0,sizeof(gap));
    memset(level,0,sizeof(level));
        gap[0]=tot;
    int u=s;
    double ans=0.,inf=INF;
    for(int i=0; i<tot; i++)//初始化当前弧为第一条弧
        cur[i]=head[i];
    while(level[s]<tot)
    {
        bool flag=0;
        for(int &i=cur[u];i!=-1;i=edge[i].next) //读取边，并更新当前弧
        {
            int v=edge[i].to;
            if(edge[i].cost>0&&level[u]==level[v]+1)
            {
                flag=1;
                if(edge[i].cost<inf)
                    inf=edge[i].cost;
                pre[v]=u;
                u=v;
                if(u==t)
                {
                    while(u!=s) //结束后u==s
                    {
                        u=pre[u];
                        edge[cur[u]].cost-=inf;
                        edge[cur[u]^1].cost+=inf;//异或是找与其配对的边
                    }
                    ans+=inf;
                    inf=INF;
                }
                break;
            }
        }
        if(!flag)
        {
            if(--gap[level[u]]==0)//gap优化
                break;
            int minn=tot;
            for(int i=head[u];i!=-1;i=edge[i].next)
            {
                int v=edge[i].to;
                if(edge[i].cost>0&&level[v]<minn)
                {
                    minn=level[v];
                    cur[u]=i;//更新当前弧，让其从度最小的线段开始
                }
            }
            level[u]=minn+1;
            gap[minn+1]++;
            if(u!=s)
                u=pre[u];
        }
    }
    return ans;
}
bool ok(double mid)
{
    build_graph(n,m,mid);
    return (1.0*n*m-sap())/2.>=exp;
}
int vis[MAXN];
void dfs(int idx)
{
    vis[idx]=1;
    for(int i=head[idx];i!=-1;i=edge[i].next)
        if(!vis[edge[i].to]&&edge[i].cost)dfs(edge[i].to);
}
int main()
{
    while(scanf("%d %d",&n,&m)!=EOF)
    {
        mm(d,0);
        for1(i,m)read(dat[i][1],dat[i][2]);
        for1(i,m)d[dat[i][1]]++,d[dat[i][2]]++;
        if(m==0){puts("1\n1\n");continue;}
        double l=1./n,r=m;
        while(r-l>1./n/n)
        {
            double mid=(r+l)/2.;
            if(ok(mid))l=mid;
            else r=mid;
        }
        ok(l);  //根据真实的比例搜出真实的割
        mm(vis,0);
        dfs(s);
        int ans=0;
        for1(i,n)ans+=vis[i];
        cout<<ans<<endl;
        for1(i,n)if(vis[i])
            cout<<i<<endl;
    }
    return 0;
}
粽子猪zZ
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
最大密度子图
taotao 的大学墓志
02-16 2685
参考文献 * 《最小割模型在信息学竞赛中的应用》分数规划分数规划问题 分数规划中有几个重要的结论，详细的证明过程不在给出,请参考论文原问题等价于 g(λ)=min{a(X)−λb(X)}g(\lambda) = min \{a(X)-\lambda b(X)\}重要性质 g(λ′)=0⇔λ′为原分数规划的最优解g(\lambda') = 0 \Leftrightarrow \lamb
poj 3155 (最大密度子图)
菜鸟起航的专栏
09-27 2058
题意：一个公司有n个人，给出了一些有冲突的人的对数(u,v)，公司决定裁人，那么总裁现在要裁掉冲突率最高的那些人（冲突率=在这些人中存在的冲突数/人数）。就是求出一些点，这些点之间的边数/点数最大。最大密度子图。 思路：胡伯涛的论文《最小割模型在信息学竞赛中的应用》介绍了两种方法： 第一种：转换为最大权闭合图的模型来求解： 设max g = f(x)= |E‘|/|V’| ，找一个子图的
无向图的极大团、最大团(Bron-Kerbosch算法)
热门推荐
yo_bc的博客
08-21 2万+
描述：团就是最大完全子图。(极大团) 给定无向图G=(V,E)。如果U包含于V，且对任意u，v属于U且有(u，v)属于E，则称U是G的完全子图。 G的完全子图U是G的团当且仅当U不包含在G的更大的完全子图中，即U就是最大完全子图。 G的最大团是指G中所含顶点数最多的团。 // 最大团： V中取K个顶点，两点间相互连接 // 最大独立集： V中取K个顶点，两点间不连接 // 最大...
POJ 3155 最大密度子图
neweryyy的博客
04-12 145
题意 传送门 POJ 3155 题解 算法思路 参考《最小割模型在信息学竞赛中的应用》。分数规划，求满足条件的 xxx 最大值 ∣E∣∣V∣≥x\frac{|E|}{|V|}\geq x∣V∣∣E∣​≥x ∣E∣|E|∣E∣ 为子图边数，∣V∣|V|∣V∣ 为子图点数，不等式左边为图的密度。不同子图密度差不小于 1/∣V∣21/|V|^21/∣V∣2，考虑到多边形点与边的关系，图密度的二分范围为 ...
poj openjudge 1111最大正向匹配
04-20
poj openjudge 1111题最大正向匹配 提交ac
poj.rar_poj
09-23
标题中的"poj.rar_poj"暗示了这是一个与POJ（Programming Online Judge）相关的压缩文件，POJ是一个在线编程竞赛平台，主要供程序员们练习和提交算法解决方案。在这个压缩包中，很可能包含了用户在POJ上参与编程挑战...
POJ.rar_poj
09-19
【标题】"POJ.rar_poj" 指的是一份与编程竞赛相关的压缩文件，主要涉及POJ（Problem Overloading Judge）平台上的问题解答。POJ是一个在线的编程竞赛平台，它为参赛者提供了大量的算法题目进行练习和比赛，以提升...
ASP外观专利图像检索平台(源代码+论文).rar
07-21
ASP外观专利图像检索平台(源代码+论文)
C++课程设计：电煤气管理系统【源码+文档】
07-21
C++课程设计：电煤气管理系统【源码+文档】 本程序是一个水电气管理信息系统，能够对高校的水电气费用进行管理, 包括了成员的基本信息，如学号、编号、姓名、成员水电气的用量等。程序的用途包括缴纳水电气费、查询一个同学水电气费用量、查看所有同学的缴费情况、增加学生信息、删除学生信息、退出系统等。在设计时也考虑到学生和教师在用水电气时的不同，学生可以免费使用一定额度的水电气，超过这个额度的以后必须付费，且付费部分水电气费的价格要高于教工的收费标准，该措施的实行是为了鼓励同学们节约资源，以免造成不必要的资源浪费。该软件主要是为了学校的管理人员提供便捷，以更快的完成水电气费用的收缴。该软件本着简洁明了，实用稳定为一体进行设计。 系统将实现以下功能： （1）实现对用户信息的录入； （2）实现水电煤气数据的录入； （3）实现计算并查询用户应缴费用，查询未缴纳费用的名单； （4）实现对人员的删除和添加；
ASP+ACCESS网上动态同学录系统（源代码+论文+系统说明+答辩PPT）.rar
07-21
ASP+ACCESS网上动态同学录系统（源代码+论文+系统说明+答辩PPT）
avaloniassification-mas笔记
07-21
avalonia
ASP+ACCESS学生排课管理系统毕业设计(源代码+论文).rar
07-21
ASP+ACCESS学生排课管理系统毕业设计(源代码+论文)
ASP+ACCESS在线教育系统设计(源代码+论文).rar
07-21
ASP+ACCESS在线教育系统设计(源代码+论文)
中国软件杯A8赛题.zip
07-21
【项目资源】：包含前端、后端、移动开发、人工智能、物联网、信息化管理、数据库、硬件开发、大数据、课程资源，毕业设计等各种技术项目的源码。包括C++、Java、python、web、C#、EDA等项目的源码。 【适用人群】：适用于希望学习不同技术领域的初学者或进阶学习者。可作为毕设项目、课程设计、大作业、工程实训或初期项目立项。 【附加价值】：项目具有较高的学习借鉴价值，也可直接拿来修改复刻。对于有一定基础或热衷于研究的人来说，可以在这些基础代码上进行修改和扩展，实现其他功能。 【沟通交流】：有任何使用上的问题，欢迎随时与博主沟通，博主会及时解答。鼓励下载和使用，并欢迎大家互相学习，共同进步。
配色识色器软件，易语言编写，非常好用
最新发布
07-21
配色识色器软件版本 编写语言：易语言 数据来源：网络 软件功能： 1、识取颜色并同时输出RGB、CMYK、16进制三种值，也可以RGB、CMYK互转； 2、可添加、导入EXCEL文件（XLS文件，配合附件中的导入工具用）的配色方案（表头：见图）； 3、可点选或搜索配色值进行配色。
java毕设&课设-ASP网络商城的设计与实现(源代码+论文).zip
07-21
计算机毕业设计资源包含（项目部署视频+源码+LW+开题报告等等），所有项目经过助教老师跑通，有问题可以私信博主解决，可以免费帮部署。
每日报告回放.pdf
07-21
每日报告回放
掌握CrowdStrike Falcon：端点检测与响应的终极指南
07-21
CrowdStrike是一家提供在线安全解决方案的公司，专注于提供基于云计算的端点保护平台。以下是关于CrowdStrike的一些关键信息： 1. **成立时间与总部**：CrowdStrike成立于2011年，总部位于美国加利福尼亚州。 2. **主要产品**：CrowdStrike的主要产品是Falcon平台，它利用人工智能和机器学习技术来检测、预防和响应网络威胁。 3. **市场地位**：CrowdStrike在价值86亿美元的“端点检测和响应”(EDR)软件市场中，所占份额约为18%，仅次于微软。 4. **技术优势**：CrowdStrike因其在检测和防御高级网络攻击方面的能力而闻名，其软件被包括微软、亚马逊AWS在内的一些最大的云服务公司提供商所使用，也包括主要的全球银行、医疗保健和能源公司。 5. **产品能力**：CrowdStrike的产品能力包括终端防病毒、威胁情报、机器学习等新技术新方向的布局，且在这些领域处于行业领先位置。 6. **SaaS模式**：CrowdStrike的Falcon终端安全平台完全基于SaaS模式，具备产品的敏捷性、易用性、可拓展性、定
POJ2135 Farm Tour —— 最小费用最大流
05-25
输出从 1 到 N 的最小费用最大流，如果无法到达，则输出 -1。 数据范围 1≤N≤1000 1≤M≤10000 1≤a,b≤N 1≤w≤10000 输入样例 4 5 1 2 1 2 3 2 3 4 3 1 3 4 2 4 3 输出样例 6 算法1 (最小费用最大流) $O...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值