room 二分图最大匹配KM
#include<bits/stdc++.h>
#define fi first
#define se second
#define INF 0x3f3f3f3f
#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
#define pqueue priority_queue
#define NEW(a,b) memset(a,b,sizeof(a))
const double pi=4.0*atan(1.0);
const double e=exp(1.0);
const int maxn=1e5+;
typedef long long LL;
typedef unsigned long long ULL;
const LL mod=1e9+;
const ULL base=1e7+;
using namespace std;
struct room{
int s[];
}a[],b[];
inline int get(room x,room y){
int sum=;
for(int i=;i<;i++){
for(int j=;j<;j++){
if(x.s[i]==y.s[j]){
sum++;
}
}
}
return sum;
}
int cx[],cy[];
int ch[][];
bool visx[],visy[];
int lx[],ly[];
int n;
int inc;
inline bool dfs(int u){
visx[u]=;
for(int i=;i<n;i++){
if(!visy[i]){
int tmp=lx[u]+ly[i]-ch[u][i];
if(tmp==){
visy[i]=;
if(cy[i]==-||dfs(cy[i])){
cx[u]=i;
cy[i]=u;
return ;
}
}
else if(tmp>){
inc=min(inc,tmp);
}
}
}
return ;
}
inline void KM(){
memset(cy,-,sizeof(cy));
memset(cx,-,sizeof(cx));
for(int i=;i<n;i++){
lx[i]=-INF;
ly[i]=;
for(int j=;j<n;j++){
lx[i]=max(lx[i],ch[i][j]);
}
}
for(int u=;u<n;u++){
while(){
memset(visx,,sizeof(visx));
memset(visy,,sizeof(visy));
inc=INF;
if(dfs(u)) break;
for(int i=;i<n;i++){
if(visx[i]){
lx[i]-=inc;
}
}
for(int i=;i<n;i++){
if(visy[i]){
ly[i]+=inc;
}
}
}
}
}
int main(){
scanf("%d",&n);
for(int i=;i<n;i++){
for(int j=;j<;j++){
scanf("%d",&a[i].s[j]);
}
}
for(int i=;i<n;i++){
for(int j=;j<;j++){
scanf("%d",&b[i].s[j]);
}
}
for(int i=;i<n;i++){
for(int j=;j<n;j++){
ch[i][j]=get(a[i],b[j]);
}
}
KM();
int sum=;
for(int i=;i<n;i++){
sum+=ch[i][cx[i]];
}
sum=*n-sum;
printf("%d\n",sum);
}
room 二分图最大匹配KM的更多相关文章
- 带权二分图最大匹配KM算法
二分图的判定 如果一个图是连通的,可以用如下的染色法判定是否二分图: 我们把X部的结点颜色设为0,Y部的颜色设为1. 从某个未染色的结点u开始,做BFS或者DFS .把u染为0,枚举u的儿子v.如果v ...
- UVA1349(带权二分图最大匹配 --> KM算法模板)
UVA1349 题意:给定一些有向带权边,求出把这些边构造成一个个环,总权值最小 解法: 对于带权的二分图的匹配问题可以用通过KM算法求解. 要求最大权匹配就是初始化g[i][j]为0,直接跑就可以: ...
- SPOJ 4206 Fast Maximum Matching (二分图最大匹配 Hopcroft-Carp 算法 模板)
题目大意: 有n1头公牛和n2头母牛,给出公母之间的m对配对关系,求最大匹配数.数据范围: 1 <= n1, n2 <= 50000, m <= 150000 算法讨论: 第一反应 ...
- HDU 2255 奔小康赚大钱(带权二分图最大匹配)
HDU 2255 奔小康赚大钱(带权二分图最大匹配) Description 传说在遥远的地方有一个非常富裕的村落,有一天,村长决定进行制度改革:重新分配房子. 这可是一件大事,关系到人民的住房问题啊 ...
- Luogu 1559 运动员最佳匹配问题(带权二分图最大匹配)
Luogu 1559 运动员最佳匹配问题(带权二分图最大匹配) Description 羽毛球队有男女运动员各n人.给定2 个n×n矩阵P和Q.P[i][j]是男运动员i和女运动员j配对组成混合双打的 ...
- 【算法笔记】二分图与KM算法(当你试图只看蓝书学算法
前言 呜,好久没写博客了,DDL 也有好多,一不留神就轮到我了呜. 看了一眼其它同学写的博客,什么数模啊,什么 CTF 啊,什么 Python 爬虫啊,感觉自己真是越来越菜了呜. 然后在我一愁莫展之际 ...
- POJ 2226二分图最大匹配
匈牙利算法是由匈牙利数学家Edmonds于1965年提出,因而得名.匈牙利算法是基于Hall定理中充分性证明的思想,它是二部图匹配最常见的算法,该算法的核心就是寻找增广路径,它是一种用增广路径求二分图 ...
- POJ2239 Selecting Courses(二分图最大匹配)
题目链接 N节课,每节课在一个星期中的某一节,求最多能选几节课 好吧,想了半天没想出来,最后看了题解是二分图最大匹配,好弱 建图: 每节课 与 时间有一条边 #include <iostream ...
- poj 2239 二分图最大匹配,基础题
1.poj 2239 Selecting Courses 二分图最大匹配问题 2.总结:看到一个题解,直接用三维数组做的,很巧妙,很暴力.. 题意:N种课,给出时间,每种课在星期几的第几节课上 ...
随机推荐
- java重写LinkedList
LinkedList重写类LinkList.java: import java.util.LinkedList;import java.util.List; public class LinkList ...
- 20165205 2017-2018-2 《Java程序设计》实验二 Java面向对象程序设计
20165205 2016-2017-2 <Java程序设计>实验二 Java面向对象程序设计 实验内容 初步掌握单元测试和TDD 理解并掌握面向对象三要素:封装.继承.多态 初步掌握UM ...
- [UGUI]图文混排(七):动态表情
帧动画脚本: http://www.cnblogs.com/lyh916/p/9194823.html 这里的动态表情,我使用的是固定间隔去刷新Image上的Sprite来实现的,即帧动画.这里可以将 ...
- pyqt 8行内就可以跑一个浏览器
pyqt 8行内就可以跑一个浏览器 from PyQt4.QtCore import * from PyQt4.QtGui import * from PyQt4.QtWebKit import * ...
- IntelliJ IDEA 注释模板设置
1.idea类注释 打开:file->setting->Editor->Filr and Code Templates->Includes->File Header 类注 ...
- 【Social Listening实战】当数据分析遭遇心理动力学:用户深层次的情感需求浮出水面
本文转自知乎 作者:苏格兰折耳喵 ----------------------------------------------------- 本文篇幅较长,分为五部分,在中间部分有关于心理分析工具的介 ...
- 《算法》第五章部分程序 part 6
▶ 书中第五章部分程序,包括在加上自己补充的代码,非确定性有穷自动机(NFA),grep 命令(利用 NFA 匹配) ● 非确定性有穷自动机(NFA) package package01; impor ...
- RDD中的cache() persist() checkpoint()
cache只有一个默认的缓存级别MEMORY_ONLY ,而persist可以根据StorageLevel设置其它的缓存级别. cache以及persist都不是action. 被重复使用的(但是)不 ...
- 2. apache整合tomcat部署集群
apache只有处理静态事物的能力, 而tomcat的强项就是处理动态的请求,所以apache和tomcat整合相互取长补短,由apache作为入口,如果是请求静态页面或者是静态文件,由apache直 ...
- log4js_Node.js中的日志管理模块使用
{ "appenders": [ // 下面一行应该是用于跟express配合输出web请求url日志的 {"type": "console" ...