【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
各种细节><
1 #include<cstdio>
2 #include<cstdlib>
3 #include<cstring>
4 #include<iostream>
5 #include<algorithm>
6 using namespace std;
7
8 const int N=200010,INF=(int)1e9+100;
9 int n,pl,mx,c[N],cntx[N],cnty[N],sumx[N],sumy[N],sx[N][2],sy[N][2],u[N],d[N],l[N],r[N],a1[N],a2[N];
10 bool num[N];
11 struct node{
12 int x,y;
13 }a[N];
14 struct nd{
15 int d,id,tmp;
16 }p[2*N];
17
18 bool cmp_num(int x,int y){return x<y;}
19 bool cmp_d(nd x,nd y){return x.d<y.d;}
20 bool cmp_a(node x,node y)
21 {
22 if(x.x==y.x) return x.y<y.y;
23 return x.x<y.x;
24 }
25 int maxx(int x,int y){return x>y ? x:y;}
26
27 void clear()
28 {
29 memset(cntx,0,sizeof(cntx));
30 memset(cnty,0,sizeof(cnty));
31 memset(sumx,0,sizeof(sumx));
32 memset(sumy,0,sizeof(sumy));
33 memset(c,0,sizeof(c));
34 memset(a1,63,sizeof(a1));
35 memset(a2,-1,sizeof(a2));
36 }
37
38 void add(int x)
39 {
40 for(int i=x;i<=mx;i+=(i&(-i))) c[i]++;
41 }
42 int getsum(int x)
43 {
44 int ans=0;
45 for(int i=x;i>=1;i-=(i&(-i))) ans+=c[i];
46 return ans;
47 }
48
49 int main()
50 {
51 freopen("a.in","r",stdin);
52 // freopen("me.out","w",stdout);
53 while(1)
54 {
55 scanf("%d",&n);
56 if(n==0) break;
57 pl=0;clear();
58 for(int i=1;i<=n;i++)
59 {
60 scanf("%d%d",&a[i].x,&a[i].y);
61 p[++pl].d=a[i].x;p[pl].id=i;p[pl].tmp=0;
62 p[++pl].d=a[i].y;p[pl].id=i;p[pl].tmp=1;
63 }
64 sort(p+1,p+1+pl,cmp_d);
65 mx=0;p[0].d=INF;
66 for(int i=1;i<=pl;i++)
67 {
68 if(p[i].d!=p[i-1].d) mx++;
69 if(p[i].tmp==0) a[p[i].id].x=mx;
70 else a[p[i].id].y=mx;
71 }
72
73 sort(a+1,a+1+n,cmp_a);
74 // for(int i=1;i<=n;i++)
75 // printf("%d %d\n",a[i].x,a[i].y);
76 for(int i=1;i<=n;i++)
77 {
78 d[i]=cntx[a[i].x];cntx[a[i].x]++;
79 l[i]=cnty[a[i].y];cnty[a[i].y]++;
80 }
81 // for(int i=1;i<=mx;i++)
82 // printf("i = %d %d %d\n",i,cntx[i],cnty[i]);
83 for(int i=1;i<=n;i++)
84 {
85 u[i]=cntx[a[i].x]-d[i]-1;
86 r[i]=cnty[a[i].y]-l[i]-1;
87 }
88 for(int i=1;i<=mx;i++)
89 {
90 sumx[i]=sumx[i-1]+cntx[i];
91 sumy[i]=sumy[i-1]+cnty[i];
92 }
93 for(int i=1;i<=n;i++)
94 {
95 int x=a[i].x,y=a[i].y;
96 sx[i][0]=sumx[x-1];
97 sx[i][1]=sumx[mx]-sumx[x];
98 sy[i][0]=sumy[y-1];
99 sy[i][1]=sumy[mx]-sumy[y];
100 }
101 // for(int i=1;i<=n;i++)
102 // {
103 // 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]);
104 // }
105 for(int i=1;i<=n;i++)
106 {
107 int x=a[i].x,y=a[i].y;
108 int BL=getsum(a[i].y-1)-d[i];
109 int TL=sx[i][0]-BL-l[i];
110 int TR=sy[i][1]-TL-u[i];
111 int BR=sy[i][0]-BL-d[i];
112 add(y);
113 if(TR+BL<a1[x]) a1[x]=TR+BL,a2[x]=TL+BR;
114 else if(TR+BL==a1[x]) a2[x]=maxx(a2[x],TL+BR);
115 // printf("%d BL = %d BR = %d TR = %d TL = %d\n",i,BL,BR,TR,TL);
116 }
117 int ans=0,nl=0;;
118 for(int i=1;i<=mx;i++)
119 {
120 if(a1[i]<INF) ans=maxx(ans,a1[i]);
121 }
122 printf("Stan: %d; Ollie:",ans);
123 memset(num,0,sizeof(num));
124 for(int i=0;i<=n;i++)
125 if(a1[i]==ans) num[a2[i]]=1;
126 for(int i=0;i<=n;i++)
127 if(num[i]) printf(" %d",i);
128 printf(";\n");
129 }
130 return 0;
131 }
