bzoj 3262: 陌上花開 -- CDQ分治
阿新 • • 發佈：2017-08-30
time line int 屬性。 inf des for blog logs
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
1
3
0
1
0
1
0
0
1
3262: 陌上花開
Time Limit: 20 Sec Memory Limit: 256 MBDescription
有n朵花，每朵花有三個屬性：花形(s)、顏色(c)、氣味(m)，又三個整數表示。現要對每朵花評級，一朵花的級別是它擁有的美麗能超過的花的數量。定義一朵花A比另一朵花B要美麗，當且僅當Sa>=Sb,Ca>=Cb,Ma>=Mb。顯然，兩朵花可能有同樣的屬性。需要統計出評出每個等級的花的數量。Input
第一行為N,K (1 <= N <= 100,000, 1 <= K <= 200,000 ), 分別表示花的數量和最大屬性值。Output
包含N行，分別表示評級為0...N-1的每級花的數量。Sample Input
10 33 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
31
3
0
1
0
1
0
0
1
HINT
1 <= N <= 100,000, 1 <= K <= 200,000
Source
#include<map> #include<cmath> #include<queue> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> using namespace std; #define inf 1000000007 #define ll long long #define N 200010 inline int rd() { int x=0,f=1;char ch=getchar();while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} return x*f; } struct qaz{int a,b,c,s,ans;}a[N],p[N]; bool cmp(qaz a,qaz b){return a.b==b.b&&a.a==b.a?a.c<b.c:(a.a==b.a?a.b<b.b:a.a<b.a);} bool cmp2(qaz a,qaz b){return a.b==b.b?a.c<b.c:a.b<b.b;} int n,k,nn,sum[N]; void add(int x,int v){for(int i=x;i<=k;i+=i&(-i)) sum[i]+=v;} int fd(int x){int tt=0;for(int i=x;i;i-=i&(-i)) tt+=sum[i];return tt;} void sol(int l,int r) { if(l==r) return; int mid=l+r>>1,i=l,j=mid+1; sol(l,mid);sol(mid+1,r); sort(p+l,p+mid+1,cmp2); sort(p+mid+1,p+r+1,cmp2); while(j<=r) { while(i<=mid&&p[i].b<=p[j].b) add(p[i].c,p[i].s),i++; p[j].ans+=fd(p[j].c);j++; } for(j=l;j<i;j++) add(p[j].c,-p[j].s); } int ans[N]; int main() { nn=rd();k=rd(); for(int i=1;i<=nn;i++) { a[i].a=rd();a[i].b=rd();a[i].c=rd(); } sort(a+1,a+nn+1,cmp); for(int i=1,ji=0;i<=nn;i++) { ji++; if(a[i].a!=a[i+1].a||a[i].b!=a[i+1].b||a[i].c!=a[i+1].c) { p[++n]=a[i]; p[n].s=ji;ji=0; } } sol(1,n); for(int i=1;i<=n;i++) ans[p[i].ans+p[i].s-1]+=p[i].s; for(int i=0;i<nn;i++) printf("%d\n",ans[i]); return 0; }
bzoj 3262: 陌上花開 -- CDQ分治