记录编号 607404 评测结果 AAAAAAAAAAAAAAAAAAAAAAAAA
题目名称 [统一省选 2025]追忆 最终得分 100
用户昵称 Gravatar李金泽 是否通过 通过
代码语言 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;
}