POJ 2464 Brownie Points II（树状数组||线段树）

转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents    by---cxlove

很久以前，用线段树做过。

题意有点悬乎：有一些点，Stan选择某个点，经过这个点画一条竖线，Ollie选择一个经过这条直接的点画了条横线。

Stan选一，三象限的点，Ollie选二、四象限的点。

Stan的策略是，自己画一条竖线之后，Ollie有很多种选择，在所有选择中，Stan能获得数目的最小值的最大值，而Ollie的选择便 是让自己越多越好。

http://poj.org/problem?id=2464  

首先对于某个方向离散化，我是按x排序，将y离散化。

建立两个树状数组，表示竖线左边的点个数，右边的点个数

然后按x递增扫描一遍，预先先把所有点插入到右树中。

然后逐步更新左树和右树。

#include<iostream>  #include<cstdio>  #include<map>  #include<cstring>  #include<cmath>  #include<vector>  #include<algorithm>  #include<set>  #include<string>  #include<queue>  #define inf 1600005  #define M 40  #define N 200000 #define maxn 300005  #define eps 1e-12#define zero(a) fabs(a)<eps  #define Min(a,b) ((a)<(b)?(a):(b))  #define Max(a,b) ((a)>(b)?(a):(b))  #define pb(a) push_back(a)  #define mp(a,b) make_pair(a,b)  #define mem(a,b) memset(a,b,sizeof(a))  #define LL unsigned long long  #define MOD 1000000007#define lson step<<1#define rson step<<1|1#define sqr(a) ((a)*(a))  #define Key_value ch[ch[root][1]][0]  #define test puts("OK");  #define pi acos(-1.0)#define lowbit(x) ((-(x))&(x))#pragma comment(linker, "/STACK:1024000000,1024000000")  using namespace std;struct Point{    int x,y;    bool operator<(const Point p)const{        return x<p.x;    }}p[N+5];int y[N],cnt;int l[N+5],r[N+5];int n;void Update(int *s,int x,int val){    for(int i=x;i<=cnt;i+=lowbit(i))        s[i]+=val;}int sum(int *s,int x){    int ret=0;    for(int i=x;i>0;i-=lowbit(i))        ret+=s[i];    return ret;}int main(){    //freopen("in.txt","r",stdin);    while(scanf("%d",&n)!=EOF&&n){        for(int i=0;i<n;i++){            scanf("%d%d",&p[i].x,&p[i].y);            y[i]=p[i].y;        }        sort(p,p+n);        sort(y,y+n);        cnt=unique(y,y+n)-y;        mem(l,0);mem(r,0);        for(int i=0;i<n;i++)            Update(r,lower_bound(y,y+cnt,p[i].y)-y+1,1);        int Stan=-1,st=0;        vector<int>Ollie;        for(int i=1;i<=n;i++){            if(i==n||p[i].x!=p[i-1].x){                for(int j=st;j<i;j++)                    Update(r,lower_bound(y,y+cnt,p[j].y)-y+1,-1);                int stan=-1,ollie=-1;                for(int j=st;j<i;j++){                    int pos=lower_bound(y,y+cnt,p[j].y)-y+1;                    int a=sum(r,cnt)-sum(r,pos)+sum(l,pos-1);                    int b=sum(l,cnt)-sum(l,pos)+sum(r,pos-1);                    if(b==ollie) stan=min(stan,a);                    else if(b>ollie) stan=a,ollie=b;                  }                if(stan>Stan){                    Stan=stan;                    Ollie.clear();                }                if(stan==Stan) Ollie.pb(ollie);                for(int j=st;j<i;j++)                    Update(l,lower_bound(y,y+cnt,p[j].y)-y+1,1);                st=i;            }        }        printf("Stan: %d; Ollie:",Stan);        sort(Ollie.begin(),Ollie.end());        cnt=unique(Ollie.begin(),Ollie.end())-Ollie.begin();        for(int i=0;i<cnt;i++) printf(" %d",Ollie[i]);        puts(";");    }    return 0;}