可用倍增LCA解题
http://codevs.cn/problem/2370/
#include<bits/stdc++.h>
using namespace std;
const int M=5e4+;
const int N=;
struct node{
int v,w;
node(int vv=,int ww=):v(vv),w(ww){}
};
vector<node>e[M];
int n,grand[M][N],dis[M][N],deep[M],book[M],root,s;
void dfs(int u){
for(int i=;i<=s;i++){
grand[u][i]=grand[grand[u][i-]][i-];
dis[u][i]=dis[u][i-]+dis[grand[u][i-]][i-];
if(!grand[u][i])
break;
}
for(int i=;i<e[u].size();i++){
int v=e[u][i].v;
if(v!=grand[u][]){
grand[v][]=u;
deep[v]=deep[u]+;
dis[v][]=e[u][i].w;
dfs(v);
}
}
}
void init(){
s=floor(log(1.0*n)/log(2.0));
deep[]=-;
dfs(root);
}
int LCA(int a,int b){
if(deep[a]>deep[b])
swap(a,b);
int ans=;
for(int i=s;i>=;i--){
if(deep[a]<deep[b]&&deep[a]<=deep[grand[b][i]])
ans+=dis[b][i],b=grand[b][i];
}
for(int i=s;i>=;i--){
if(grand[a][i]!=grand[b][i])
ans+=dis[a][i],ans+=dis[b][i],a=grand[a][i],b=grand[b][i];
}
if(a!=b)
ans+=dis[a][]+dis[b][];
return ans;
}
int main(){ scanf("%d",&n);
for(int i=;i<n;i++){
int u,v,w;
scanf("%d%d%d",&u,&v,&w);
u++;
v++;///fu
book[u]=;
e[u].push_back(node(v,w));
e[v].push_back(node(u,w));
}
for(int i=;i<=n;i++){
if(!book[i]){
root=i;
break;
}
}
init();
int m;
scanf("%d",&m);
while(m--){
int u,v;
scanf("%d%d",&u,&v);
u++,v++;
printf("%d\n",LCA(u,v));
}
return ;
}
http://codevs.cn/problem/1036/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;
const int M=3e4+;
const int N=;
struct node{
int v,w;
node(int vv=,int ww=):v(vv),w(ww){}
};
vector<node>e[M];
int s,grand[M][N],dis[M][N],book[M],deep[M],root,n;
void dfs(int u){
for(int i=;i<=s;i++){
grand[u][i]=grand[grand[u][i-]][i-];
dis[u][i]=dis[u][i-]+dis[grand[u][i-]][i-];
if(!grand[u][i])
break;
}
for(int i=;i<e[u].size();i++){
int v=e[u][i].v;
if(v!=grand[u][]){
grand[v][]=u;
deep[v]=deep[u]+;
dis[v][]=e[u][i].w;
dfs(v);
} // cout<<"!!"<<endl;
}
}
void init(){
s=floor(log(1.0*n)/log(2.0));
deep[]=-;
dfs(root);
}
int LCA(int a,int b){
if(deep[a]>deep[b])
swap(a,b);
int ans=;
for(int i=s;i>=;i--){
if(deep[b]>deep[a]&&deep[a]<=deep[grand[b][i]])
ans+=dis[b][i],b=grand[b][i];
}
for(int i=s;i>=;i--)
if(grand[a][i]!=grand[b][i])
ans+=dis[a][i]+dis[b][i],b=grand[b][i],a=grand[a][i];
if(a!=b)
ans+=dis[a][]+dis[b][];
return ans;
}
int main(){ scanf("%d",&n);
for(int i=;i<n;i++){
int u,v;
scanf("%d%d",&u,&v);
book[v]=;
e[u].push_back(node(v,));
e[v].push_back(node(u,));
}
for(int i=;i<=n;i++)
if(!book[i]){
root=i;
break;
} init();
int m,u,v;
scanf("%d",&m);
scanf("%d",&u);
int ans=;
for(int i=;i<m;i++){
scanf("%d",&v);
ans+=LCA(u,v);
u=v;
}
printf("%d\n",ans);
return ;
}
http://oj.fjaxyz.com:3389/problem.php?id=223
当LCA(a,b）==a时才能参赛
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;
typedef long long ll;
const int M=1e4+;
struct node{
int v;
ll w;
node(int vv=,ll ww=):v(vv),w(ww){}
};
vector<node>e[M];
int grand[M][],s,n,root,deep[M];
ll dis[M][],num,numt;
void dfs(int u){
for(int i=;i<=s;i++){
grand[u][i]=grand[grand[u][i-]][i-];
dis[u][i]=dis[u][i-]+dis[grand[u][i-]][i-];
if(!grand[u][i])
break;
}
for(int i=;i<e[u].size();i++){
int v=e[u][i].v;
if(v!=grand[u][]){
grand[v][]=u;
deep[v]=deep[u]+;
dis[v][]=e[u][i].w;
dfs(v);
}
}
}
void init(){
s=floor(log(1.0*n)/log(2.0));
deep[]=-;
dfs(root);
}
void LCA(int u,int v){
int a=u,b=v;
if(deep[a]>deep[b])
swap(a,b) ;
ll ans=;
for(int i=s;i>=;i--){
if(deep[b]>deep[a]&&deep[a]<=deep[grand[b][i]])
ans+=dis[b][i],b=grand[b][i];
}
for(int i=s;i>=;i--){
if(grand[b][i]!=grand[a][i])
ans+=dis[b][i]+dis[a][i],a=grand[a][i],b=grand[b][i];
}
if(a!=b)
ans+=dis[a][]+dis[b][],a=grand[a][],b=grand[b][];
if(a!=u)
return ;
// cout<<"!!"<<endl; numt+=ans;
num++;
}
int main(){
int m;
scanf("%d%d",&n,&m);
for(int i=;i<n;i++){
int u,v,w;
scanf("%d%d%d",&u,&v,&w);
e[u].push_back(node(v,w));
e[v].push_back(node(u,w));
}
root=;
init();
while(m--){
int u,v;
scanf("%d%d",&u,&v); LCA(u,v);
}
printf("%d\n",num);
printf("%d\n",numt); return ;
}
zoj 3198
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=3320
吐槽一下输出格式，言不表意！！
题意：求用最小权值的路径连接起3个节点，问这些路径的和；
分析，答案就是俩俩之间路径和/2；
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;
const int M=5e4+;
struct node{
int v,w;
node(int vv,int ww):v(vv),w(ww){}
};
vector<node>e[M];
int grand[M][],dis[M][],deep[M],root,s,n,book[M];
void dfs(int u){
for(int i=;i<=s;i++){
grand[u][i]=grand[grand[u][i-]][i-];
dis[u][i]=dis[grand[u][i-]][i-]+dis[u][i-];
if(!grand[u][i])
break;
}
for(int i=;i<e[u].size();i++){
int v=e[u][i].v;
if(v!=grand[u][]){
grand[v][]=u;
deep[v]=deep[u]+;
dis[v][]=e[u][i].w;
dfs(v);
}
}
}
void init(){
s=floor(log(1.0*n)/log(2.0));
deep[]=-;
dfs(root);
}
int LCA(int a,int b){
if(deep[a]>deep[b])
swap(a,b);
int ans=;
for(int i=s;i>=;i--){
if(deep[a]<deep[b]&&deep[a]<=deep[grand[b][i]])
ans+=dis[b][i],b=grand[b][i];
}
for(int i=s;i>=;i--){
if(grand[a][i]!=grand[b][i])
ans+=dis[a][i]+dis[b][i],a=grand[a][i],b=grand[b][i];
}
if(a!=b)
ans+=dis[a][]+dis[b][];
return ans;
}
int main(){
int flag=;
while(~scanf("%d",&n)&&n){
if(flag)
flag=;
else
printf("\n");
for(int i=;i<=n;i++){
e[i].clear();
deep[i]=;
book[i]=;
for(int j=;j<=;j++)
grand[i][j]=dis[i][j]=; }
for(int i=;i<n;i++){
int u,v,w;
scanf("%d%d%d",&u,&v,&w);
u++;
v++;
book[v]=;
e[u].push_back(node(v,w));
e[v].push_back(node(u,w));
}
for(int i=;i<=n;i++){
if(!book[i]){
root=i;
break;
}
}
init();
int m;
scanf("%d",&m);
while(m--){
int u,v,k;
scanf("%d%d%d",&u,&v,&k);
u++,v++,k++;
int ans=;
ans+=LCA(u,v);
// cout<<LCA(u,v)<<"~~~~~~~~1"<<endl;
ans+=LCA(u,k);
// cout<<LCA(u,k)<<"~~~~~~~~2"<<endl;
ans+=LCA(v,k);
// cout<<LCA(k,v)<<"~~~~~~~~3"<<endl;
printf("%d\n",ans/);
} }
return ;
}
可用倍增LCA解题的更多相关文章
- [板子]倍增LCA
倍增LCA板子,没有压行,可读性应该还可以.转载请随意. #include <cstdio> #include <cstring> #include <algorithm ...
- 洛谷P3128 [USACO15DEC]最大流Max Flow [倍增LCA]
题目描述 Farmer John has installed a new system of pipes to transport milk between the stalls in his b ...
- Gym100685G Gadget Hackwrench（倍增LCA）
题目大概说一棵边有方向的树,q个询问,每次询问结点u是否能走到v. 倍增LCA搞即可: 除了par[k][u]表示u结点往上走2k步到达的结点, 再加上upp[k][u]表示u结点往上走2k步经过边的 ...
- Codeforces 418d Big Problems for Organizers [树形dp][倍增lca]
题意: 给你一棵有n个节点的树,树的边权都是1. 有m次询问,每次询问输出树上所有节点离其较近结点距离的最大值. 思路: 1.首先是按照常规树形dp的思路维护一个子树节点中距离该点的最大值son_di ...
- hdu 4674 Trip Advisor（缩点+倍增lca）
花了一天半的时间,才把这道题ac= = 确实是道好题,好久没敲这么长的code了,尤其是最后的判定,各种销魂啊~ 题目中给出的条件最值得关注的就是:每个点最多只能在一个环内->原图是由一个个边连 ...
- Tsinsen A1505. 树（张闻涛） 倍增LCA,可持久化线段树,DFS序
题目:http://www.tsinsen.com/A1505 A1505. 树(张闻涛) 时间限制:1.0s 内存限制:512.0MB 总提交次数:196 AC次数:65 平均分: ...
- codevs 1036 商务旅行 （倍增LCA）
/* 在我还不知道LCA之前 暴力跑的SPFA 70分 三个点TLE */ #include<iostream> #include<cstdio> #include<cs ...
- hdu 2586 How far away ？倍增LCA
hdu 2586 How far away ?倍增LCA 题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=2586 思路: 针对询问次数多的时候,采取倍增 ...
- 洛谷P4180 [Beijing2010组队]次小生成树Tree（最小生成树,LCT,主席树,倍增LCA,倍增,树链剖分）
洛谷题目传送门 %%%TPLY巨佬和ysner巨佬%%% 他们的题解 思路分析 具体思路都在各位巨佬的题解中.这题做法挺多的,我就不对每个都详细讲了,泛泛而谈吧. 大多数算法都要用kruskal把最小 ...
随机推荐
- 吴裕雄--天生自然 PHP开发学习：While 循环
<html> <body> <?php $i=1; while($i<=5) { echo "The number is " . $i . &q ...
- CSS3 box-shadow 效果大全（内阴影，外阴影，三边阴影，双边阴影，单边阴影，细线描边…）
/* offset-x | offset-y | color */ box-shadow: 60px -16px teal; /* offset-x | offset-y | blur-radius ...
- 2020/1/28 PHP代码审计之代码执行漏洞
0x00代码执行原理 应用程序在调用一些能够将字符串转换为代码的函数(如PHP中的eval)时,没有考虑用户是否控制这个字符串,将造成代码执行漏洞. 该漏洞主要存在于eval().assert().p ...
- vue 动画原理 part1
Vue动画原理 增加和删除css增加样式实现一个过渡效果也就是动画效果 1.需要动画效果的标签外包裹一个transition标签 会被自动分析css样式,然后自动构建一个动画流程 transition ...
- python format输出
http://www.cnblogs.com/nulige/p/6115793.html 2.Format 方式 [[fill]align][sign][#][0][width][,][.precis ...
- 16. docker 网络 端口映射
一.本地操作 1.如何将 nginx 暴露给外界 创建 nginx 服务器 docker run --name web -d nginx 查看 nginx 的 ip地址 docker network ...
- Spring Cloud服务间调用鉴权
学习使用Spring Cloud 微服务间的调用都是RestFul风格,如何保证调用之间的安全性,这是一个很重要的问题. 通过查阅资料http://wiselyman.iteye.com/blog/2 ...
- MyBatis从入门到精通(第4章)：MyBatis动态SQL【foreach、bind、OGNL用法】
(第4章):MyBatis动态SQL[foreach.bind.OGNL用法] 4.4 foreach 用法 SQL 语句中有时会使用 IN 关键字,例如 id in (1,2,3).可以使用 ${i ...
- CertUtil: -hashfile 失败: 0xd00000bb (-805306181)
使用CertUtil验证Python安装文件的时候出现了这个错误. CertUtil: -hashfile 失败: 0xd00000bb (-805306181) 代码是这样 certutil -ha ...
- drf三大认证：认证组件-权限组件-权限六表-自定义认证组件的使用
三大认证工作原理简介 认证.权限.频率 源码分析: from rest_framework.views import APIView 源码分析入口: 内部的三大认证方法封装: 三大组件的原理分析: 权 ...