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[*N],D,maxn;//开始我写成的是4*N所以报错T T

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

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

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

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

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

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

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

代码如下：

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

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

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

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