洛谷 P7520 - [省选联考 2021 A 卷] 支配(支配树)
真·支配树不 sb 的题。
首先题面已经疯狂暗示咱们建出支配树对吧,那咱就老老实实建呗。由于这题数据范围允许 \(n^2\) 算法通过,因此可以考虑 \(\mathcal O(n^2)\) 地建立支配树,具体来说我们枚举每个点 \(x\),将这个点暂时地从图中删除,如果对于图中另一个点 \(y\) 满足删除 \(x\) 后 \(1\) 不能到达 \(y\),那么 \(x\) 就在 \(y\) 的支配集中,这样我们再对整个 DAG DFS 一遍求出每个点的 DFS 序,然后取 DFS 序最大的点作为每个点在支配树上的父亲即可。
接下来考虑怎样计算答案。首先显然的一件事情是,我们加入一条边后最多只会让某些点的支配集大小变小,而不会使支配集大小变大,因此我们只需考虑有哪些点在加入这条边后,存在某个点原来能支配它而现在不能即可。注意到一个性质,就是对于一个点 \(x\),如果 \(fa_x\)(当然有些人喜欢称这个东西为 \(idom_x\),反正能看懂就行了吧)的支配集改变,那么 \(x\) 的支配集也会改变。因为根据支配树的性质,每个点的支配集一定是该点到 \(1\) 路径上所有点组成的集合,因此如果 \(fa_x\) 的支配集改变,就必然存在某个 \(fa_x\) 的祖先 \(y\),满足存在路径 \(1\to fa_x\) 且不经过 \(y\),这样就存在路径 \(1\to fa_x\to x\) 不经过 \(x\),\(y\) 就从 \(x\) 的支配集中消失了。同理,如果 \(fa_x\) 的支配集没变,但 \(x\) 的支配集改变,必然是因为 \(fa_x\) 无法支配 \(x\),因为如果存在某个 \(x\) 的祖先 \(y\ne fa_x\),满足 \(y\) 不再支配 \(x\) 且 \(y\) 能支配 \(fa_x\),那就能推出这个 \(1\to x\) 且不经过 \(y\) 的路径肯定不经过 \(fa_x\),从而 \(fa_x\) 不支配 \(x\)。
因此我们考虑每次询问对整棵树进行 DFS,如果走到一个点发现 \(fa_x\) 不支配 \(x\),答案就加上 \(x\) 子树的大小并 return
,那么怎么判断 \(fa_x\) 是否支配 \(x\) 呢?显然如果 \(fa_x\) 不支配 \(x\) 那么必然存在路径 \(1\to u\to v\to x\) 满足这条路径不经过 \(fa_x\),而这又 obviously 等价于 \(1\to u,v\to x\) 均不经过 \(fa_x\),前者可以通过建支配树时预处理出的“删掉点 \(x\) 后是否存在 \(1\to y\) 的路径的数组 \(ban_{x,y}\)”求出,而关于后者我们发现都是形如”删掉 \(fa_x\) 后 \(x\) 能否在反图上到达 \(y\)“,因此我们再建一个 \(ban\_fa_{x,y}\) 维护这个东西即可。时间复杂度 \(\mathcal O(n^2+nq)\)
卡常技巧:交换 \(ban\) 和 \(ban\_fa\) 的两维后,效率大约能快 25%,具体原理见这儿。
const int MAXN=3e3;
const int MAXM=MAXN<<1;
int n,m,qu;
struct graph{
int hd[MAXN+5],nxt[MAXM+5],to[MAXM+5],ec=0;
void adde(int u,int v){to[++ec]=v;nxt[ec]=hd[u];hd[u]=ec;}
} g,rv_g,dt;
int dfn[MAXN+5],rid[MAXN+5],tim=0,fa[MAXN+5];
bool ban[MAXN+5][MAXN+5],ban_fa[MAXN+5][MAXN+5];
void dfs(int x){
rid[dfn[x]=++tim]=x;
for(int e=g.hd[x];e;e=g.nxt[e]){
int y=g.to[e];if(!dfn[y]) dfs(y);
}
}
void dfs_ban(int x,int ban_id){
if(x==ban_id) return;ban[x][ban_id]=1;
// printf("{%d,%d}\n",ban_id,x);
for(int e=g.hd[x];e;e=g.nxt[e]){
int y=g.to[e];
if(!ban[y][ban_id]) dfs_ban(y,ban_id);
}
}
void dfs_ban_fa(int x,int ban_id){
if(x==fa[ban_id]) return;ban_fa[x][ban_id]=1;
for(int e=rv_g.hd[x];e;e=rv_g.nxt[e]){
int y=rv_g.to[e];
if(!ban_fa[y][ban_id]) dfs_ban_fa(y,ban_id);
}
}
int siz[MAXN+5];
void dfssiz(int x){
siz[x]=1;
for(int e=dt.hd[x];e;e=dt.nxt[e]){
int y=dt.to[e];dfssiz(y);
siz[x]+=siz[y];
}
}
int X,Y,res=0;
void dfscalc(int x){
if(x^1){
if(ban[X][fa[x]]&&ban_fa[Y][x]){
res+=siz[x];return;
}
} for(int e=dt.hd[x];e;e=dt.nxt[e]){
int y=dt.to[e];dfscalc(y);
}
}
int main(){
scanf("%d%d%d",&n,&m,&qu);
for(int i=1,u,v;i<=m;i++){
scanf("%d%d",&u,&v);
g.adde(u,v);rv_g.adde(v,u);
} dfs(1);for(int i=1;i<=n;i++) dfs_ban(1,i);
// for(int i=1;i<=n;i++) printf("%d %d\n",dfn[i],rid[i]);
for(int i=1;i<=n;i++) for(int j=1;j<=n;j++)
if(!ban[j][i]&&(i^j)) chkmax(fa[j],dfn[i]);
for(int i=2;i<=n;i++) fa[i]=rid[fa[i]];fa[1]=0;
for(int i=2;i<=n;i++) dt.adde(fa[i],i),dfs_ban_fa(i,i);
// for(int i=2;i<=n;i++) printf("%d\n",fa[i]);
dfssiz(1);
while(qu--){scanf("%d%d",&X,&Y);res=0;dfscalc(1);printf("%d\n",res);}
return 0;
}
洛谷 P7520 - [省选联考 2021 A 卷] 支配(支配树)的更多相关文章
- 洛谷 P7515 - [省选联考 2021 A 卷] 矩阵游戏(差分约束)
题面传送门 emmm--怎么评价这个题呢,赛后学完差分约束之后看题解感觉没那么 dl,可是现场为啥就因为种种原因想不到呢?显然是 wtcl( 先不考虑"非负"及" \(\ ...
- 洛谷 P7516 - [省选联考 2021 A/B 卷] 图函数(Floyd)
洛谷题面传送门 一道需要发现一些简单的性质的中档题(不过可能这道题放在省选 D1T3 中偏简单了?) u1s1 现在已经是 \(1\text{s}\) \(10^9\) 的时代了吗?落伍了落伍了/ ...
- 洛谷P6623——[省选联考 2020 A 卷] 树
传送门:QAQQAQ 题意:自己看 思路:正解应该是线段树/trie树合并? 但是本蒟蒻啥也不会,就用了树上二次差分 (思路来源于https://www.luogu.com.cn/blog/dengy ...
- 洛谷 P6624 - [省选联考 2020 A 卷] 作业题(矩阵树定理+简单数论)
题面传送门 u1s1 这种题目还是相当套路的罢 首先看到 \(\gcd\) 可以套路地往数论方向想,我们记 \(f_i\) 为满足边权的 \(\gcd\) 为 \(i\) 的倍数的所有生成树的权值之和 ...
- [题解] LOJ 3300 洛谷 P6620 [省选联考 2020 A 卷] 组合数问题 数学,第二类斯特林数,下降幂
题目 题目里要求的是: \[\sum_{k=0}^n f(k) \times X^k \times \binom nk \] 这里面出现了给定的多项式,还有组合数,这种题目的套路就是先把给定的普通多项 ...
- [省选联考 2021 A 卷] 矩阵游戏
很巧妙的一个构造. 我是没有想到的. 自己的思维能力可能还是不足. 考虑先满足\(b\)对\(a\)的限制,把\(a\)的第一行和第一列设\(0\),推出这个\(a\). 接下来考虑对这个\(a\), ...
- [省选联考 2021 A/B 卷] 卡牌游戏
垃圾福建垫底选手来看看这题. 大家怎么都写带 \(log\) 的. 我来说一个线性做法好了. 那么我们考虑枚举 \(k\) 作为翻转完的最小值. 那么构造出一个满足条件的操作,我们在 \(a_i\) ...
- [省选联考 2020 A 卷] 组合数问题
题意 [省选联考 2020 A 卷] 组合数问题 想法 自己在多项式和数论方面还是太差了,最近写这些题都没多少思路,看完题解才会 首先有这两个柿子 \(k*\dbinom{n}{k} = n*\dbi ...
- luoguP6623 [省选联考 2020 A 卷] 树(trie树)
luoguP6623 [省选联考 2020 A 卷] 树(trie树) Luogu 题外话: ...想不出来啥好说的了. 我认识的人基本都切这道题了. 就我只会10分暴力. 我是傻逼. 题解时间 先不 ...
随机推荐
- 力扣 - 剑指 Offer 53 - I. 在排序数组中查找数字 I
题目 剑指 Offer 53 - I. 在排序数组中查找数字 I 思路1 一般来说,首先想到的是使用一个变量,从头开始遍历整个数组,记录target数组出现的次数,但是这样的时间复杂度是O(n),还是 ...
- 使用Servlet前Tomcat介绍
虚拟目录的映射方式:让tomcat服务器自动映射tomcat服务器会自动管理webapps目录下的所有web应用,并把它映射成虚似目录.换句话说,tomcat服务器webapps目录中的web应用,外 ...
- spring源码分析(二)- 容器基础
1.基本用法 用过Spring的都知道,bean是Spring中最基础也是最核心的.首先看一个简单的例子. 一个类和一个配置文件 package bean; public class MyBean { ...
- Scrum Meeting 16
第16次例会报告 日期:2021年06月11日 会议主要内容概述: 最后一次例会 一.进度情况 我们采用日报的形式记录每个人的具体进度,链接Home · Wiki,如下记录仅为保证公开性: 组员 负责 ...
- BUAA-OO-UML
BUAA-OO-UML 作业架构设计分析 第一次作业 类图如下: 这个架构十分简明,就是在底层数据和调用者之间建立起一层隔离层.但其实可以将转换过程延迟到调用阶段. 第二次作业 类图如下: 架构基本同 ...
- Noip模拟47 2021.8.25
期望得分:55+24+53 实际得分:0+0+3 乐死 累加变量清零了吗? 打出更高的部分分暴力删了吗? 样例解释换行你看见了吗? T1 Prime 打出55分做法没删原来的暴力,结果就轻松挂55分 ...
- Photoshop cc 绿色版 最新版 下载
Photoshop cc 绿色版 下载 Photoshop cc 绿色版 最新版下载百度网盘下载 Photoshop 下载提取码: dh6z 作为一个程序员, 不懂点基本的作图都不配"新时代 ...
- Python ValueError: Attempted relative import in non-package Relative import相对引用 错误
包含相对路径import的python脚本不能直接运行,只能作为module被引用. 例如 from . import mod1 有这样代码的文件只能最为moulule为不能直接运行.相对路径就是相对 ...
- 网关服务spring cloud zuul
Zuul例子配置文件 spring.application.name=switch-gateway server.port=5555 请求路由 传统路由方式 zuul.routes.api-a-url ...
- zabbix 报警发送企业威信
1.组册企业微信,创建应用 2.下载脚本文件: https://raw.githubusercontent.com/OneOaaS/weixin-alert/master/weixin_linux_a ...