可用倍增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把最小 ...
随机推荐
- ES6 之 Math对象的扩展
1.ES5 http://www.w3school.com.cn/jsref/jsref_obj_math.asp 2.ES6 Math.trunc() - 取整,去掉一个数的小数部分 console ...
- import datetime
import datetimenow = datetime.datetime.now()print('当前时间:',now) 当前时间: 2019-11-21 11:11:58.093122
- template.js的介绍
artTemplate 介绍 artTemplate 是新一代 javascript 模板引擎,它采用预编译方式让性能有了质的飞跃,并且充分利用 javascript 引擎特性,使得其性能无论在前端还 ...
- 关于Java编码规范
一.尽量使用卫语句 卫语句概念 条件表达式通常有两种表现形式,第一种形式是:所有分支都属于正常行为:第二种形式则是:条件表达式提供的答案中只有一种是正常行为,其他都是不常见的情况.这两类条件表达式有不 ...
- dp--背包--开心的金明
题目描述 金明今天很开心,家里购置的新房就要领钥匙了,新房里有一间他自己专用的很宽敞的房间.更让他高兴的是,妈妈昨天对他说:“你的房间需要购买哪些物品,怎么布置,你说了算,只要不超过N元钱就行”.今天 ...
- h5集成环信在线客服自定义窗口
自定义客服窗口从底部弹出 <!DOCTYPE html> <html> <head> <meta charset="UTF-8"> ...
- redis的过期策略
1.了解redis 什么是Redis,为啥用缓存? Redis是用内存当缓存的.Redis主要是基于内存来进行高性能.高并发的读写操作的. 内存是有限的,比如Redis就只能用10个G,你一直往里面写 ...
- aop 实现原理
aop 底层采用代理机制实现 接口 + 实现类 :spring 采用 jdk 的 动态代理 只有实现类:spring 采用 cglib 字节码增强 aop专业术语 1.target(目标) 需要被代理 ...
- linux 中的.so和.a文件
Linux下的.so是基于Linux下的动态链接,其功能和作用类似与windows下.dll文件. 下面是关于.so的介绍: 一.引言 通常情况下,对函数库的链接是放在编译时期(compile tim ...
- 基于JSP+Servlet开发在线租车系统 java 源码
运行环境: 最好是java jdk 1.8,我们在这个平台上运行的.其他版本理论上也可以.IDE环境: Eclipse,Myeclipse,IDEA都可以tomcat环境: Tomcat 7.x,8. ...