BZOJ 3671 NOI2014 随机数生成器 贪心+暴力

题目大意：。。。。。有点长自己看吧

首先既然是排序后的序列字典序最小，那么一定要选尽量小的数字走

然后T是1~m*n的序列 所以不存在重复（一开始我居然把这个条件看漏了）

好的这题贪心
每次选择没有被标记的最小点，然后把左下方和右上方都标记掉（记得标记重复时break，不然就挂了）

注意5000*5000的数组开两个int就是极限了 开多了妥妥MLE 所以T数组记得重复利用

暴力跑了38秒 不过这题每一行能够选择的区域一定是连续的 可以对于每一行维护一个l和r 每次更新取最值即可 这个代码跑了28秒 还是好慢QAQ

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 5010
using namespace std;
typedef long long ll;
int T[M*M],map[M][M],l[M],r[M];
ll seed,a,b,c,d;
int n,m,q;
inline int Rand()
{
	return seed=(a*seed*seed%d+b*seed%d+c)%d;
}
void kill(int x,int y)
{
	int i;
	for(i=1;i<x;i++)
		r[i]=min(r[i],y);
	for(i=x+1;i<=m;i++)
		l[i]=max(l[i],y);
}
int main()
{
	int i,j,k,x,y,flag=0;
	cin>>seed>>a>>b>>c>>d;
	cin>>m>>n>>q;
	for(i=1;i<=m;i++)
		r[i]=n;
	for(i=1;i<=m*n;i++)
		T[i]=i,swap(T[i],T[Rand()%i+1]);
	for(i=1;i<=q;i++)
	{
		scanf("%d%d",&x,&y);
		swap(T[x],T[y]);
	}
	k=0;
	for(i=1;i<=m;i++)
		for(j=1;j<=n;j++)
			map[i][j]=T[++k];
	for(i=1;i<=m;i++)
		for(j=1;j<=n;j++)
			T[map[i][j]]=i<<16|j;
	for(i=1;i<=m*n;i++)
	{
		x=T[i]>>16;
		y=T[i]&65535;
		if(y>=l[x]&&y<=r[x])
		{
			if(flag)
				putchar(' ');
			flag=1;
			printf("%d",map[x][y]);
			kill(x,y);
		}
	}
}