hdu6601 2019杭电多校第二场 Keen On Everything But Triangle

看到可以离线询问区间，瞬间写莫队，后来偷听到学弟说只要记录最大的30个，我思考了一蛤，发现应该是最大的50个，然后继续莫队，结果一直超时，还想了好多办法优化常数，竟然没想到主席树区间求第k大枚举区间最大的50个。。。。卧槽我今天讲课内容才给学弟们讲过主席树啊。。。。

其实最难的想法还是枚举三角形时，对边排序，然后从长到端枚举，假设有 a<=b<=c<=d.如果b+c>d 那么成立，否则a<=b

#include
#define maxl 200010
using namespace std;

int n,q,cnt,tot;
int a[maxl],num[maxl],rt[maxl];
long long b[110];
struct node
{
	int ls,rs,sum;
}tree[maxl*40];

inline void insert(int num,int &x,int l,int r)
{
	tree[++tot]=tree[x];x=tot;
	++tree[x].sum;
	if(l==r) return;
	int mid=(l+r)>>1;
	if(num<=mid)
		insert(num,tree[x].ls,l,mid);
	else
		insert(num,tree[x].rs,mid+1,r);
}

inline int query(int i,int j,int k,int l,int r)
{
	if(l==r)
		return l;
	int tp=tree[tree[j].ls].sum-tree[tree[i].ls].sum;
	int mid=(l+r)>>1;
	if(k<=tp) 
		return query(tree[i].ls,tree[j].ls,k,l,mid);
	else
		return query(tree[i].rs,tree[j].rs,k-tp,mid+1,r);
}

inline void prework()
{
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]),num[i]=a[i];
	sort(num+1,num+1+n);
	cnt=unique(num+1,num+1+n)-num-1;
	for(int i=1;i<=n;i++)
		a[i]=lower_bound(num+1,num+1+cnt,a[i])-num;
	for(int i=0;i<=tot;i++)
		tree[i].ls=tree[i].rs=tree[i].sum=0;
	tot=0;
	for(int i=1;i<=n;i++)
		rt[i]=0;
	for(int i=1;i<=n;i++)
	{
		rt[i]=rt[i-1];
		insert(a[i],rt[i],1,cnt);
	}
}

inline void mainwork()
{
	int l,r,up;
	long long ans=0;
	for(int i=1;i<=q;i++)
	{
		scanf("%d%d",&l,&r);
		up=min(50,r-l+1);
		for(int j=1;j<=up;j++)
			b[j]=num[query(rt[l-1],rt[r],r-l+1-j+1,1,cnt)];
		ans=-1;
		for(int j=1;j<=up-2;j++)
		if(b[j]