題目鏈接

https://atcoder.jp/contests/agc034/tasks/agc034_d

題解

降智了沒(méi)想出來(lái)建圖……(不過(guò)這場(chǎng)是真的DE放反了)
注意到曼哈頓距離有一個(gè)重要的性質(zhì): 合法的最大，也就是對(duì)于任意維數(shù)兩個(gè)點(diǎn)的曼哈頓距離而言，假設(shè)我們把每個(gè)絕對(duì)值符號(hào)任意地改變成正號(hào)或者負(fù)號(hào)，那么合法的(實(shí)際的)曼哈頓距離是這所有值里最大的。證明顯然。
那么既然我們求的是最大距離和，就可以忽略合不合法的問(wèn)題，四種拆絕對(duì)值的方法不對(duì)連接的兩個(gè)點(diǎn)的坐標(biāo)有要求。
考慮費(fèi)用流建圖，建兩排點(diǎn)分別表示兩種球，然后在中間建四個(gè)點(diǎn)代表四種匹配方式，兩邊的每個(gè)點(diǎn)和中間的每個(gè)點(diǎn)連對(duì)應(yīng)匹配方式的權(quán)值，求最大費(fèi)用最大流即可。
時(shí)間復(fù)雜度\(O(MFMC(n,10n))\), 似乎費(fèi)用流算法的復(fù)雜度可以被估計(jì)為\(O(mC\log m)\)其中\(C\)為最大流量，而這里最大流量是\(O(S)\)的，故復(fù)雜度\(O(nS\log n)\).
題解里說(shuō)可以模擬費(fèi)用流做到\(O(S\log n)\)... 瑟瑟發(fā)抖

代碼

#include<bits/stdc++.h> #define llong long long #define mkpr make_pair #define riterator reverse_iterator using namespace std;inline int read() {int x = 0,f = 1; char ch = getchar();for(;!isdigit(ch);ch=getchar()) {if(ch=='-') f = -1;}for(; isdigit(ch);ch=getchar()) {x = x*10+ch-48;}return x*f; }const llong INF = 1e12;namespace NetFlow {const int N = 2006;const int M = 10000;struct AEdge{int u,v,wl,wr; llong c;} ae[M+3];struct Edge{int u,v,nxt,w; llong c;} e[(M<<1)+3];int fe[N+3];llong dis[N+3];int que[N+5];bool inq[N+3];int lst[N+3];int n,m,en,s,t; llong mf,mc;void addedge(int u,int v,int w,llong c){ // printf("addedge %d %d %d %lld\n",u,v,w,c);en++; e[en].u = u,e[en].v = v,e[en].w = w,e[en].c = c;e[en].nxt = fe[u]; fe[u] = en;en++; e[en].u = v,e[en].v = u,e[en].w = 0,e[en].c = -c;e[en].nxt = fe[v]; fe[v] = en;}bool spfa(){for(int i=1; i<=n; i++) dis[i] = -INF;int hd = 1,tl = 2; que[1] = s; dis[1] = 0;while(hd!=tl){int u = que[hd]; hd++; if(hd>n+1) hd-=n+1;for(int i=fe[u]; i; i=e[i].nxt){int v = e[i].v;if(e[i].w>0&&dis[e[i].v]<dis[u]+e[i].c){dis[e[i].v] = dis[u]+e[i].c; lst[e[i].v] = i;if(!inq[e[i].v]){inq[e[i].v] = true;que[tl] = e[i].v; tl++; if(tl>n+1) tl-=n+1;}}}inq[u] = false;}return dis[t]!=-INF;}void calcflow(){int flow = 1e5;for(int u=t; u!=s; u=e[lst[u]].u){flow = min(flow,e[lst[u]].w);}for(int u=t; u!=s; u=e[lst[u]].u){e[lst[u]].w -= flow; e[lst[u]^1].w += flow;}mf += flow; mc += 1ll*flow*dis[t];}llong mfmc(int _n,int _s,int _t){n = _n,s = _s,t = _t; mf = 0,mc = 0ll;while(spfa()) {calcflow();} return mc;} } using NetFlow::addedge;const int N = 1000; struct Point {int x,y,cnt; } a[N+3],b[N+3]; int n;int main() {NetFlow::en = 1;scanf("%d",&n);for(int i=1; i<=n; i++) scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].cnt);for(int i=1; i<=n; i++) scanf("%d%d%d",&b[i].x,&b[i].y,&b[i].cnt);for(int i=1; i<=n; i++){addedge(1,i+6,a[i].cnt,0);addedge(i+6,3,a[i].cnt,-a[i].x-a[i].y);addedge(i+6,4,a[i].cnt,-a[i].x+a[i].y);addedge(i+6,5,a[i].cnt,a[i].x-a[i].y);addedge(i+6,6,a[i].cnt,a[i].x+a[i].y);}for(int i=1; i<=n; i++){addedge(i+n+6,2,b[i].cnt,0);addedge(3,i+n+6,b[i].cnt,b[i].x+b[i].y);addedge(4,i+n+6,b[i].cnt,b[i].x-b[i].y);addedge(5,i+n+6,b[i].cnt,-b[i].x+b[i].y);addedge(6,i+n+6,b[i].cnt,-b[i].x-b[i].y);}llong ans = NetFlow::mfmc(n+n+6,1,2);printf("%lld\n",ans);return 0; }

總結(jié)

以上是生活随笔為你收集整理的AtCoder AGC034D Manhattan Max Matching (费用流)的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。