概述
正题
题目链接:https://www.luogu.com.cn/problem/P2354
解题思路
以 1 ∼ n ∗ m 1sim n*m 1∼n∗m的数字组成的 n ∗ m n*m n∗m的矩阵,求一条路径使得路径上的数排序后字典序最小。
解题思路
考虑从小到大依次枚举来判断是否可行。
每次判断一个点是否可以经过然后加入一个必经点。
考虑这个必经点的限制条件,一个 ( x , y ) (x,y) (x,y)的必经点使得你无法走到任何一个 ( a , b ) ( a > x , b < y ) (a,b)(a>x,b<y) (a,b)(a>x,b<y)或者 ( a , b ) ( a < x , b > y ) (a,b)(a<x,b>y) (a,b)(a<x,b>y)
因为必经点最多只有 n + m + 1 n+m+1 n+m+1个,所以每行维护一个 m i n n minn minn和 m a x n maxn maxn来表示该行对 y y y的限制即可,然后加入必经点时暴力修改即可。
时间复杂度 : O ( n m + ( n + m + 1 ) n ) :O(nm+(n+m+1)n) :O(nm+(n+m+1)n)
c o d e code 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);
}
}
}
话说
考场上傻逼逼的写了个 60 p t s 60pts 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(' ');
}
}
}
最后
以上就是简单西装为你收集整理的P2354,jzoj3757-[NOI2014]随机数生成器【贪心,暴力】正题的全部内容,希望文章能够帮你解决P2354,jzoj3757-[NOI2014]随机数生成器【贪心,暴力】正题所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复