三维偏序（陌上花开）

思路：

首先我们可以以x作为第一关键字进行排序，（y和z作为次要关键字），这样我们就可以保证从左到右是以x单调递增的。

然后我们将这个区间一分为二，分别对于每一个序列以y作为第一关键字排序，这样一来我们得到的两个序列拥有以下的性质：

左序列的任意x值都小于右序列的任意x值

每一个序列里的y值都是单调递增的

最后我们分别在序列头和序列中点设定两个指针j和i，对于每一个i，j都从上一个不满足条件的地方开始枚举，将j对应的y小于i对应的y的j加入树状数组中，查询z值满足条件的数量。

注意去重

#include <bits/stdc++.h>
using namespace std;

#define ull unsigned long long
#define ll long long
#define edl '\n'

const int N = 1e6 + 10;
const int M = 1e3 + 10;

const int mod = 1e9 + 7;
ll n, m;
string s;

struct QWQ
{
	ll x, y, z;
	ll sum,ans;
} a[N];

bool cmp(QWQ A, QWQ B)
{
	if(A.x == B.x)
	{
		if(A.y == B.y) return A.z < B.z;
		return A.y < B.y;
	}
	return A.x < B.x;
}
bool cmp1(QWQ A,QWQ B)
{
//	if(A.y == B.y) return A.z < B.z;
	return A.y < B.y;
}

ll AN[N],cnt[N];

template <typename T>
struct Fenwick {
	int n;
	std::vector<T> a;
	
	Fenwick(int n_ = 0) {
		init(n_);
	}
	
	void init(int n_) {
		n = n_;
		a.assign(n, T{});
	}
	
	void add(int x, const T &v) {
		for (int i = x + 1; i <= n; i += i & -i) {
			a[i - 1] = a[i - 1] + v;
		}
	}
	//[0,x)
	T sum(int x) {
		T ans{};
		for (int i = x; i > 0; i -= i & -i) {
			ans = ans + a[i - 1];
		}
		return ans;
	}
	//左闭右开
	T rangeSum(int l, int r) {
		return sum(r) - sum(l);
	}
	
	//return 第一个sum(x) > k的x
	int select(const T &k) {
		int x = 0;
		T cur{};
		for (int i = 1 << __lg(n); i; i /= 2) {
			if (x + i <= n && cur + a[x + i - 1] <= k) {
				x += i;
				cur = cur + a[x - 1];
			}
		}
		return x;
	}
};

Fenwick<int>f(N);

void cdq(int l,int r)
{
	if(l == r) return;
	ll mid = (l + r) >> 1;
	cdq(l,mid);
	cdq(mid + 1,r);
	sort(a + l,a + mid + 1,cmp1);
	sort(a + mid + 1,a + r + 1,cmp1);
	int i = l;
	int j = mid + 1;
	while(j <= r)
	{
		while(a[i].y <= a[j].y && i <= mid)
		{
			f.add(a[i].z,a[i].sum);
			i ++;
		}
		a[j].ans += f.sum(a[j].z + 1);
		j ++;
	}
	for(j = l; j < i; j ++ ) f.add(a[j].z,-a[j].sum);
}

void solve()
{
	cin >> n >> m;
	for (int i = 1; i <= n; i ++ ) cin >> a[i].x >> a[i].y >> a[i].z;
	sort(a + 1, a + n + 1, cmp);
	int nn = 0;
	for (int i = 1; i <= n; i ++ )
	{
		if (a[i].x == a[nn].x && a[i].y == a[nn].y && a[i].z == a[nn].z) a[nn].sum ++;
		else
		{
			a[++nn] = a[i];
			a[nn].sum = 1;
		}
	}
	cdq(1,nn);
	for(int i = 1; i <= nn; i ++ ) cnt[a[i].ans + a[i].sum - 1] += a[i].sum;
	for(int i = 0; i < n; i ++ ) cout << cnt[i] << '\n';
}

int main()
{
	std::ios::sync_with_stdio(false);
	std::cin.tie(nullptr);
	int t = 1;
//	cin >> t;
	while (t -- ) solve();
	return 0;
}