hdu6601 Keen On Everything But Triangle (主席树+斐波那契和三角形组成的关系)

题意：

给长度为n的数列，m个询问
每次询问给区间（L，R）
从区间内选择三个数组成三角形
问组成的三角形的最大周长
如果不能组成三角形输出-1

Sample Input

5 3
2 5 6 5 2
1 3
2 4
2 5

Sample Output

13
16
16

分析:

显然边长要最长优先选择长的边
用主席树实现选取区间第k大

假设前4长的边依次为abcd
根据两边之和大于第三边
如果b+c>a则可以组成三角形，得到答案退出
如果b+c<=a则不能组成三角形
且除a外其他任意两条边的和都不可能大于a
因此接下来考虑bcd，
直到取得答案或者没有答案

至于枚举是否会超时，
先看这个：
斐波那契和三角形组成的关系
因为最坏的情况是数列中的数全部都是斐波那契数列中的数，不能组成三角形
题目的数据范围是1e9,斐波那契到达这个数只有44项，
因此我们最多判断不超过50次
所以枚举不会超时

code：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
typedef long long ll;
using namespace std;
const int maxm=1e5+5;
struct Node{
    int lc,rc,cnt;
}a[maxm*40];
int r[maxm],cnt;
int b[maxm];
int xx[maxm],num;//离散化
int n,m;
void init(){
    cnt=0;
}
int build(int l,int r){
    int k=cnt++;
    if(l!=r){
        int mid=(l+r)/2;
        a[k].lc=build(l,mid);
        a[k].rc=build(mid+1,r);
    }
    return k;
}
int update(int x,int l,int r,int last){
    int k=cnt++;
    a[k]=a[last];
    a[k].cnt++;
    if(l!=r){
        int mid=(l+r)/2;
        if(x<=mid)a[k].lc=update(x,l,mid,a[last].lc);
        else a[k].rc=update(x,mid+1,r,a[last].rc);
    }
    return k;
}
int kth(int k,int l,int r,int x,int y){
    if(l==r)return l;
    int mid=(l+r)/2;
    int res=a[a[y].lc].cnt-a[a[x].lc].cnt;
    if(res>=k)return kth(k,l,mid,a[x].lc,a[y].lc);
    return kth(k-res,mid+1,r,a[x].rc,a[y].rc);
}
signed main(){
    while(scanf("%d%d",&n,&m)!=EOF){
        init();
        for(int i=1;i<=n;i++){
            scanf("%d",&b[i]);
            xx[i]=b[i];
        }
        sort(xx+1,xx+1+n);
        num=unique(xx+1,xx+1+n)-xx-1;
        r[0]=build(1,num);
        for(int i=1;i<=n;i++){
            int pos=lower_bound(xx+1,xx+1+num,b[i])-xx;
            r[i]=update(pos,1,num,r[i-1]);
        }
        for(int i=1;i<=m;i++){
            int L,R;
            scanf("%d%d",&L,&R);
            int len=R-L+1;//区间内数字个数
            if(len<3){
                puts("-1");
            }else{
                ll ans=-1;
                ll x,y,z;//三条边长度
                x=xx[kth(len,1,num,r[L-1],r[R])];
                y=xx[kth(len-1,1,num,r[L-1],r[R])];
                for(int i=len-2;i>=1;i--){
                    z=xx[kth(i,1,num,r[L-1],r[R])];
                    if(x-y<z){
                        ans=x+y+z;
                        break;
                    }else{
                        x=y;
                        y=z;
                    }
                }
                printf("%I64d\n",ans);
            }
        }
    }
    return 0;
}