poj 2481 cows 树状数组 老少皆宜版解答

请看几段区间

1      --------------

2    -------------

3       --------

4           -----

先按y排序，再按x排序；

排好序后就跟stars那题一样了

一次遍历排序后的数组，由于比当前遍历元素strong的区间只可能存在于已经遍历过的元素中；所以我们可以放心大胆的直接对树状数组求和

然后向后更新树状数组中统治x的点上的值，因为那些位置都多了一个比  以他们为左端点的区间  强壮  的区间；

然后注意一下相同区间的处理即可

跑了1157ms的代码。其实感觉也蛮快的了*-*

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define MAX 100010
int c[MAX];
int val[MAX];
struct point{
    int x,y;
    int id;
}a[MAX];
int cmp(point a,point b)
{
    if(a.y==b.y)
        return a.x<b.x;
    return a.y>b.y; 
}
int lowbit(int x)
{
    return x&(-x);
}
void update(int x,int d)
{
    while(x<=MAX)
    {
        c[x]+=d;
        x+=lowbit(x);
    }
}
int sum(int x)
{
    int ans=0;
    while(x>0)
    {
        ans+=c[x];
        x-=lowbit(x);
    }
    return ans;
}
int main()
{
    int n,i,x,y;
    while(scanf("%d",&n),n)
    {
        memset(c,0,sizeof(c));
        memset(val,0,sizeof(val));
        for(i=0;i<n;i++)
        {
            scanf("%d%d",&a[i].x,&a[i].y);
            a[i].id=i;
            a[i].x++;
        }
        sort(a,a+n,cmp);
        val[a[0].id]=sum(a[0].x);
        update(a[0].x,1);
        for(i=1;i<n;i++)
        {
            if(a[i].x==a[i-1].x&&a[i].y==a[i-1].y)
             {
               val[a[i].id]=val[a[i-1].id];
               update(a[i].x,1);
             }
               else 
               {
                                 val[a[i].id]=sum(a[i].x); 
                      update(a[i].x,1);
               }
        }
        for(i=0;i<n;i++)
        {
            if(!i)
                printf("%d",val[i]);
            else printf(" %d",val[i]);
        }
        printf("\n");
    }
    return 0;
}