Guessing Camels 【Gym - 101485G】【CDQ+樹狀數組求解三維偏序】

題目鏈接

題面


  有三個人在賭駱駝,駱駝的標號是1~N,但是這三人無聊,他們想知道他們賭的駱駝中有幾對駱駝是相對位置相同的,什麼是相對位置相同?就是譬如駱駝A和駱駝B,在三個人的紙條上都是A在B的前面,那麼就代表了相對位置是相同的了。

  所以,我們可以記錄每個駱駝的三個座標,然後求一個最長上生子序列,用三維偏序的方法,去求解答案。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
#include <unordered_map>
#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define Big_INF 0x3f3f3f3f3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
#define MP(a, b) make_pair(a, b)
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 2e5 + 7;
int N;
struct node
{
    int a[3]; bool op;
    friend bool operator < (node e1, node e2) { return e1.a[0] < e2.a[0]; }
} t[maxN], cop[maxN];
ll ans, tree[maxN];
void update(int x, ll val) { while(x <= N) { tree[x] += val; x += lowbit(x); } }
ll query(int x) { ll sum = 0; while(x) { sum += tree[x]; x -= lowbit(x); } return sum; }
void cdq(int l, int r)
{
    if(l == r) return;
    int mid = HalF;
    cdq(l, mid); cdq(mid + 1, r);
    int ql = l, qr = mid + 1, pos = l;
    while(ql <= mid || qr <= r)
    {
        if(qr > r || (ql <= mid && t[ql].a[1] < t[qr].a[1]))
        {
            cop[pos++] = t[ql];
            update(t[ql].a[2], 1);
            ql++;
        }
        else
        {
            cop[pos++] = t[qr];
            ans += query(t[qr].a[2]);
            qr++;
        }
    }
    for(int i=l; i<=mid; i++) update(t[i].a[2], -1);
    for(int i=l; i<=r; i++) t[i] = cop[i];
}
int main()
{
    scanf("%d", &N);
    for(int i=0, v; i<3; i++)
    {
        for(int j=1; j<=N; j++)
        {
            scanf("%d", &v);
            t[v].a[i] = j;
        }
    }
    sort(t + 1, t + N + 1);
    cdq(1, N);
    printf("%lld\n", ans);
    return 0;
}

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章