#	User	Rating
1	tourist	3985
2	jiangly	3814
3	jqdai0815	3682
4	Benq	3529
5	orzdevinwang	3526
6	ksun48	3517
7	Radewoosh	3410
8	hos.lyric	3399
9	ecnerwala	3392
9	Um_nik	3392

#	User	Contrib.
1	cry	169
2	maomao90	162
2	Um_nik	162
4	atcoder_official	161
5	djm03178	158
6	-is-this-fft-	157
7	adamant	155
8	awoo	154
8	Dominater069	154
10	luogu_official	150

It is not uncommon to have interactive tree problems where you are allowed to query some connected component of the tree and use the return value to determine whether the answer is in the connected component or outside the connected component (Link to example problem). The general approach for these kinds of problems is to always choose a connected component of size $$$\frac{n}{2}$$$. However, there are also problems where the allowed queries are more restricted, preventing $$$\frac{n}{2}$$$ halving from being possible. This blog covers one of those types of problems.

Definitions

Subtree: $$$S(r, u)$$$ contains the set of vertices in the subtree of vertex $$$u$$$ if the tree was rooted at vertex $$$r$$$.
Neighbour: $$$N(u)$$$ contains the set of vertices that are directly adjacent to vertex $$$u$$$.
Extended subtree: $$$ES(r, V) = \bigcup_{v\in V} S(r, v)\text{ if } v\in N(r)$$$. In other words, an extended subtree is a combination of the subtrees of a chosen set of vertices that are directly adjacent to the root.

Problem Structure

There is a hidden special vertex in a tree with $$$n$$$ vertices. Find the special vertex using at most $$$\lceil\log_{1.5}n\rceil$$$ of the following query:

Choose an extended subtree of the tree. The grader will return whether the special vertex is in the chosen extended subtree. More formally, choose any vertex $$$r$$$ and a subset of neighbours $$$V \subseteq N(r)$$$, then the grader will return whether the special vertex $$$x \in ES(r, V)$$$.

Solution

Let us denote the size of the chosen extended subtree as $$$x$$$. If the grader returns that the special vertex is in the extended subtree, we can combine the root with the extended subtree to form a new tree with $$$x + 1$$$ vertices. Otherwise, the vertices outside of the extended subtree form a tree with $$$n - x$$$ vertices. This means that with each query, we can reduce the size of the tree to at least $$$\max(x + 1, n - x)$$$.

Since we are allowed to use $$$\lceil\log_{1.5}n\rceil$$$ queries, we will be able to solve this problem if the size of the tree reduces to at least $$$\frac{2n}{3}$$$ after each query. Solving the inequality $$$\max(x + 1, n - x) \le \frac{2n}{3}$$$, we get $$$\frac{n}{3} \le x \le \frac{2n}{3} - 1$$$. Now, let us see whether it is always possible to find an extended subtree of size between $$$\frac{n}{3}$$$ and $$$\frac{2n}{3} - 1$$$.

Let us root the tree at its centroid. Recall that all the children of the centroid have subtree sizes less than or equal to $$$\frac{n}{2}$$$. If there exists a child with subtree size more than or equal to $$$\frac{n}{3}$$$, then that lone subtree can be the chosen extended subtree since $$$\frac{n}{3} \le x \le \frac{n}{2} \le \frac{2n}{3} - 1$$$ $$$^\dagger$$$. From this point onwards, we will assume that all children of the centroid have subtree sizes less than $$$\frac{n}{3}$$$.

We will add the children of the centroid one by one into the extended subtree until the first time that the size of the extended subtree becomes more than or equal to $$$\frac{n}{3}$$$. Note that this is always possible as taking all children will result in a size of $$$n - 1$$$ which is more than or equal to $$$\frac{n}{3}$$$. Let $$$x$$$ be the size of the resultant extended subtree and let $$$p$$$ be the size of the last child that was added into the extended subtree. We can obtain the following inequalities: $$$x \ge \frac{n}{3}$$$, $$$x - p < \frac{n}{3}$$$, and $$$p < \frac{n}{3}$$$. Combining the second and third inequality, we obtain $$$x < \frac{2n}{3}$$$. Hence, $$$\frac{n}{3} \le x \le \frac{2n}{3} - 1$$$ and the desired extended subtree is found.

$$$^\dagger$$$ Note that $$$\frac{n}{2} \le \frac{2n}{3} - 1$$$ might not be true when $$$n \le 5$$$, however, this is not very important as $$$n - 1 \le \lceil \log_{1.5} n \rceil$$$ when $$$n \le 5$$$, so we can just query every subtree except for vertex 1 to get the answer.

Optimality

It might feel like $$$\lceil\log_{1.5}n\rceil$$$ is a loose bound and we can do better. However, we can prove that it is not possible to do any better. Special thanks to pavement and jamessngg for helping me to come up with this proof. Consider the tree below.

Jellyfish tree

There are $$$n = 3k + 1$$$ vertices in the tree. The closest we can get to an extended subtree of size $$$\frac{n}{2}$$$ is either a single chain of size $$$k$$$ or two chains having size $$$2k$$$. Both of those options will result in a tree of size $$$2k + 1 \approx \frac{2n}{3}$$$ in the worst case after a query. This shows that for all $$$n = 3k + 1$$$, there exists a tree where it is not possible to reduce the size of the tree by less than $$$\frac{2n}{3}$$$ in a single query. We can use similar trees but with different chain lengths to come up with a similar proof for $$$n = 3k + 2$$$ and $$$n = 3k + 3$$$ to obtain the proof for all $$$n$$$.

An issue with this proof is that it does not account for amortisation. For the tree that was shown above, even though it is not possible to reduce the size of the tree by less than $$$\frac{2n}{3}$$$ in the first query, the next query will be done on a single chain of length $$$2k + 1$$$, which makes it easy for the remaining queries to reduce the size of the chain by half each time. Hence, this proof is not quite complete yet and I would appreciate it if anyone could tell me if they come up with a new idea.

Example

Link to problem

Abridged statement

You are given a tree consisting of $$$n$$$ nodes and there are two distinct nodes $$$a$$$ and $$$b$$$ that you have to determine. You can ask up to 40 queries of the form: $$$k\ \ x_1\ \ x_2\ \ldots\ x_k$$$ and the grader will return the distance between the lowest common ancestor of $$$x_1, x_2, \ldots, x_k$$$ if the tree was rooted at $$$a$$$ and the lowest common ancestor of $$$x_1, x_2, \ldots, x_k$$$ if the tree was rooted at $$$b$$$.

Solution

If we query all the vertices, we will be able to determine the distance between $$$a$$$ and $$$b$$$. Let this distance be $$$d$$$. Suppose we query an extended subtree together with its root. If the return value is $$$0$$$, it means that both $$$a$$$ and $$$b$$$ are outside of the extended subtree. If the return value is $$$d$$$, it means that both $$$a$$$ and $$$b$$$ are in the extended subtree or the root. Lastly, if the return value is between $$$1$$$ and $$$d - 1$$$, we know the distance that $$$a$$$ and $$$b$$$ are away from the root, so we can do two separate binary searches to find each of them.

Using the above-mentioned algorithm to choose the optimal extended subtree, we can solve this problem using $$$1 + \log_{1.5} n + 2\log_2 n \le 40$$$ queries. Note that there is a solution that uses only $$$2\log_2 n$$$ queries in the editorial of the problem, however, I feel like my solution is easier to come up with as it was the solution that I came up with during the contest.

Code

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

const int MAXN = 1005;

int n;
vector<int> adj[MAXN];
bool bad[MAXN];
vector<int> ans;

int ask(vector<int> qv) {
    cout << "? " << qv.size();
    for (int i : qv) {
        cout << ' ' << i;
    }
    cout << endl;
    int d; cin >> d;
    return d;
}
void answer(int x, int y) {
    cout << "! " << x << ' ' << y << endl;
    exit(0);
}

int sub[MAXN];
void dfs(int u, int p) {
    sub[u] = 1;
    for (int v : adj[u]) {
        if (v == p || bad[v]) {
            continue;
        }
        dfs(v, u);
        sub[u] += sub[v];
    }
}
int getc(int u, int p, int s) {
    for (int v : adj[u]) {
        if (v == p || bad[v]) {
            continue;
        }
        if (sub[v] > s / 2) {
            return getc(v, u, s);
        }
    }
    return u;
}
vector<int> stk;
void pushsub(int u, int p) {
    stk.push_back(u);
    for (int v : adj[u]) {
        if (v == p || bad[v]) {
            continue;
        }
        pushsub(v, u);
    }
}
void getl(int u, int p, int tl) {
    if (tl == 0) {
        stk.push_back(u);
    }
    for (int v : adj[u]) {
        if (v == p || bad[v]) {
            continue;
        }
        getl(v, u, tl - 1);
    }
}
int search(int r, vector<int> stk, int d) {
    int lo = 0, hi = stk.size() - 1;
    while (lo < hi) {
        int mid = lo + hi >> 1;
        vector<int> qv;
        qv.push_back(r);
        for (int i = lo; i <= mid; i++) {
            qv.push_back(stk[i]);
        }
        if (ask(qv) == d) {
            hi = mid;
        } else {
            lo = mid + 1;
        }
    }
    return stk[hi];
}

int main() {
    cin >> n;
    for (int i = 1; i < n; i++) {
        int u, v; cin >> u >> v;
        adj[u].push_back(v);
        adj[v].push_back(u);
    }
    vector<int> qv(n);
    iota(qv.begin(), qv.end(), 1);
    int d = ask(qv);
    int r = 1;
    while (1) {
        dfs(r, -1);
        if (sub[r] == 2) {
            assert(d == 1);
            for (int v : adj[r]) {
                if (!bad[v]) {
                    answer(r, v);
                    return 0;
                }
            }
            assert(0);
        }
        assert(sub[r] > 2);
        r = getc(r, -1, sub[r]);
        dfs(r, -1);
        int sm = 1;
        pair<int, int> rg = {-1, -1};
        for (int i = 0; i < adj[r].size(); i++) {
            int v = adj[r][i];
            if (bad[v]) {
                continue;
            }
            if (sub[v] >= sub[r] / 3) {
                rg = {i, i};
                break;
            }
            sm += sub[v];
            if (sm >= sub[r] / 3) {
                rg = {0, i};
                break;
            }
        }
        auto [s, e] = rg;
        stk.clear();
        stk.push_back(r);
        for (int i = s; i <= e; i++) {
            int v = adj[r][i];
            if (bad[v]) {
                continue;
            }
            pushsub(v, r);
        }
        int td = ask(stk);
        if (td == d) {
            for (int i = 0; i < adj[r].size(); i++) {
                if (s <= i && i <= e) {
                    continue;
                }
                int v = adj[r][i];
                bad[v] = 1;
            }
        } else if (td == 0) {
            for (int i = s; i <= e; i++) {
                int v = adj[r][i];
                bad[v] = 1;
            }
        } else {
            stk.clear();
            for (int i = s; i <= e; i++) {
                int v = adj[r][i];
                if (bad[v]) {
                    continue;
                }
                getl(v, r, td - 1);
            }
            int x = search(r, stk, td);
            stk.clear();
            for (int i = 0; i < adj[r].size(); i++) {
                if (s <= i && i <= e) {
                    continue;
                }
                int v = adj[r][i];
                if (bad[v]) {
                    continue;
                }
                getl(v, r, d - td - 1);
            }
            int y = search(r, stk, d - td);
            answer(x, y);
            return 0;
        }
    }
    return 0;
}

Conclusion

I hope this will be interesting and useful to some people 🥺. Please let me know in the comments if you can come up with a better proof that handles the amortisation or maybe even come up with a more optimal algorithm that can perform better than $$$\log_{1.5} n$$$ queries by making use of the amortisation. If you are aware of other problems that can be solved using this technique, feel free to let me know in the comments as well. Thank you for your support 😁.

Comments (6)

Write comment?

Vladithur

15 months ago, # |

← Rev. 3 →

+11

Nice blog, but I think there's a small mistake in this proof:

We can obtain the following inequalities: $$$x \ge \frac{n}{3}$$$, $$$x - p < \frac{n}{3}$$$, and $$$p < \frac{n}{3}$$$. Combining the second and third inequality, we obtain $$$x < \frac{2n}{3}$$$. Hence, $$$\frac{n}{3} \le x \le \frac{2n}{3}−1$$$ and the desired extended subtree is found.

But $$$\frac{2n}{3}$$$ is not an integer, so $$$x < \frac{2n}{3}$$$ does not imply $$$x \le \frac{2n}{3}−1$$$. In fact, it's impossible to get $$$\frac{n}{3} \le x \le \frac{2n}{3}−1$$$ on the tree you give just below.

My attempt at fixing the proof

→ Reply

Sacharlemagne

The first question is the problem "search" from this EJOI isn't it?
What I did was everytime i "eliminate" one side of the tree if the root has two children, and then i find the centroid in the tree after i "cut off" one of the sides, then i query on that centroid
If its in it, you have reduced by 1/2 because you will eliminate one of the sides of its children and you will be left with one children which has at most n/2 nodes, if its not you have reduced by 1/2 because the rest of the tree is at most n/2.
This gives i think 2log2(n) queries which passed to 100 for me (I participated in the mirror)

15 months ago, # ^ |

Also i know this might be the same solution its just that the point of view here is different

FatihCihan

he says the it is sufficient to do log2(n) operations in some problems , and that ejoi problem is one of them. His technique is good for different problems

lrvideckis

9 months ago, # |

← Rev. 2 →

An issue with this proof is that it does not account for amortisation. For the tree that was shown above, even though it is not possible to reduce the size of the tree by less than 2n3 in the first query, the next query will be done on a single chain of length 2k+1, which makes it easy for the remaining queries to reduce the size of the chain by half each time. Hence, this proof is not quite complete yet and I would appreciate it if anyone could tell me if they come up with a new idea.

I'm thinking of a case like this (perfect ternary tree), where it should have 1/3, 2/3 splits all the way down. But there is one doubt, which is certain decompositions have some edges upwards, which means it's not perfectly 1/3,2/3

9 months ago, # ^ |

So I tested this case (by doing knapsack-like DP on each node to get best split; then recursing on best split, and finding max depth of edge-centroid-decomp-tree 255249875)

so a perfect ternary tree with max depth l has 1+3+3^2+...+3^l = (3^(l+1)-1)/2 nodes. and the above code gave:

n: 1 1
n: 4 3
n: 13 5
n: 40 7
n: 121 9
n: 364 11
n: 1093 13
n: 3280 15
n: 9841 17
n: 29524 19

So assuming this pattern is correct: (big assumption)

a tree with n=(3^(l+1)-1)/2 nodes has max edge-centroid-decomp-tree depth = 2*l+1, whereas log1.5((3^(l+1)-1)/2) has an oblique asymptote of l*log1.5(3) + 1

The oblique asymptote of n=l*log1.5(3) + 1 checks out because we can solve for l, then plug it into (3^(l+1)-1)/2: (3^((n-1)/log1.5(3)+1)-1)/2, and it equals log1.5(n) exactly

Doing a similar thing but instead with 2*l+1, and noting 2 = log(3) / log(sqrt(3)), we get that max edge-centroid-decomp-tree depth in this case = log(n)/log(sqrt(3)) + log(2)/log(sqrt(3)) = log(2n/sqrt(3)) / log(sqrt(3))

basically showing amortized complexity is at least as bad as log with base = sqrt(3)

maomao90's blog

Definitions

Problem Structure

Solution

Optimality

Example

Abridged statement

Solution

Conclusion