Submission #747893

# Submission time Handle Problem Language Result Execution time Memory
747893 2023-05-25T06:39:20 Z FulopMate Janjetina (COCI21_janjetina) C++17
15 / 110
1500 ms 11728 KB
#include <bits/stdc++.h>

// #define MULTI_TEST_CASE
// #define TEST_TEXT

using namespace std;

#define ll long long
#define MAX(a, b) (a) = max((a), (b))
#define MIN(a, b) (a) = min((a), (b))
#define all(a) (a).begin(), (a).end()
#define sortedpair(a, b) {min((a), (b)), max((a), (b))}

const ll MOD = 1e9+7;

struct edge {
    int to, w;
};

int n, k;
vector<vector<edge>> v;
ll ans = 0;
vector<bool> voltcentroid;
vector<int> siz;
const int st_size = 262144;
vector<ll> st_ut, st_ut_max;

void add(vector<ll> &st, int d, int dif = 1, int ind = 1, int l = 0, int r = st_size-1){
    if(l == r){
        st[ind]+=dif;
        return;
    }
    int mid = (l+r)/2;
    if(d <= mid) add(st, d, dif, ind*2, l, mid);
    else add(st, d, dif, ind*2+1, mid+1, r);
    st[ind] = st[ind*2] + st[ind*2+1];
}

ll query(vector<ll> &st, int ql, int qr, int ind = 1, int l = 0, int r = st_size-1){
    if(qr < ql)return 0;
    if(r < ql || l > qr)return 0;
    if(ql <= l && r <= qr) return st[ind];
    int mid = (l+r)/2;
    return query(st, ql, qr, ind*2, l, mid) + query(st, ql, qr, ind*2+1, mid+1, r);
}

int c_szamoz(int x, int p = -1){
    if(voltcentroid[x])return 0;
    int cur = 1;
    for(edge&i : v[x]){
        if(i.to == p)continue;
        cur += c_szamoz(i.to, x);
    }
    return siz[x] = cur;
}

int c_keres(int x, int meret, int p = -1){
    for(edge&i : v[x]){
        if(i.to == p)continue;
        if(!voltcentroid[i.to]){
            if(siz[i.to] > meret/2){
                return c_keres(i.to, meret, x);
            }
        }
    }
    return x;
}

struct c_ut {
    int w, x, d;
    bool operator<(const c_ut & other) const {
        return -w < -other.w;
    }
};

vector<int> inserts;

void reset_st(){
    while(inserts.size()){
        add(st_ut, inserts.back(), -1);
        inserts.pop_back();
    }
}

void megold(int x = 0){
    int meret = c_szamoz(x);
    int centroid = c_keres(x, meret);

    st_ut.assign(st_size*2, 0);
    // st_ut_max.assign(st_size*2, 0);

    add(st_ut, 0); inserts.push_back(0);

    vector<bool> seen(n, 0);
    seen[centroid] = 1;
    priority_queue<c_ut> q;
    for(edge&i : v[centroid]){
        if(!voltcentroid[i.to]){
            q.push({i.w, i.to, 1});
        }
    }
    while(!q.empty()){
        c_ut c = q.top();
        q.pop();
        if(seen[c.x])continue;
        seen[c.x] = 1;
        ans += query(st_ut, 0, c.w-c.d-k);
        add(st_ut, c.d); inserts.push_back(c.d);
        for(edge&i : v[c.x]){
            if(!seen[i.to] && !voltcentroid[i.to]){
                q.push({max(c.w, i.w), i.to, c.d+1});
            }
        }
    }

    reset_st();
    seen.assign(n, 0);
    for(edge&i : v[centroid]){
        if(!voltcentroid[i.to]){
            seen[centroid] = 1;
            priority_queue<c_ut> q;
            q.push({i.w, i.to, 1});
            while(!q.empty()){
                c_ut c = q.top();
                q.pop();
                if(seen[c.x])continue;
                seen[c.x] = 1;
                ans -= query(st_ut, 0, c.w-c.d-k);
                add(st_ut, c.d); inserts.push_back(c.d);
                for(edge&i : v[c.x]){
                    if(!seen[i.to] && !voltcentroid[i.to]){
                        q.push({max(c.w, i.w), i.to, c.d+1});
                    }
                }
            }
            reset_st();
        }
    }


    // add(st_ut, 0);

    // for(edge&i : v[centroid]){
    //     if(!voltcentroid[i.to]){
    //         c_szamol(i.to, i.w, 1, centroid);
    //         c_frissit(i.to, i.w, 1, centroid);
    //     }
    // }


    voltcentroid[centroid] = true;
    for(edge&i : v[centroid]){
        if(!voltcentroid[i.to]){
            megold(i.to);
        }
    }
}

void solve(){
    cin>>n>>k;
    v.assign(n, {});
    voltcentroid.assign(n, 0);
    siz.assign(n, 0);
    for(int i = 0; i < n-1; i++){
        int a, b, c; cin>>a>>b>>c; a--; b--;
        v[a].push_back({b, c});
        v[b].push_back({a, c});
    }
    megold();
    cout<<ans*2<<endl;
}

int main(){
    ios_base::sync_with_stdio(0); cin.tie(0);
    int _t = 1;
#ifdef MULTI_TEST_CASE
    cin >> _t;
#endif
    for(int _i = 0; _i < _t; _i++){
        #ifdef TEST_TEXT
        cout<<"Case #"<<_i+1<<": ";
        #endif
        solve();
    }
    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 2 ms 4436 KB Output is correct
2 Correct 3 ms 4416 KB Output is correct
3 Correct 12 ms 4448 KB Output is correct
4 Correct 134 ms 4564 KB Output is correct
5 Correct 105 ms 4564 KB Output is correct
6 Correct 121 ms 4564 KB Output is correct
7 Correct 126 ms 4544 KB Output is correct
8 Correct 121 ms 4552 KB Output is correct
9 Correct 115 ms 4436 KB Output is correct
10 Correct 102 ms 4436 KB Output is correct
11 Correct 104 ms 4500 KB Output is correct
12 Correct 113 ms 4436 KB Output is correct
13 Correct 136 ms 4492 KB Output is correct
14 Correct 144 ms 4492 KB Output is correct
15 Correct 118 ms 4492 KB Output is correct
16 Correct 109 ms 4500 KB Output is correct
17 Correct 111 ms 4436 KB Output is correct
18 Correct 117 ms 4496 KB Output is correct
19 Correct 108 ms 4496 KB Output is correct
20 Correct 137 ms 4500 KB Output is correct
21 Correct 130 ms 4496 KB Output is correct
22 Correct 133 ms 4504 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 2 ms 4308 KB Output is correct
2 Correct 2 ms 4308 KB Output is correct
3 Correct 13 ms 4444 KB Output is correct
4 Correct 112 ms 4544 KB Output is correct
5 Correct 1185 ms 5736 KB Output is correct
6 Execution timed out 1572 ms 11728 KB Time limit exceeded
7 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 2 ms 4436 KB Output is correct
2 Correct 3 ms 4416 KB Output is correct
3 Correct 12 ms 4448 KB Output is correct
4 Correct 134 ms 4564 KB Output is correct
5 Correct 105 ms 4564 KB Output is correct
6 Correct 121 ms 4564 KB Output is correct
7 Correct 126 ms 4544 KB Output is correct
8 Correct 121 ms 4552 KB Output is correct
9 Correct 115 ms 4436 KB Output is correct
10 Correct 102 ms 4436 KB Output is correct
11 Correct 104 ms 4500 KB Output is correct
12 Correct 113 ms 4436 KB Output is correct
13 Correct 136 ms 4492 KB Output is correct
14 Correct 144 ms 4492 KB Output is correct
15 Correct 118 ms 4492 KB Output is correct
16 Correct 109 ms 4500 KB Output is correct
17 Correct 111 ms 4436 KB Output is correct
18 Correct 117 ms 4496 KB Output is correct
19 Correct 108 ms 4496 KB Output is correct
20 Correct 137 ms 4500 KB Output is correct
21 Correct 130 ms 4496 KB Output is correct
22 Correct 133 ms 4504 KB Output is correct
23 Correct 2 ms 4308 KB Output is correct
24 Correct 2 ms 4308 KB Output is correct
25 Correct 13 ms 4444 KB Output is correct
26 Correct 112 ms 4544 KB Output is correct
27 Correct 1185 ms 5736 KB Output is correct
28 Execution timed out 1572 ms 11728 KB Time limit exceeded
29 Halted 0 ms 0 KB -