正题

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

---

解题思路

以$1\sim n\ast m$的数字组成的$n\ast 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(' ');}}
}