记录编号 |
607404 |
评测结果 |
AAAAAAAAAAAAAAAAAAAAAAAAA |
题目名称 |
[统一省选 2025]追忆 |
最终得分 |
100 |
用户昵称 |
李金泽 |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
78.379 s |
提交时间 |
2025-10-14 18:31:39 |
内存使用 |
753.63 MiB |
显示代码纯文本
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#define N 100005
#define M 317
#define S 1563
#define ul unsigned long long
#define ri register int
using namespace std;
int T,n,ns,m,q,s=350,w=64,wh[N+80],h[N],cnt,in[N],a[N],b[N],ra[N],rb[N],op,x,y,k,l,r,ans;
struct edge{int v,n;}e[N<<1];
void ad(int u,int v){e[++cnt]={v,h[u]};h[u]=cnt;in[v]++;}
void swap(int &x,int &y){int t=x;x=y;y=t;}
int min(int x,int y){return x<y?x:y;}
struct bitset{
ul t[S];
void clear(){for(ri i=0;i<S;i++)t[i]=0;}
void ud(int i,bool k)
{
int x=wh[i];i-=x*w;
if(k)t[x]|=((ul)1<<i);
else t[x]&=~((ul)1<<i);
}
ul gs(int i)
{
int x=wh[i];i-=x*w;
return t[x]&((ul)1<<i);
}
void operator=(bitset y)
{
for(ri i=0;i<ns;i++)t[i]=y.t[i];
}
void operator|=(bitset y)
{
for(ri i=0;i<ns;i++)t[i]|=y.t[i];
}
void operator&=(bitset y)
{
for(ri i=0;i<ns;i++)t[i]&=y.t[i];
}
void operator^=(bitset y)
{
for(ri i=0;i<ns;i++)t[i]^=y.t[i];
}
void output()
{
for(ri i=0;i<n;i++)printf("%d",gs(i)?1:0);
}
}vis[N],v[M],c[M],R;
void topu()
{
queue<int>q;
for(int i=0;i<n;i++)if(!in[i])q.push(i);
while(q.size())
{
int u=q.front();q.pop();
for(ri i=h[u];i;i=e[i].n)
{
int v=e[i].v;
vis[v]|=vis[u];
if(!(--in[v]))q.push(v);
}
}
}
int read(){
int sum=0;bool f=0;char c=getchar();
for(;c<48||c>57;c=getchar())if(c==45)f=1;
for(;c>=48&&c<=57;c=getchar())sum=sum*10+(c&15);
return f?-sum:sum;
}
int main(){
freopen("recall.in","r",stdin);freopen("recall.out","w",stdout);
for(ri i=0;i<N+80;i++)wh[i]=i/w;
read();T=read();
while(T--)
{
memset(h,0,sizeof(h));memset(in,0,sizeof(in));
cnt=0;for(ri i=0;i<M;i++)v[i].clear(),c[i].clear();
n=read();m=read();q=read();ns=wh[(n+63)];
while(m--)x=read()-1,y=read()-1,ad(y,x);
for(ri i=0;i<n;i++)a[i]=read()-1,ra[a[i]]=i;
for(ri i=0;i<n;i++)b[i]=read()-1,rb[b[i]]=i;
for(ri i=0;i<n;i++)vis[i].ud(i,1);
topu();
for(ri i=0;i<n;i++)
for(ri j=0;j*s<n;j++)
{
if(a[i]>=j*s)v[j].ud(i,1);
if(b[i]>=j*s)c[j].ud(i,1);
}
while(q--)
{
op=read();
switch(op)
{
case 1:{
x=read()-1;y=read()-1;
if(a[x]>a[y])swap(x,y);
for(ri i=0;i*s<n;i++)
if(a[x]<i*s&&i*s<=a[y])
v[i].ud(x,1),v[i].ud(y,0);
swap(ra[a[x]],ra[a[y]]);
swap(a[x],a[y]);
break;
}
case 2:{
x=read()-1;y=read()-1;
if(b[x]>b[y])swap(x,y);
for(ri i=0;i*s<n;i++)
if(b[x]<i*s&&i*s<=b[y])
c[i].ud(x,1),c[i].ud(y,0);
swap(rb[b[x]],rb[b[y]]);
swap(b[x],b[y]);
break;
}
case 3:{
k=read()-1;l=read()-1;r=read()-1;R.clear();
//计算l<=ay<=r的y集合
//计算l
x=(l+s-1)/s;
if(x*s<n)R=v[x];
for(ri i=l;i<min(n,x*s);i++)R.ud(ra[i],1);
//计算r
x=r/s+1;
if(x*s<n)R^=v[x];
for(ri i=r+1;i<min(n,x*s);i++)R.ud(ra[i],0);
//合并
R&=vis[k];
//计算max(by)
int p=0,q=0;ans=-1;
for(;(p+1)*s<n;p++)
{
for(;q<ns&&!(R.t[q]&c[p+1].t[q]);q++);
if(q>=ns)break;
}
for(int i=p*s;i<min(n,(p+1)*s);i++)
if(R.gs(rb[i]))
ans=i;
printf("%d\n",ans+1);
break;
}
}
}
for(ri i=0;i<n;i++)vis[i].clear();
}
return 0;
}