可用倍增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把最小 ...
随机推荐
- 201612-2 工资计算 Java
思路: 税+税后所得A=税前工资S. 因为工资是整百的数,每次减100来判断.好理解但是超时. import java.util.Scanner; //只有90分,超时了 public class M ...
- html+css 通信课上 2019。3.22
数据通信 http协议:无状态.无连接.单向的应用层协议:采用请求/响应模型:通信请求只能由客户端发起,服务端对请求做出应答处理 服务器推送数据的解决方案:轮询( ajax) :让浏览器几秒就发送一次 ...
- 如何向女朋友介绍MySQL索引
目录 一.前言 二.正文 三.索引的类型 四.动态查找树 五.B-Tree 1.B-Tree特征 2.B-Tree的查找(select) 3.B-Tree的插入(insert) 4.B-Tree的删除 ...
- Nginx和php交互的两种方式
Unix socket 也叫IPC socket 也就是进程间通信套接字用于同一台主机上的不同进程间交换数据 TCP socket IP socket要利用主机的传输层(tcp),可以用于同一台主机 ...
- C++代码质量度量工具大阅兵
姊妹篇:Java代码质量度量工具大阅兵: https://www.cnblogs.com/jiangxinnju/p/10010177.html cppcheck cppcheck: https:// ...
- VC调用VB写的COM
VB. 步骤: 1.创建类库:类库的创建必须分为接口和实现类:给外面提供的是COM接口: 创建了接口和类之后还要创建"Guid",这个在"工具->创建GUID&qu ...
- error_reporting() 设置 PHP 的报错级别并返回当前级别
error_reporting() 设置 PHP 的报错级别并返回当前级别. 语法 error_reporting(report_level) 如果参数 level 未指定,当前报错级别将被返回.下面 ...
- liunx 常用操作(自用)
Centos7解压文件 tar -zxvf 文件名[test.tar.gz] Centos7安装vim yum -y install vim* Centos7安装ifconfig yum instal ...
- c#之初识结构(Struct)
C# 结构(Struct) 首先结构是值类型数据结构.它使得一个单一变量可以存储各种数据类型的相关数据.struct 关键字用于创建结构.通俗说:结构就是一个可以包含不同数据类型的集合.它是一种可以自 ...
- springboot的http监控接口启动器的配置
基于SpringBoot框架企业级应用系统开发全面实战()->03.07_http监控_recv.mp4 监控接口启动器 自定义监控接口启动器的配置 ====================== ...