[BZOJ2125]最短路(圆方树DP)
题意:仙人掌图最短路。
算法:圆方树DP,$O(n\log n+Q\log n)$
首先建出仙人掌圆方树(与点双圆方树的区别在于直接连割边,也就是存在圆圆边),然后考虑点u-v的最短路径,显然就是:在圆方树上u-v的路径上的所有边权之和,加上每个环(方点)中连出去的两个点的最短距离。
现在问题就是:如何求出环上两个点的最短路径。考虑这样设定边权,首先显然圆圆边的边权就是原图的边权,然后设一个环在搜索树中深度最小的点为这个环的根,则方圆边的边权是环的根到这个点的最短距离,这个可以在Tarjan的时候直接求出。
但是圆方树问题通常需要在LCA处分圆方点讨论。首先如果LCA是圆点,那么直接做即可。如果是方点,就需要决定要不要走环的另一侧,这个同样直接讨论即可。
具体见代码,感觉思路还是比较清晰的。
#include<cstdio>
#include<algorithm>
#define rep(i,l,r) for (int i=l; i<=r; i++)
using namespace std; const int N=;
int n,m,Q,u,v,w,tot,tim,top,dep[N],len[N],type[N],stk[N];
int dfn[N],low[N],dis[N],lst[N],fa[N][],sm[N][]; struct E{
int cnt,h[N],to[N<<],nxt[N<<],val[N<<];
void add(int u,int v,int w){ to[++cnt]=v; val[cnt]=w; nxt[cnt]=h[u]; h[u]=cnt; }
}G,G1; void work(int x,int k){
tot++; int t; len[tot]=dis[stk[top]]-dis[x]+lst[stk[top]];
do{
t=stk[top--];
int A=dis[t]-dis[x],B=len[tot]-A;
G1.add(tot,t,min(A,B)); type[t]=(A<=B);
}while (t!=k);
G1.add(x,tot,);
} void Tarjan(int x,int pre){
//printf("%d\n",x);
dfn[x]=low[x]=++tim; stk[++top]=x;
for (int i=G.h[x],k; i; i=G.nxt[i]){
if ((k=G.to[i])==pre) continue;
if (!dfn[k]){
dis[k]=dis[x]+G.val[i]; Tarjan(k,x);
//printf("%d %d %d %d\n",x,k,dfn[x],low[k]);
if (low[k]>dfn[x]) top--,G1.add(x,k,G.val[i]);
else if (low[k]==dfn[x]) work(x,k);
low[x]=min(low[x],low[k]);
}else low[x]=min(low[x],dfn[k]),lst[x]=G.val[i];
}
} void dfs(int x,int pre){
for (int i=G1.h[x],k; i; i=G1.nxt[i])
fa[k=G1.to[i]][]=x,dep[k]=dep[x]+,sm[k][]=G1.val[i],dfs(k,x);
} int lca(int u,int v){
if (dep[u]<dep[v]) swap(u,v);
int t=dep[u]-dep[v],res=;
for (int i=; ~i; i--) if (t&(<<i)) res+=sm[u][i],u=fa[u][i];
if (u==v) return res;
for (int i=; ~i; i--) if (fa[u][i]!=fa[v][i])
res+=sm[u][i]+sm[v][i],u=fa[u][i],v=fa[v][i];
if (fa[u][]<=n) return sm[u][]+sm[v][]+res;
int A=sm[u][],B=sm[v][],mn;
if (type[u]==type[v]) mn=min(abs(A-B),len[fa[u][]]-abs(A-B));
else mn=min(A+B,len[fa[u][]]-A-B);
return res+mn;
} int main(){
freopen("bzoj2125.in","r",stdin);
freopen("bzoj2125.out","w",stdout);
scanf("%d%d%d",&n,&m,&Q); tot=n;
rep(i,,m) scanf("%d%d%d",&u,&v,&w),G.add(u,v,w),G.add(v,u,w);
Tarjan(,); dfs(,);
//rep(i,1,tot) printf("%d ",low[i]); puts("");
rep(j,,) rep(i,,tot)
fa[i][j]=fa[fa[i][j-]][j-],sm[i][j]=sm[i][j-]+sm[fa[i][j-]][j-];
rep(i,,Q) scanf("%d%d",&u,&v),printf("%d\n",lca(u,v));
return ;
}
[BZOJ2125]最短路(圆方树DP)的更多相关文章
- BZOJ2125 最短路 圆方树、倍增
传送门 对仙人掌建立圆方树,然后对边定权 对于圆点和圆点之间的边,是原来仙人掌上的桥,边权保持不变 对于圆点和方点之间的边,将圆方树看做以一个圆点为根的有根树之后,一个方点的父亲一定是一个圆点.对于这 ...
- [BZOJ2125]最短路[圆方树]
题意 给定仙人掌,多次询问两点之间的最短路径. \(n\le 10000, Q\le 10000\) 分析 建出圆方树,分路径 lca 是圆点还是方点讨论. 预处理出根圆点到每个圆点的最短距离 \( ...
- BZOJ1023:[SHOI2008]cactus仙人掌图(圆方树,DP,单调队列)
Description 如果某个无向连通图的任意一条边至多只出现在一条简单回路(simple cycle)里,我们就称这张图为仙人掌图(cactus). 所谓简单回路就是指在图上不重复经过任何一个顶点 ...
- 【BZOJ】2125: 最短路 圆方树(静态仙人掌)
[题意]给定带边权仙人掌图,Q次询问两点间最短距离.n,m,Q<=10000 [算法]圆方树处理仙人掌问题 [题解]树上的两点间最短路问题,常用倍增求LCA解决,考虑扩展到仙人掌图. 先对仙人掌 ...
- [BZOJ4316]小C的独立集(圆方树DP)
题意:求仙人掌图直径. 算法:建出仙人掌圆方树,对于圆点直接做普通的树上DP(忽略方点儿子),方点做环上DP并将值直接赋给父亲. 建图时有一个很好的性质,就是一个方点在邻接表里的点的顺序正好就是从环的 ...
- [BZOJ5463][APIO2018]铁人两项(圆方树DP)
题意:给出一张图,求满足存在一条从u到v的长度大于3的简单路径的有序点对(u,v)个数. 做了上一题[HDU5739]Fantasia(点双连通分量+DP),这个题就是一个NOIP题了. 一开始考虑了 ...
- [HDU5739]Fantasia(圆方树DP)
题意:给一张无向点带有权无向图.定义连通图的权值为图中各点权的乘积,图的权值为其包含的各连通图的权和.设z_i为删除i点后图的权值,求$S = (\sum\limits_{i=1}^{n}i\cdot ...
- 洛谷4630APIO2018铁人两项(圆方树+dp)
QWQ神仙题啊(据说是今年第一次出现圆方树的地方) 首先根据题目,我们就是求对于每一个路径\((s,t)\)他的贡献就是两个点之间的点数,但是图上问题我并没有办法很好的解决... 这时候考虑圆方树,我 ...
- 仙人掌&圆方树学习笔记
仙人掌&圆方树学习笔记 1.仙人掌 圆方树用来干啥? --处理仙人掌的问题. 仙人掌是啥? (图片来自于\(BZOJ1023\)) --也就是任意一条边只会出现在一个环里面. 当然,如果你的图 ...
随机推荐
- C++——设计模式说明
一.设计模式6大原则 名称 解释0.单一职责原则(SRP) 就一个类而言,应该仅有一个引起它变化的原因.一."开放-封闭"原则(OCP) 在软件设计模式中,这种不能修改,但可以扩展 ...
- 【luogu 1439 最长公共子序列】
题目描述 给出1-n的两个排列P1和P2,求它们的最长公共子序列. 输入输出格式 输入格式: 第一行是一个数n, 接下来两行,每行为n个数,为自然数1-n的一个排列. 输出格式: 一个数,即最长公共子 ...
- wait for it
- CSS3边框会动的信封
<!DOCTYPE html><html> <head> <title>酷炫的CSS3</title> <meta charset=& ...
- Java中文乱码问题(转)
解决JSP中文乱码问题 大家在JSP的开发过程中,经常出现中文乱码的问题,可能一至困扰着大家,现把JSP开发中遇到的中文乱码的问题及解决办法写出来供大家参考.首先了解一下Java中文问题的由来: Ja ...
- 用eval转化对象
var str = '{"name": "tom","age": 12,"sex": "man"}' ...
- uoj198【CTSC2016】时空旅行
传送门:http://uoj.ac/problem/198 [题解] 首先y.z是没有用的.. 然后式子就是w = (x0-xi)^2+ci的最小值,化出来可以变成一个直线的形式. 然后我们可以用线段 ...
- 【AtCoder】ARC082 F - Sandglass
[链接]F - Sandglass [题意]给定沙漏A和B,分别装着a和X-a的沙子,开始时A在上B在下,每秒漏1,漏完不再漏.给定n,有n个时刻ai沙漏倒转.给定m个询问,每次询问给定初值a和时刻t ...
- vue入门介绍
vue:解决前端大型应用的开发,将之前几十个.几百个.更多的HTML页面集成1个html页面(单页面应用)jquery:前端方法库bootstrap:UI组件库 angular/react这两个框架都 ...
- UART硬件流控制信号的使用(图)DTR 【转】
UART硬件流控制信号的使用(图) 转自:http://blog.163.com/zy_tommy/blog/static/86926777201321925451164/ 2013-03-19 14 ...