HDU6601 Keen On Everything But Triangle（主席树模板题）

题目

N sticks are arranged in a row, and their lengths are a1,a2,…,aN.

There are Q querys. For i-th of them, you can only use sticks between li-th to ri-th. Please output the maximum circumference of all the triangles that you can make with these sticks, or print −1 denoting no triangles you can make.

题目大意

给你一排棍，q次询问，问第l个棍到第r个棍可以组成的最大周长的三角形

解析

主席树的模板题，需要查询区间第k大。

首先需要知道这样一件事：在一个特定的序列中，想找最大周长的三角形一定是，他的三条边首先可能为 （第一大，第二大，第三大），如果这样无法构成三角形，那么第一大边就可以舍弃，考虑（第二大，第三大，第四大）。因为第一大与之后的任意边都一定无法组成三角形。

那这样扫下去，循环次数会很大吗？答案是不会的，循环次数最多的构造也就类似于斐波那契数列这样的，对于本题的范围来说，循环次数都不会超过50次

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long LL;
const int MAXN=400005;
inline int read(){
	int x=0; char c=getchar();
	while(c<'0'||c>'9') c=getchar();
	while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
	return x;
}

int cnt;
int lc[MAXN*20],rc[MAXN*20],sum[MAXN*20],root[MAXN];
int build(int l,int r){
	int x=++cnt,mid=l+r>>1;
	sum[x]=0;
	if(l<r){
		lc[x]=build(l,mid);
		rc[x]=build(mid+1,r);
	}
	return x;
}
int update(int pre,int l,int r,int val){
	int x=++cnt;
	lc[x]=lc[pre]; rc[x]=rc[pre]; sum[x]=sum[pre]+1;
	if(l<r){
		int mid=l+r>>1;
		if(val<=mid) lc[x]=update(lc[pre],l,mid,val);
		else rc[x]=update(rc[pre],mid+1,r,val);
	}
	return x;
}
inline int query(int u,int v,int l,int r,int k){
	if(l==r) return l;
	int x=sum[rc[v]]-sum[rc[u]];
	int mid=l+r>>1;
	if(x>=k) return query(rc[u],rc[v],mid+1,r,k);
	else return query(lc[u],lc[v],l,mid,k-x);
}

int a[MAXN],b[MAXN];
int main(){
	int n,T;
	while(scanf("%d%d",&n,&T)==2){
		cnt=0;
		for(int i=1;i<=n;i++)
			b[i]=a[i]=read();

		sort(b+1,b+n+1);
		int m=unique(b+1,b+n+1)-(b+1);
		root[0]=build(1,m);
		for(int i=1;i<=n;i++){
			a[i]=lower_bound(b+1,b+m+1,a[i])-b;
			root[i]=update(root[i-1],1,m,a[i]);
		}

		while(T--){
			int l=read(),r=read();
			if(r-l<2) printf("-1\n");
			else{
				LL q1=query(root[l-1],root[r],1,m,1),
					q2=query(root[l-1],root[r],1,m,2);
				bool tag=false;
				for(int i=3;i<=r-l+1;i++){
					LL q3=query(root[l-1],root[r],1,m,i);
					if(b[q3]+b[q2]>b[q1]){
						printf("%lld\n",b[q1]+b[q2]+b[q3]);
						tag=true; break;
					}
					q1=q2; q2=q3;
				}
				if(!tag) printf("-1\n");
			}
		}
	}
	return 0;
}