cdq分治的入门题目,学习一下这个很实用的玩意儿。
大概思想就是把在线变为离线,去掉时间的限制,然后每次对于询问区间分治处理,块内之间的影响可以递归处理,然后剩下的就是块与块之间的影响,这个其实就很简单了,随便搞搞就好,反正是静态查询,双指针啊啥的随便上。。
这题的话第一维排序,第二维cdq,第三维树状数组。
直接做就好啦。
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=2e5+5;
int m;
typedef long long ll;
struct node
{
int a,b,c,cnt,ans;
bool operator < (const node &x) const
{
if (x.a!=a)return a<x.a;
if (x.b!=b)return b<x.b;
return c<x.c;
}
}a[N],tmp[N];
bool cmp(const node &x,const node &y)
{
if (x.b!=y.b)return x.b<y.b;
if (x.c!=y.c)return x.c<y.c;
return x.a<y.a;
}
int n,K,tot,c[N+10],tim[N],cnt[N],T;
inline int lowbit(int x)
{
return x&(-x);
}
inline void add(int x,int val)
{
while (x<=K)
{
if (tim[x]!=T)c[x]=0;
tim[x]=T;
c[x]+=val;
x+=lowbit(x);
}
}
inline int query(int x)
{
int ret=0;
while (x>0)
{
if (tim[x]==T)
ret+=c[x];
x-=lowbit(x);
}
return ret;
}
inline void cdq(int l,int r)
{
if (l==r)
{
a[l].ans+=a[l].cnt-1;
return;
}
int mid=(l+r)>>1;
int l1,l2;
l1=l,l2=mid+1;
fo(i,l,r)
{
if (a[i].a<=mid)tmp[l1++]=a[i];
else tmp[l2++]=a[i];
}
fo(i,l,r)a[i]=tmp[i];
cdq(l,mid);
int j=l;T++;
fo(i,mid+1,r)
{
for(;j<=mid&&a[j].b<=a[i].b;j++)add(a[j].c,a[j].cnt);
a[i].ans+=query(a[i].c);
}
cdq(mid+1,r);
l1=l,l2=mid+1;
fo(i,l,r)
if (l1<=mid&&(l2>r||cmp(a[l1],a[l2])))tmp[i]=a[l1++];
else tmp[i]=a[l2++];
fo(i,l,r)a[i]=tmp[i];
}
int main()
{
scanf("%d%d",&n,&K);
fo(i,1,n)
{
scanf("%d%d%d",&a[i].a,&a[i].b,&a[i].c);
a[i].cnt++;
}
sort(a+1,a+1+n);
fo(i,1,n)
{
if (i==1||a[i-1]<a[i])a[++tot]=a[i];
else a[tot].cnt++;
}
fo(i,1,tot)a[i].a=i;
sort(a+1,a+1+tot,cmp);
cdq(1,tot);
fo(i,1,tot)cnt[a[i].ans]+=a[i].cnt;
fo(i,0,n-1)printf("%d\n",cnt[i]);
}