@description@

考虑一个二维平面,执行共 2*N 次操作:

前 N 次,第 i 次在 (rx, ry) 处放置 rc 个红色球;

后 N 次,第 i 次在 (bx, by) 处放置 bc 个蓝色球。

保证放置的红色球总数 = 放置的蓝色球总数。

请将这些球两两配对,使得所有配对中 (bx, by) 与 (rx, ry) 的 |rx - bx| + |ry - by| 之和最大。

Constraints

1≤N≤1000, 0≤RXi,RYi,BXi,BYi≤10^9, 1≤RCi,BCi≤10

∑RCi=∑BCi

Input

输入形式如下:

N

RX1 RY1 RC1

RX2 RY2 RC2



RXN RYN RCN

BX1 BY1 BC1

BX2 BY2 BC2



BXN BYN BCN

Output

输出配对后曼哈顿距离之和的最大值。

Sample Input 1

2

0 0 1

3 2 1

2 2 1

5 0 1

Sample Output 1

8

Sample Input 2

Copy

3

0 0 1

2 2 1

0 0 2

1 1 1

1 1 1

3 3 2

Sample Output 2

16

@solution@

假如不考虑数据范围,我们可以将所有红球与所有蓝球连边,边权为匹配的曼哈顿距离。

然后跑一个二分图的最大权完美匹配。

这看起来非常好,但是 N <= 1000,边建出来一共有 N^2 条,会炸。

我们考虑怎么才能优化建图。

注意我们为什么要建 N^2 条边:哈密顿距离中带有绝对值,使得边权由两点共同决定。

在最大化问题中,有一个去掉绝对值的常用方法(当然不是零点分段):|x| = max{x, -x}。

这道题同理,|rx - bx| + |ry - by| 可以表示成以下几种情况的最大值:

(rx - bx) + (ry - by) = (+ rx + ry) + (- bx - by)

(rx - bx) + (by - ry) = (+ rx - ry) + (- bx + by)

(bx - rx) + (ry - by) = (- rx + ry) + (+ bx - by)

(bx - rx) + (by - ry) = (- rx - ry) + (+ bx + by)

这样的转化使得边权可以由两点分别的权之和决定,就可以拆开,不用再建 O(n^2) 条边了。

具体来说,我们先 s 向 N 个红球连容量为 rc,费用为 0 的边;再 N 个蓝球向 t 连容量为 bc,费用为 0 的边。

中间另建 4 个点,表示以上 4 种情况。N 个红球连过去 4 个点,4 个点连向 N 个蓝球,费用如上面所谈论的。

注意 N 个红球与 N 个蓝球要交错着连(正正 对应 负负)。

@accepted code@

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MAXV = 4000;
const int MAXE = 20000;
const int INF = (1<<30);
const ll LINF = (1LL<<60);
struct FlowGraph{
struct edge{
int to, cap, flow; ll dis;
edge *nxt, *rev;
}edges[MAXE + 5], *adj[MAXV + 5], *cur[MAXV + 5], *ecnt;
FlowGraph() {ecnt = &edges[0];}
int s, t;
void addedge(int u, int v, int c, ll w) {
edge *p = (++ecnt), *q = (++ecnt);
p->to = v, p->cap = c, p->flow = 0, p->dis = -w;
p->nxt = adj[u], adj[u] = p;
q->to = u, q->cap = 0, q->flow = 0, q->dis = w;
q->nxt = adj[v], adj[v] = q;
p->rev = q, q->rev = p;
// printf("%d %d %d %lld\n", u, v, c, w);
}
ll h[MAXV + 5], d[MAXV + 5], f[MAXV + 5];
int hp[MAXV + 5];
void update(int x, ll k) {
f[x] = k;
while( x ) {
hp[x] = x;
if( (x<<1) <= t && f[hp[x]] > f[hp[x<<1]] )
hp[x] = hp[x<<1];
if( (x<<1|1) <= t && f[hp[x]] > f[hp[x<<1|1]] )
hp[x] = hp[x<<1|1];
x >>= 1;
}
}
bool relabel() {
for(int i=s;i<=t;i++)
h[i] += d[i], d[i] = f[i] = LINF, hp[i] = i, cur[i] = adj[i];
update(s, d[s] = 0);
while( f[hp[1]] != LINF ) {
int x = hp[1]; update(x, LINF);
for(edge *p=adj[x];p;p=p->nxt) {
if( p->cap > p->flow ) {
ll w = p->dis + h[x] - h[p->to];
if( d[x] + w < d[p->to] )
update(p->to, d[p->to] = d[x] + w);
}
}
}
return !(d[t] == LINF);
}
bool vis[MAXV + 5];
int aug(int x, int tot) {
if( x == t ) return tot;
int sum = 0; vis[x] = true;
for(edge *&p=cur[x];p;p=p->nxt) {
ll w = p->dis + h[x] - h[p->to];
if( !vis[p->to] && p->cap > p->flow && d[x] + w == d[p->to] ) {
int del = aug(p->to, min(tot - sum, p->cap - p->flow));
p->flow += del, p->rev->flow -= del, sum += del;
if( sum == tot ) break;
}
}
vis[x] = false;
return sum;
}
ll min_cost_max_flow(int _s, int _t) {
s = _s, t = _t; ll cost = 0;
while( relabel() ) {
int del = aug(s, INF);
cost += (d[t] + h[t]) * del;
}
return -cost;
}
}G;
int N, s, t;
int main() {
scanf("%d", &N), s = 1, t = 2*N+6;
for(int i=1;i<=N;i++) {
int x, y, c; scanf("%d%d%d", &x, &y, &c);
G.addedge(s, i+1, c, 0);
G.addedge(i+1, 2*N+2, INF, + x + y);
G.addedge(i+1, 2*N+3, INF, + x - y);
G.addedge(i+1, 2*N+4, INF, - x + y);
G.addedge(i+1, 2*N+5, INF, - x - y);
}
for(int i=1;i<=N;i++) {
int x, y, c; scanf("%d%d%d", &x, &y, &c);
G.addedge(i+N+1, t, c, 0);
G.addedge(2*N+5, i+N+1, INF, + x + y);
G.addedge(2*N+4, i+N+1, INF, + x - y);
G.addedge(2*N+3, i+N+1, INF, - x + y);
G.addedge(2*N+2, i+N+1, INF, - x - y);
}
printf("%lld\n", G.min_cost_max_flow(s, t));
}
/*
(+ rx + ry) + (- bx - by)
(+ rx - ry) + (- bx + by)
(- rx + ry) + (+ bx - by)
(- rx - ry) + (+ bx + by)
*/

@details@

为什么我会傻到以为图匹配是 NP 问题。。。

特别是这个图还是个二分图。。。

想了半天的贪心,然后不停地叉掉,然后继续贪。。。

@atcoder - AGC034D@ Manhattan Max Matching的更多相关文章

  1. 【杂题】[AGC034D] Manhattan Max Matching【费用流】

    Description 有一个无限大的平面,有2N个位置上面有若干个球(可能重复),其中N个位置是红球,N个位置是蓝球,红球与蓝球的总数均为S. 给出2N个位置和上面的球数,现要将红球与蓝球完美匹配, ...

  2. [AGC034D]Manhattan Max Matching:费用流

    前置姿势 \(k\)维空间内两点曼哈顿距离中绝对值的处理 戳这里:[CF1093G]Multidimensional Queries 多路增广的费用流 据说这个东西叫做ZKW费用流? 流程其实很简单, ...

  3. 「AGC034D」 Manhattan Max Matching

    「AGC034D」 Manhattan Max Matching 传送门 不知道这个结论啊... (其实就是菜嘛) 首先 \(O(n^2)\) 的建边显然不太行. 曼哈顿距离有这样一个性质,如果将绝对 ...

  4. AtCoder ABC 126F XOR Matching

    题目链接:https://atcoder.jp/contests/abc126/tasks/abc126_f 题目大意 给定两个整数 M 和 K ,用小于 2M 的的所有自然数,每个两个,用这些数排成 ...

  5. 【AtCoder】AGC034

    AGC034 刷了那么久AtCoder我发现自己还是只会ABCE(手动再见 A - Kenken Race 大意是一个横列,每个点可以跳一步或者跳两步,每个格子是空地或者石头,要求每一步不能走到石头或 ...

  6. [2019多校联考(Round 6 T3)]脱单计划 (费用流)

    [2019多校联考(Round 6 T3)]脱单计划 (费用流) 题面 你是一家相亲机构的策划总监,在一次相亲活动中,有 n 个小区的若干男士和 n个小区的若干女士报名了这次活动,你需要将这些参与者两 ...

  7. CSS 3学习——transition 过渡

    以下内容根据官方规范翻译以及自己的理解整理. 1.介绍 这篇文档介绍能够实现隐式过渡的CSS新特性.文档中介绍的CSS新特性描述了CSS属性的值如何在给定的时间内平滑地从一个值变为另一个值. 2.过渡 ...

  8. 《算法》第六章部分程序 part 6

    ▶ 书中第六章部分程序,包括在加上自己补充的代码,包括二分图最大匹配(最小顶点覆盖)的交替路径算法和 HopcroftKarp 算法 ● 二分图最大匹配(最小顶点覆盖)的交替路径算法 package ...

  9. AtCoder Beginner Contest 247 E - Max Min // 容斥原理

    原题链接:E - Max Min (atcoder.jp) 题意: 给定一个数组,求满足最大值为X且最小值为Y的区间个数. 思路:容斥原理 因为必须要包含端点,直接求是不容易的.因此考虑去求不一定包含 ...

随机推荐

  1. php连接数据库查询方法(还少一种pdo方法)

    <?php header("content-type:text/html;charset=utf-8"); $conn = mysql_connect("local ...

  2. HTTP协议②缓存

    首先介绍一下HTTP的报文信息,主要分成两个部分 1.包含属性的首部(header)---------------附加信息(cookie,缓存信息等)与缓存相关的规则信息,均包含在header中 2. ...

  3. Java的split()方法

    在java中,split方法很常用,但是有些时候我们发现并不管用,不管用的情况:split("."),split("^"), 管用的情况:split(" ...

  4. Win7下设置WiFi热点

    Win7下设置WiFi热点 今天研究了下Win7设置WIFI热点,Connectify神马的都是浮云~亲測可用,现拿出来分享下~ 1.点击"開始",再点击"执行" ...

  5. 2017-09-23-ArchData

    layout: post title: ArchData category: Technical tags: [分布式,区块链,AI,大数据] ArchData 技术峰会 神经网络和函数式编程 杨博: ...

  6. OSGi教程:Class Space Consistency

    此教程基于OSGi Core Release 7 OSGi类空间的一致性 详细内容上面英文教程有详细解答 下面主要是一些个人见解,若有不当之处,欢迎指出: "Class space cons ...

  7. 荷畔微风 - 在函数计算FunctionCompute中使用WebAssembly

    WebAssembly 是一种新的W3C规范,无需插件可以在所有现代浏览器中实现近乎原生代码的性能.同时由于 WebAssembly 运行在轻量级的沙箱虚拟机上,在安全.可移植性上比原生进程更加具备优 ...

  8. 模拟19 题解(waiting)

    T1,千万别转化成链了!! 直接数就可以,dfs搜索每种情况,对于搜到的点,如果子树大小过大,直接return,相等说明可以,小的话向上累加, 优化是先预处理子树大小,若子树小,不用搜了直接加上就行 ...

  9. DataIntputStream / DataOutputStream 类

    1. DataInputStream类(熟悉)   (1)基本概念 java.io.DataInputStream类用于读取java中的基本数据类型.   (2)常用的方法 DataInputStre ...

  10. I / O流 类

    一.概述 1 基本概念 I/O就是Input/Output的简写,也就是输入/输出的含义. I/O流就是指像流水一样源源不断地进行读写的过程.   2 基本分类   根据读写数据的单元分为:字节流 和 ...