题目背景
spj来源：loj-Robin。已获得授权。
https://www.luogu.org/paste/dxytr6gc
附spj，一些修改部分未按照代码规范，请各位谅解。
题目描述
魔法国度 Magic Land 里最近出现了一个大盗 Frank，他在 Magic Land 四处作案，专门窃取政府机关的机密文件（因而有人怀疑 Frank 是敌国派来的间谍）。
为了捉住 Frank，Magic Land 的安全局重拳出击！
Magic Land 由 N 个城市组成，并且这 N 个城市又由恰好 N-1 条公路彼此连接起来，使得任意两个城市间都可以通过若干条公路互达。从数据结构的角度我们也可以说，这 N 个城市和 N-1 条公路形成了一棵树。
例如，下图就是 Magic Land 的一个可能格局（4 个城市用数字编号，3 条公路用字母编号）：
大盗 Frank 能够在公路上以任意速度移动。
比方说，对于上图给出的格局，在 0.00001 秒钟内（或者任意短的一段时间内），Frank 就可以从城市 1 经过城市 2 到达城市 4，中间经过了两条公路。
想要生擒 Frank 困难重重，所以安全局派出了经验丰富的警探，这些警探具有非凡的追捕才能：
-
只要有警探和 Frank 同处一个城市，那么就能够立刻察觉到Frank，并且将其逮捕。
-
虽然 Frank 可以在公路上以任意快的速度移动，但是如果有警探和 Frank 在同一条公路上相遇，那么警探也可以立刻察觉到 Frank 并将其逮捕。
安全局完全不知道 Frank 躲在哪个城市，或者正在哪条公路上移动，所以需要制定一个周密的抓捕计划，计划由若干 步骤组成。在每一步中，可以做如下几件事中的一个：
-
在某个城市空降一位警探。警探可以直接从指挥部空降到 Magic Land 的任意一个城市里。此操作记为“L x”，表示在编号为 x 的城市里空降一位警探。耗时 1 秒。
-
把留在某个城市里的一位警探直接召回指挥部。以备在以后的步骤中再度空降到某个城市里。此操作记为“B x”。表示把编号为 x 的城市里的一位警探召回指挥部。耗时 1 秒。
-
让待在城市 x 的一位警探沿着公路移动到城市 y，此操作记为“M x y”。耗时 1 秒。当然，前提是城市 x 和城市 y 之间有公路。如果在警探移动的过程中，大盗 Frank 也在同一条公路上，那么警探就抓捕到了Frank。
现在，由你来制定一套追捕计划，也就是给出若干个步骤，需要保证：无论大盗 Frank 一开始躲在哪儿，也无论 Frank 在整个过程中如何狡猾地移动（Frank大盗可能会窃取到追捕行动的计划书，所以他一定会想尽办法逃避），他一定会被缉拿归案。
希望参与的警探越少越好，因为经验丰富的警探毕竟不多。
例如对于前面所给的那个图示格局，一个可行的计划如下：
-
L 2 在城市 2 空降一位警探。注意这一步完成之后，城市 2 里不会有 Frank，否则他将被捉住。
-
L 2 再在城市 2 空降一位警探。
-
M 2 1 让城市 2 的一位警探移动到城市 1。注意城市 2 里还留有另一位警探。这一步完成之后，城市 1 里不会有 Frank，公路 A 上也不会有 Frank。也就是说，假如 Frank 还没有被逮捕，那么他只能是在城市 3 或城市 4 里，或者公路 B 或公路 C 上。
-
B 1 召回城市 1 的一位警探。注意虽然召回了这位警探，但是由于我们始终留了一位警探在城市 2 把守，所以 Frank 仍然不可能跑到城市 1 或者是公路 A 上。
-
L 3 在城市 3 空降一位警探。注意这一步可以空降在此之前被召回的那位警探。这一步完成之后，城市 3 里不会有 Frank，否则他会被捉住。
-
M 3 2 让城市 3 里的一位警探移动到城市 2。这一步完成之后，如果 Frank 还没有被捉住，那他只能是在公路 C 上或者城市 4 里。注意这一步之后，城市 2 里有两位警探。
-
M 2 4 让城市 2 里的一位警探移动到城市 4。这一步完成之后，Frank 一定会被捉住，除非他根本就没来 Magic Land。
这个计划总共需要 2 位警探的参与。可以证明：如果自始至终只有 1 名或者更少的警探参与，则 Frank 就会逍遥法外。
你的任务很简单：对于一个输入的 Magic Land 的格局，计算 S，也就是为了追捕 Frank 至少需要投入多少位警探，并且给出相应的追捕计划步骤。
输入格式
输入文件给出了 Magic Land 的格局。
第一行一个整数 N，代表有 N 个城市，城市的编号是 1~N。
接下来 N-1 行，每行有两个用空格分开的整数 x i ，y i ，代表城市 x i ，y i之间有公路相连。保证 1≤x i ,y i ≤N
输出格式
向输出文件输出你所给出的追捕计划。
第一行请输出一个整数 S，代表追捕计划需要多少位警探。
第二行请输出一个整数 T，代表追捕计划总共有多少步。
接下来请输出 T 行，依次描述了追捕计划的每一步。每行必须是以下三种形式之一：
”L x”，其中 L 是大写字母，接着是一个空格，再接着是整数 x，代表在城市 x 空降一位警探。你必须保证 1≤x≤N。
“B x”，其中 B 是大写字母，接着是一个空格，再接着是整数 x，代表召回城市 x 的一位警探。你必须保证 1≤x≤N，且你的计划执行到这一步之前，城市 x 里面确实至少有一位警探。
“M x y”，其中 M 是大写字母，接着是一个空格，再接着是整数 x，再跟一个空格，最后一个是整数 y。代表让城市 x 的一位警探沿着公路移动到城市 y。你必须保证 1≤x, y≤N，且你的计划执行到这一步之前，城市 x 里面确实至少有一位警探，且城市 x, y 之前确实有公路。
必须保证输出的 S 确实等于追捕计划中所需要的警探数目。
输入输出样例
输入 #1复制
4 1 2 3 2 2 4
输出 #1复制
2 7 L 2 L 2 M 2 1 B 1 L 3 M 3 2 M 2 4
说明/提示
对于任何一个测试点：
如果输出的追捕计划不合法，或者整个追捕计划的步骤数 T 超过了 20000，或者追捕计划结束之后，不能保证捉住 Frank，则不能得分。
否则，用你输出的 S 和我们已知的标准答案 S * 相比较：
1. 若 S<S * ，则得到 120%的分。
2. 若 S=S * ，则得到 100%的分。
3. 若 S * <S≤S * +2，则得到 60%的分。
4. 若 S * +2<S≤S * +4，则得到 40%的分。
5. 若 S * +4<S≤S * +8，则得到 20%的分。
6. 若 S>S * +8，则得到 10%的分。
输入保证描述了一棵连通的 N 结点树，1≤N≤1 000。
下面设 ��→�Tu→v 表示树 �T 以 �u 为根时 �v 的子树内所有点的生成子树，�(�)S(T) 表示树 �T 的答案，一条路径的“旁枝”为把路径上的点和边去掉后留下的子树的集合（这里和路径上的点相连的边仍然保留）。
一定存在一个最优方案，这其中能找出一个“警长”，它走过的是一条简单路径，称为“大道”，设“大道”内“警长”走到的第 �i 个节点为 ��vi，且“警长”在 ��vi 时其他警察会把和 ��vi 相连的所有旁枝查完，答案即为所有“大道”的所有“旁枝”的答案最大值加一。
有一个比较显然的性质是 �⊆�A⊆B 时 �(�)≤�(�)S(A)≤S(B)。
还有一个重要性质是若存在一个度数大于等于三的点 �u 和三条边 (�,�),(�,�),(�,�)(u,a),(u,b),(u,c)，则 �(�)>min(�(��→�),�(��→�),�(��→�))S(T)>min(S(Tu→a),S(Tu→b),S(Tu→c))，如果因为假设 �(��→�)=�(��→�)=�(��→�)=�(�)S(Tu→a)=S(Tu→b)=S(Tu→c)=S(T)，且先全部清查 ��→�Tu→a 后再查 ��→�Tu→b 最后再查 ��→�Tu→c，那么必然有一个时刻要把全部警力放到 ��→�Tu→b 里，那么此时 �u 无人值守，��→�Tu→c 又没查完，��→�Tu→a 就白查了。
如果把“大道”设为 �T 的重心，那么可以得到一个 log∣�∣log∣T∣ 级别的构造方案，所以 �(�)S(T) 最多是 log∣�∣log∣T∣ 级别的。
考虑找到这么样的一个“大道”使得答案最优，仍然设第 �i 个节点为 ��vi，点数为 �l，那么不难得到下面的充要条件：
∀1≤�≤�,(��,�)∈�,�≠��−1,�≠��+1,�(���→�)<�(�)∀1≤i≤l,(vi,s)∈T,s=vi−1,s=vi+1,S(Tvi→s)<S(T)
称使得 ∀(�,�)∈�,�(��→�)<�(�)∀(u,v)∈T,S(Tu→v)<S(T) 的 �v 为 �T 的 Hub，此时 Hub 单独构成一个“大道”。
设一个边集 �={(�,�)∣(�,�)∈�,�(��→�)=�(��→�)=�(�)}X={(u,v)∣(u,v)∈T,S(Tu→v)=S(Tv→u)=S(T)}，当 �T 内没有 Hub 时 �X 必然非空，因为对于一个非 Hub 点 �v 设 �(�)={(�,�)∣(�,�)∈�,�(��→�)=�(�)}f(v)={(u,v)∣(u,v)∈T,S(Tv→u)=S(T)}，那么没有 Hub 时根据抽屉原理必然有至少一组 �(�)=�(�)f(u)=f(v)，那么 (�,�)∈�(u,v)∈X。
且因为上面的重要性质不会有三条 �X 边连在一个点上。
且 �X 联通，否则存在一条 (�,�)∉�(x,y)∈/X 使得把 (�,�)(x,y) 两边均有 �X 内的边，那么必然存在 ��→�⊆��→�Ta→b⊆Tx→y 且 (�,�)∈�(a,b)∈X，那么 �(��→�)≤�(��→�)≤�(�)S(Ta→b)≤S(Tx→y)≤S(T)，则 �(��→�)=�(��→�)=�(�)S(Ta→b)=S(Tx→y)=S(T)，同理 �(��→�)=�(�)S(Ty→x)=S(T)，得出矛盾。
所以 �X 里的边构成了一条路径，会发现这条路径满足了上面的充要条件，构成了一条“大道”。
好像这个结论就是一个对这个做法第一步的证明？
对于一棵树 �T 和一个节点 �r，考虑把节点分成下面四种类型：
-
H：�r 是 �T 的Hub。 -
E：�r 是 �T 的“大道”的端点或者“大道”的端点对应的旁枝中的点。 -
I：�r 是 �T 的“大道”中的点但不是端点。 -
M：其他情况。
注意 Hub 也算“大道”的端点。
主要思想随便选一个根，然后树形 dp，dp 完子树后要支持把一条边挂到根上去、根从儿子换到父亲的换根操作和两个子树信息的合并操作，维护的信息要记录答案，“大道”端点，根的类型，以及根如果是 M 型的话根所在旁枝的递归信息，如果不是 H 那么这里的“大道”是上面 �X 构成的那条，“把一条边挂上去”这个操作可以视为合并操作完成。
换根时前两个信息不变，此时新根度数为一，分为这么一些情况：
-
若原根为
E或H则新根必然为E。 -
若原根为
I，答案为一时这棵子树就是条链，新根为E，否则新根可以直接设为M，递归信息要处理的子树内只有一条边容易直接计算。 -
若原根为
M，那么新根也为M，递归信息需要递归调用换根来完成，由于递归信息内答案严格递减所以递归次数的级别最多为 �(log�)O(logn)。
合并时根在两边的度数可能不止为一了，可以分为两边答案相不相等来讨论：
两边答案相等时：
-
若两边均为
H则根为H答案不变。 -
若一边为
H一边为I则根为I答案不变。 -
若一边为
H一边为E则根为E答案不变，“大道”的端点可以设为根和E中的那条“大道”中离根更远的那个端点，因为“警长”多走不会使这边的答案更劣，且走到根后就能把另一边给扫了，而走没走到根的话是会导致答案变大的。 -
若两边均为
E则根为I答案不变，“大道”的端点可以设为两边的的“大道”中离根更远的那个端点，理由一样。 -
若一边为
M则根为H答案加一，因为设M那边为 �1T1 另一边为 �2T2 合并后为 �T，��vi 在 �1T1 的“大道”上且根在 ��→�Ta→b 内，那么 �(��→�)≥�(�2)=�(�1)S(Ta→b)≥S(T2)=S(T1)，�(���→��−1)=�(���→��+1)=�(�1)S(Tvi→vi−1)=S(Tvi→vi+1)=S(T1)，根据上面的重要性质所以 �(�)>�(�1)S(T)>S(T1)，此时根设为Hub就能取到 �(�)=�(�1)+1S(T)=S(T1)+1 的下界。 -
其他情况（其实只有两边均为
I和一边I一边E）和上一种一样。
两边答案不相等时，不妨设答案大的那边为 �1T1 小的那边为 �2T2：
-
若 �1T1 中根为
H或I则答案为 �(�1)S(T1)，“大道”是 �1T1 的“大道”，根自然也是 �1T1 中的类型。 -
若 �1T1 中根为
E则情况类似两边相等时的第三种讨论，根为E答案为 �(�1)S(T1)，“大道”的端点可以设为根和 �1T1 中的那条“大道”中离根更远的那个端点。 -
若 �1T1 中根为
M，设根所在 �1T1 的大道旁枝和 �2T2 的并为 �3T3，则又要根据记录的递归信息来递归计算 �3T3 的信息，和上面一样递归次数的级别最多为 �(log�)O(logn)，然后要合并 �1T1 和 �3T3 又要讨论： -
若 �(�1)>�(�3)S(T1)>S(T3)，则根为
M答案为 �(�1)S(T1)，“大道”是 �1T1 的“大道”，递归信息为 �3T3 的信息。 -
若 �(�1)=�(�3)S(T1)=S(T3)，则情况类似两边相等时的第五种讨论，根为
H答案加一。 -
�(�1)<�(�3)S(T1)<S(T3) 是不可能的，因为把根设为“大道”后 �(�3)S(T3) 的答案为“根所在 �1T1 的大道旁枝”和 �2T2 这两个部分的较大值加一，这是不会大于 �(�1)S(T1) 的
至此我们终于完成了这个分讨。
求出答案后要构造方案，这个就直接算完求出“大道”以后递归求旁枝方案就行了，单次求答案是单 log 的，由于递归层数是 log 的所以时间复杂度双 log，比较难卡满。
因为同一条边不会走两次，所以构造方案的步数最多为 3�3n
代码如下:
#include<bits/stdc++.h>
#define ll long long
#define INF 2147483647
inline int inp(){
char c = getchar();
while(c < '0' || c > '9')
c = getchar();
int sum = 0;
while(c >= '0' && c <= '9'){
sum = sum * 10 + c - '0';
c = getchar();
}
return sum;
}
int head[100010];
int nxt[20010];
int end[20010];
char type[100000];
int num1[100000], num2[100000];
int cnt = 0;
int cou = 0;
int f[20010];
void link(int a, int b){
nxt[++cou] = head[a];
head[a] = cou;
end[cou] = b;
}
void dfs(int cur, int last){
int max = 0;
bool mt = true;
for(int x = head[cur]; x != -1; x = nxt[x]){
if(end[x] != last){
dfs(end[x], cur);
if(f[end[x]] > max){
max = f[end[x]];
mt = false;
} else if(f[end[x]] == max)
mt = true;
}
}
if(mt)
f[cur] = max + 1;
else
f[cur] = max;
}
void dfs2(int cur, int last){
int max = 0;
bool mt = true;
int degree = 0;
int pos;
for(int x = head[cur]; x != -1; x = nxt[x]){
if(end[x] != last){
dfs(end[x], cur);
if(f[end[x]] > max){
max = f[end[x]];
mt = false;
pos = end[x];
} else if(f[end[x]] == max)
mt = true;
degree++;
}
}
// printf("%d %d pos = %d\n", cur, last, pos);
if(degree == 0){
type[++cnt] = 'B';
num1[cnt] = cur;
return ;
}
for(int x = head[cur]; x != -1; x = nxt[x]){
if(end[x] != pos && end[x] != last){
type[++cnt] = 'L';
num1[cnt] = cur;
type[++cnt] = 'M';
num1[cnt] = cur;
num2[cnt] = end[x];
dfs2(end[x], cur);
}
}
type[++cnt] = 'M';
num1[cnt] = cur;
num2[cnt] = pos;
dfs2(pos, cur);
if(mt)
f[cur] = max + 1;
else
f[cur] = max;
}
int main(){
memset(head, -1, sizeof(head));
int n = inp();
for(int i = 1; i < n; i++){
int u = inp();
int v = inp();
link(u, v);
link(v, u);
}
int root = 0;
int min = INF;
for(int i = 1; i <= n; i++){
dfs(i, 0);
if(f[i] < min){
min = f[i];
root = i;
}
}
dfs2(root, 0);
printf("%d\n%d\n", min, cnt + 1);
printf("L %d\n", root);
for(int i = 1; i <= cnt; i++){
putchar(type[i]);
printf(" %d", num1[i]);
if(type[i] == 'M')
printf(" %d", num2[i]);
putchar('\n');
}
}拜拜!
749
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
