二分图 （ 未完待更...

最新推荐文章于 2021-05-31 21:04:15 发布
最新推荐文章于 2021-05-31 21:04:15 发布
阅读量99 收藏
点赞数
分类专栏： 图论--网络流
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lizhiwei2017/article/details/80137235
版权

过山车 （ 二分图最大匹配

#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#include<algorithm>
using namespace std;
const int N = 100000;
const int INF = 100000000; 
struct Node{
	int ne;
	int to;
	int w;
}e[N<<1];
int head[N],deth[N];
int n,m,k,x,y,st,en,val,cnt;
void init()
{
	memset(head,-1,sizeof(head));
	memset(deth,0,sizeof(deth));
	cnt = st = 0;
	en = n + m + 1;
}
void add(int u,int v,int val)
{
	e[cnt].to = v;
	e[cnt].w = val;
	e[cnt].ne = head[u];
	head[u] = cnt ++;
}
bool bfs()
{
	memset(deth,0,sizeof(deth));
	queue<int>q;
	deth[st] = 1;
	q.push(st);
	while(!q.empty())
	{
		int now = q.front();
		q.pop();
		for(int i=head[now];~i;i=e[i].ne)
		{
			int to = e[i].to;
			if(deth[to] == 0 && e[i].w)
			{
				deth[to] = deth[now] + 1;
				q.push(to);
			}
		}
	}
	if(deth[en] == 0)
		return 0;
	return 1;
}
int dfs(int u,int dist)
{
	if(u == en)
	return dist;
	int temp = dist;
	for(int i=head[u];~i;i=e[i].ne)
	{
		int to = e[i].to;
		if(deth[to] == deth[u] + 1 && e[i].w)
		{
			int di = dfs(to,min(e[i].w,temp));
			if(di>0)
			{
				temp -= di;
				e[i].w -= di;
				e[i^1].w += di;
				if(temp <= 0)
				return dist;
			}
		}
	}
	return dist - temp;
}
int dinic()
{
	int ans = 0;
	while(bfs())
	{
		while(int di = dfs(st,INF))
		ans += di;
	}
	return ans;
}
int main()
{
	while(scanf("%d",&k) && k)
	{
		scanf("%d%d",&n,&m);
		init();
		for(int i=1;i<=n;i++)
		{
			add(st,i,1);
			add(i,st,0);
		}
		for(int i=n+1;i<=n+m;i++)
		{
			add(i,en,1);
			add(en,i,0);
		}
		for(int i=1;i<=k;i++)
		{
			scanf("%d%d",&x,&y);
			add(x,n+y,1);
			add(n+y,x,0);
		}
		printf("%d\n",dinic());
	}
	return 0;
}

匈牙利做法

#include<iostream>
#include<cstring> 
#include<cstdio>
using namespace std;
const int maxn = 555;
const int maxm = 5555;
struct Node{
	int to;
	int ne;
}e[maxn];
int head[maxn],match[maxn];
int cnt,ans;
bool vis[maxn];
int k,m,n;
void init()
{
	memset(head,-1,sizeof(head));
	memset(match,-1,sizeof(match));
	cnt = ans = 0;
}
void add(int u,int v)
{
	e[cnt].to = v;
	e[cnt].ne = head[u];
	head[u] = cnt ++;
}
bool findpath(int u)
{
	for(int i=head[u];~i;i=e[i].ne)
	{
		int to = e[i].to;
		if(!vis[to])
		{
			vis[to] = 1;
			if(match[to] == -1 || findpath(match[to]))
			{
				match[to] = u;
				return true;
			}
		}
	}
	return false;
}
void Hungary()
{
	for(int i=1;i<=m;i++)
	{
		memset(vis,0,sizeof(vis));
		if(findpath(i))
		ans ++;
	}
}
int main()
{
	int u,v;
	while(scanf("%d",&k) && k)
	{
		init();
		scanf("%d%d",&m,&n);
		for(int i=1;i<=k;i++)
		{
			scanf("%d%d",&u,&v);
			add(u,v);
		}
		Hungary();
		printf("%d\n",ans);
	}
	return 0;
}


二分图多重匹配 

全文阅读请移步 二分图多重匹配 
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
const int N = 10000;
const int INF = 100000000;
struct Node{
	int to;
	int ne;
	int w;
}e[N<<1];
int head[N],deth[N];
int T,x,y,st,k,en,n,m,val,cnt;
void init()
{
	memset(head,-1,sizeof(head));
	cnt = st = 0;
	en = n + m + 1;
}
void add(int u,int v,int val)
{
	e[cnt].to = v;
	e[cnt].w = val;
	e[cnt].ne = head[u];
	head[u] = cnt ++;
}
int bfs()
{
	memset(deth,0,sizeof(deth));
	queue<int>q;
	deth[st] = 1;
	q.push(st);
	while(!q.empty())
	{
		int now = q.front();
		q.pop();
		for(int i=head[now];~i;i=e[i].ne)
		{
			int to = e[i].to;
			if(deth[to] == 0 && e[i].w>0)
			{
				deth[to] = deth[now] + 1;
				q.push(to);
			}
		}
	}
	if(deth[en] == 0)
	return 0;
	return 1;
}
int dfs(int u,int dist)
{
	if(u == en)
	return dist;
	int temp = dist;
	for(int i =head[u];~i;i=e[i].ne)
	{
		int to = e[i].to;
		if(deth[to] == deth[u] + 1 && e[i].w > 0)
		{
			int di = dfs(to,min(e[i].w,temp));
			if(di>0)
			{
				temp -= di;
				e[i].w -= di;
				e[i^1].w += di;
				if(temp <= 0)
				return dist;
			}
		}
	}
	return dist - temp;
}
int dinic()
{
	int ans = 0;
	while(bfs())
	{
		while(int di = dfs(st,INF))
		ans += di;
	}
	return ans;
} 
int main()
{
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d",&n,&m);
		init();
		int sum = 0;
		for(int i=1;i<=m;i++)
		{
			scanf("%d",&x);
			sum += x;
			add(n+i,en,x);
			add(en,n+i,0);
		}
		for(int i=1;i<=n;i++)
		{
			scanf("%d%d",&x,&y);
			add(st,i,x);
			add(i,st,0);
			while(y--)
			{
				scanf("%d",&k);
				add(i,n+k,1);
				add(n+k,i,0);
			}
		}
		if(sum == dinic())
		puts("Yes");
		else puts("No");
	}
	return 0;
}

假期的宿舍 ( 二分图完美匹配？？？

1. 有床的向 汇点 连边 
2.  源点 向需要床的连边，边权都为1 
3. 跑最大流判断maxflow 是否和需要住宿的人 sum 相等即可 

#include<iostream>
#include<queue>
#include<cstdio>
#include<cstring>
using namespace std;
const int N = 1000;
const int INF = 100000000;
struct Node{
	int ne;
	int to;
	int w;
}e[N<<1];
int head[N],deth[N];
int bed[N],home[N];
int st,en,n,m,x,y,val,cnt,k,T;
void init()
{
	memset(head,-1,sizeof(head));
	memset(deth,0,sizeof(deth));
	st = cnt = 0;
	en = n * 2 + 1;
}
void add(int u,int v,int val)
{
	e[cnt].to = v;
	e[cnt].w = val;
	e[cnt].ne = head[u];
	head[u] = cnt ++;
}
bool bfs()
{
	memset(deth,0,sizeof(deth));
	queue<int>q;
	deth[st] = 1;
	q.push(st);
	while(!q.empty())
	{
		int now = q.front();
		q.pop();
		for(int i=head[now];~i;i=e[i].ne)
		{
			int to = e[i].to;
			if(deth[to] == 0 && e[i].w>0)
			{
				deth[to] = deth[now] + 1;
				q.push(to);
			}
		}
	}
	if(deth[en] == 0)
	return 0;
	return 1;
}
int dfs(int u,int dist)
{
	if(u == en)
	return dist;
	int temp = dist;
	for(int i=head[u];~i;i=e[i].ne)
	{
		int to = e[i].to;
		if(deth[to] == deth[u] + 1 && e[i].w)
		{
			int di = dfs(to,min(e[i].w,temp));
			if(di>0)
			{
				temp -= di;
				e[i].w -= di;
				e[i^1].w += di;
				if(temp <= 0)
				return dist;
			}
		}
	}
	return dist - temp;
}
int dinic()
{
	int ans = 0;
	while(bfs())
	{
		while(int di = dfs(st,INF))
		ans += di; 
	}
	return ans;
}
int main()
{
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&n);
		init();
		int sum = 0;
		for(int i=1;i<=n;i++)
		{
			scanf("%d",&bed[i]);//有床的向汇点连边
			if(bed[i])
			{
				add(i+n,en,1);
				add(en,i+n,0);
			}
		}
		for(int i=1;i<=n;i++)
		{
			scanf("%d",&home[i]);
			if(!bed[i] || (bed[i] && !home[i])) //源点向需要床的连边
			{
				sum ++;
				add(st,i,1);
				add(i,st,0);
			}
		}
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=n;j++)
			{
				scanf("%d",&x);
				if(x ||i == j)
				{
					add(i,j+n,1);
					add(j+n,i,0);
				}
			}
		}
		if(sum == dinic())
		puts("^_^");
		else puts("T_T");
	}
	return 0;
}

二分图带权匹配

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<cmath>
using namespace std;
#define ABSs(x) ((x)<0?-(x):(x))
const int N = 105;
struct Node{
	int ne;
	int to;
	int w;
	int co;
}e[N*(N+2)<<1+5];
int head[N<<1+6],hh[N*N],mm[N*N];
int path[N<<1+5],v[N<<1+5],dis[N<<1+5];
int n,m,len,st,en,cnt;
char s[N];
void init()
{
	memset(head,-1,sizeof(head));
	cnt =  0;
}
void add(int u,int v,int co,int val)
{
	e[cnt].to = v;
	e[cnt].w = val;
	e[cnt].co = co;
	e[cnt].ne = head[u];
	head[u] = cnt ++;
}
bool spfa(int st,int en)
{
	memset(path,-1,sizeof(path));
	memset(v,0,sizeof(v));
	memset(dis,127,sizeof(dis));
	queue<int>q;
	q.push(st);
	v[st] = 1;
	dis[st] = 0;
	while(!q.empty())
	{
		int now = q.front();
		q.pop();
		for(int i=head[now];~i;i=e[i].ne)
		{
			int to = e[i].to;
			if(e[i].w > 0 && dis[to] > dis[now] + e[i].co)
			{
				dis[to] = dis[now] + e[i].co;
				path[to] = i;
				if(!v[to])
				{
					v[to] = 1;
					q.push(to);
				}
			}
		}
		v[now] = 0;
	}
	return dis[en] != dis[N<<1+4];
}
int MCMF(int st,int en)
{
	int ans = 0;
	while(spfa(st,en))
	{
		int minn = 1e9;
		for(int i=path[en];~i;i=path[e[i^1].to])
		{
			if(minn > e[i].w)
			minn = e[i].w;
		}
		ans += minn*dis[en];
		for(int i=path[en];~i;i=path[e[i^1].to])
		{
			e[i].w -= minn;
			e[i^1].w += minn;
		}
	}
	return ans;
}
int main()
{
	while(scanf("%d%d",&n,&m) && n+m)
	{
		init();
		int hlen = 0;
		int mlen = 0;
		for(int i=0;i<n;i++)
		{
			scanf("%s",s);
			for(int j=0;j<m;j++)
			{
				if(s[j] == 'H')
				{
					hh[hlen++] = i*m+j;
				}
				else if(s[j] == 'm')
				{
					mm[mlen ++] = i*m+j;
				}
			}
		}
		int x0,x1,y0,y1;
		for(int i=0;i<mlen;i++)
		{
			for(int j=0;j<hlen;j++)
			{
				x0 = mm[j]/m;
				y0 = mm[j]%m;
				x1 = hh[i]/m;
				y1 = hh[i]%m;
				add(i+1,j+mlen+1,ABSs(x0-x1) + ABSs(y0-y1),1);
				add(j+mlen+1,i+1,- ABSs(x0-x1) - ABSs(y0-y1),0);
			}
		}
		n = mlen;
		m = hlen;
		st = 0;
		en = n + m + 1;
		for(int i=0;i<n;i++)
		{
			add(st,i+1,0,1);
			add(i+1,st,0,0);
		}
		for(int j=0;j<m;j++)
		{
			add(j+n+1,en,0,1);
			add(en,j+n+1,0,0);
		}
		printf("%d\n",MCMF(st,en));
	}
	return 0;
}

Xiao-wei1998
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
图的基础总结
2301_80841939的博客
02-21 1016
邻接矩阵是一个二维数组，其中每个元素表示相应顶点对之间的边的权重。对于无向图，邻接矩阵是对称的；对于有向图，邻接矩阵不一定对称。邻接表则是一种链式存储结构，对于每个顶点，都有一个与之相邻的顶点的链表。除了邻接矩阵和邻接表之外，还有一些其他的图的表示方式，如边集数组、邻接多重表等。图的基本操作。
kingshard配置文件ks.yaml示例
03-12
kingshard配置文件ks.yaml示例
二分图匹配问题（未完待续）
Royen_的博客
05-31 193
匹配问题（学习笔记与例题） 文章目录匹配问题（学习笔记与例题）前置知识一、二分图1. 定义2. 最大匹配增广路径3. 匈牙利算法HK算法（Hopcroft Karp算法）4. 二分图的判断实战应用[kuangbin带你飞]专题十 匹配问题A - Fire NetB - The Accomodation of StudentsC - CoursesD - 棋盘游戏总结 前置知识 一、二分图 1. 定义 二分图又称作二部图，是图论中的一种特殊模型。 设 G =（V,E）是一个无向图，如果顶点V可分割为两个
CodeForces Round #403 (Div.2) A-F
dezhen7015的博客
03-06 115
精神不佳，选择了在场外同步划水 没想到实际做起来手感还好，早知道就报名了…… 该打 未完待续233 A. Andryusha and Socks 模拟，模拟大法好。注意每次是先判断完能不能收进柜子，再算桌子上的袜子数量的。 1 /*by SilverN*/ 2 #include<iostream> 3 #include<algo...
【知识点总结】【CSP考前复习】图论大杂烩【未完】
Cyan_rose
10-23 889
OI生涯大汇总.jpg
【复习】算法模版<二>
weixin_30507481的博客
05-12 115
从现在就开始复习算法，都忘干干静静了。 一、最短路问题： 1.dijstra无优化版 1 #include<iostream> 3 #define INF 0x7f 4 #define maxn 10010 5 #define maxm 50010 6 using namespace std; 7 int num,k,minl,n,m,s; ...
多目标跟踪算法SORT
weixin_43508499的博客
04-16 5238
目录一. 目标检测简介1.1 目标检测1.2 IOU指标二. 目标跟踪简介三. 卡尔曼滤波器3.1 什么是卡尔曼滤波器3.2 具体计算过程3.3 换个角度看四. 匈牙利算法五. SORT 一. 目标检测简介 给定一张图像，找出我们需要的类别位置，并给定一个检测框，检测框一般包含位置坐标，检测框置信度，每个类别的分数。 1.1 目标检测 1.2 IOU指标 IOU指标就是Intersection Over Union，也叫交并比。 如上图所示，红色框表示groundTruth，蓝色框和绿色框是两个检测框
读书笔记_算法第四版（二）
pifutan的博客
08-18 632
算法第四版（谢路云译） 官方网站：http://algs4.cs.princeton.edu/home/有部分源代码和部分课后习题答案。 个人练习代码：https://github.com/morefans/AlgorithmsFourthEdition   第4章 图 4.1 无向图 l  图是由一组顶点和一组能够将两个顶点相连的边组成的。 l  自环：即一条连接一个顶点和其他
一句话题解（20170801~20170125）
weixin_30332705的博客
09-23 553
8.1 bzoj 4720 noip2016 换教室 floyd预处理+期望（薛定谔的猫） bzoj 4318 OSU！ 三次函数期望值 从一次、二次推得 8.2 bzoj 1076 状压+期望DP 逆拓扑序（贪心常用手段防止现在过度影响未来）lim边界问题曾WA多次 bzoj 1012 树状数组（线段树）第二 暴力更新后缀第一（明显可卡） 暴力查询要TLE 单调栈没调出来 //bzo...
ACM模版-f_zyj v 2.0——更新通知
逐梦者
09-03 4790
ACM模版距离 ACM模版-f_zyj v 1.1\text{ACM模版-f_zyj v 1.1} 版成工已经一年整了，这一年，我每次发现其中有不足时，都会在我在博客 ACM在线模版－f-zyj\text{ACM在线模版－f-zyj} 中对其进行更新，稀稀拉拉的一年过去了，我发现增删改的地方实在不少，所以总是有朋友问我什么时候会将这些更新整理到 PDFPDF 格式中……这个嘛，我总是一拖再拖，因为整
11.综合布线系统图（未完）.dwg
06-08
11.综合布线系统图（未完）.dwg
大数据Java面试总结（未完待续）.docx
11-12
大数据Java面试总结（未完待续）.docx
openresty-1.13.6.1-win32
01-22
openresty 是一个基于 Nginx 与Lua 的高性能 Web 平台,其内部集成了大量精良的 Lua 库、第三方模块以及大多数的依赖项。用于方便地搭建能够处理超高并发
2、大数据篇(未完待续...).pdf
05-04
Hadoop、Samba、Spark、MySQL、HIve、Scala等在Centos7下的安装详细教程！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！
Java_Spring Boot 3主分支2其他分支和Spring Cloud微服务的分布式配置演示Spring Cl.zip
05-22
Java_Spring Boot 3主分支2其他分支和Spring Cloud微服务的分布式配置演示Spring Cl
ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计
05-22
ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)+编程项目+毕业设计ERP客户关系系统设计(含源代码+毕业设计文档)
基于MATLAB实现的V两幅图像中有重叠部分，通过数字图像相关算法可以找到两幅图像相同的点+使用说明文档.rar
05-22
CSDN IT狂飙上传的代码均可运行，功能ok的情况下才上传的，直接替换数据即可使用，小白也能轻松上手 【资源说明】 基于MATLAB实现的V两幅图像中有重叠部分，通过数字图像相关算法可以找到两幅图像相同的点+使用说明文档.rar 1、代码压缩包内容 主函数：main.m； 调用函数：其他m文件；无需运行 运行结果效果图； 2、代码运行版本 Matlab 2020b；若运行有误，根据提示GPT修改；若不会，私信博主（问题描述要详细）； 3、运行操作步骤 步骤一：将所有文件放到Matlab的当前文件夹中； 步骤二：双击打开main.m文件； 步骤三：点击运行，等程序运行完得到结果； 4、仿真咨询 如需其他服务，可后台私信博主； 4.1 期刊或参考文献复现 4.2 Matlab程序定制 4.3 科研合作 功率谱估计： 故障诊断分析： 雷达通信：雷达LFM、MIMO、成像、定位、干扰、检测、信号分析、脉冲压缩 滤波估计：SOC估计 目标定位：WSN定位、滤波跟踪、目标定位 生物电信号：肌电信号EMG、脑电信号EEG、心电信号ECG 通信系统：DOA估计、编码译码、变分模态分解、管道泄漏、滤波器、数字信号处理+传输+分析+去噪、数字信号调制、误码率、信号估计、DTMF、信号检测识别融合、LEACH协议、信号检测、水声通信 5、欢迎下载，沟通交流，互相学习，共同进步！
全球国家列表和国家代码最详细版本
最新发布
05-22
全球国家列表和国家代码最全最详细版本，国家country，code
cholmod-config.cmake
08-25
如果您需要更多关于cholmod-config.cmake文件的信息，请参考引用中提供的链接。<span class="em">1</span><span class="em">2</span><span class="em">3 #### 引用[.reference_title] - *1* *2* *3* [CMake知识补充...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值