【HDOJ6601】Keen On Everything But Triangle（主席树）

题意：给定一个长为n的序列，有q次询问，每次询问[l,r]这段区间内挑三个数，能组成的三角形的最大周长，无解输出-1

n,q<=1e5，a[i]<=1e9

思路：题解写法和我的不太一样

先说题解做法，显然最坏情况下是斐波那契数列的形式，大概是log2(1e9)项就没有-1了，所以维护一个可以取某一段中前50大的数字的数据结构，取出来之后sort看相邻3个取max即可

我的做法是用主席树把区间权值抠出来，树上二分权值，如果当前结点的右儿子的size>=3就必定是最优解，否则把右儿子中的所有数取出来再递归左儿子，最后同sort处理

原理应该差不多

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 typedef unsigned int uint;
 typedef unsigned long long ull;
 typedef pair<int,int> PII;
 typedef pair<ll,ll> Pll;
 typedef vector<int> VI;
 #define N  110000
 #define M  1100000
 #define fi first
 #define se second
 #define MP make_pair
 #define pi acos(-1)
 #define mem(a,b) memset(a,b,sizeof(a))
 #define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
 #define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
 #define lowbit(x) x&(-x)
 #define Rand (rand()*(1<<16)+rand())
 #define id(x) ((x)<=B?(x):m-n/(x)+1)
 #define ls p<<1
 #define rs p<<1|1

 const ll MOD=,inv2=(MOD+)/;
       double eps=1e-;
       ll INF=1e14;

 struct arr
 {
     int l,r,s;
 }t[N*];

 int c[N],root[N],flag,m,cnt;
 ll ans;

 int read()
 {
    int v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 void update(int l,int r,int x,int &p)
 {
     t[++cnt]=t[p];
     p=cnt;
     t[p].s++;
     if(l==r) return;
     int mid=(l+r)>>;
     if(x<=mid) update(l,mid,x,t[p].l);
      else update(mid+,r,x,t[p].r);
 }

 int find(int l,int r,int x,int p1,int p2)
 {
     //printf("find l=%d r=%d\n",l,r);
     int tmp=t[p2].s-t[p1].s;
     if(x<=) return ;
     if(!tmp) return ;
     if(l==r)
     {
         rep(i,,min(,tmp)) c[++m]=l;
         return min(,tmp);
     }
     int mid=(l+r)>>;
     x-=find(mid+,r,x,t[p1].r,t[p2].r);
     find(l,mid,x,t[p1].l,t[p2].l);
     return tmp;
 }

 void query(int l,int r,int p1,int p2)
 {
     if(l>r) return;
     if(l==r)
     {
         rep(i,,min(,t[p2].s-t[p1].s)) c[++m]=l;
         return;
     }
     int mid=(l+r)>>;
     int tmp=t[t[p2].r].s-t[t[p1].r].s;
     if(tmp>=)
     {
         find(mid+,r,,t[p1].r,t[p2].r);
         return;
     }
      else
      {
          if(tmp) find(mid+,r,tmp,t[p1].r,t[p2].r);
          query(l,mid,t[p1].l,t[p2].l);
      }
 }

 int main()
 {
     int n,q;
     while(scanf("%d%d",&n,&q)!=EOF)
     {
         rep(i,,cnt) t[i].l=t[i].r=t[i].s=;
         cnt=; root[]=;
         rep(i,,n)
         {
             int x=read();
             root[i]=root[i-];
             update(,1e9,x,root[i]);
         }
         while(q--)
         {
             int x=read(),y=read();
             ans=;
             m=;
             flag=;
             query(,1e9,root[x-],root[y]);
             //printf("flag=%d\n",flag);
             sort(c+,c+m+);
             //printf("m=%d\n",m);
             //rep(i,1,m) printf("%d\n",c[i]);
             rep(i,,m-)
              if(c[i]+c[i-]>c[i+]) ans=max(ans,0ll+c[i-]+c[i]+c[i+]);

              if(ans==) printf("-1\n");
               else printf("%I64d\n",ans);
              rep(i,,m) c[i]=;
         }

     }
     return ;
 }