首页 > 代码库 > (树状数组) poj 2481

(树状数组) poj 2481

Cows
Time Limit: 3000MS Memory Limit: 65536K
Total Submissions: 13421 Accepted: 4442

Description

Farmer John‘s cows have discovered that the clover growing along the ridge of the hill (which we can think of as a one-dimensional number line) in his field is particularly good. 

Farmer John has N cows (we number the cows from 1 to N). Each of Farmer John‘s N cows has a range of clover that she particularly likes (these ranges might overlap). The ranges are defined by a closed interval [S,E]. 

But some cows are strong and some are weak. Given two cows: cowi and cowj, their favourite clover range is [Si, Ei] and [Sj, Ej]. If Si <= Sj and Ej <= Ei and Ei - Si > Ej - Sj, we say that cowi is stronger than cowj

For each cow, how many cows are stronger than her? Farmer John needs your help!

Input

The input contains multiple test cases. 
For each test case, the first line is an integer N (1 <= N <= 105), which is the number of cows. Then come N lines, the i-th of which contains two integers: S and E(0 <= S < E <= 105) specifying the start end location respectively of a range preferred by some cow. Locations are given as distance from the start of the ridge. 

The end of the input contains a single 0.

Output

For each test case, output one line containing n space-separated integers, the i-th of which specifying the number of cows that are stronger than cowi

Sample Input

31 20 33 40

Sample Output

1 0 0

Hint

Huge input and output,scanf and printf is recommended.
#include<iostream>#include<cstdio>#include<cstring>#include<string>#include<cmath>#include<algorithm>#include<cstdlib>using namespace std;#define N 100010int c[N],val[N];struct node{      int x,y,id;}e[N];bool cmp(node a,node b){      if(a.y!=b.y) return a.y>b.y;      return a.x<b.x;}int lowbit(int x){      return x&(-x);}void update(int pos,int m){      while(pos<N)      {           c[pos]+=m;           pos+=lowbit(pos);      }}int sum(int x){      int sum=0;      while(x>0)      {            sum+=c[x];            x=x-lowbit(x);      }      return sum;}int main(){      int n;      while(scanf("%d",&n)!=EOF)      {            if(n==0)                  break;            memset(c,0,sizeof(c));            for(int i=1;i<=n;i++)            {                 scanf("%d%d",&e[i].x,&e[i].y);                 e[i].x++,e[i].y++;                 e[i].id=i;            }            sort(e+1,e+1+n,cmp);            val[e[1].id]=sum(e[1].x);            update(e[1].x,1);            for(int i=2;i<=n;i++)            {                if(e[i].x==e[i-1].x&&e[i].y==e[i-1].y)                        val[e[i].id]=val[e[i-1].id];                else val[e[i].id]=sum(e[i].x);                update(e[i].x,1);            }            printf("%d",val[1]);            for(int i=2;i<=n;i++)                  printf(" %d",val[i]);            printf("\n");      }      return 0;}

  

(树状数组) poj 2481