LCA

题意:LCA模板题,输入n和m,表示n个点m条边,下面m行是边的信息,两端点和权,后面的那个字母无视掉,没用的。接着k,下面k个询问lca,输出即可

有人说要考虑不连通的情况,我没考虑AC了,另外可能有u,u这样的询问,不过这不影响,照样是写模板,没有特判,一样能过

还是Tarjan快一些

LCA转RMQ在线算法

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;
#define N 40010
#define M 25int tot;
int __pow[M];
int head[N];
struct edge{int u,v,w,next;
}e[2*N];
int ver[2*N],R[2*N],first[N],dir[N];
int dp[2*N][25];
bool vis[N];inline void add(int u , int v ,int w ,int &k)
{e[k].u = u; e[k].v = v; e[k].w = w;e[k].next = head[u]; head[u] = k++;u = u^v; v = u^v; u = u^v;e[k].u = u; e[k].v = v; e[k].w = w;e[k].next = head[u]; head[u] = k++;
}void dfs(int u ,int dep)
{vis[u] = true; first[u] = ++tot; ver[tot] = u; R[tot] = dep;for(int k=head[u]; k!=-1; k=e[k].next)if( !vis[e[k].v] ){int v = e[k].v , w = e[k].w;dir[v] = dir[u] + w;dfs(v,dep+1);ver[++tot] = u; R[tot] = dep;}
}void ST(int len)
{int K = (int)(log((double)(len)) / log(2.0));for(int i=1; i<=len; i++) dp[i][0] = i;for(int j=1; j<=K; j++)for(int i=1; i+__pow[j]-1 <= len; i++){int a = dp[i][j-1] , b = dp[i+__pow[j-1]][j-1];if(R[a] < R[b]) dp[i][j] = a;else            dp[i][j] = b;}
}int RMQ(int x ,int y)
{int K = (int)(log((double)(y-x+1)) / log(2.0));int a = dp[x][K] , b = dp[y-__pow[K]+1][K];if(R[a] < R[b]) return a;else            return b;
}int LCA(int u ,int v)
{int x = first[u] , y = first[v];if(x > y) swap(x,y);int index = RMQ(x,y);return ver[index];
}int main()
{for(int i=0; i<M; i++) __pow[i] = (1<<i);int n,m,k,str[10];while(scanf("%d%d",&n,&m)!=EOF){k = 0;memset(head,-1,sizeof(head));memset(vis,false,sizeof(vis));while(m--){int u,v,w;scanf("%d%d%d%s",&u,&v,&w,str);add(u,v,w,k);}tot = dir[1] = 0;dfs(1,1);ST(tot);int q;scanf("%d",&q);while(q--){int u,v,lca;scanf("%d%d",&u,&v);lca = LCA(u,v);printf("%d\n",dir[u] + dir[v] - 2*dir[lca]);}}return 0;
}

Tarjan离线算法

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
#define N 40010
#define M 20010int head[N];
struct edge{int u,v,w,next;
}e[2*N];
int __head[N];
struct ask{int u,v,lca,next;
}ea[M];
int fa[N],ance[N],dir[N];
bool vis[N];inline void add_edge(int u ,int v ,int w ,int &k)
{e[k].u = u; e[k].v = v; e[k].w = w;e[k].next = head[u]; head[u] = k++;u = u^v; v = u^v; u = u^v;e[k].u = u; e[k].v = v; e[k].w = w;e[k].next = head[u]; head[u] = k++;
}inline void add_ask(int u ,int v ,int &k)
{ea[k].u = u; ea[k].v = v; ea[k].lca = -1;ea[k].next = __head[u]; __head[u] = k++;u = u^v; v = u^v; u = u^v;ea[k].u = u; ea[k].v = v; ea[k].lca = -1;ea[k].next = __head[u]; __head[u] = k++;
}int find(int x){return x == fa[x] ? x : fa[x] = find(fa[x]);
}void Tarjan(int u)
{vis[u] = true;ance[u] = fa[u] = u;for(int k=head[u]; k!=-1; k=e[k].next)if( !vis[e[k].v] ){int v = e[k].v , w = e[k].w;dir[v] = dir[u] + w;Tarjan(v);fa[v] = u;}for(int k=__head[u]; k!=-1; k=ea[k].next)if( vis[ea[k].v] )ea[k].lca = ea[k^1].lca = ance[find(ea[k].v)];
}int main()
{int n,m,q,k; char str[10];while(scanf("%d%d",&n,&m)!=EOF){memset(head,-1,sizeof(head));memset(__head,-1,sizeof(__head));memset(vis,false,sizeof(vis));k = 0;while(m--){int u,v,w;scanf("%d%d%d%s",&u,&v,&w,str);add_edge(u,v,w,k);}scanf("%d",&q);k = 0;for(int i=0; i<q; i++){int u ,v;scanf("%d%d",&u,&v);add_ask(u,v,k);}dir[1] = 0;Tarjan(1);for(int i=0; i<q; i++){int s = i*2 , u = ea[s].u , v = ea[s].v , lca = ea[s].lca;printf("%d\n",dir[u] + dir[v] - 2*dir[lca]);}}return 0;
}

转载于:https://www.cnblogs.com/scau20110726/archive/2013/05/27/3102068.html

poj 1986 Distance Queries相关推荐

  1. poj 1986 Distance Queries LCA

    题目链接:http://poj.org/problem?id=1986 Farmer John's cows refused to run in his marathon since he chose ...

  2. POJ 1986 Distance Queries(LCA)

    [题目链接] http://poj.org/problem?id=1986 [题目大意] 给出一棵树,问任意两点间距离. [题解] u,v之间距离为dis[u]+dis[v]-2*dis[LCA(u, ...

  3. POJ - 1986 Distance Queries 倍增求LCA

    题意还是很有病的,说了半天后边的方向都没用的,意思就是跟样例一样求两点间的距离,裸的LCA 两个点x,y的距离=dis[x]+dis[y]-2*dis[LCA(x,y)]; 三个点x,y,z的距离=d ...

  4. 【POJ No. 1986】 距离查询 Distance Queries

    [POJ No. 1986] 距离查询 Distance Queries 北大OJ 题目地址 [题意] 约翰有N 个农场,标记为1-N .有M 条垂直和水平的道路连接农场,每条道路的长度各不相同.每个 ...

  5. Open3D Distance Queries 距离查询

    Distance Queries 距离查询 Open3D 中的类RaycastingScene提供了一组距离查询,可用于将三角形网格转换为隐式函数.查询到表面的距离或确定点是否位于网格内. 在本教程中 ...

  6. POJ 1986:Distance Queries(倍增求LCA)

    http://poj.org/problem?id=1986 题意:给出一棵n个点m条边的树,还有q个询问,求树上两点的距离. 思路:这次学了一下倍增算法求LCA.模板. dp[i][j]代表第i个点 ...

  7. poj1986 Distance Queries(lca又是一道模版题)

    题目链接:http://poj.org/problem?id=1986 题意:就是老问题求val[u]+val[v]-2*val[root]就行.还有这题没有给出不联通怎么输出那么题目给出的数据一定 ...

  8. poj pku图论、网络流入门题总结、汇总

    poj pku图论.网络流入门题总结.汇总 分类: acm图论 2010-08-25 18:49 243人阅读 评论(0) 收藏 举报 网络算法networkgraphconstructioninte ...

  9. POJ 图论分类 + DP(较全 自己又加了点)

    DP -----------动态规划 状态压缩DP 2411 (棋盘规模较大)状态压缩DP+DFS+滚动数组 2664 (棋盘规模较小)直接递推即可(DP) 2506 (棋盘规模较小)直接递推即可(D ...

最新文章

  1. 使用Caffe基于cifar10进行物体识别
  2. ICML2020 | 一行代码就能实现的测试集上分技巧
  3. Opencv之通过url抓取图片并通过opencv可视化
  4. R count函数_【函数分享】PHP函数str_word_count()分享(202098)
  5. c语言酒店管理系统,基于C#的酒店管理系统(V3.1)最新版
  6. IDEA 之搭建spring-boot maven报错Project ‘org.springframework.boot:spring-boot-starter-parent:2.2.0.RELEAS
  7. 海底捞成功的全套培训体系(收藏)
  8. 本人GitHub网址https://github.com/LH34128
  9. edvac是商用计算机吗,计算机的基础作业1
  10. BCH编码(15,5)1bit串行,5bit并行
  11. 总有人问 IT 行业有发展前景吗?那就用我的经历来告诉你
  12. SQL Server中的六种约束:主键约束,外键约束,唯一约束,非空约束,检查约束,默认约束
  13. python 验证码图片 模拟登录_【python】带图片验证码的登录自动化实战
  14. IBM Thinkpad 笔记本型号和保修(包括水货)全揭秘
  15. JAVA学习56_用UltraEdit代替“笨重”的IDE,实现轻巧编程!
  16. 如何将浏览器设置对应的国家?
  17. 【翻译】TypeScript: Documentation - Narrowing
  18. java.lang.NumberFormatException 数字转换异常的错误及解决办法
  19. Window2012 R2下 IIS的详细配置
  20. 理论计算机科学逻辑博导,软件学院研究生论文导师一览表

热门文章

  1. 几个用于序列化的代码片段
  2. 世界领先的界面设计公司:The Skins Factory
  3. 华为服务器如何登录修改密码,如何修改云服务器的登录密码
  4. 对计算机网络用户而言 掌握网络,计算机网络的特点
  5. android新拟态实现方法,Android 新拟态UI (Neumorphism)
  6. 贵州大学java期末考试_CONTRIBUTING.md
  7. 华为鸿蒙系统技术细节盘点
  8. Vim和Vi的常用命令
  9. 常用 SQL Server 规范集锦
  10. mysql开方_MySQL数学函数的实际用法