【poj2464】树状数组

这道题。。太特么多细节了。。

题意：在平面直角坐标系中给你N个点，stan和ollie玩一个游戏，首先stan在竖直方向上画一条直线，该直线必须要过其中的某个点，然后ollie在水平方向上画一条直线，该直线的要求是要经过一个stan之前画过的点。 这时候平面就被分割成了四块，两个人这时候会有一个得分，stan的得分是平面上第1、3象限内的点的个数，ollie的得分是平面上第2、4象限内的点的个数，在统计的时候在所画线上的点都不计算在内。求最终stan使得自己的最差得分最高，并且输出此时ollie的得分。

题解：

我们可以枚举哪颗星星是中心点，然后就可以知道他们所确定的直线。

线上可以维护四个值：up，down，left，right，分别表示线上四个方位有多少颗星星。

然后我们只要求BL,就可以知道其它：

TL=横坐标比x小的星星总数-BL-left

TR=y坐标比x大的星星总数-TL-up

BR=y坐标比x小的星星总数-BL-down

各种细节><

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 using namespace std;

 const int N=,INF=(int)1e9+;
 int n,pl,mx,c[N],cntx[N],cnty[N],sumx[N],sumy[N],sx[N][],sy[N][],u[N],d[N],l[N],r[N],a1[N],a2[N];
 bool num[N];
 struct node{
     int x,y;
 }a[N];
 struct nd{
     int d,id,tmp;
 }p[*N];

 bool cmp_num(int x,int y){return x<y;}
 bool cmp_d(nd x,nd y){return x.d<y.d;}
 bool cmp_a(node x,node y)
 {
     if(x.x==y.x) return x.y<y.y;
     return x.x<y.x;
 }
 int maxx(int x,int y){return x>y ? x:y;}

 void clear()
 {
     memset(cntx,,sizeof(cntx));
     memset(cnty,,sizeof(cnty));
     memset(sumx,,sizeof(sumx));
     memset(sumy,,sizeof(sumy));
     memset(c,,sizeof(c));
     memset(a1,,sizeof(a1));
     memset(a2,-,sizeof(a2));
 }

 void add(int x)
 {
     for(int i=x;i<=mx;i+=(i&(-i))) c[i]++;
 }
 int getsum(int x)
 {
     int ans=;
     for(int i=x;i>=;i-=(i&(-i))) ans+=c[i];
     return ans;
 }

 int main()
 {
     freopen("a.in","r",stdin);
     // freopen("me.out","w",stdout);
     while()
     {
         scanf("%d",&n);
         if(n==) break;
         pl=;clear();
         for(int i=;i<=n;i++)
         {
             scanf("%d%d",&a[i].x,&a[i].y);
             p[++pl].d=a[i].x;p[pl].id=i;p[pl].tmp=;
             p[++pl].d=a[i].y;p[pl].id=i;p[pl].tmp=;
         }
         sort(p+,p++pl,cmp_d);
         mx=;p[].d=INF;
         for(int i=;i<=pl;i++)
         {
             if(p[i].d!=p[i-].d) mx++;
             if(p[i].tmp==) a[p[i].id].x=mx;
             else a[p[i].id].y=mx;
         }

         sort(a+,a++n,cmp_a);
         // for(int i=1;i<=n;i++)
             // printf("%d %d\n",a[i].x,a[i].y);
         for(int i=;i<=n;i++)
         {
             d[i]=cntx[a[i].x];cntx[a[i].x]++;
             l[i]=cnty[a[i].y];cnty[a[i].y]++;
         }
         // for(int i=1;i<=mx;i++)
             // printf("i = %d  %d %d\n",i,cntx[i],cnty[i]);
         for(int i=;i<=n;i++)
         {
             u[i]=cntx[a[i].x]-d[i]-;
             r[i]=cnty[a[i].y]-l[i]-;
         }
         for(int i=;i<=mx;i++)
         {
             sumx[i]=sumx[i-]+cntx[i];
             sumy[i]=sumy[i-]+cnty[i];
         }
         for(int i=;i<=n;i++)
         {
             int x=a[i].x,y=a[i].y;
             sx[i][]=sumx[x-];
             sx[i][]=sumx[mx]-sumx[x];
             sy[i][]=sumy[y-];
             sy[i][]=sumy[mx]-sumy[y];
         }
         // for(int i=1;i<=n;i++)
         // {
             // printf("%d  sx0 %d sx1 %d sy0 %d sy1 %d d %d u %d l %d r %d\n",i,sx[i][0],sx[i][1],sy[i][0],sy[i][1],d[i],u[i],l[i],r[i]);
         // }
         for(int i=;i<=n;i++)
         {
             int x=a[i].x,y=a[i].y;
             int BL=getsum(a[i].y-)-d[i];
             int TL=sx[i][]-BL-l[i];
             int TR=sy[i][]-TL-u[i];
             int BR=sy[i][]-BL-d[i];
             add(y);
             if(TR+BL<a1[x]) a1[x]=TR+BL,a2[x]=TL+BR;
             else if(TR+BL==a1[x]) a2[x]=maxx(a2[x],TL+BR);
             // printf("%d  BL = %d  BR = %d  TR = %d TL = %d\n",i,BL,BR,TR,TL);
         }
         int ans=,nl=;;
         for(int i=;i<=mx;i++)
         {
             if(a1[i]<INF) ans=maxx(ans,a1[i]);
         }
         printf("Stan: %d; Ollie:",ans);
         memset(num,,sizeof(num));
         for(int i=;i<=n;i++)
             if(a1[i]==ans) num[a2[i]]=;
         for(int i=;i<=n;i++)
             if(num[i]) printf(" %d",i);
         printf(";\n");
     }
     return ;
 }