Keen On Everything But Triangle HDU 6601（主席树+斐波那契增速）

大致题意

给一个长度为n的数组，m次询问，求区间[l,r]范围内的ai不重复选取可以组成的最大的三角形的周长，不能构成则输出-1。(n<=100000,ai<=1e9)

思路

考虑枚举三角形最大的一条边，那么较小的两条边要大于该边。考虑将边的大小升序排列，那么最坏情况即a[i-2]+a[i-1]=a[i] ，增长速度为斐波那契数列，斐波那契数列第49=1e9，所以只要枚举区间内前49大的值，从大到小去凑三角形，凑出即可。

代码

贴一份主席树的代码，要离散化的，好久不用，今天写的时候漏洞百出。

#include<bits/stdc++.h>
using namespace std;
#define maxn 100005
#define maxm 1000006
#define ll long long int
#define INF 0x3f3f3f3f
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
#define mem(a) memset(a,0,sizeof(a))
#define sqr(x) (x*x)
#define inf (ll)2e18+1
#define PI acos(-1)
int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
     while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
     return f*x;
}
int n,m,root[maxn],a[maxn],b[maxn],cnt,tot,len;
struct node{int l,r,sum;}T[maxn*40];
vector<int>v;
int getid(int x){return lower_bound(v.begin(),v.begin()+len,x)-v.begin()+1;}
void update(int l,int r,int &x,int y,int pos)
{
    T[++cnt]=T[y],T[cnt].sum++,x=cnt;
    if(l==r)return ;
    int mid=(l+r)/2;
    if(mid>=pos)update(l,mid,T[x].l,T[y].l,pos);
    else update(mid+1,r,T[x].r,T[y].r,pos);
}
int query(int l,int r,int x,int y,int k)
{
    if(l==r)return l;
    int mid=(l+r)/2;
    int sum=T[T[y].l].sum-T[T[x].l].sum;
    if(sum>=k)return query(l,mid,T[x].l,T[y].l,k);
    else return query(mid+1,r,T[x].r,T[y].r,k-sum);
}
int main()
{
    while(~scanf("%d%d",&n,&m)){
        cnt=0;
        v.clear();
        inc(i,1,n){
            a[i]=read();
            v.push_back(a[i]);
        }
        sort(v.begin(),v.end());
        len=unique(v.begin(),v.end())-v.begin();
        inc(i,1,n)update(1,len,root[i],root[i-1],getid(a[i]));
        int x,y;
        inc(i,1,m){
            x=read();y=read();
            ll ans=-1;
            if(y-x+1<3)printf("-1\n");
            else{
                int pp=0;
                inc(j,1,y-x+1){
                    b[j]=v[query(1,len,root[x-1],root[y],y-x+1-j+1)-1];
                    pp++;
                    if(pp>=50)break;
                }
                inc(j,1,pp-2)if(b[j]<b[j+1]+b[j+2]){ans=1ll*b[j]+b[j+1]+b[j+2];break;}
                printf("%lld\n",ans);
            }
        }
        inc(i,0,cnt)T[i].l=T[i].r=T[i].sum=0;
    }
    return 0;
}