G. Unusual Entertainment

G. Unusual Entertainment

A tree is a connected graph without cycles.

A permutation is an array consisting of n distinct integers from 1 to n in any order. For example, [5,1,3,2,4] is a permutation, but [2,1,1] is not a permutation (as 1 appears twice in the array) and [1,3,2,5] is also not a permutation (as n=4, but 5 is present in the array).

After a failed shoot in the BrMeast video, Alex fell into depression. Even his birthday did not make him happy. However, after receiving a gift from Timofey, Alex's mood suddenly improved. Now he spent days playing with the gifted constructor. Recently, he came up with an unusual entertainment.

Alex builds a tree from his constructor, consisting of n vertices numbered from 1 to n, with the root at vertex 1. Then he writes down each integer from 1 to n in some order, obtaining a permutation p. After that, Alex comes up with q triples of integers l,r,x. For each triple, he tries to determine if there is at least one descendant of vertex x among the vertices pl,pl+1,,pr.

A vertex u is a descendant of vertex v if and only if dist(1,v)+dist(v,u)=dist(1,u), where dist(a,b) is the distance between vertices a and b. In other words, vertex v must be on the path from the root to vertex u.

Alex told Zakhar about this entertainment. Now Alex tells his friend q triples as described above, hoping that Zakhar can check for the presence of a descendant. Zakhar is very sleepy, so he turned to you for help. Help Zakhar answer all of Alex's questions and finally go to sleep.

Input

The first line of the input contains a single integer t (1t104) — the number of test cases.

The first line of each test case contains two integers n,q (1n,q105) — the number of vertices in the tree and the number of questions, respectively.

Each of the next n1 lines contains two integers ui and vi (1ui,vin), indicating that there is an edge between vertices ui and vi (it is guaranteed that the resulting graph is a tree).

The next line contains n integers p1,p2,,pn (1pin) — the permutation p (it is guaranteed that each integer from 1 to n appears exactly once).

Then follow q lines describing Alex's questions. The i-th line contains three integers l,r,x (1lrn, 1xn), as described in the statement.

It is guaranteed that the sum of n and the sum of q over all test cases do not exceed 105.

Output

For each of Alex's questions, print "Yes" (without quotes) if the described descendant exists, otherwise print "No" (without quotes).

You can output the answer in any case (for example, the strings "yEs", "yes", "Yes", and "YES" will be recognized as a positive answer).

Example

input

3
3 5
1 2
2 3
1 2 3
1 2 2
1 2 3
2 3 1
1 2 3
2 3 3
10 10
2 6
2 7
2 4
1 7
2 8
10 6
8 5
9 4
3 4
10 2 5 9 1 7 6 4 3 8
8 9 8
7 8 1
7 10 6
4 8 9
5 5 10
7 10 1
9 9 2
9 10 6
6 6 2
10 10 6
1 1
1
1 1 1

output

YES
NO
YES
NO
YES

NO
YES
YES
YES
NO
YES
YES
NO
NO
NO

YES

 

解题思路

  问 plpr 中是否存在一个节点是 x 的子节点,等价于问 plpr 中是否存在一个节点出现在以 x 为根的子树中。

  为此可以考虑 dfs 序,如果一个节点 i 在以 x 为根的子树中,那么必然有 tinxtinpitoutpitoutx。所以问题就变成了是否存在一个 i 同时满足 {lirtinxtinpitoutx

  这就是一个二维数点问题,可以参考这篇博客。把第一、二个约束分别看作是横坐标和纵坐标的某个区间范围,那么点的坐标就是 (i,tinpi),每次询问就相当于求左下角为 (l,tinx) 和右上角为 (r,toutx) 这个矩形内点的数量是否大于 0

  AC 代码如下,时间复杂度为 O(qlogq+(n+q)logn)

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;

const int N = 1e5 + 10, M = N * 2;

int n, m;
int head[N], e[M], ne[M], idx;
int tin[N], tout[N], sz;
int p[N];
struct Node {
    int x, y, c, idx;
}q[N * 4];
int tr[N];
int ans[N];

void add(int u, int v) {
    e[idx] = v, ne[idx] = head[u], head[u] = idx++;
}

void dfs(int u, int pre) {
    tin[u] = ++sz;
    for (int i = head[u]; i != -1; i = ne[i]) {
        if (e[i] != pre) {
            dfs(e[i], u);
        }
    }
    tout[u] = sz;
}

int lowbit(int x) {
    return x & -x;
}

void modify(int x, int c) {
    for (int i = x; i <= n; i += lowbit(i)) {
        tr[i] += c;
    }
}

int query(int x) {
    int ret = 0;
    for (int i = x; i; i -= lowbit(i)) {
        ret += tr[i];
    }
    return ret;
}

void solve() {
    scanf("%d %d", &n, &m);
    memset(head, -1, n + 10 << 2);
    idx = 0;
    for (int i = 0; i < n - 1; i++) {
        int u, v;
        scanf("%d %d", &u, &v);
        add(u, v), add(v, u);
    }
    for (int i = 1; i <= n; i++) {
        scanf("%d", p + i);
    }
    sz = 0;
    dfs(1, -1);
    for (int i = 0, j = 0; i < m; i++) {
        int l, r, x;
        scanf("%d %d %d", &l, &r, &x);
        int x1 = l, y1 = tin[x], x2 = r, y2 = tout[x];
        q[j++] = {x2, y2, 1, i};
        q[j++] = {x1 - 1, y1 - 1, 1, i};
        q[j++] = {x1 - 1, y2, -1, i};
        q[j++] = {x2, y1 - 1, -1, i};
    }
    sort(q, q + 4 * m, [&](Node &a, Node &b) {
        return a.x < b.x;
    });
    memset(tr, 0, n + 10 << 2);
    memset(ans, 0, n + 10 << 2);
    for (int i = 0, j = 1; i < m << 2; i++) {
        while (j <= n && j <= q[i].x) {
            modify(tin[p[j++]], 1);
        }
        ans[q[i].idx] += q[i].c * query(q[i].y);
    }
    for (int i = 0; i < m; i++) {
        printf("%s\n", ans[i] ? "YES" : "NO");
    }
    printf("\n");
}

int main() {
    int t;
    scanf("%d", &t);
    while (t--) {
        solve();
    }
    
    return 0;
}

  补充用主席树的做法。

  用权值线段树来维护点 (i,tinpi) 的纵坐标的出现次数。依次从 1n 遍历 i,进行单点修改(tinpi 的出现次数加 1),维护出 n 个版本的线段树,root[i] 表示只维护前 i 个点所得到线段树的根节点。那么对于询问 (l,r,x),由于满足可减性,root[r]root[l1] 的版本中相同值域区间对应的节点的信息差就是 [l,r] 中的点的信息,只需根据这两个版本的信息差来查询值域 [tinx,toutx] 内是否至少有一次记录即可。

  AC 代码如下,时间复杂度为 O((n+q)logn)

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;

const int N = 1e5 + 10, M = N * 2;

int head[N], e[M], ne[M], idx;
int tin[N], tout[N], sz;
int p[N];
struct Node {
    int l, r, s;
}tr[N * 20];
int root[N], tot;

void add(int u, int v) {
    e[idx] = v, ne[idx] = head[u], head[u] = idx++;
}

void dfs(int u, int pre) {
    tin[u] = ++sz;
    for (int i = head[u]; i != -1; i = ne[i]) {
        if (e[i] != pre) dfs(e[i], u);
    }
    tout[u] = sz;
}

int build(int l, int r) {
    int u = ++tot;
    if (l == r) return u;
    int mid = l + r >> 1;
    tr[u].l = build(l, mid);
    tr[u].r = build(mid + 1, r);
    tr[u].s = 0;
    return u;
}

int insert(int u, int l, int r, int x) {
    int v = ++tot;
    tr[v] = tr[u];
    if (l == r) {
        tr[v].s++;
        return v;
    }
    int mid = l + r >> 1;
    if (x <= mid) tr[v].l = insert(tr[u].l, l, mid, x);
    else tr[v].r = insert(tr[u].r, mid + 1, r, x);
    tr[v].s = tr[tr[v].l].s + tr[tr[v].r].s;
    return v;
}

int query(int u, int v, int l, int r, int ql, int qr) {
    if (l >= ql && r <= qr) return tr[v].s - tr[u].s;
    int mid = l + r >> 1, ret = 0;
    if (ql <= mid) ret = query(tr[u].l, tr[v].l, l, mid, ql, qr);
    if (qr >= mid + 1) ret += query(tr[u].r, tr[v].r, mid + 1, r, ql, qr);
    return ret;
}

void solve() {
    int n, m;
    scanf("%d %d", &n, &m);
    memset(head, -1, n + 10 << 2);
    idx = 0;
    for (int i = 0; i < n - 1; i++) {
        int u, v;
        scanf("%d %d", &u, &v);
        add(u, v), add(v, u);
    }
    for (int i = 1; i <= n; i++) {
        scanf("%d", p + i);
    }
    sz = 0;
    dfs(1, -1);
    tot = 0;
    root[0] = build(1, n);
    for (int i = 1; i <= n; i++) {
        root[i] = insert(root[i - 1], 1, n, tin[p[i]]);
    }
    while (m--) {
        int l, r, x;
        scanf("%d %d %d", &l, &r, &x);
        printf("%s\n", query(root[l - 1], root[r], 1, n, tin[x], tout[x]) ? "YES" : "NO");
    }
    printf("\n");
}

int main() {
    int t;
    scanf("%d", &t);
    while (t--) {
        solve();
    }
    
    return 0;
}

 

参考资料

  Разбор Codeforces Round 909 (Div. 3):https://codeforces.com/blog/entry/122407

  Codeforces Round 909 (Div. 3) A-G 讲解:https://www.bilibili.com/video/BV1Uu4y1b7Dc/

posted @   onlyblues  阅读(45)  评论(1编辑  收藏  举报
相关博文:
阅读排行:
· 单线程的Redis速度为什么快?
· 展开说说关于C#中ORM框架的用法!
· Pantheons:用 TypeScript 打造主流大模型对话的一站式集成库
· SQL Server 2025 AI相关能力初探
· 为什么 退出登录 或 修改密码 无法使 token 失效
历史上的今天:
2022-11-21 B. Elimination of a Ring
Web Analytics
点击右上角即可分享
微信分享提示