【BZOJ3672】【UOJ#6】【NOI2014】随机数生成器

时间:2024-08-22 20:34:26

暴力出奇迹

原题:

【BZOJ3672】【UOJ#6】【NOI2014】随机数生成器

2≤N,M≤5000

0≤Q≤50000

0≤a≤300

0≤b,c≤108

0≤x0<d≤108

1≤ui,vi≤N×M

恩首先容易看出来这个棋盘直接模拟搞出来就行了,不用反演矩阵乘之类的奇怪的东西

然后又容易发现只需要遍历从1~n*m的数尽量答案里塞就是最优答案 = =|||

然后贪心搞一下,从1~n*m遍历,对于每一个n记录一个top和一个bottom表示第i行能取第bottom[i]到top[i]列的数

容易发现暴力维护的复杂度是资瓷的 = =|||

易证不会存在因为之前钦定选择某些数导致后面选不够n+m-1个数的情况

然后暴力搞就可以辣

因为2.5e7的数组只能开两个所以有一个要重复使用……

注意答案的长度是2n…………

注意一开始递推的姿势,比如两个int和一个longlong相乘要先int乘longlong再乘int之类的

代码:

 #include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
#define ll long long
int rd(){int z=; char ch=getchar();
while(ch<''||ch>'') ch=getchar();
while(ch>=''&&ch<=''){z=(z<<)+(z<<)+ch-''; ch=getchar();}
return z;
}
ll a,b,c,d,n,m,o; int x[]; int nm;
int T[],tp[],bt[];
int ans[],ast=;
int main(){//freopen("ddd.in","r",stdin);
int l,r;
cin>>x[]>>a>>b>>c>>d>>n>>m>>o; nm=n*m;
for(int i=;i<=nm;++i) x[i]=((((x[i-]*a)%d)*x[i-])%d+(x[i-]*b)%d+c)%d;
for(int i=;i<=nm;++i) T[i]=i;
for(int i=;i<=nm;++i) swap(T[i],T[x[i]%i+]);
while(o--) l=rd(),r=rd(),swap(T[l],T[r]);
for(int i=;i<=nm;++i) x[T[i]]=i-;
for(int i=;i<=n;++i) tp[i]=m,bt[i]=;
for(int i=;i<=nm;++i){
l=x[i]/m+,r=x[i]%m+;
if((r>=bt[l])&(r<=tp[l])){
for(int j=;j<l;++j)if(r<tp[j]) tp[j]=r;
for(int j=l+;j<=n;++j)if(r>bt[j]) bt[j]=r;
ans[++ast]=i;
if(ast==n+m-) break;
}
}
for(int i=;i<ast;++i) printf("%d ",ans[i]);
cout<<ans[ast];
return ;
}