poj 2464 线段树统计区间..两棵树

仔细分析的话，不用说了，我只是尽自己的力量去写一个尽量好看的代码。。。

思路说一下：

      就是先统计这个点的上区间有多少点，下区间有多少点，然后x排序再求上下区间，这个时候便可以得知已经更新进去的点一定在该点的左下方、左上方、下方。。。

      由此统计出每个点被作为标准点画出，每个人的得分情况。。

     

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<algorithm>
#include<string.h>
using namespace std;
using std::sort;
using std::unique;
const int N = 200022;
const int INF = 0x7fffffff;
struct point
{
    int x,y,k;
}P[N];
int lt[N],lb[N],rt[N],rb[N],Ollie[N];
int Up[N],Down[N];
int X[N];
int sum[N<<2];
bool cmp1(point a,point b)
{
     if(a.y==b.y)
     return a.x<b.x;
     return a.y<b.y;
}
bool cmp2(point a,point b)
{
    if(a.x==b.x)
    return a.y<b.y;
    return a.x<b.x;
}
void update(int t,int l,int r,int i)
{
     sum[t]++;
     if(l==r)return ;
     int m=(l+r)>>1;
     if(i<=m)update(t<<1,l,m,i);
     else update(t<<1|1,m+1,r,i);
}
int query(int t,int l,int r,int L,int R)
{
     if(L>R)return 0;
     if(L<=l&&r<=R)return sum[t];
     int m=(l+r)>>1;
     int ans=0;
     if(L<=m)ans+=query(t<<1,l,m,L,R);
     if(R>m)ans+=query(t<<1|1,m+1,r,L,R);
     return ans;
}
int Max(int a,int b)
{
    return a>b?a:b;
}
int main()
{
    int n;
    while(~scanf("%d",&n)&&n)
    {
        for(int i=1;i<=n;i++)
        {
            scanf("%d%d",&P[i].x,&P[i].y);
            P[i].k=i;
        }
        sort(P+1,P+n+1,cmp1);
        int topy=0;
        for(int i=1;i<=n;i++)
        {
            if(i==n){P[i].y=topy;break;}
            if(P[i].y==P[i+1].y)P[i].y=topy;
            else P[i].y=topy++;
        }
        memset(sum,0,sizeof(sum));
        for(int i=1;i<=n;i++)update(1,0,topy,P[i].y);
        for(int i=1;i<=n;i++)
        {
            Up[P[i].k]=query(1,0,topy,P[i].y+1,topy);
            Down[P[i].k]=query(1,0,topy,0,P[i].y-1);
        }
        sort(P+1,P+n+1,cmp2);
        int topx=0;
        for(int i=1;i<=n;i++)
        {
            if(i==n){P[i].x=topx;break;}
            if(P[i].x==P[i+1].x)P[i].x=topx;
            else P[i].x=topx++;
        }
        memset(sum,0,sizeof(sum));
        int Stan=0,top=0;
        for(int i=0;i<=topx;i++)X[i]=INF;
        for(int i=1;i<=n;i++)
        {
            lt[i]=query(1,0,topy,P[i].y+1,topy);
            lb[i]=query(1,0,topy,0,P[i].y-1);
            rb[i]=Down[P[i].k]-lb[i];
            for(int j=i-1;j>=1&&P[i].x==P[j].x;j--)lb[i]--;
            rt[i]=Up[P[i].k]-lt[i];
            for(int j=i+1;j<=n&&P[i].x==P[j].x;j++)rt[i]--;
            if(lb[i]+rt[i]<X[P[i].x])
            {
                X[P[i].x]=lb[i]+rt[i];
            }
            update(1,0,topy,P[i].y);
        }
        for(int i=0;i<=topx;i++)Stan=Max(Stan,X[i]);
        for(int i=1;i<=n;i++)
        {
            if(Stan==X[P[i].x])
            Ollie[top++]=lt[i]+rb[i];
        }
        sort(Ollie,Ollie+top);
        top=unique(Ollie,Ollie+top)-Ollie;
        printf("Stan: %d; Ollie:",Stan);
        for(int i=0;i<top;i++)printf(" %d",Ollie[i]);
        printf(";\n");
    }
    return 0;
}