P7520 [省选联考 2021 A 卷] 支配

题目

题目

思路

先对原图建支配树（题目名提示），建议使用nlogn，对常数大，厌氧选手留一点空间余地（当然要是太难打n²也可以，就是要卡常（开了O2都不叫事））。
剩下我们慢慢考虑。
由数据范围得，查询单次复杂度最高O(n)（预处理O(xn²+ynlogn)（其中x，y为常数）情况下，为什么n²一会讲），所以我们需要O(n)查询。
接下来我们研究单次查询。
每次查询给的边为(u,v)，由于是O(n)，所以显然是每个点都扫一遍，然后O(1)判断是否有贡献。
要想O(1)判断，显然在不预处理的情况下不可能（除非宁用随机数等算法并且宁是欧皇）。
所以我们要预处理。
一个点的受支配集只会往小变（显然）。
一个点（设为x）的受支配集，显然只有在1到x的新路径跳过了一个支配x的点的情况下才会变。
1到x的新路径跳过了一个支配x的点，我们设这个支配点为k，那么这个表述可以变为当且仅当有3条路径1→u，u→v，v→x均不经过k且存在，这个点x才会做出贡献。
1→u可以预处理成$n\ast n$的方阵A，$A_{i,j}=1$表示1可以在不经过j的情况下到i，时间复杂度O(n²)。
u→v显然成立。
v→x可以预处理成$n\ast n$的方阵B，$B_{i,j}=1$表示i可以在不经过j支配树上的父节点的情况下到j，时间复杂度O(n²)（这里常数会比上面的O(n²)小一些）。
这样我们就做到了O(n²)预处理，O(n)单次查询，卡卡常可在不开O2的情况下AC本题（有可能是因为本人常数大）。
code:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
//When I wrote this code,God and I unterstood what was I doing 
inline int read()
{
    int ret,c,f=1;
    while (((c=getchar())> '9'||c< '0')&&c!='-');
    if (c=='-') f=-1,ret=0;
    else ret=c-'0';
    while ((c=getchar())>='0'&&c<='9') ret=(ret<<1)+(ret<<3)+(c^48);
    return ret*f;
}
inline void write(int x)
{
	if (x<10)
	{
		putchar(x+'0');
		return; 
	}
	write(x/10);
	putchar(x%10+'0');
	return; 
}
struct fP{
	int to[6010],nxt[6010];
	int head[3003],tot;
	void reset()
	{
		memset(head,0,sizeof(head));
		tot=0;
		return;
	}
	void add(int x,int y)
	{
		to[++tot]=y;
		nxt[tot]=head[x];
		head[x]=tot;
		return;
	}
} wj,jw,kyx;
int n,q,m,id[3003],dfn[3003],dep[3003],fa[3003],f[3003],semi[3003],idom[3003],st[20][3003],lyw[3003],bq[3003],ans[3003],x,y,tot;
bool zp[3003][3003],zp2[3003][3003],wjj[3003];
void Tarjan(int x)
{
	dfn[x]=++tot;
	lyw[tot]=x;
	for (int i=wj.head[x];i;i=wj.nxt[i])
	{
		if (!dfn[wj.to[i]])
		{
			fa[wj.to[i]]=x;
			Tarjan(wj.to[i]);
		}
	}
	return;
}
int find(int x)
{
	if (x==f[x]) return x;
	int y=find(f[x]);
	if (dfn[semi[bq[f[x]]]]<dfn[semi[bq[x]]]) bq[x]=bq[f[x]];
	return f[x]=y;
}
void dfs(int x,int y)
{
	if (idom[y]==x) return;
	zp[y][x]=1;
	for (int i=jw.head[x];i;i=jw.nxt[i])
	{
		if (!zp[y][jw.to[i]]) dfs(jw.to[i],y);
	}
	return;
}
void df2(int x,int y)
{
	zp2[y][x]=1;
	if (x==y)
	{
		return;
	}
	for (int i=wj.head[x];i;i=wj.nxt[i])
	{
		if (!zp2[y][wj.to[i]]) df2(wj.to[i],y);
	}
	return;
}
int main()
{
	n=read(),m=read(),q=read();
	for (int i=1;i<=m;i++)
	{
		x=read(),y=read();
		wj.add(x,y),jw.add(y,x);
	}
	Tarjan(1);
	for (int i=1;i<=n;i++) bq[i]=f[i]=semi[i]=i;
	df2(1,0);
	for (int i=1;i<=n;i++) df2(1,i);
	for (int i=tot;i>=2;i--)
	{
		for (int j=jw.head[lyw[i]];j;j=jw.nxt[j])
		{
			if (dfn[jw.to[j]]==0) continue;
			find(jw.to[j]);
			if (dfn[semi[bq[jw.to[j]]]]<dfn[semi[lyw[i]]]) semi[lyw[i]]=semi[bq[jw.to[j]]];
		}
		f[lyw[i]]=fa[lyw[i]];
		kyx.add(semi[lyw[i]],lyw[i]);
		for (int j=kyx.head[fa[lyw[i]]];j;j=kyx.nxt[j])
		{
			find(kyx.to[j]);
			if (semi[bq[kyx.to[j]]]==fa[lyw[i]]) idom[kyx.to[j]]=fa[lyw[i]];
			else idom[kyx.to[j]]=bq[kyx.to[j]];
		}
		kyx.head[fa[lyw[i]]]=0;
	}
	for (int i=2;i<=tot;i++)
	{
		if (idom[lyw[i]]!=semi[lyw[i]]) idom[lyw[i]]=idom[idom[lyw[i]]];
	}
	for (int i=2;i<=n;i++)
	{
		dfs(i,i);
	}
	while (q--)
	{
		x=read(),y=read();
		int s=0;
		for (int i=1;i<=n;i++)
		{
			if (idom[i]!=1&&zp2[idom[i]][x]&&x!=idom[i]&&zp[i][y]) wjj[i]=1;
		}
		for (int i=2;i<=tot;i++)
		{
			if (wjj[idom[lyw[i]]]) wjj[lyw[i]]=1;
		}
		for (int i=1;i<=n;i++)
		{
			if (wjj[i])
			{
				s++;
				wjj[i]=0;
			}
		}
		write(s);
		puts("");
	}
	return 0;
}
//Now,only God know