「BZOJ3262」陌上花开

有n朵花,每朵花有三个属性:花形(s)、颜色(c)、气味(m),又三个整数表示。现要对每朵花评级,一朵花的级别是它拥有的美丽能超过的花的数量。定义一朵花A比另一朵花B要美丽,当且仅当Sa>=Sb,Ca>=Cb,Ma>=Mb。显然,两朵花可能有同样的属性。需要统计出评出每个等级的花的数量。

Input

第一行为N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分别表示花的数量和最大属性值。
以下N行,每行三个整数si, ci, mi (1 <= si, ci, mi <= K),表示第i朵花的属性

Output

包含N行,分别表示评级为0…N-1的每级花的数量。

Sample Input

10 3
3 3 3
2 3 3
2 3 1
3 1 1
3 1 2
1 3 1
1 1 2
1 2 2
1 3 2
1 2 1

Sample Output

3
1
3
0
1
0
1
0
0
1

HINT

1 <= N <= 100,000, 1 <= K <= 200,000

这是一道CDQ分治经典题

BZOJ1176:还是这道题,加一加容斥

#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
template<typename T> void read(T &x){
	x = 0;char ch = getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}
}
struct fl{
	int a,b,c,ans = 0;
} a[100010];
bool cmp1(fl a,fl b){
	if(a.a!=b.a)return a.a<b.a;
	if(a.b!=b.b)return a.b<b.b;
	if(a.c!=b.c)return a.c<b.c;
	return a.ans < b.ans;
}
bool cmp2(fl a,fl b){
	return a.b<b.b;
}
int n,k,level[100010],buc[200020],szsz[200020];
void add(int p,int v){while(p<200020){szsz[p]+=v;p+=p&(-p);}}
int query(int p){int ret = 0;while(p){ret+=szsz[p];p-=p&(-p);}return ret;}
void solve(int st,int sz){
	if(sz == 1)return;
	int mid = sz/2;
	solve(st,mid);solve(st+mid,sz-mid);
	sort(a+st,a+st+mid,cmp2);sort(a+st+mid,a+st+sz,cmp2);
	int p = 0;
	for(int i=st+mid;i<st+sz;i++){
		while(p<mid && a[st+p].b<=a[i].b){
			add(a[st+p].c,1);
			++p;
		}
		a[i].ans+=query(a[i].c);
	}
	for(int i=0;i<p;i++)add(a[st+i].c,-1);
}
int main() {
	read(n);read(k);
	for(int i=0;i<n;i++){
		read(a[i].a);read(a[i].b);read(a[i].c);
	}
	sort(a,a+n,cmp1);
	solve(0,n);
	sort(a,a+n,cmp1);
	for(int i=n-1;i>=0;i--){
		buc[a[i].ans]++;
		if(i&&a[i].a == a[i-1].a && a[i].b == a[i-1].b && a[i].c == a[i-1].c)a[i-1].ans = a[i].ans;
	}
	for(int i=0;i<n;i++)cout<<buc[i]<<'\n';
	return 0;
}

发表评论

电子邮件地址不会被公开。 必填项已用*标注