「WC2018即时战略」
题目描述
小 M 在玩一个即时战略 (Real Time Strategy) 游戏。不同于大多数同类游戏，这个游戏的地图是树形的。也就是说，地图可以用一个由 (n) 个结点，(n - 1) 条边构成的连通图来表示。这些结点被编号为 (1 sim n)。
每个结点有两种可能的状态：“已知的”或“未知的”。游戏开始时，只有 (1) 号结点是已知的。
在游戏的过程中，小 M 可以尝试探索更多的结点。具体来说，小 M 每次操作时需要选择一个已知的结点 (x)，和一个不同于 (x) 的任意结点 (y)（结点 (y) 可以是未知的）。然后游戏的自动寻路系统会给出 (x) 到 (y) 的最短路径上的第二个结点 (z)，也就是从 (x) 走到 (y) 的最短路径上与 (x) 相邻的结点。此时，如果结点 (z) 是未知的，小 M 会将它标记为已知的。这个游戏的目标是：利用至多 (T) 次探索操作，让所有结点的状态都成为已知的。然而小 M 还是这个游戏的新手，她希望得到你的帮助。
解题思路 :
首先有一个比较直观的暴力，random_shuffle一个询问顺序，同时维护一棵“已知树”。
每次从根节点开始询问，回答要么是当前点的儿子，要么是一个未知节点，如果是当前点的儿子就进入儿子节点，否则就把未知节点添加进树。
这样子做复杂度和询问次数都是 (O(n^2))，加上一条链的暴力可以得到 (65) 分。
实际上每次如果是已知节点的话，只需要进入儿子对应的子树询问即可，所以很容易想到用点分树维护这个“已知树”，每次直接找到这个儿子对应的点分中心进行询问，树高变成 (logn) 。
但是加点操作会破坏点分树的性质，使得树高会大于 (logn) 以至于退化到平方级别的复杂度，在这里可以用替罪羊树的思想，每次加完点后暴力向上检查子树的平衡性暴力重构，(alpha) 这里一般设 (0.7) 。
不过由于我维护点分树信息的时候用 (map) 存了每个儿子对应的点分中心是什么，所以我的复杂度是 $O(nlog^2n) $ ，有点卡常数，不保证所有地方都能过。
#include "rts.h"
/*program by mangoyang*/
#include<bits/stdc++.h>
#define inf (0x7f7f7f7f)
#define Max(a, b) ((a) > (b) ? (a) : (b))
#define Min(a, b) ((a) < (b) ? (a) : (b))
typedef long long ll;
using namespace std;
const int N = 1000005;
map<int, int> mp[N], ss[N];
int vis[N], vi[N], sz[N], in[N], cc[N], mx[N], fa[N];
int size[N], id[N], ps[N], all, mn, rt, Root = 1;
map<int, int>::iterator it;
namespace Line{
int s[2] = {1, 1};
inline void solve(int n){
for(int i = 1, x; i < n; i++){
int now = id[i], tw = 1;
while(!vis[now]){
x = explore(s[tw], now);
if(vis[x]) x = explore(s[tw^=1], now);
vis[x] = 1, s[tw] = x;
}
}
}
}
inline void cleartag(int u){
in[u] = 1;
for(map<int,int>::iterator it = mp[u].begin(); it != mp[u].end(); it++)
if(it->second) cleartag(it->second);
}
inline void getsize(int u, int ff){
int now = 0; sz[u] = 1;
for(map<int, int>::iterator it = mp[u].begin(); it != mp[u].end(); it++){
int v = it->first;
if(v == ff || !in[v]) continue;
getsize(v, u), sz[u] += sz[v];
if(sz[v] >= now) now = sz[v];
}
now = max(now, all - sz[u]);
if(now <= mn) mn = now, rt = u;
}
inline void rebuild(int u){
int last = all; in[u] = 0;
for(map<int, int>::iterator it = mp[u].begin(); it != mp[u].end(); it++){
int v = it->first;
if(!in[v]) continue;
mn = all = sz[v] >= sz[u] ? last - sz[u] : sz[v];
getsize(v, u);
size[rt] = all, fa[rt] = u, mx[u] = Max(mx[u], all);
mp[u][v] = rt, ss[u][rt] = v, mp[v][u] = 0, rebuild(rt);
}
}
inline void update(int u){
int ned = 0;
for(int x = u; x != Root; x = fa[x]){
size[fa[x]]++;
if(size[x] > mx[fa[x]]) mx[fa[x]] = size[x];
if(mx[fa[x]] >= size[fa[x]] * 0.735) ned = fa[x];
}
if(!ned) return;
if(ned){
cleartag(ned);
all = mn = size[ned], getsize(ned, fa[ned]);
if(ned == Root) Root = rt;
size[rt] = all, fa[rt] = fa[ned];
int k = ss[fa[rt]][ned];
ss[fa[rt]][rt] = k, mp[fa[rt]][k] = rt, rebuild(rt);
}
}
inline void addnode(int pos){
for(register int u = Root; ; ){
int x = explore(u, pos);
if(!vis[x]){
fa[x] = u, size[x] = vis[x] = 1;
mp[u][x] = ss[u][x] = x, mp[x][u] = 0, update(x);
break;
} u = mp[u][x];
}
}
void play(int n, int T, int datatype){
srand(19262333);
for(int i = 1; i < n; i++) id[i] = i + 1;
random_shuffle(id + 1, id + n), vis[1] = 1;
if(datatype == 3) return (void) (Line::solve(n));
for(int i = 1; i < n; i++) while(!vis[id[i]]) addnode(id[i]);
}