P4630-[APIO2018]Duathlon铁人两项【圆方树】

1|0正题

题目链接:https://www.luogu.com.cn/problem/P4630


1|1题目大意

n个点m条边的一张无向图,求有多少对三元组(s,c,f)满足sft且存在一条从sf的简单路径经过c


1|2解题思路

一个比较显然的结论是在一个点双中的三个点(a,b,c)那么必然存在一条ab的简单路径经过c。因为一定存在两条不交的a>cc>b的路径,那么如果一条a>cc>b的路径交了,那么另一条就一定不交。

然后从一个点双出来后就不能再回到这个点双了,所以我们可以考虑在圆方树上做这个问题。

设定义圆点的权值为1,方点的权值为连接的圆点数量,这样我们在圆方树上走的时候就可以固定经过进入和离开这个点双的点了。

然后问题就变为了求每条圆点之间路径的点权和的和。

用树形dp搞就好了,时间复杂度O(n)


1|3code

#include<cstdio> #include<cstring> #include<algorithm> #include<vector> #include<stack> using namespace std; const int N=2e5+10; int n,m,num,cnt,dfc,w[N]; int low[N],dfn[N],siz[N]; vector<int> G[N],T[N]; stack<int> s; long long ans; void tarjan(int x){ dfn[x]=low[x]=++dfc; w[x]=-1;s.push(x);num++; for(int y:T[x]) if(!dfn[y]){ tarjan(y); low[x]=min(low[x],low[y]); if(dfn[x]==low[y]){ ++cnt;int k; do{ k=s.top(); G[cnt].push_back(k); G[k].push_back(cnt); w[cnt]++;s.pop(); }while(k!=y); G[cnt].push_back(x); G[x].push_back(cnt); w[cnt]++; } } else low[x]=min(low[x],dfn[y]); return; } void solve(int x,int fa){ siz[x]=(x<=n); for(int y:G[x]){ if(y==fa)continue; solve(y,x); ans+=2ll*siz[y]*siz[x]*w[x]; siz[x]+=siz[y]; } ans+=2ll*siz[x]*(num-siz[x])*w[x]; return; } int main() { scanf("%d%d",&n,&m); for(int i=1;i<=m;i++){ int x,y;scanf("%d%d",&x,&y); T[x].push_back(y); T[y].push_back(x); } cnt=n; for(int i=1;i<=n;i++) if(!dfn[i]){ num=0; tarjan(i); solve(i,0); } printf("%lld\n",ans); return 0; }

__EOF__

本文作者QuantAsk
本文链接https://www.cnblogs.com/QuantAsk/p/14273423.html
关于博主:退役OIer,GD划水选手
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角推荐一下。您的鼓励是博主的最大动力!
posted @   QuantAsk  阅读(70)  评论(0编辑  收藏  举报
编辑推荐:
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
· Linux系列:如何用heaptrack跟踪.NET程序的非托管内存泄露
阅读排行:
· TypeScript + Deepseek 打造卜卦网站:技术与玄学的结合
· 阿里巴巴 QwQ-32B真的超越了 DeepSeek R-1吗?
· 【译】Visual Studio 中新的强大生产力特性
· 张高兴的大模型开发实战:(一)使用 Selenium 进行网页爬虫
· 【设计模式】告别冗长if-else语句:使用策略模式优化代码结构
点击右上角即可分享
微信分享提示