poj 1986 Distance Queries
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相关推荐
- poj 1986 Distance Queries LCA
题目链接:http://poj.org/problem?id=1986 Farmer John's cows refused to run in his marathon since he chose ...
- POJ 1986 Distance Queries(LCA)
[题目链接] http://poj.org/problem?id=1986 [题目大意] 给出一棵树,问任意两点间距离. [题解] u,v之间距离为dis[u]+dis[v]-2*dis[LCA(u, ...
- POJ - 1986 Distance Queries 倍增求LCA
题意还是很有病的,说了半天后边的方向都没用的,意思就是跟样例一样求两点间的距离,裸的LCA 两个点x,y的距离=dis[x]+dis[y]-2*dis[LCA(x,y)]; 三个点x,y,z的距离=d ...
- 【POJ No. 1986】 距离查询 Distance Queries
[POJ No. 1986] 距离查询 Distance Queries 北大OJ 题目地址 [题意] 约翰有N 个农场,标记为1-N .有M 条垂直和水平的道路连接农场,每条道路的长度各不相同.每个 ...
- Open3D Distance Queries 距离查询
Distance Queries 距离查询 Open3D 中的类RaycastingScene提供了一组距离查询,可用于将三角形网格转换为隐式函数.查询到表面的距离或确定点是否位于网格内. 在本教程中 ...
- POJ 1986:Distance Queries(倍增求LCA)
http://poj.org/problem?id=1986 题意:给出一棵n个点m条边的树,还有q个询问,求树上两点的距离. 思路:这次学了一下倍增算法求LCA.模板. dp[i][j]代表第i个点 ...
- poj1986 Distance Queries(lca又是一道模版题)
题目链接:http://poj.org/problem?id=1986 题意:就是老问题求val[u]+val[v]-2*val[root]就行.还有这题没有给出不联通怎么输出那么题目给出的数据一定 ...
- poj pku图论、网络流入门题总结、汇总
poj pku图论.网络流入门题总结.汇总 分类: acm图论 2010-08-25 18:49 243人阅读 评论(0) 收藏 举报 网络算法networkgraphconstructioninte ...
- POJ 图论分类 + DP(较全 自己又加了点)
DP -----------动态规划 状态压缩DP 2411 (棋盘规模较大)状态压缩DP+DFS+滚动数组 2664 (棋盘规模较小)直接递推即可(DP) 2506 (棋盘规模较小)直接递推即可(D ...
最新文章
- 使用Caffe基于cifar10进行物体识别
- ICML2020 | 一行代码就能实现的测试集上分技巧
- Opencv之通过url抓取图片并通过opencv可视化
- R count函数_【函数分享】PHP函数str_word_count()分享(202098)
- c语言酒店管理系统,基于C#的酒店管理系统(V3.1)最新版
- IDEA 之搭建spring-boot maven报错Project ‘org.springframework.boot:spring-boot-starter-parent:2.2.0.RELEAS
- 海底捞成功的全套培训体系(收藏)
- 本人GitHub网址https://github.com/LH34128
- edvac是商用计算机吗,计算机的基础作业1
- BCH编码(15,5)1bit串行,5bit并行
- 总有人问 IT 行业有发展前景吗?那就用我的经历来告诉你
- SQL Server中的六种约束:主键约束,外键约束,唯一约束,非空约束,检查约束,默认约束
- python 验证码图片 模拟登录_【python】带图片验证码的登录自动化实战
- IBM Thinkpad 笔记本型号和保修(包括水货)全揭秘
- JAVA学习56_用UltraEdit代替“笨重”的IDE,实现轻巧编程!
- 如何将浏览器设置对应的国家?
- 【翻译】TypeScript: Documentation - Narrowing
- java.lang.NumberFormatException 数字转换异常的错误及解决办法
- Window2012 R2下 IIS的详细配置
- 理论计算机科学逻辑博导,软件学院研究生论文导师一览表