P5327 [ZJOI2019]语言

 

 P5327 [ZJOI2019]语言

题目描述

详见:P5327 [ZJOI2019]语言

简要题意:给定一棵树和一些链,问树上处于同一条链的不同点对数。

Solution

对于每一个点i,考虑以它为端点的可行路径有哪些。

我们可以发现,i可以到达的节点会组成一个斯坦纳树,这棵斯坦纳树包含了i即经过i链。

我们进一步可知,这棵斯坦纳树就是以i和经过i的所有链的端点u,v为关键点的最小斯坦纳树。

所以我们考虑对于每一个点,维护经过它的链组成的斯坦纳树是什么,顺便维护答案(斯坦纳树的边数)。

 

事实上,斯坦纳树的边数就是其中所有点的深度-dfs序相邻的点的LCA的深度,如下:

设斯坦纳树的点集为V,记i的dfs序为  dfn[i],深度为  dep[i]  。

若满足  \forall_{i<j}\;\;dfn[V[i]]<dfn[V[j]]  ,即dfs序按升序排序,则

ans_u=\sum_{i}{dep[V[i]]}-\sum_{i}{dep[lca(V[i],V[i+1])]}-dep[lca(V[k],V[1])]

Ans=\sum ans_i

 

可以用线段树合并(启发式合并也可)支持上述信息的维护。

线段树的下标为  dfn,并记录  [l,r]  中dfn最大的点ladfn 最小的点fi,以及 ans  。

合并时   ans_{l,r}=ans_{l,mid}+ans_{mid+1,r}-dep[lca(la_{l,mid},fi_{mid+1,r})]

统计答案时  Ans+=ans_{l,r}-dep[lca(fi_{l,r},la_{l,r})] 

 

因此,对于每一条路径  (u,v)  ,我们可以用树上差分把它拆成  (1,u),(1,v),(1,lca),(1,fa[lca])  这四条代价分别为+1,+1,-1,-1的链,分别放入权值线段树,并不断向上合并到根,沿路统计答案即可。

st表实现LCA,时间复杂度  O(nlgn)

如果写倍增LCA的会多一个log,卡卡常应该也能过。

 

PS:这道题写了1h,调了5h,原因是线段树动态开点要开log倍的内存。但我一直以为是dfs爆栈了(因为今早有一题先例),由于要用到dfs序,所以甚至在try coding 手工栈,However,WA声依旧,改得代码奇丑无比,自闭了一个下午,最后柳暗花明又一村,有一种想右转直行的冲动。

Code

最后贴一个丑陋的代码(早把罪恶的手工栈删了)。

#include<bits/stdc++.h>
#define RG register
typedef long long ll;
const int MAXN=4e5+50;
std::vector<int> e[MAXN],Del[MAXN];
int DFN=0,t=0,n,m;
int dfn[MAXN],st[MAXN][21],fa[MAXN];
int pre[MAXN],rt[MAXN],Log[MAXN];
ll dep[MAXN],ans=0;
inline int read()
{RG int x=0,f=1; char c=getchar();while (c<'0'||c>'9') { if(c=='-') f=-1; c=getchar(); }while (c>='0'&&c<='9') { x=(x<<3)+(x<<1)+(c^48); c=getchar(); }return x*f;
}
inline void get_dfn(int x,int father)
{dfn[x]=++DFN;pre[x]=++t; st[t][0]=x;fa[x]=father; dep[x]=dep[father]+1;for (RG int i=0;i<e[x].size();i++)if (e[x][i]!=fa[x]) get_dfn(e[x][i],x),st[++t][0]=x;
}
inline void get_st()
{Log[1]=0;for (RG int i=2;i<=t;i++) Log[i]=Log[i>>1]+1;for (RG int i=1;i<=Log[t];i++)for (RG int j=1;j<=t-(1<<i)+1;j++) st[j][i]=dep[st[j][i-1]]<dep[st[j+(1<<(i-1))][i-1]]?st[j][i-1]:st[j+(1<<(i-1))][i-1];
}
inline int get_lca(int x,int y)
{x=pre[x],y=pre[y];if (x>y) std::swap(x,y);RG int k=Log[y-x+1];return std::max(1,dep[st[x][k]]<dep[st[y-(1<<k)+1][k]]?st[x][k]:st[y-(1<<k)+1][k]);
}
struct Segment_Tree
{int segnum=0;struct treenode{int lson,rson,fi,la; ll sum,ans; } tree[MAXN<<4];#define lft tree[x].lson#define rht tree[x].rsoninline void up(int x){RG int lca=get_lca(tree[lft].la,tree[rht].fi);//cout<<lft<<" "<<rht<<" "<<lca<<endl;tree[x].ans=tree[lft].ans+tree[rht].ans-dep[lca];tree[x].fi=tree[lft].fi?tree[lft].fi:tree[rht].fi;tree[x].la=tree[rht].la?tree[rht].la:tree[lft].la;}inline void change(int &x,int p,int v,int L,int R){if (!x) x=++segnum;if (L==R){tree[x].sum+=v;tree[x].ans=tree[x].sum?dep[p]:0;tree[x].la=tree[x].fi=tree[x].sum?p:0;return;}RG int mid=(L+R)>>1;if (dfn[p]<=mid) change(lft,p,v,L,mid);else change(rht,p,v,mid+1,R);up(x);}inline ll query(int x){ RG int lca=get_lca(tree[x].la,tree[x].fi);return tree[x].ans-dep[lca];}inline void merge(int &x,int v,int L,int R){if (!x||!v) { x|=v; return; }if (L==R){tree[x].sum+=tree[v].sum;tree[x].ans|=tree[v].ans;tree[x].la|=tree[v].la;tree[x].fi|=tree[v].fi;return;}RG int mid=(L+R)>>1;merge(lft,tree[v].lson,L,mid);merge(rht,tree[v].rson,mid+1,R);up(x);}
} segment;
inline void solve(int x)
{for (RG int i=0;i<e[x].size();i++)if (e[x][i]!=fa[x]) solve(e[x][i]);for (RG int i=0;i<Del[x].size();i++) segment.change(rt[x],Del[x][i],-1,1,n);ans+=segment.query(rt[x]);//cout<<ans<<endl;if (fa[x]) segment.merge(rt[fa[x]],rt[x],1,n);
}
int main()
{//freopen("a.in","r",stdin);n=read(),m=read();for (RG int i=1;i<n;i++){int u=read(),v=read();e[u].push_back(v);e[v].push_back(u);}dep[0]=-1;get_dfn(1,0);get_st();for (RG int i=1;i<=m;i++) {RG int u=read(),v=read(),lca=get_lca(u,v);//cout<<lca<<endl;segment.change(rt[u],v,1,1,n); segment.change(rt[v],u,1,1,n); segment.change(rt[v],v,1,1,n); segment.change(rt[u],u,1,1,n); //segment.change(rt[lca],u,-1,1,n); segment.change(rt[lca],v,-1,1,n);//if (fa[lca]) segment.change(rt[fa[lca]],u,-1,1,n),segment.change(rt[fa[lca]],v,-1,1,n);//cout<<lca<<endl;Del[lca].push_back(u),Del[lca].push_back(v);if (fa[lca]) Del[fa[lca]].push_back(u),Del[fa[lca]].push_back(v);}solve(1);printf("%lld\n",ans>>1);return 0;
}
/*
12 4
1 2
2 3
2 4
4 5
1 6
6 7
6 8
6 9
9 10
9 11
1 124 12
8 10
7 11
3 10
*/

还附加了一个小样例,然而样例输出咕了。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315831.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

EasyOffice-.NetCore一行代码导入导出Excel,生成Word

Excel和Word操作在开发过程中经常需要使用&#xff0c;这类工作不涉及到核心业务&#xff0c;但又往往不可缺少。以往的开发方式在业务代码中直接引入NPOI、Aspose或者其他第三方库&#xff0c;工作繁琐&#xff0c;耗时多&#xff0c;扩展性差——比如基础库由NPOI修改为EPPlu…

CF1572B. Xor of 3

CF1572B. Xor of 3 题意&#xff1a; 给你个01序列&#xff0c;你有一种操作&#xff1a;每次选位置x&#xff0c;然后位置x&#xff0c;x1&#xff0c;x2的值变为三者的异或值。 现在要让所有的数都等于0&#xff0c;请输出存在的合法操作序列 题解&#xff1a; 首先如果有…

P4103 [HEOI2014]大工程

P4103 [HEOI2014]大工程 题目描述 详见&#xff1a;P4103 [HEOI2014]大工程 Solution 显然是虚树的板子题啊&#xff08;我也不造我为啥调了1h&#xff09; 直接建虚树&#xff0c;DP。 两两路径和很好求&#xff0c;只要考虑每条边的贡献&#xff0c;设为询问的关键点个…

【中部武汉】理想离家并不遥远

▼我叫Mike, 是一名“高端”程序猿。毕业于武汉某男生比例占据 70%的高校。我的家乡是湖北襄阳&#xff0c;一个历史悠久的文化名城。7年前&#xff0c;我还是一个刚毕业的毛头小子&#xff0c;在学校也算是半个风云人物&#xff0c;年轻气盛&#xff0c;心比天高&#xff0c;放…

P2495 [SDOI2011]消耗战

P2495 [SDOI2011]消耗战 题目描述 详见&#xff1a;P2495 [SDOI2011]消耗战 Solution 此题是虚树的经典版子题吧qwq。 是不是直接贴代码就行了啊&#xff08;反正不是luogu题解&#xff0c;没有人查&#xff09;。 大概就是先建出虚树&#xff08;一般给定一堆关键点的题…

使用Try.NET创建可交互.NET文档

原文地址&#xff1a;Create Interactive .NET Documentation with Try .NET[1] 原文作者&#xff1a;Maria 译文地址&#xff1a;https://www.cnblogs.com/lwqlun/p/10894497.html 译者&#xff1a;Lamond Lu背景当我们编写开发人员使用的文档时&#xff0c;我们需要捕捉他们的…

Link-Cut Tree

Link-Cut Tree 概述. LCT是一种支持动态维护树上路径信息的数据结构&#xff0c;其本质是实链剖分&#xff0c;通过其他数据结构维护实链的信息达到维护路径及一些子树信息的效果&#xff08;通常为splay&#xff09; 刚开始学的时候感觉很难&#xff0c;但做过几道题了解套…

2021-2022ACM赛季小总结

最近一个月比赛月&#xff0c;做题量和频率太高&#xff0c;写博客浪费时间&#xff0c;所以一直没写。。 当前总结&#xff1a; CCPC桂林铜&#xff0c;银牌题线段树差一行调出&#xff0c;可惜又无奈&#xff0c;缺少运气 ICPC沈阳铜&#xff0c;银牌题的tarjan没想出来&…

ASP.NET Core模块化前后端分离快速开发框架介绍

源码地址GitHub&#xff1a;https://github.com/iamoldli/NetModular演示地址地址&#xff1a;http://118.24.75.170:6220/账户&#xff1a;admin密码&#xff1a;admin前端框架演示地址(临时)地址&#xff1a;http://progqx5cu.bkt.clouddn.com/skins/index.html#/账户&#x…

[POJ2888] Magic Bracelet

[POJ2888] Magic Bracelet 题目描述 简要题意&#xff1a;给圆上个点染色&#xff0c;颜色有种&#xff0c;其中对颜色不能相邻&#xff0c;循环同构&#xff0c;多组数据&#xff0c;询问染色方案数。 Solution 大概就是一道挺显然的Burnside题&#xff08;一般染色&#x…

手写编译器

假期没事想搞搞编译器&#xff0c;但学校课程未教过编译原理&#xff0c;自己只好在b站看了stanford的CS143《编译原理》&#xff0c;并结合相关资料&#xff0c;打算寒假期间自己手写一个简易的编译器 ReactjsbootStrap开发自制编译语言Monkey的编译器&#xff1a;创建简易的页…

使用命令行导出 SQL Server 数据层应用程序

点击上方蓝字关注“汪宇杰博客”我们可以使用 SSMS 导出 SQL Server 数据库的数据层应用程序。我在本地机器上使用这种方法已经有好几年了。如果不知道什么是 DAC&#xff0c;您可以参考 Microsoft 文档&#xff1a;https://docs.microsoft.com/en-us/sql/relational-databases…

Codeforces 1065 E. Side Transmutations

Codeforces 1065 E. Side Transmutations 题目描述 题意&#xff1a;一个字符串长度为&#xff0c;字符集大小为 &#xff0c;再给出种交换方式&#xff0c;表示可以交换。如果两个不同的串能够通过任意的交换方式得到相同的串&#xff0c;则视为相同&#xff0c;求方案数。…

P1 用reactjs和bootstrap创建页面IDE

文章目录页面IDE展示&#xff1a;要求&#xff1a;安装好Node.jsreactjs组件化开发页面IDE展示&#xff1a; 要求&#xff1a; 希望能够达到eclipse辅助编程的功能 关键字提醒&#xff0c;高亮 1&#xff0c; 有类似C语言的语法结构 2&#xff0c; 支持变量绑定 3&#xff0c;…

所以Apache基金会不受美国法律约束?

开源中国昨天发布的一篇文章《开源界也要注意&#xff0c;Apache 基金会与 GitHub 都受美国法律约束》引发了开源届乃至整个 IT 行业的热烈讨论&#xff0c;其中有个别声音认为文章的说法有误&#xff0c;甚至制造了恐慌。文中我们引用了两处内容&#xff0c;分别是 Apache 软件…

Loj#6485. LJJ 学二项式定理

Loj#6485. LJJ 学二项式定理&#xff08;单位根反演&#xff09; 题目描述 题目描述 题意&#xff1a;求下面式子的答案QAQ。 [∑((ni)⋅si⋅aimod4)]mod998244353[\sum(\tbinom{n}{i}\cdot s^i \cdot a_{i\;\;mod\;\;4}) ]\;mod\;\;998244353 [∑((in​)⋅si⋅aimod4​)]mod…

团体程序设计天梯赛-练习集

团体程序设计天梯赛-练习集 L2-001 紧急救援 (25 分) 知识点&#xff1a;最短路dij #include<bits/stdc.h> #define debug(x,y) printf("%s %d\n",x,y); typedef long long ll; using namespace std; inline int read(){int s0,w1;char chgetchar();while(…

C#中的Explicit和Implicit了解一下吧

今天在研究公司项目框架的时候看到了下面的用法&#xff0c;public static implicit operator JsonData(int data);。貌似很久没用过这种隐式转换的写法了&#xff0c;因此重新温习一下C#中转换相关的知识。作者&#xff1a;依乐祝原文地址&#xff1a;https://www.cnblogs.com…

bzoj#2125. 最短路

bzoj#2125. 最短路 题目描述 Description 给一个N个点M条边的连通无向图&#xff0c;满足每条边最多属于一个环&#xff0c;有Q组询问&#xff0c;每次询问两点之间的最短路径。 Input 输入的第一行包含三个整数&#xff0c;分别表示N和M和Q 下接M行&#xff0c;每行三个整数v…

Dapper的正确使用姿势

本文demo适用于MySQL优点高性能、易排查、易运维、灵活可控缺点和EF相比&#xff0c;手写sql当修改表结构不易发现bug。习惯了EF后再来使用Dapper&#xff0c;会很难适应那种没有了强类型的安全感。不过可以用单元测和心细来避免。问题&#xff1a;IDbConnection需不需要手动Op…