P2354,jzoj3757-[NOI2014]随机数生成器【贪心,暴力】

正题

题目链接:https://www.luogu.com.cn/problem/P2354

---

解题思路

以 $1\sim n*m$的数字组成的 $n*m$的矩阵，求一条路径使得路径上的数排序后字典序最小。

---

解题思路

考虑从小到大依次枚举来判断是否可行。

每次判断一个点是否可以经过然后加入一个必经点。

考虑这个必经点的限制条件，一个 $(x,y)$的必经点使得你无法走到任何一个 $(a,b)(a>x,b<y)$或者 $(a,b)(a<x,b>y)$

因为必经点最多只有 $n+m+1$个，所以每行维护一个 $minn$和 $maxn$来表示该行对 $y$的限制即可，然后加入必经点时暴力修改即可。

时间复杂度 $:O(nm+(n+m+1)n)$

---

$code$

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#define lowbit(x) (x&-x)
#define t(i,j) (a[(i-1)*m+j])
using namespace std;
const int N=5100;
long long x0,a0,b0,c0,d0;
int n,m,q,k,a[N*N],num[N*N],minn[N],maxn[N];
int read() {
	int x=0,f=1; char c=getchar();
	while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
	while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();
	return x*f;
}
void print(int x){
	if (x>9) print(x/10); putchar(x%10+48); return;
}
int main()
{
	freopen("data.out","w",stdout);
	scanf("%lld%lld%lld%lld%lld",&x0,&a0,&b0,&c0,&d0);
	scanf("%d%d%d",&n,&m,&q);k=n*m;
	for(int i=1;i<=k;i++){
		x0=(a0*x0*x0+b0*x0+c0)%d0;
		a[i]=i;swap(a[i],a[x0%i+1]);
	}
	for(int i=1;i<=q;i++){
		int u,v;
		u=read();v=read();
		swap(a[u],a[v]);
	}
	for(int i=1;i<=n;i++)
		minn[i]=1,maxn[i]=m;
	for(int i=1;i<=k;i++)
		num[a[i]]=i;
	memset(a,0,sizeof(a));
	int z=0;
	for(int i=1;i<=k;i++){
		int x=(num[i]-1)/m+1,y=(num[i]-1)%m+1;
		if(y>=minn[x]&&y<=maxn[x]){
			printf("%d ",i);
			for(int j=x+1;j<=n;j++)minn[j]=max(minn[j],y);
			for(int j=1;j<x;j++)maxn[j]=min(maxn[j],y);
		}
	}
}

---

话说

考场上傻逼逼的写了个 $60pts$的二维树状数组代码就丢这了

#pragma GCC optimize(2)
%:pragma GCC optimize(3)
%:pragma GCC optimize("Ofast")
%:pragma GCC optimize("inline")
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#define lowbit(x) (x&-x)
#define t(i,j) (a[(i-1)*m+j])
using namespace std;
const int N=5100;
long long x0,a0,b0,c0,d0;
int n,m,q,k,a[N*N],num[N*N];
int read() {
	int x=0,f=1; char c=getchar();
	while(!isdigit(c)) {if(c=='-')f=-f;c=getchar();}
	while(isdigit(c)) x=(x<<1)+(x<<3)+c-48,c=getchar();
	return x*f;
}
void print(int x){
	if (x>9) print(x/10); putchar(x%10+48); return;
}
void change(int x,int y){
	for(int i=x;i<=n;i+=lowbit(i))
		for(int j=y;j<=m;j+=lowbit(j))
			t(i,j)++;
	return;
}
int ask(int x,int y){
	int ans=0;
	if(!x||!y)return 0;
	for(int i=x;i;i-=lowbit(i))
		for(int j=y;j;j-=lowbit(j))
			ans+=t(i,j);
	return ans;
}
int main()
{
	scanf("%lld%lld%lld%lld%lld",&x0,&a0,&b0,&c0,&d0);
	scanf("%d%d%d",&n,&m,&q);k=n*m;
	for(int i=1;i<=k;i++){
		x0=(a0*x0*x0+b0*x0+c0)%d0;
		a[i]=i;swap(a[i],a[x0%i+1]);
	}
	for(int i=1;i<=q;i++){
		int u,v;
		u=read();v=read();
		swap(a[u],a[v]);
	}
	for(int i=1;i<=k;i++)
		num[a[i]]=i;
	memset(a,0,sizeof(a));
	int z=0;
	for(int i=1;i<=k;i++){
		int x=(num[i]-1)/m+1,y=(num[i]-1)%m+1;
		y=m-y+1;
		if(!ask(x-1,y-1)&&(ask(n,y)+ask(x,m)-ask(x,y)==z)){
			change(x,y);z++;
			print(i);putchar(' ');
		}
	}
}