树状数组

例题

「POJ2352」Stars
题目描述

天文学家常常检查星星地图，星星都有它的x,y坐标，星星的等级的是左下方星星数量决定，包含正左和正下。
例如,看看上面的星图。星星5的等级为3 (由星星1、2和4决定的)。星星2的等级为1（由星星1决定的）。在这张地图上0级的星星有一颗，1级的星星有两颗，2级的星星有一颗，3级的星星有一颗，

你要编写一个程序,计算每个等级的星星的数量。

输入

第一行为星星的数量N（(1<=N<=15000)

接下来N行，每行为星星的x,y坐标，用空格来分开(0<=X,Y<=32000)，每一个点上只有一个星星，按Y的升序给出坐标，如果Y相同，则按照X的升序给出。

输出

输出应该包含N行,每行一个数字。第一行0级星星的数量,第二行1级星星的数量等等,最后一行包含n – 1星星的数量。

样例输入
5
1 1
5 1
7 1
3 3
5 5
样例输出
1
2
1
1
0

分析

由于是按y的升序输入，使用只需要找x比Xi小的，这里用树状数组实现。

代码

#include<bits/stdc++.h>
using namespace std;
#define N 15000+5
#define MAXN 32000+5
#define INF 0x3f3f3f3f
#define lowbit(x) x&(-x)
#define full(a,b) memset(a,b,sizeof a)
int n,tree[MAXN],ji[N];
int sum(int x){//求[1~x]中的和 
	int ret=0;
	while(x>0){
		ret+=tree[x];
		x-=lowbit(x);
	}
	return ret;
}
void update(int x,int y){//把tree[x]+y 
	while(x<=MAXN){
		tree[x]+=y;
		x+=lowbit(x);
	}
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		int x,y;scanf("%d%d",&x,&y);
		++x;//0的lowbit值为0，会死循环 
		ji[sum(x)]++;//sum(x)是第i颗星的级数 
		update(x,1);//坐标为x的星又增加了一颗 
	}
	for(int i=0;i<n;i++)
	printf("%d\n",ji[i]);
	return 0;
}