2017 多校训练第一场 HDU 6040 Hints of sd0061

题意：

给你一个长度为n的序列，序列由题面给的函数生成。然后m次询问，询问这个序列上第bi小的数。

分析：

n最大达到1e7，单纯的快排是不行的。但是可以利用快排的思想，每次找到一个划分，序列中比划分值小的数都在划分的左侧，反之则在右侧，这样每阶段都可以省一段时间。题面给的函数返回类型是unsigned，这个也需要注意。

这就需要一个利器的STL了，C++中的nth_element(arr,arr+k,arr+n)，用法是将长度为n的数组arr进行划分，第k位置上就是第k+1大的数（下标从0开始算，也即有k个数比arr[k]小），最重要的是这个函数近似线性！！！

注意：
因为输入保证任意两个小的之和小于第三个
所以查询数列的间隔一定大于等于斐波那契
也就是从大到小查询的话，每次至少能去掉一半的区间

#include <bits/stdc++.h>
#define test TEST 
using namespace std;

typedef long long ll;
const int maxn=1e7+7;
int b[105],id[105];
unsigned ans[105];

bool vis[maxn];
unsigned a[maxn];
unsigned x,y,z;
unsigned rng61() {
	unsigned t;
	x ^= x << 16;
	x ^= x >> 5;
	x ^= x << 1;
	t = x;
	x = y;
	y = z;
	z = t ^ x ^ y;
	return z;
}

bool cmp(const int &xx,const int &yy){
	return b[xx]>b[yy];
}
int main(int argc, char const *argv[])
{
	ll n,m;
	int cas=1;
	while(~scanf("%lld %lld %u %u %u",&n,&m,&x,&y,&z)){
		memset(vis,0,sizeof(vis));
		for(int i=0;i<m;i++){
			scanf("%d",&b[i]);
			id[i]=i;
		}
		for(int i=0;i<n;i++){
			a[i]=rng61();
		}
		sort(id,id+m,cmp);//挑选询问大的编号放在前面,可以降低复杂度
		int p,last=n;
		for(int i=0;i<m;i++){
			p=id[i];
			nth_element(a,a+b[p],a+last);
			ans[p]=a[b[p]];
			last=b[p];
		}
		//输出答案
		printf("Case #%d:",cas++);
		for(int i=0;i<m;i++){
			printf(" %u",ans[i]);
		}printf("\n");
	}
	return 0;
}