Vijos P1066 弱弱的战壕

Vijos P1066 弱弱的战壕

---

题目

描述

永恒和mx正在玩一个即时战略游戏，名字嘛~~恕本人记性不好，忘了-_-b。
mx在他的基地附近建立了n个战壕，每个战壕都是一个独立的作战单位，射程可以达到无限（“mx不赢定了？！？”永恒ftING…@_@）。
但是，战壕有一个弱点，就是只能攻击它的左下方，说白了就是横纵坐标都不大于它的点（mx：“我的战壕为什么这么菜”ToT）。这样，永恒就可以从别的地方进攻摧毁战壕，从而消灭mx的部队。
战壕都有一个保护范围，同它的攻击范围一样，它可以保护处在它左下方的战壕。所有处于它保护范围的战壕都叫做它的保护对象。这样，永恒就必须找到mx的战壕中保护对象最多的点，从而优先消灭它。
现在，由于永恒没有时间来计算，所以拜托你来完成这个任务：
给出这n个战壕的坐标xi、yi，要你求出保护对象个数为0，1，2……n-1的战壕的个数。
格式

输入格式

第一行，一个正整数n（1<=n<=15000）
接下来n行，每行两个数xi,yi，代表第i个点的坐标
（1<=xi,yi<=32000）
注意：可能包含多重战壕的情况（即有数个点在同一坐标）
输出格式

输出n行，分别代表保护对象为0，1，2……n-1的战壕的个数。
样例1

样例输入1

5
1 1
5 1
7 1
3 3
5 5

样例输出1

1
2
1
1
0

限制

各点2s（算是宽限吧^_^）
来源

URAL1028战役版

---

题解

因为最后输出的答案与输入的顺序无关，只要输出保护对象为0~n-1的战壕的个数即可，所以我们很容易想到以x为第一关键字，y为第二关键字排序，然后按顺序算出第y个点时，线段1到y[i]中的数和为x（即为y点保护对象的个数）把ans[x]++，然后把y[i]加入线段树，最后输出1~n-1的所有ans即可

---

代码

#include<cstdio>
#include<cstdlib>
using namespace std;

int n;
int x[15005],y[15005],ans[15005];
struct data{
    int l,r,sum;
} tr[100005];

int readln()
{
    int x=0;
    char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while ('0'<=ch&&ch<='9') x=x*10+ch-48,ch=getchar();
    return x;
}

void build(int l,int r,int rt)
{
    tr[rt].l=l;tr[rt].r=r;
    if (l==r) return;
    int mid=(l+r)>>1;
    build(l,mid,rt<<1);build(mid+1,r,rt<<1|1);
}

void qsort(int l,int r)
{
    int i=l,j=r,t=rand()%(r-l+1)+l,midx=x[t],midy=y[t];
    do {
        while (x[i]<midx||(x[i]==midx&&y[i]<midy)) i++;
        while (x[j]>midx||(x[j]==midx&&y[j]>midy)) j--;
        if (i<=j) {
            t=x[i];x[i]=x[j];x[j]=t;
            t=y[i];y[i]=y[j];y[j]=t;
            i++;j--;
        }
    } while (i<=j);
    if (i<r) qsort(i,r);
    if (l<j) qsort(l,j);
}

void updata(int x,int rt)
{
    int l=tr[rt].l,r=tr[rt].r,mid=(tr[rt].l+tr[rt].r)>>1;
    if (l==r) {tr[rt].sum++;return;}
    if (x<=mid) updata(x,rt<<1);
    else updata(x,rt<<1|1);
    tr[rt].sum=tr[rt<<1].sum+tr[rt<<1|1].sum;
}

int query(int l,int r,int rt)
{
    int ll=tr[rt].l,rr=tr[rt].r;
    if (ll==rr) return tr[rt].sum;
    int mid=(ll+rr)>>1;
    if (r<=mid) return query(l,r,rt<<1);
    else if (mid<l) return query(l,r,rt<<1|1);
    else return query(l,mid,rt<<1)+query(mid+1,r,rt<<1|1);
}

int main()
{
    n=readln();
    build(1,32000,1);
    for (int i=1;i<=n;i++) x[i]=readln(),y[i]=readln();
    qsort(1,n);
    for (int i=1;i<=n;i++)
    {
        ans[query(1,y[i],1)]++;
        updata(y[i],1);
    }
    for (int i=0;i<n;i++) printf("%d\n",ans[i]);
    return 0;
}