【问题描述】
卡普地公司举办了「世界第一的猛汉王」全球大会，来自世界各地的猛汉为了争夺「猛汉王」的名号前来一决高下。现在举行的是弓箭组选拔赛。卡普地公司为比赛新建了一张PVP地图——「猛汉竞技场」。有许多使用弓的猛汉在这里互相较量。他们中的一些装填了「接击瓶」，这使得他们在接近战中会占有一定优势，但是在远程战中会相当劣势。具体来说如下：
假设q装填了「接击瓶」而p没有，则当他们的曼哈顿距离大于D时，p压制q，反之q压制p。如果p和q都装填了「接击瓶」或者都没有，则他们之间仍然会存在一个客观上的单向压制关系，但是在比赛刚开始时无法得知。
竞技场上一共有n+m个猛汉，其中n个装填了「接击瓶」，另外m个没有。每个猛汉降临到竞技场时有一个坐标(x, y)。Mark Douglas作为上一届的猛汉王正在观看这场比赛，他希望得知场上有多少个「猛汉三角」。「猛汉三角」是指三个人u、v、w满足u压制v，v压制w，w压制u，且三人中至少有一人装填了「接击瓶」且至少有一人没有。由于场上尚存在一些不明了的压制关系，所以Mark希望知道可能的「猛汉三角」数量的最小值和最大值。
【输入格式】
输入文件名为mhw.in。
输入第一行为三个正整数n m D。
接下来n行每行两个正整数x y，表示装填了「接击瓶」的猛汉的坐标。
接下来m行每行两个正整数x y，表示没有装填「接击瓶」的猛汉的坐标。
可能有多个猛汉站在同一个位置。
【输出格式】
输出文件名为mhw.out。
输出两个数min max，表示答案的最小值和最大值。
【样例输入与输出】
example_mhw1.in | example_mhw1.out |
2 2 1 1 2 1 1 3 1 2 2 | 0 2 |
难度确实大！
首先我们将猛汉们分为两个阵营（记为黑色与白色）。我们一个猛汉三角就是形如a->c->b（a，b同阵营），然后在确定a,b之间的压制关系得到最小或者最大值。
我们以白色阵营为例。我们设cover[v]为v压制的对方阵营的数量，cover[v,u]为v,u共同压制的对方阵营的数量。maxans同理。
。对于cover[v,u]，我们就算另一方阵营对
的贡献就行了。如果一个猛汉被对面的k个人压制，那么他对
的贡献就是
。
求cover数组可以用扫描+线段树。
代码：
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<ctime>
#include<queue>
#define ll long long
#define N 100005 using namespace std;
inline ll Get() {ll x=,f=;char ch=getchar();while(ch<''||ch>'') {if(ch=='-') f=-;ch=getchar();}while(''<=ch&&ch<='') {x=(x<<)+(x<<)+ch-'';ch=getchar();}return x*f;} ll n,m,lim;
ll mn,mx;
ll cc;
ll d[N*];
ll cnt;
ll c[][N]; struct line {
ll x,op,col,id;
ll l,r;
bool operator <(const line &a)const {
if(x!=a.x) return x<a.x;
return op<a.op;
}
}st[N*]; struct tree {
ll l,r;
ll sum[];
}tr[N*]; void build(ll v,ll l,ll r) {
tr[v].l=l,tr[v].r=r;
if(l==r) return ;
ll mid=l+r>>;
build(v<<,l,mid),build(v<<|,mid+,r);
} void Modify(ll v,ll l,ll r,ll tag,ll f) {
if(tr[v].l>r||tr[v].r<l) return ;
if(l<=tr[v].l&&tr[v].r<=r) {
tr[v].sum[tag]+=f;
return ;
}
Modify(v<<,l,r,tag,f),Modify(v<<|,l,r,tag,f);
} ll query(ll v,ll pos,ll tag) {
if(tr[v].l>pos||tr[v].r<pos) return ;
if(tr[v].l==tr[v].r) return tr[v].sum[tag];
return query(v<<,pos,tag)+query(v<<|,pos,tag)+tr[v].sum[tag];
} void work() {
for(ll i=;i<=cnt;i++) {
if(st[i].op==-) {
Modify(,st[i].l,st[i].r,st[i].col,);
} else if(st[i].op==) {
c[st[i].col][st[i].id]=query(,st[i].l,st[i].col^);
} else {
Modify(,st[i].l,st[i].r,st[i].col,-);
}
}
} int main() {
n=Get(),m=Get(),lim=Get();
ll a,b,x,y;
for(ll i=;i<=n;i++) {
a=Get(),b=Get();
x=a+b,y=a-b;
st[++cnt]=(line) {x-lim,-,,i,y-lim,y+lim};
st[++cnt]=(line) {x,,,i,y};
st[++cnt]=(line) {x+lim,,,i,y-lim,y+lim};
d[++d[]]=y;
d[++d[]]=y+lim;
d[++d[]]=y-lim;
}
for(ll i=;i<=m;i++) {
a=Get(),b=Get();
x=a+b,y=a-b;
st[++cnt]=(line) {x-lim,-,,i,y-lim,y+lim};
st[++cnt]=(line) {x,,,i,y};
st[++cnt]=(line) {x+lim,,,i,y-lim,y+lim};
d[++d[]]=y;
d[++d[]]=y+lim;
d[++d[]]=y-lim;
}
sort(st+,st++cnt);
sort(d+,d++d[]);
cc=unique(d+,d++d[])-d;
for(ll i=;i<=cnt;i++) {
st[i].l=lower_bound(d+,d+cc,st[i].l)-d;
st[i].r=lower_bound(d+,d+cc,st[i].r)-d;
}
build(,,cc); work();
for(ll i=;i<=m;i++) c[][i]=n-c[][i];
sort(c[]+,c[]++n);
sort(c[]+,c[]++m); for(ll i=;i<=n;i++) {
mn+=(n-i)*c[][i];
mx+=(i-)*c[][i];
ll x=1ll*(m-c[][i])*(m-c[][i]-)/;
mn-=x,mx-=x;
}
for(ll i=;i<=m;i++) {
mn+=1ll*(m-i)*c[][i];
mx+=1ll*(i-)*c[][i];
ll x=1ll*(n-c[][i])*(n-c[][i]-)/;
mn-=x,mx-=x;
}
cout<<mn<<" "<<mx;
return ;
}