[kuangbin带你飞]专题十 匹配问题 二分图最大权匹配
二分图最大权匹配有km算法和网络流算法
km算法模板默认解决最大权匹配的问题 而使用最小费用最大流 是解决最小权匹配问题
这两种办法都可以求最大最小权 需要两次取反
TAT 感觉讲km会很难的样子...
P hdu2255
km的模板题
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<map>
#include<string>
#include<vector>
#include<queue>
#include<iostream>
using namespace std;
#define L long long
bool visx[305] , visy[305] ;
int sla[305] ;
int linker[305];
int lx[305] , ly[305];
int n ;
int c[305][305]; bool fin(int u ){
visx[u] = false ;
for(int v = 1; v <= n ; v ++ ){
if(visy[v]){
int tmp = lx[u] + ly[v] - c[u][v] ;
if(tmp == 0){
visy[v] = false ;
if(linker[v] == -1 || fin(linker[v])){
linker[v] = u ;
return true ;
}
}
else {
sla[v] = min (sla[v] , tmp ) ;
}
}
}
return false ;
}
int km(){
memset(linker, -1 ,sizeof(linker));
memset(ly , 0 ,sizeof(ly)) ;
for(int i = 1; i <= n ; i ++ ){
lx[i] = 0 ;
for(int j = 1; j <= n; j ++ ){
lx[i] = max(lx[i] , c[i][j]) ;
}
}
for(int i = 1; i<= n ; i ++ ){
for(int j = 1 ; j <= n ; j ++ )sla[j] = 999999999;
while(true){
memset(visx , true , sizeof(visx)) ;
memset(visy , true , sizeof(visy)) ;
if(fin(i))break ;
int d = 999999999 ;
for(int j = 1; j<= n ; j ++ ){
if(visy[j]){
d = min(d , sla[j]) ;
}
}
for(int j = 1; j<= n ; j ++ ){
if(visx[j] == false){
lx[j] -= d;
}
if(visy[j] == false){
ly[j] += d ;
}
else {
sla[j] -= d;
}
}
}
}
int res = 0 ;
for(int j = 1; j <= n ; j ++ ){
if(linker[j] != -1){
res += c[linker[j]][j];
}
}
return res ;
}
int main(){
while(scanf("%d",&n)!=EOF){
for(int i = 1; i <= n ;i ++) {
for(int j = 1; j <= n ; j ++ )
scanf("%d",&c[i][j]) ;
}
int ans = km() ;
printf("%d\n",ans );
}
}
Q hdu3488
是一个要求将所有的点连为一些环 并且边权的加和最小
环的合并可以用最小费用最大流来做
可以转化为 每个点上都连着两条边 使总边和最小 即 求最小权匹配
建图 是左边n个点 右边n个点 那么 最后可以求出一个完备匹配 左边的x点和右边的x点上都有一条边相连
km算法和网络列解法都是可以的 km100+ms 费用流900+ms
km :
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<map>
#include<string>
#include<vector>
#include<queue>
#include<iostream>
using namespace std;
#define L long long
bool visx[305] , visy[305] ;
int sla[305] ;
int linker[305];
int lx[305] , ly[305];
int n , m ;
int c[305][305];
bool fin(int u ){
visx[u] = false ;
for(int v = 1; v <= n ; v ++ ){
if(visy[v]){
int tmp = lx[u] + ly[v] - c[u][v] ;
if(tmp == 0){
visy[v] = false ;
if(linker[v] == -1 || fin(linker[v])){
linker[v] = u ;
return true ;
}
}
else {
sla[v] = min (sla[v] , tmp ) ;
}
}
}
return false ;
}
int km(){
memset(linker, -1 ,sizeof(linker));
memset(ly , 0 ,sizeof(ly)) ;
for(int i = 1; i <= n ; i ++ ){
lx[i] = -999999999 ;
for(int j = 1; j <= n; j ++ ){
lx[i] = max(lx[i] , c[i][j]) ;
}
}
for(int i = 1; i<= n ; i ++ ){
for(int j = 1 ; j <= n ; j ++ )sla[j] = 999999999;
while(true){
memset(visx , true , sizeof(visx)) ;
memset(visy , true , sizeof(visy)) ;
if(fin(i))break ;
int d = 999999999 ;
for(int j = 1; j<= n ; j ++ ){
if(visy[j]){
d = min(d , sla[j]) ;
}
}
for(int j = 1; j<= n ; j ++ ){
if(visx[j] == false){
lx[j] -= d;
}
if(visy[j] == false){
ly[j] += d ;
}
else {
sla[j] -= d;
}
}
}
}
int res = 0 ;
for(int j = 1; j <= n ; j ++ ){
if(linker[j] != -1){
res += c[linker[j]][j];
}
}
return res ;
}
int main(){
int t ;
scanf("%d",&t);
while(t--){
scanf("%d%d",&n,&m) ;
for(int i = 1; i<= n ; i ++ ){
for(int j = 1; j <= n ; j ++)
c[i][j] = - 999999999;
}
for(int i = 1 ; i <= m ; i ++ ){
int u , v, w;
scanf("%d%d%d",&u,&v,&w);
c[u][v] = max(c[u][v] , -w );
}
int ans = km() ;
printf("%d\n", -ans);
}
}
费用流 :
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<map>
#include<string>
#include<vector>
#include<queue>
#include<iostream>
using namespace std;
#define L long long
int n ;
int m ;
int cnt ;
struct ed{
int v , nex ,cap , flow ,cost ;
}e[405 * 405] ;
int head[405];
int tol;
int pre[405];
int dis[405];
bool vis[405];
void init(){
cnt = 0;
memset(head, -1 ,sizeof(head)) ;
}
void add(int u , int v ,int cap ,int cost){
e[cnt].v = v;
e[cnt].nex = head[u];
e[cnt].cap = cap ;
e[cnt].cost = cost ;
e[cnt].flow = 0 ;
head[u] = cnt ;
cnt ++ ;
e[cnt].v = u ;
e[cnt].nex = head[v] ;
e[cnt].cap = 0 ;
e[cnt].cost = -cost ;
e[cnt].flow = 0 ;
head[v] = cnt ;
cnt ++ ;
}
bool spfa(int s, int t){
queue<int >que;
for(int i = 0; i<= 2*n + 1; i ++ ){
dis[i] = 999999999 ;
vis[i] = false ;
pre[i] = -1 ;
}
dis[s] = 0;
vis[s] = true ;
que.push(s) ;
while(!que.empty()){
int u = que.front() ; que.pop () ;
vis[u] = false ;
for(int i = head[u] ; i != -1 ; i = e[i].nex) {
int v = e[i].v ;
if(e[i].cap > e[i].flow && dis[v] > dis[u] + e[i].cost ){
dis[v] = dis[u] + e[i].cost ;
pre[v] = i ;
if(vis[v] == false ){
vis[v] = true ;
que.push(v) ;
}
}
}
}
if(pre[t] == -1)return false ;
else return true ;
}
int fyl(int s , int t , int &cost){
int flow = 0;
cost = 0;
while(spfa(s,t)){
int minn = 999999999;
for(int i = pre[t] ; i != -1 ; i = pre[e[i^1].v]){
if(minn > e[i].cap - e[i].flow){
minn = e[i].cap - e[i].flow ;
}
}
for(int i = pre[t] ; i != -1 ; i = pre[e[i^1].v]){
e[i].flow += minn ;
e[i^1].flow -= minn ;
cost += e[i].cost * minn ;
}
flow += minn ;
}
return flow ;
}
int main(){
int t ;
scanf("%d",&t);
while(t--){
init();
scanf("%d%d",&n,&m);
for(int i = 1; i <= m ; i ++ ){
int u , v , w;
scanf("%d%d%d",&u,&v,&w);
add(u,v+n,1,w);
}
for(int i = 1; i <= n ; i ++ ){
add(0,i,1,0);
add(i+n, n*2 + 1 , 1 , 0) ;
}
int cost ;
int flow = fyl(0,n*2+1,cost);
printf("%d\n",cost);
}
}
[kuangbin带你飞]专题十 匹配问题 二分图最大权匹配的更多相关文章
- [kuangbin带你飞]专题十 匹配问题
A-L 二分匹配 M-O 二分图多重匹配 P-Q 二分图最大权匹配 R-S 一般图匹配带花树 模板请自己找 ID Origin Title 61 / 72 Problem A HD ...
- [kuangbin带你飞]专题十 匹配问题 一般图匹配
过去做的都是二分图匹配 即 同一个集合里的点 互相不联通 但是如果延伸到一般图上去 求一个一般图的最大匹配 就要用带花树来解决 带花树模板 用来处理一个无向图上的最大匹配 看了一会还是不懂 抄了一遍 ...
- [kuangbin带你飞]专题十 匹配问题 二分图多重匹配
二分图的多重匹配问题不同于普通的最大匹配中的"每个点只能有最多一条边" 而是"每个点连接的边数不超过自己的限定数量" 最大匹配所解决的问题一般是"每个 ...
- [kuangbin带你飞]专题十 匹配问题 二分匹配部分
刚回到家 开了二分匹配专题 手握xyl模板 奋力写写写 终于写完了一群模板题 A hdu1045 对这个图进行 行列的重写 给每个位置赋予新的行列 使不能相互打到的位置 拥有不同的行与列 然后左行右列 ...
- [kuangbin带你飞]专题十五 数位DP
ID Origin Title 62 / 175 Problem A CodeForces 55D Beautiful numbers 30 / 84 Problem B HD ...
- [kuangbin带你飞]专题十六 KMP & 扩展KMP & Manacher 题解报告
来刷kuangbin字符串了,字符串处理在ACM中是很重要的,一般比赛都会都1——2道有关字符串处理的题目,而且不会很难的那种,大多数时候都是用到一些KMP的性质或者找规律. 点击标题可跳转至VJ比赛 ...
- [kuangbin带你飞]专题十四 数论基础
ID Origin Title 111 / 423 Problem A LightOJ 1370 Bi-shoe and Phi-shoe 21 / 74 Problem B ...
- 【算法系列学习】DP和滚动数组 [kuangbin带你飞]专题十二 基础DP1 A - Max Sum Plus Plus
A - Max Sum Plus Plus https://vjudge.net/contest/68966#problem/A http://www.cnblogs.com/kuangbin/arc ...
- [kuangbin带你飞]专题十二 基础DP1
ID Origin Title 167 / 465 Problem A HDU 1024 Max Sum Plus Plus 234 / 372 Problem B HDU 1 ...
随机推荐
- HiveServer连接优化
引言 数据平台目前通过Hive SQL的方式提供数据分析服务,系统使用多台HiveServer(JDBCServer)接收客户端连接请求,实际使用场景中频频出现HiveServer内存消耗过多导致 ...
- 【3】JAVA---地址App小软件(AddPanel.class)(表现层)
添加地址信息界面. 年龄和地址必须是数字,否则会弹出窗口提示. 地址信息不能为空. /* * AddPanel.java * * Created on __DATE__, __TIME__ */ pa ...
- 关于Cookie和Session【转载】
当你第一次访问一个网站的时候,网站服务器会在响应头内加上Set-Cookie:PHPSESSID=nj1tvkclp3jh83olcn3191sjq3(php服务器),或Set-Cookie JSES ...
- Sublime Text3 配置markdown插件
sublime是一个亮骚的文本编辑器,而且是跨三大平台,而markdown是一门标记语法,对于记录真是神器,具体语法百度很多,下面教你在sublime上配置markdown. 这两个神器结合起来简直好 ...
- Android SimpleAdapter GridView (网格图片点击放大显示)
GridView网格视图 GridView网格视图是按照行,列分布的方式来显示多个组件,通常用于显示图片或是图标等,在使用网格视图时,首先需要要在屏幕上添加GridView组件. 常用属性: 1. a ...
- Eclipse自动插件依赖的一种配置解决方式
Eclipse的插件具有以下特点: (1)每一个插件有自己独立的classloader (2)插件资源的交互通过MENIFEST.MF中"Export-Package, Require-Bu ...
- MySQL主从复制(Master-Slave)与读写分离(MySQL-Proxy)
MySQL主从复制(Master-Slave)即同步与读写分离(MySQL-Proxy)即集群
- 怎样在 SSASserver之间做同步
简单介绍: 从SQL Server 2005開始,分析服务就支持了同步的功能.本文将介绍怎样在SQL Server 2012下同步Adventureworks的分析服务数据库.通过同步的功能,我 ...
- 选课 树形dp+路径输出
#include<iostream> #include<cstdio> #include<cstring> #define maxn 2010 using name ...
- spring源码分析
编译问题 spring-4.0.5.release编译是用jdk8编译的,为啥可以运行在jdk7的环境? 源码分析 spring源码分析,由一个点各个击破,比如依赖注入,autowired. spri ...