[CF468D] Tree

Little X has a tree consisting of n nodes (they are numbered from 1 to n). Each edge of the tree has a positive length. Let's define the distance between two nodes v and u (we'll denote it d(v, u)) as the sum of the lengths of edges in the shortest path between v and u.

A permutation p is a sequence of n distinct integers p1, p2, ..., pn(1 ≤ pi ≤ n). Little X wants to find a permutation p such that sum is maximal possible. If there are multiple optimal permutations, he wants to find the lexicographically smallest one. Help him with the task!

Input

The first line contains an integer n (1 ≤ n ≤ 105).

Each of the next n - 1 lines contains three space separated integers ui,  vi, wi (1 ≤  ui,  vi ≤  n; 1 ≤  wi ≤  105), denoting an edge between nodes ui and vi with length equal to wi.

It is guaranteed that these edges form a tree.

Output

In the first line print the maximum possible value of the described sum. In the second line print n integers, representing the lexicographically smallest permutation.

Examples
Input
Copy
2
1 2 3
Output
Copy
6
2 1
Input
Copy
5
1 2 2
1 3 3
2 4 4
2 5 5
Output
Copy
32
2 1 4 5 3




写了一上午,我太辣鸡了。
dis(i,pi)=(depi+deppi2×deplca(i,pi))=2×depi2×deplca(i,pi)
前面一定是不变的, 后面我们可以把它降低到0, 我们选择树的重心作为根,可以把后边变成0.
第一问解决。把一个点看成两个, 一个是入点,一个是出点.
记为in[x], out[x];要满足子树x可以匹配到答案, 我们必须让 $\large in[x] + out[x] <= \sum in[i]$

$in[x]$
$ f[i, j] = $
$\LARGE f[i, j]=\frac{\sum_{to}^{ } f[nxt[nxt[i,j],j],to] + f[nxt[nxt[i,j],j],j]}{deg[j]+1} + 1$
latex炸了不写了




#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <set>
using namespace std;
inline int read() {
    int res = 0;char ch=getchar();
    while(!isdigit(ch)) ch=getchar();
    while(isdigit(ch)) res=(res<<3)+(res<<1)+(ch^48), ch=getchar();
    return res;
}
#define reg register
#define N 100005
#define ll long long

int n;
struct edge {
    int nxt, to, val;
}ed[N*2];int head[N], cnt;
inline void add(int x, int y, int z) {
    ed[++cnt] = (edge) {head[x], y, z};
    head[x] = cnt;
}
ll ans;
int siz[N], root, mrt = 1e9;

void dfs(int x, int fa) 
{
    siz[x] = 1;
    int tmp = 0;
    for (reg int i = head[x] ; i ; i = ed[i].nxt)
    {
        int to = ed[i].to;
        if (to == fa) continue;
        dfs(to, x);
        siz[x] += siz[to];
        tmp = max(tmp, siz[to]);
    }
    tmp = max(tmp, n - siz[x]);
    if (tmp < mrt) mrt = tmp, root = x;
}

void gfs(int x, int fa, ll d) {
    ans += d;
    for (reg int i = head[x] ; i ; i = ed[i].nxt) 
    {
        int to = ed[i].to;
        if (to == fa) continue;
        gfs(to, x, d + ed[i].val);
    }
}

set <int> all;
set <int> tr[N];
set < pair <int, int> > inout;
#define makp(a, b) make_pair(a, b)

int belong[N], sum[N];
void ffs(int x, int fa, int bel) {
    belong[x] = bel;
    tr[bel].insert(x);
    siz[x] = 1;
    for (reg int i = head[x] ; i ; i = ed[i].nxt) 
    {
        int to = ed[i].to;
        if (to == fa) continue;
        ffs(to, x, bel);
        siz[x] += siz[to];
    }
}

inline int Find(int x) 
{
    set < pair <int, int> > :: iterator it = inout.lower_bound(makp(n - x + 1, 0));
    bool must = 0;
    if (it != inout.end() and (*it).first == n - x + 1 and (*it).second != belong[x])  must = 1;
    if (belong[x]) {
        inout.erase(makp(sum[belong[x]], belong[x]));
        sum[belong[x]]--;
        inout.insert(makp(sum[belong[x]], belong[x]));
    }
    if (must) {
        int k = *tr[(*it).second].begin();
        all.erase(k);
        tr[(*it).second].erase(k);
        if (tr[(*it).second].size()) all.insert(*tr[(*it).second].begin());
        inout.erase(makp(sum[belong[k]], belong[k]));
        inout.insert(makp(--sum[belong[k]], belong[k]));
        return k;        
    } else {
        set <int> :: iterator itt = all.begin();
        if(belong[*itt] == belong[x] and x != root) itt++;
        all.erase(itt);
        int k = *itt;
        if (belong[k]) {
            tr[belong[k]].erase(k);
            if (tr[belong[k]].size()) all.insert(*tr[belong[k]].begin());
            inout.erase(makp(sum[belong[k]], belong[k]));
            inout.insert(makp(--sum[belong[k]], belong[k]));        
        }
        return k;
    }
    return 0;
}

int main()
{
    n = read();
    for (reg int i = 1 ; i < n ; i ++) {
        int x = read(), y = read(), z = read();
        add(x, y, z), add(y, x, z);
    }
    
    root = 1;
    dfs(1, 0);
    ans = 0;
    gfs(root, 0, 0);
    printf("%lld\n", ans * 2);
    all.insert(root);
    
    for (reg int i = head[root] ; i ; i = ed[i].nxt) 
    {
        int to = ed[i].to;
        ffs(to, root, to);
        sum[to] = siz[to] * 2;
        inout.insert(makp(sum[to], to));
        siz[root] += siz[to];
        all.insert(*tr[to].begin());
    }
    
    for (reg int i = 1 ; i <= n ; i ++)
        printf("%d ", Find(i));
    return 0;
}
 
    

 

 

转载于:https://www.cnblogs.com/BriMon/p/9547292.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值