POJ 2352 star level

题目链接：

http://poj.org/problem?id=2352

题目大意：
对于每一颗星星来说，都有一个属于自己的level，这个值为其他星星x，y坐标均不大于本星星的个数。输入时按先y由小到大，再x由小到大排列，无相同位置的星星

仔细一想其实这道题目根本不需要用到y，我是反向来思考的，构建一个保存x坐标由0到maxn的线段树，因为最后一个星星肯定不可能x，y同时不大于其他星星，所

以从后面开始看，只计算比它x小或相等的星星的个数，用cnt[ans]++，来记录成为ans水平的个数，再删去这个点，继续找上一个点对应的level

在这里要注意的是数组的范围，我就因为这个改了半天错，第一次发现数组范围定的不够大也是WA

他这里星星个数为15000，但是我们保存的是x的坐标，所以要看x最大达32000，所以用4*32000的大小保存tree[]，我的代码里懒得修改了就写成了

#define N 15005
int cnt[N],num[2*N],tree[8*N],D,maxn;//开始我写成的是4*N所以报错T T

#include <cstdio>
#include <cstring>
using namespace std;
#define N 15005
int cnt[N],num[2*N],tree[8*N],D,maxn;

int max(int a,int b)
{
    return a>b?a:b;
}

void update(int i)
{
    for(;i^1;i>>=1)
        tree[i>>1]=tree[i]+tree[i^1];
}

int query(int a,int b)
{
    int i=D+a-1,j=D+b+1,ans=0;
    for(;i^j^1;i>>=1,j>>=1)
    {
        if(~i&1) ans+=tree[i^1];
        if(j&1)  ans+=tree[j^1];
    }
    return ans;
}

int main()
{
    int n,x[N],y;
    while(scanf("%d",&n)!=EOF){
        maxn=0;
        memset(num,0,sizeof(num));
        //memset(tree,0,sizeof(tree));
        for(int i=1;i<=n;i++){
            scanf("%d%d",&x[i],&y);
            maxn=max(maxn,x[i]);
            num[x[i]]++;
        }
        for(D=1;D<maxn+1+2;D<<=1);
        for(int i=1;i<=maxn+1;i++) tree[D+i]=num[i-1];
        for(int i=D-1;i>=1;i--) tree[i]=tree[i<<1]+tree[i<<1|1];
        for(int i=0;i<=n-1;i++) cnt[i]=0;
        //for(int i=1;i<=2*D-1;i++) printf("%d
",tree[i]);
        for(int i=n;i>=1;i--){
           // printf("%d",query(1,x[i]+1));
            tree[D+x[i]+1]--;
            update(D+x[i]+1);
            cnt[query(1,x[i]+1)]++;
        }
        for(int i=0;i<n;i++) printf("%d
",cnt[i]);
    }
    return 0;
}

当然从前往后思考也是一样的，那样的话是不断将x添入线段树，这里有一份是学长的代码，果然比我反向思考的要简洁好多，而且从他的代码也得到了其实tree的底层在没有

得到顶点的必要时，也是可以不必强求定位1<<n这种2的几次方的大小的，而在这里正好适用，因为每次找个数和只求到i^j^1即可（这个代表左右子树）；

代码如下：

#include <cstdio>
#include <cstring>
const int maxn = 32010;
int tree[maxn<<2],a[maxn];
int N,D;

void update(int i){
    for(;i^1;i >>= 1){
        tree[i>>1] = tree[i]+tree[i^1];
    }
}

int query(int x,int y){
    int i = D+x-1,j = D+y+1,ans = 0;
    for(;i^j^1;i >>= 1,j >>= 1){
        if(~i&1)    ans += tree[i^1];
        if(j&1)     ans += tree[j^1];
    }
    return ans;
}

int main(){
    while(scanf("%d",&N) != EOF){
        memset(tree,0,sizeof(tree));
        memset(a,0,sizeof(a));
        D = 32010;
        for(int i = 0;i < N;i++){
            int x,y;
            scanf("%d%d",&x,&y);
            a[query(1,++x)]++;
            tree[D+x]++;
            update(D+x);
        }
        for(int i = 0;i < N;i++)    printf("%d
",a[i]);
    }
    return 0;
}