【BZOJ3671】【NOI2014】随机数生成器

暴力即正解。

#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 5010
#define NN 25001000
using namespace std;
struct Fiona
{
	short x,y;
	int w;
}s[NN];
int acl[N],acr[N];
long long seed,a,b,c,p,n,m,q;
int rdm()
{
	return seed=(a*seed*seed+b*seed+c)%p;
}

bool cmp(Fiona A,Fiona B)
{
	return A.w<B.w;
}

int main()
{
//	freopen("test.in","r",stdin);
	int i,j,k,cnt;
	int x,y;
	scanf("%lld%lld%lld%lld%lld%lld%lld%lld",&seed,&a,&b,&c,&p,&n,&m,&q);
	for(cnt=i=0;i<n;i++)
	{
		for(j=0;j<m;j++)
		{
			s[++cnt].w=cnt;
			s[cnt].x=i;
			s[cnt].y=j;
		}
		acr[i]=m-1;
	}
	for(i=1;i<=cnt;i++)
	{
		swap(s[i].w,s[rdm()%i+1].w);
	}
	for(i=1;i<=q;i++)
	{
		scanf("%d%d",&x,&y);
		if(x==y)continue;
		swap(s[x].w,s[y].w);
	}
	sort(s+1,s+cnt+1,cmp);
	for(cnt=i=1;i<n+m-1;i++,cnt++)
	{
		x=s[cnt].x;
		y=s[cnt].y;
		while(y<acl[x]||y>acr[x])
		{
			cnt++;
			x=s[cnt].x;
			y=s[cnt].y;
		}
		printf("%d ",s[cnt].w);
		for(j=0;j<x;j++)
		{
			acr[j]=min(acr[j],y);
		}
		for(j=x+1;j<n;j++)
		{
			acl[j]=max(acl[j],y);
		}
	}
		x=s[cnt].x;
		y=s[cnt].y;
		while(y<acl[x]||y>acr[x])
		{
			cnt++;
			x=s[cnt].x;
			y=s[cnt].y;
		}
		printf("%d",s[cnt].w);
	return 0;
}