How far away

&题解:

和上篇是一样的题,这用的是lca方法做的, 不知道为什么,把数组开到80000 就a了 >_<

哈 我现在知道为什么了,因为我的rmq数组没有乘2,rmq比较的是depth数组,但是depth数组的范围是maxn*2,所以rmq的数组乘2就好了

&代码:

#include <cstdio>
#include <bitset>
#include <iostream>
#include <set>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <map>
#include <queue>
#include <vector>
using namespace std;
#define INF 0x3f3f3f3f
#define ll long long
#define fo(i,a,b) for(int i=(a);i<=(b);i++)
#define fd(i,a,b) for(int i=(a);i>=(b);i--) const int max_v = 40000 + 7;
int rq[max_v*2][20], mm[max_v*2], n, m; void initRMQ(int n, int b[]) {
mm[0] = -1;
for(int i = 1; i <= n; i++) {
mm[i] = (i & (i - 1)) ? mm[i - 1] : mm[i - 1] + 1;
rq[i][0] = i;
}
for(int j = 1; j <= mm[n]; j++)
for(int i = 1; i + (1 << j) - 1 <= n; i++) {
int x = rq[i][j - 1], y = rq[i + (1 << (j - 1))][j - 1];
rq[i][j] = b[x] < b[y] ? x : y;
}
} int RMQ(int x, int y, int b[]) {
if(x > y) swap(x, y);
int k = mm[y - x + 1];
int x2 = rq[x][k], y2 = rq[y - (1 << k) + 1][k];
return b[x2] < b[y2] ? x2 : y2;
} int dis[max_v];
struct Edge {
int u, v, w, nx;
} G[2 * max_v];
int root;
int vs[max_v * 2];
int depth[max_v * 2];
int id[max_v]; int cnt, ft[max_v * 2], vis[max_v]; void Add(int u, int v, int w) {
G[cnt].u = u; G[cnt].v = v; G[cnt].w = w;
G[cnt].nx = ft[u];
ft[u] = cnt++;
} void dfs(int v, int p, int d, int &k) {
id[v] = k;
vs[k] = v;
depth[k++] = d;
// vis[v] = 1;
for(int i = ft[v]; ~i; i = G[i].nx) {
// printf("%d %d-- %d\n", G[i].u, G[i].v, p);
if(G[i].v != p) {
dis[G[i].v] = dis[G[i].u] + G[i].w;
dfs(G[i].v, G[i].u, d + 1, k);
vs[k] = v;
depth[k++] = d;
}
}
} void initLCA() {
int k = 1;
memset(vis, 0, sizeof(vis));
dfs(root, -1, 0, k);
initRMQ(k - 1, depth);
} int LCA(int u, int v) {
return vs[RMQ(min(id[u], id[v]), max(id[u], id[v]) + 1, depth)];
} int main() {
//("E:1.in", "r", stdin);
int T; scanf("%d", &T);
while(T--) {
scanf("%d%d", &n, &m);
cnt = 0;
memset(ft, -1, sizeof(ft));
fo(i, 1, n - 1) {
int x, y, val;
scanf("%d%d%d", &x, &y, &val);
Add(x, y, val);
Add(y, x, val);
}
root = 1;
dis[root] = 0;
initLCA();
fo(i, 1, m) {
int u, v;
scanf("%d%d", &u, &v);
printf("%d\n", dis[u] + dis[v] - 2 * dis[LCA(u, v)]);
}
}
return 0;
}

这个是用临街矩阵写的,RT了,但有注释

#include <cstdio>
#include <bitset>
#include <iostream>
#include <set>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <map>
#include <queue>
#include <vector>
using namespace std;
#define INF 0x3f3f3f3f
#define ll long long
#define fo(i,a,b) for(int i=(a);i<=(b);i++)
#define fd(i,a,b) for(int i=(a);i>=(b);i--)
#define pii pair<int,int> //RMQ 相关
const int maxn = 40000 + 7, max_v = maxn;
int dp[maxn][20], mm[maxn], n, m; void initRMQ(int n, int b[]) {
mm[0] = -1;
for(int i = 1; i <= n; i++) {
mm[i] = (i & (i - 1)) ? mm[i - 1] : mm[i - 1] + 1;
dp[i][0] = i;
}
for(int j = 1; j <= mm[n]; j++)
for(int i = 1; i + (1 << j) - 1 <= n; i++) {
int x = dp[i][j - 1], y = dp[i + (1 << (j - 1))][j - 1];
dp[i][j] = b[x] < b[y] ? x : y;
// dp[i][j] = min(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1]);
}
} int RMQ(int x, int y, int b[]) {
if(x > y) swap(x, y);
int k = mm[y - x + 1];
int x2 = dp[x][k], y2 = dp[y - (1 << k) + 1][k];
return b[x2] < b[y2] ? x2 : y2;
} int dis[maxn];
//LCA 相关
vector<pii> G[max_v]; //图的邻接表
int root;
int vs[max_v * 2]; //dfs访问的顺序
int depth[max_v * 2]; //节点深度
int id[max_v]; //各顶点在vs中首次出现的下标 //v:现节点 p:父节点 d:深度 k:编号
void dfs(int v, int p, int d, int &k) {
id[v] = k;
vs[k] = v;
depth[k++] = d;
for(int i = 0; i < G[v].size(); i++) {
if(G[v][i].first != p) {
dis[G[v][i].first] = dis[v] + G[v][i].second;
dfs(G[v][i].first, v, d + 1, k);
vs[k] = v;
depth[k++] = d;
}
}
} void initLCA(int V) {
int k = 1;
dfs(root, -1, 0, k);
initRMQ(k - 1, depth);
} int LCA(int u, int v) {
return vs[RMQ(min(id[u], id[v]), max(id[u], id[v]) + 1, depth)];
} int main() {
//("E:1.in", "r", stdin);
int T; scanf("%d", &T);
while(T--) {
scanf("%d%d", &n, &m);
fo(i, 1, n) G[i].clear();
fo(i, 1, n - 1) {
int x, y, val;
scanf("%d%d%d", &x, &y, &val);
G[x].push_back(make_pair(y, val));
G[y].push_back(make_pair(x, val));
}
root = 1;
dis[root] = 0;
initLCA(n);
fo(i, 1, m) {
int u, v;
scanf("%d%d", &u, &v);
printf("%d\n", dis[u] + dis[v] - 2 * dis[LCA(u, v)]);
}
// for(int i = 0; i < 8; i++) {
// printf("%3d", i);
// } printf("\n");
// for(int i = 0; i < 8; i++) {
// printf("%3d", dis[i]);
// } printf("\n");
}
return 0;
}

HDU 2586 How far away(LCA+邻接表)的更多相关文章

  1. HDU 2586 How far away(dfs+邻接表)

    How far away [题目链接]How far away [题目类型]dfs+邻接表 &题意: 题目大意:一个村子里有n个房子,这n个房子用n-1条路连接起来,接下了有m次询问,每次询问 ...

  2. hdu 2586(最近公共祖先LCA)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2586 思路:在求解最近公共祖先的问题上,用到的是Tarjan的思想,从根结点开始形成一棵深搜树,非常好 ...

  3. HDU 2586 How far away? LCA 转化成RMQ

    链接:http://acm.hdu.edu.cn/showproblem.php?pid=2586 [题意] 给出一个N 个和N-1条边的连通图,询问任意两点间的距离.N<=40000 . [分 ...

  4. HDU 2586 倍增法求lca

    How far away ? Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)To ...

  5. TTTTTTTTTTTTTTTTT HDU 2586 How far away LCA的离线算法 Tarjan

    链接: How far away ? Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Other ...

  6. 【HDU 2586 How far away?】LCA问题 Tarjan算法

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2586 题意:给出一棵n个节点的无根树,每条边有各自的权值.给出m个查询,对于每条查询返回节点u到v的最 ...

  7. HDU 2586 How far away ?【LCA】

    任意门:http://acm.hdu.edu.cn/showproblem.php?pid=2586 How far away ? Time Limit: 2000/1000 MS (Java/Oth ...

  8. LCA(最近公共祖先)--tarjan离线算法 hdu 2586

    HDU 2586 How far away ? Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/ ...

  9. HDU 4858 项目管理(邻接表 暴力模拟)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4858 我们建造了一个大项目!这个项目有n个节点,用很多边连接起来,并且这个项目是连通的! 两个节点间可 ...

随机推荐

  1. windows下IDEA的terminal配置bash命令

    使用git-bash.exe会单独打开一个窗口,而我们希望是在终端内置的命令行.这里我使用bash.exe 在IDEA中,打开settings,设置相应的bash路径 settings–>Too ...

  2. [译] 理解 LSTM(Long Short-Term Memory, LSTM) 网络

    本文译自 Christopher Olah 的博文 Recurrent Neural Networks 人类并不是每时每刻都从一片空白的大脑开始他们的思考.在你阅读这篇文章时候,你都是基于自己已经拥有 ...

  3. DirectX using C++_error X3539:ps1_x is no longer supported...解决方案

    问题来源 在研究HLSL时编译一个demo出现了error X3539的问题 解决方案 将代码中的ps_1_1 改为ps_2_0 PixelShader = compile ps_1_1 PS(); ...

  4. php -- 类对象调用静态方法

    以前一直以为 静态方法的调用:类名::静态方法 非静态方法的调用:类对象->非静态方法 最近研究一个类,发现一个比较奇怪的问题,用“类对象->静态方法”这种方式居然成功的调用了静态方法.很 ...

  5. Mybatis sql映射文件浅析 Mybatis简介(三) 简介

    Mybatis sql映射文件浅析 Mybatis简介(三)   简介 除了配置相关之外,另一个核心就是SQL映射,MyBatis 的真正强大也在于它的映射语句. Mybatis创建了一套规则以XML ...

  6. ansible运维工具(一)

    运维工具介绍 OS Provisioning: PXE, Cobbler(repository, distritution,profile) PXE: dhcp, tftp, (http, ftp) ...

  7. 我的预约订单页面List

    <%@ page language="java" contentType="text/html;charset=UTF-8"%> <%@ ta ...

  8. numpy和tensorflow中的广播机制

    广播的引出 numpy两个数组的相加.相减以及相乘都是对应元素之间的操作. import numpy as np x = np.array([[2,2,3],[1,2,3]]) y = np.arra ...

  9. LeetCode - 804. Unique Morse Code Words

    International Morse Code defines a standard encoding where each letter is mapped to a series of dots ...

  10. Nest.js 守卫

    Docs: https://docs.nestjs.com/guards 当调用者具有足够的权限时,特定路由才可用 // app.guard.ts import { CanActivate, Exec ...