emmmmmmmmmmmmmm
首先可以发现,对于每个点,只包含y值大于它的矩形。故用y值从大到小排序u。
然后到了重头戏:计算。我的第一次想法是:边插入边计算。可是,我们可以找到反例:
1 2
3 4
按顺序计算4点的时候,我们发现它同排的3点已经算过(1,2,3,4)了!
仔细思考一下,将算法改进为同排一次插入,再只算与前面一位的点的差*后面的空。可以借助线段树或树状数组实现。
#include<cstdio>
#include<algorithm>
using namespace std;
const int N = 200002;
struct node {
int x, y;
}s[N];
int n, tmp, cnt;
long long ans, t[N << 2];
bool cmp(const node a, const node b) {
if(a.y < b.y || (a.y == b.y && a.x > b.x))
return 0;
return 1;
}
bool cmp1(const node a, const node b) {
return a.x < b.x;
}
bool cmp2(const node a, const node b) {
return a.y < b.y;
}
void add(const int l, const int r, const int id, const int num) {
if(l == r) {
t[id] = 1;
return;
}
int mid = l + r >> 1;
if(num <= mid)
add(l, mid, id << 1, num);
else
add(mid + 1, r, id << 1 | 1, num);
t[id] = t[id << 1] + t[id << 1 | 1];
}
long long ask(const int l, const int r, const int id, const int L, const int R) {
if(L > r || R < l)
return 0;
if(L <= l && R >= r)
return t[id];
int mid = l + r >> 1;
return ask(l, mid, id << 1, L, R) + ask(mid + 1, r, id << 1 | 1, L, R);
}
int main() {
scanf("%d", &n);
for(int i = 1; i <= n; i ++) {
scanf("%d %d", &s[i].x, &s[i].y);
}
sort(s + 1, s + n + 1, cmp1);
for(int i = 1; i <= n; i ++)
if(tmp == s[i].x)
s[i].x = cnt;
else {
tmp = s[i].x;
s[i].x = ++ cnt;
}
tmp = cnt = 0;
sort(s + 1, s + n + 1, cmp2);
for(int i = 1; i <= n; i ++)
if(tmp == s[i].y)
s[i].y = cnt;
else {
tmp = s[i].y;
s[i].y = ++ cnt;
};
sort(s + 1, s + n + 1, cmp);
for(int i = 1; i <= n; i ++) {
tmp = s[i].y;
int j = i;
while(tmp == s[j].y) {
add(1, n, 1, s[j].x);
j ++;
}
ans += ask(1, n, 1, 1, s[i].x) * ask(1, n, 1, s[i].x, n);
i ++;
while(tmp == s[i].y) {
ans += ask(1, n, 1, s[i - 1].x + 1, s[i].x) * ask(1, n, 1, s[i].x, n);
i ++;
}
i --;
}
printf("%lld\n", ans);
return 0;
}
如果有错误的话,请大佬指出。