一道很好的圆方树入门题

感谢PinkRabbit巨佬的博客,讲的太好啦

首先是构建圆方树的代码,也比较好想好记

void tarjan(int u) {
dfn[u] = low[u] = ++dfn_clk; //初始化dfn和low数组
stk[++tp] = u; //把u加入栈中
for(int i = head[u]; i; i = e[i].next) {
int v = e[i].to;
if(!dfn[v]) { //v还未访问
tarjan(v); //先访问
low[u] = min(low[u], low[v]); //然后更新u的信息
if(low[v] == dfn[u]) { //找到一个以u为顶点的点双
++tot; //新建一个方点
for(int x = 0; x != v; --tp) { //把栈中在v及其之前的点都向方点连边并弹出
x = stk[tp];
G[tot].push_back(x);
G[x].push_back(tot);
}
G[tot].push_back(u); //注意不能把u弹出
G[u].push_back(tot); //因为u可能在多个点双中
}
}
else low[u] = min(low[u], dfn[v]);
}
}

注释写的还算详细\(QWQ\)

考虑这一题怎么做

题目大意

给你一张无向图,让你求这样的有序三元组\(<s,c,f>\)的个数,使得存在一条简单路径依次经过\(s,c,f\)

Solution

首先我们把圆方树建出来

考虑如下性质,对于在同一个点双中的两点\(s,t\),还有一个给定的也在这个点双中的点\(c\),一定存在一条简单路径依次经过\(s,c,t\),貌似挺显然的

在这题中,假设钦定了路径的两个端点\(s,t\),由上面的性质,那么能作为中间点的点集就是在圆方树上\(s\)到\(t\)的路径所经过的方点所代表的点双的并集(不包括\(s,t\))。这句话是本题的突破点,虽然有点拗口

然后是一个很套路的处理,把方点的点权设为点双的大小,圆点的点权设为\(-1\),这样的话上面要求的值就转化为\(s\)到\(t\)路径上的点权之和了

直接枚举\(s\)和\(t\)显然不行,考虑枚举每个点对答案的贡献,即

\[w[u]=val[u]*经过u的路径条数
\]

然后用树形\(dp\)就可以\(O(n)\)的统计答案了

另外,注意图不一定联通,所以需要单独统计每个联通块中的答案

#include <bits/stdc++.h>

using namespace std;

#define N 500000
#define M 200000 int n, m, tot;
int head[N+5], eid;
int dfn[N+5], low[N+5], dfn_clk;
int stk[N+5], tp, val[N+M+5], vis[N+M+5], cnt[N+M+5], sz[N+M+5], S;
long long ans;
vector<int> G[N+5]; struct Edge {
int next, to;
}e[2*M+5]; void addEdge(int u, int v) {
e[++eid].next = head[u];
e[eid].to = v;
head[u] = eid;
} void tarjan(int u) {
dfn[u] = low[u] = ++dfn_clk; //初始化dfn和low数组
stk[++tp] = u; //把u加入栈中
for(int i = head[u]; i; i = e[i].next) {
int v = e[i].to;
if(!dfn[v]) { //v还未访问
tarjan(v); //先访问
low[u] = min(low[u], low[v]); //然后更新u的信息
if(low[v] == dfn[u]) { //找到一个以u为顶点的点双
++tot; //新建一个方点
for(int x = 0; x != v; --tp) { //把栈中在u之前的点都向方点连边并弹出
x = stk[tp];
G[tot].push_back(x);
G[x].push_back(tot);
}
G[tot].push_back(u); //注意不能把u弹出
G[u].push_back(tot); //因为u可能在多个点双中
}
}
else low[u] = min(low[u], dfn[v]);
}
} void getcnt(int u, int fa) {
if(u <= n) cnt[u] = 1;
for(int i = 0; i < G[u].size(); ++i) {
int v = G[u][i];
if(v == fa) continue;
getcnt(v, u);
cnt[u] += cnt[v];
}
} void dp(int u, int fa) {
vis[u] = 1;
if(u <= n) sz[u] = 1;
long long sum = 1LL*cnt[u]*(S-cnt[u]);
for(int i = 0; i < G[u].size(); ++i) {
int v = G[u][i];
if(v == fa) continue;
dp(v, u);
sum += 1LL*sz[u]*sz[v];
sz[u] += sz[v];
}
ans += 1LL*val[u]*sum;
} int main() {
scanf("%d%d", &n, &m);
tot = n;
for(int i = 1, x, y; i <= m; ++i) {
scanf("%d%d", &x, &y);
addEdge(x, y), addEdge(y, x);
}
for(int i = 1; i <= n; ++i)
if(!dfn[i]) tarjan(i);
for(int i = 1; i <= n; ++i) val[i] = -1;
for(int i = n+1; i <= tot; ++i) val[i] = G[i].size();
for(int i = 1; i <= tot; ++i)
if(!vis[i]) {
getcnt(i, 0);
S = cnt[i];
dp(i, 0);
}
printf("%lld\n", ans*2);
return 0;
}

洛谷P4630 铁人两项--圆方树的更多相关文章

  1. [APIO2018] Duathlon 铁人两项 圆方树,DP

    [APIO2018] Duathlon 铁人两项 LG传送门 圆方树+简单DP. 不会圆方树的话可以看看我的另一篇文章. 考虑暴力怎么写,枚举两个点,答案加上两个点之间的点的个数. 看到题面中的一句话 ...

  2. [APIO2018]铁人两项 --- 圆方树

     [APIO2018] 铁人两项 题目大意: 给定一张图,问有多少三元组(a,b,c)(a,b,c 互不相等)满足存在一条点不重复的以a为起点,经过b,终点为c的路径 如果你不会圆方树 ------- ...

  3. [APIO2018]铁人两项 [圆方树模板]

    把这个图缩成圆方树,把方点的权值设成-1,圆点的权值设成点双的size,算 经过这个点的路径的数量*这个点的点权 的和即是答案. #include <iostream> #include ...

  4. [APIO2018]铁人两项——圆方树+树形DP

    题目链接: [APIO2018]铁人两项 对于点双连通分量有一个性质:在同一个点双里的三个点$a,b,c$,一定存在一条从$a$到$c$的路径经过$b$且经过的点只被经过一次. 那么我们建出原图的圆方 ...

  5. 【Luogu4630】【APIO2018】 Duathlon 铁人两项 (圆方树)

    Description ​ 给你一张\(~n~\)个点\(~m~\)条边的无向图,求有多少个三元组\(~(x, ~y, ~z)~\)满足存在一条从\(~x~\)到\(~z~\)并且经过\(~y~\)的 ...

  6. LOJ 2587 「APIO2018」铁人两项——圆方树

    题目:https://loj.ac/problem/2587 先写了 47 分暴力. 对于 n<=50 的部分, n3 枚举三个点,把图的圆方树建出来,合法条件是 c 是 s -> f 路 ...

  7. loj2587 「APIO2018」铁人两项[圆方树+树形DP]

    主要卡在一个结论上..关于点双有一个常用结论,也经常作为在圆方树/简单路径上的良好性质,对于任意点双内互不相同的三点$s,c,t$,都存在简单路径$s\to c\to t$,证明不会.可以参见clz博 ...

  8. [BZOJ5463][APIO2018]铁人两项(圆方树DP)

    题意:给出一张图,求满足存在一条从u到v的长度大于3的简单路径的有序点对(u,v)个数. 做了上一题[HDU5739]Fantasia(点双连通分量+DP),这个题就是一个NOIP题了. 一开始考虑了 ...

  9. 洛谷4630APIO2018铁人两项(圆方树+dp)

    QWQ神仙题啊(据说是今年第一次出现圆方树的地方) 首先根据题目,我们就是求对于每一个路径\((s,t)\)他的贡献就是两个点之间的点数,但是图上问题我并没有办法很好的解决... 这时候考虑圆方树,我 ...

随机推荐

  1. Vmware安装CentOs7+gitlab(一)

    本篇文章主要介绍了VMware安装Centos7超详细过程(图文),具有一定的参考价值,感兴趣的小伙伴们可以参考一下 1.软硬件准备 软件:推荐使用VMwear,我用的是VMwear 12 镜像:Ce ...

  2. 005. [转] SSH端口转发

    玩转SSH端口转发 SSH有三种端口转发模式,本地端口转发(Local Port Forwarding),远程端口转发(Remote Port Forwarding)以及动态端口转发(Dynamic ...

  3. vue实例化

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  4. JavaScript局部变量变量和函数命名提升

    之前接触了一些javascript局部变量命名提升的问题但是一直没有总结今天特地好好总结一下 变量提升 一个变量的作用域是程序源代码中定义的这个变量的区域.全局变量拥有全局作用域,在javascrip ...

  5. Swiper4.x使用方法

    1.首先加载插件,需要用到的文件有swiper.min.js和swiper.min.css文件.可下载Swiper文件或使用CDN. <!DOCTYPE html> <html> ...

  6. 转:Git Submodule管理项目子模块

    使用场景 当项目越来越庞大之后,不可避免的要拆分成多个子模块,我们希望各个子模块有独立的版本管理,并且由专门的人去维护,这时候我们就要用到git的submodule功能. 常用命令 git clone ...

  7. 本地服务器硬件信息获取指令wmic

    获取BIOS序列号 wmic bios list full | find "SerialNumber" SerialNumber=P50168VB 获取CPUID(WIN32_PR ...

  8. VS Code常用快捷键大全

    常用 General 按 Press 功能 Function Ctrl + Shift + P,F1 显示命令面板 Show Command Palette Ctrl + P 快速打开 Quick O ...

  9. scrapy安装失败:error:Microsoft Visual C++ 14.0 is reuired.及同类型安装问题解决办法

    今天在安装scrapy的时候(pip install Scrapy),出现了如下错误: building 'twisted.test.raiser' extensionerror: Microsoft ...

  10. anaconda 环境新建/删除/拷贝 jupyter notebook上使用python虚拟环境 TensorFlow

    naconda修改国内镜像源 国外网络有时太慢,可以通过配置把下载源改为国内的通过 conda config 命令生成配置文件,这里使用清华的镜像: https://mirrors.tuna.tsin ...