首页 > 代码库 > [BZOJ4317]Atm的树
[BZOJ4317]Atm的树
[BZOJ4317]Atm的树
试题描述
Atm有一段时间在虐qtree的题目,于是,他满脑子都是tree,tree,tree……
于是,一天晚上他梦到自己被关在了一个有根树中,每条路径都有边权,一个神秘的声音告诉他,每个点到其他的点有一个距离(什么是距离不用说吧),他需要对于每个点回答:从这个点出发的第k小距离是多少;
如果atm不能回答出来,那么明天4019的闹钟将不会响,4019全寝可能就迟到了,所以atm希望你帮帮他。
输入
第一行,两个正整数n,k,表示树的点数,询问的是第几小距离;
第二~n行,每行三个正整数x,y,w,表示x和y之间有一条边,x为父亲,边权为w;
输出
n行, 每行一个数,第i行输出从i开始第k小距离
输入示例
5 2 1 5 2 1 2 4 2 3 6 2 4 5
输出示例
4 5 10 9 6
数据规模及约定
100% n<=15000, 边权在1~10之间,为了方便,保证1为根;K<=5000
题解
参见这道题。
#include <iostream> #include <cstdio> #include <algorithm> #include <cmath> #include <stack> #include <vector> #include <queue> #include <cstring> #include <string> #include <map> #include <set> using namespace std; const int BufferSize = 1 << 16; char buffer[BufferSize], *Head, *Tail; inline char Getchar() { if(Head == Tail) { int l = fread(buffer, 1, BufferSize, stdin); Tail = (Head = buffer) + l; } return *Head++; } int read() { int x = 0, f = 1; char c = Getchar(); while(!isdigit(c)){ if(c == ‘-‘) f = -1; c = Getchar(); } while(isdigit(c)){ x = x * 10 + c - ‘0‘; c = Getchar(); } return x * f; } #define maxn 15010 #define maxm 30010 #define maxlog 15 int n, m, head[maxn], nxt[maxm], to[maxm], dist[maxm]; void AddEdge(int a, int b, int c) { to[++m] = b; dist[m] = c; nxt[m] = head[a]; head[a] = m; swap(a, b); to[++m] = b; dist[m] = c; nxt[m] = head[a]; head[a] = m; return ; } int mnd[maxlog][maxn<<1], Log[maxn<<1], clo, dep[maxn], dfn[maxn]; void build(int u, int pa) { mnd[0][dfn[u] = ++clo] = dep[u]; for(int e = head[u]; e; e = nxt[e]) if(to[e] != pa) dep[to[e]] = dep[u] + dist[e], build(to[e], u), mnd[0][++clo] = dep[u]; return ; } void rmq_init() { Log[1] = 0; for(int i = 2; i <= clo; i++) Log[i] = Log[i>>1] + 1; for(int j = 1; (1 << j) <= clo; j++) for(int i = 1; i + (1 << j) - 1 <= clo; i++) mnd[j][i] = min(mnd[j-1][i], mnd[j-1][i+(1<<j-1)]); return ; } int cdist(int a, int b) { int l = dfn[a], r = dfn[b]; if(l > r) swap(l, r); int t = Log[r-l+1]; return dep[a] + dep[b] - (min(mnd[t][l], mnd[t][r-(1<<t)+1]) << 1); } #define maxnode 1800010 struct Node { int v, r, siz; Node() {} Node(int _, int __): v(_), r(__) {} }; struct Treap { Node ns[maxnode]; int ToT, fa[maxnode], ch[maxnode][2]; void maintain(int o) { if(!o) return ; ns[o].siz = 1; for(int i = 0; i < 2; i++) if(ch[o][i]) ns[o].siz += ns[ch[o][i]].siz; return ; } void rotate(int u) { int y = fa[u], z = fa[y], l = 0, r = 1; if(z) ch[z][ch[z][1]==y] = u; if(ch[y][1] == u) swap(l, r); fa[u] = z; fa[y] = u; fa[ch[u][r]] = y; ch[y][l] = ch[u][r]; ch[u][r] = y; maintain(y); maintain(u); return ; } void Insert(int& o, int v) { if(!o) { ns[o = ++ToT] = Node(v, rand()); return maintain(o); } bool d = v > ns[o].v; Insert(ch[o][d], v); fa[ch[o][d]] = o; if(ns[ch[o][d]].r > ns[o].r) { int t = ch[o][d]; rotate(t); o = t; } return maintain(o); } int query(int o, int v) { if(!o) return 0; int ls = ch[o][0] ? ns[ch[o][0]].siz : 0; if(v < ns[o].v) return query(ch[o][0], v); return ls + 1 + query(ch[o][1], v); } } sol; int f[maxn], siz[maxn], rt, size; bool vis[maxn]; void getroot(int u, int pa) { siz[u] = 1; f[u] = 0; for(int e = head[u]; e; e = nxt[e]) if(to[e] != pa && !vis[to[e]]) { getroot(to[e], u); siz[u] += siz[to[e]]; f[u] = max(f[u], siz[to[e]]); } f[u] = max(f[u], size - siz[u]); if(f[rt] > f[u]) rt = u; return ; } void dfs(int u, int pa) { siz[u] = 1; for(int e = head[u]; e; e = nxt[e]) if(to[e] != pa && !vis[to[e]]) dfs(to[e], u), siz[u] += siz[to[e]]; return ; } int fa[maxn]; void solve(int u) { vis[u] = 1; for(int e = head[u]; e; e = nxt[e]) if(!vis[to[e]]) { dfs(to[e], u); f[rt = 0] = size = siz[to[e]]; getroot(to[e], u); fa[rt] = u; solve(rt); } return ; } int Rt[maxn], Rtfa[maxn]; void update(int s) { sol.Insert(Rt[s], 0); for(int u = s; fa[u]; u = fa[u]) { int d = cdist(s, fa[u]); // printf("insert:: %d(fa: %d): %d\n", u, fa[u], d); sol.Insert(Rtfa[u], d); sol.Insert(Rt[fa[u]], d); } return ; } int query(int s, int x) { int ans = sol.query(Rt[s], x); // printf("base: %d\n", ans); for(int u = s; fa[u]; u = fa[u]) { int d = cdist(s, fa[u]); // printf("query:: %d(fa: %d): %d (%d - %d)\n", u, fa[u], x - d, sol.query(Rt[fa[u]], x - d), sol.query(Rtfa[u], x - d)); ans += sol.query(Rt[fa[u]], x - d) - sol.query(Rtfa[u], x - d); } return ans; } int main() { n = read(); int K = read(), sum = 0; for(int i = 1; i < n; i++) { int a = read(), b = read(), c = read(); AddEdge(a, b, c); sum += c; } build(1, 0); rmq_init(); f[rt = 0] = size = n; getroot(1, 0); solve(rt); // for(int i = 1; i <= n; i++) printf("(%d)%d%c", i, fa[i], i < n ? ‘ ‘ : ‘\n‘); for(int i = 1; i <= n; i++) update(i); for(int i = 1; i <= n; i++) { int l = 1, r = sum; while(l < r) { int mid = l + r >> 1; if(query(i, mid) - 1 < K) l = mid + 1; else r = mid; } printf("%d\n", l); } return 0; }
又 1A 啦 2333333!
[BZOJ4317]Atm的树
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。