[不完全动态凸包]SGU277

由于SDOI某题是线段树套不完全动态凸包，以前没写过，所以抱着试试的心态去写了写。

（感觉还是蛮好写的对吧？）

由于不想手写平衡树，又因为有一个二分斜率的东西，所以set要搞一个开关。

于是写啊写，搞了好久，终于对了。

然后悲剧发生了。。。

常数太大，死交不过。尼玛？根本优化不动。

一怒一下，改成手写。。。然后悲剧又发生了。。。

还是常数太大，还是死交不过。。。。尼玛？！！！！

唉。。。

我的做法是分别维护上凸和下凸（两者只是叉积符号相反）。

set里按(x,y)升序排，插入就向前向后维护一下。

set比较烦人的就是begin和end等等判越界的问题。

然后就顺手把SGU277 A了。

说明正确性和复杂度还是没问题的。

#include<cstdio>
#include<cstring>
#include<set>
using namespace std;
#define mp make_pair
#define x first
#define y second
template<class T> inline void read(T&x){bool fu=0;char c;for(c=getchar();c<=32;c=getchar());if(c=='-')fu=1,c=getchar();for(x=0;c>32;c=getchar())x=x*10+c-'0';if(fu)x=-x;};
typedef long long ll;
typedef pair<ll,ll> node;
typedef set<node>::iterator sit;

ll operator*(node a,node b){return 1LL*a.x*b.y-1LL*a.y*b.x;}
node operator-(node a,node b){return mp(a.x-b.x,a.y-b.y);}
void Read(node&a){read(a.x),read(a.y);}
const ll INF=ll(1e18);
int ca,i;char ch;int F;ll ans[2];
set<node> S[2];
node now,oo=mp(INF,INF),zero;sit it,pos;node p1,p2;
#define B S[t].begin()
#define E S[t].end()
#define IN insert
#define ER erase
#define up(x,y) x=max(x,y)
#define k0 (t==0?1:-1)
#define L(a) ((a)==S[t].begin()?oo:*(pos=a,--pos))
#define R(a) ((pos=a,++pos)==S[t].end()?oo:*pos)
#define LB S[t].lower_bound(now)
#define FI S[t].find(now)
int xj(node a,node b,node c){if(max(max(a,b),c)==oo)return 0;    return ((b-a)*(c-a)>0)*2-1;}// >0:1 <0:-1 error:0
#define C(p,K) {  p1=L(p);if(p1==oo)p1=zero;
                  p2=R(p);if(p2==oo)p2=zero;
                  ans[t]+=K*((p1*(*p)+(*p)*p2)-p1*p2);   }
void ins()
{
    for(int t=0;t<=1;t++){
        if((it=LB)!=E&&((*it)==now||xj(L(it),now,*it)*k0>0))continue;
        if((it=LB)!=E)for(;xj(now,*it,R(it))*k0>0;S[t].erase(it++))C(it,-1);
        if((it=LB)!=B)for(it--;xj(L(it),*it,now)*k0>0;S[t].erase(it--))C(it,-1);
        S[t].IN(now);C(FI,1);
    }
}

int main()
{
    for(i=0;i<=2;i++)Read(now),ins();
    for(read(ca);ca;ca--){Read(now),ins();printf("%lld
",ans[1]-ans[0]);}
    return 0;
}

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<set>
#include<iostream>
using namespace std;
#define mp make_pair
#define rep(i,x,y) for(i=(x);i<=(y);i++)
#define _rep(i,x,y) for(i=(x);i>=(y);i--)
#define REP(i,x,y) for(int i=(x);i<=(y);i++)
#define _REP(i,x,y) for(int i=(x);i>=(y);i--)
#define CL(S,x) memset(S,x,sizeof(S))
#define CP(S1,S2) memcpy(S1,S2,sizeof(S2))
#define ALL(x,S) for(x=S.begin();x!=S.end();x++)
#define max3(a,b,c) max(max(a,b),c)
typedef long long ll;
typedef double ld;

inline char getc(){char c;for(c=getchar();c<=32;c=getchar());return c;}

const int N=400010;
const ll INF=ll(1e18);
const ld eps=1e-6;
const int 
int thec=0x7fffffff;

int ca,ca0,n,m,i,j,k,l,r;char ch;int F;ll ans;

#define DC(x) x^=(ans&thec)
template<class T> inline void read(T&x){bool fu=0;char c;for(c=getchar();c<=32;c=getchar());if(c=='-')fu=1,c=getchar();for(x=0;c>32;c=getchar())x=x*10+c-'0';if(fu)x=-x;DC(x);};
template<class T> inline void read(T&x,T&y){read(x);read(y);};

struct node{ll x,y;ld k;};
typedef set<node>::iterator sit;
node u,now,oo=(node){INF,INF,0};sit it,it1,it2,pos;
#define B S[t].begin()
#define E S[t].end()
#define LB lower_bound
#define IN insert
#define ER erase
#define up(x,y) x=max(x,y)

bool operator<(node a,node b){if(F==0)return mp(a.x,a.y)<mp(b.x,b.y);else if(F==1)return a.k>b.k;else return a.k<b.k;}
bool operator==(node a,node b){return mp(a.x,a.y)==mp(b.x,b.y);}
void Read(node&a){read(a.x,a.y);}
ll operator*(node a,node b){return 1LL*a.x*b.y-1LL*a.y*b.x;}
node operator-(node a,node b){return (node){a.x-b.x,a.y-b.y,0};}
int xj(node a,node b,node c){if(max3(a,b,c)==oo)return 0;    return ((b-a)*(c-a)>0)*2-1;}// >0:1 <0:-1 error:0
ld calc(node a,node b){return ld(b.y-a.y)/(b.x-a.x+eps);}
ll mul(node a,node b){return a.x*b.x+a.y*b.y;}

void P(node a){cout<<a.x<<','<<a.y<<' ';}
void Pl(node a){cout<<a.x<<','<<a.y<<endl;}

vector<int> V[1<<20];
struct tnode
{
    set<node> S[2];
    #define k0 (t==0?1:-1)
    #define L(a) ((a)==S[t].begin()?oo:*(pos=a,--pos))
    #define R(a) ((pos=a,++pos)==S[t].end()?oo:*pos)
    void ins()
    {
        REP(t,0,1){
            it=S[t].LB(now);
            if(it!=E&&((*it)==now||xj(L(it),now,*it)*k0>0))continue;
            if(it!=E)for(;xj(now,*it,R(it))*k0>0;S[t].erase(it++));
            if(it!=B)for(it--;xj(L(it),*it,now)*k0>0;S[t].erase(it--));
            if((it1=S[t].LB(now))!=E){u=*it1;S[t].ER(it1);u.k=calc(now,u);S[t].IN(u);}
            if((it2=S[t].LB(now))!=B){u=*(--it2);now.k=calc(u,now);}else now.k=(t==0)?INF:-INF;
            S[t].IN(now);
        }
    }
    void find()
    {
        now.k=ld(-now.x)/(now.y+(now.y==0)*eps);
        REP(t,0,1){
            F=t+1;it=S[t].LB(now);
            if(it!=E)    up(ans,mul(*it,now));
            if(it!=B)    up(ans,mul(*(--it),now));
        }
        F=0;
    }
}T[1<<20];

void tins(int i,int tx,int ty,int d){
    T[i].ins();
    if(tx!=ty)if(d<=(tx+ty)/2)tins(i*2,tx,(tx+ty)/2,d);else tins(i*2+1,(tx+ty)/2+1,ty,d);
}
void find(int i,int tx,int ty,int l,int r){
    if(ty<l||tx>r)return;
    if(tx>=l&&ty<=r){T[i].find();return;}
    find(i*2,tx,(tx+ty)/2,l,r);find(i*2+1,(tx+ty)/2+1,ty,l,r);
}


int main()
{
    //freopen("vset.in","r",stdin);freopen("vset.out","w",stdout);
    read(ca);n=ca;ch=getc();if(ch=='E')thec=0;
    rep(ca0,1,ca)
    {
        ch=getc();
        if(ch=='A')Read(now),tins(1,1,n,++m);
        else if(ch=='Q')
        {
            Read(now);read(l,r);
            ans=-INF;find(1,1,n,l,r);printf("%I64d
",ans);
        }
    }
        
    scanf("
");
    return 0;
}