NC13947 contest(CDQ分治+树状数组+三维偏序)

题目链接

题意：
$n支队伍，每支参加了三场比赛$
$如果三场比赛里至少一场x队比y队排名高$
$那么x队自认为比y强$
$问有多少组(x,y)$
$x自认为比y强，y自认为比x强$
$（x，y）和（y，x）算同一组$
题解：
$首先对第一场排名进行排序，保证i+1自认为比i强$
$然后就需要找有多少队j(1<=j<i)后两场比i强$
$这样明显一看是一个n^2算法$
$但是n<=2e5，肯定是不允许n^2$

$那换一个想法，反着写$
$总共有n\ast (n-1)/2种组合，有多少种组合不可能$
$减去之后就是答案$
$不可能的组一定是其中一组三场比赛都比另一组排名高$
$所以想找到不可能的组合$
$就是找多少队j(1<=j<i)三场比赛比i排名都低$
$这样问题就很明显了，是一个裸的三维偏序问题$
$三个维分别代表三场比赛$
$三维偏序问题需要用cdq处理第二维，BIT处理第三维$
$首先对第一维进行从小到大排序，保证了第一维的有序$
$然后就行cdq分治，同时用BIT进行维护左区间第三维对应位置的个数$
$用两个指针分别看左右区间的当前位$
$由于保证了左区间的第一维一定比右区间的第一维小$
$如果左区间的当前位第二维小，把这一位的第三维加入BIT$
$否则可以直接通过BIT统计当前有多少数比右区间这一位的第三维小$
$并加入答案，然后把归并的顺序赋给原数组，往上传递$
$(讲的不是很清楚，可以到网上查阅CDQ分治或者三维偏序的讲解，那里会有更具体的解法)$

AC代码

/*
    Author:zzugzx
    Lang:C++
    Blog:blog.csdn.net/qq_43756519
*/
#include
using namespace std;
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define endl '\n'
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
//const int mod=1e9+7;
const int mod=998244353;
const double eps = 1e-10;
const double pi=acos(-1.0);
const int maxn=1e6+10;
const ll inf=0x3f3f3f3f;
const int dir[4][2]={{0,1},{1,0},{0,-1},{-1,0}};

ll n,ans,c[maxn];
struct node{
    ll a,b,c;
}p[maxn],p1[maxn];
bool cmp(node a,node b){
    return a.a<b.a;
}
int lowbit(int x){
    return x&(-x);
}
void add(int x,ll v){
    while(x<maxn){
        c[x]+=v;
        x+=lowbit(x);
    }
}
ll query(int x){
    ll res=0;
    while(x){
        res+=c[x];
        x-=lowbit(x);
    }
    return res;
}
void cdq(int l,int r){
    if(l==r)return;
    int mid=l+r>>1;
    cdq(l,mid);cdq(mid+1,r);
    int i=l,j=mid+1,k=l;
    while(i<=mid&&j<=r){
        if(p[i].b<p[j].b)add(p[i].c,1),p1[k++]=p[i++];
        else ans+=query(p[j].c),p1[k++]=p[j++];
    }
    while(i<=mid)add(p[i].c,1),p1[k++]=p[i++];
    while(j<=r)ans+=query(p[j].c),p1[k++]=p[j++];
    for(int i=l;i<=mid;i++)add(p[i].c,-1);
    for(int i=l;i<=r;i++)p[i]=p1[i];
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);cout.tie(0);
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    cin>>n;
    for(int i=1;i<=n;i++)cin>>p[i].a>>p[i].b>>p[i].c;
    sort(p+1,p+1+n,cmp);
    cdq(1,n);
    cout<<n*(n-1)/2-ans;
    return 0;
}