Submission #373163

#	Time	Username	Problem	Language	Result	Execution time	Memory
373163		BartolM	Factories (JOI14_factories)	C++17	100 / 100	4848 ms	131484 KiB

factories

#include "factories.h"
#include <bits/stdc++.h>

using namespace std;

#define X first
#define Y second
#define mp make_pair
#define pb push_back
typedef long long ll;
typedef pair <int, int> pii;
typedef pair <int, pii> pip;
typedef pair <pii, int> ppi;
typedef pair <ll, ll> pll;

const int INF=0x3f3f3f3f;
const ll MAX=(ll)INF*INF;
const int MAXN=5e5+5;
const int LOG=20;
const int OFF=(1<<19);

int n, dtime=0, ind;
ll dist[MAXN];
int par[LOG][MAXN], dep[MAXN], L[MAXN], R[MAXN];
vector <pii> ls[MAXN];
ll da[MAXN], db[MAXN];
vector <int> v;

int cmp(int a, int b) {
    return L[a]<L[b];
}

void dfs(int node, int rod, ll d) {
    par[0][node]=rod;
    dist[node]=d;
    dep[node]=dep[rod]+1;
    L[node]=dtime++;
//    printf("node: %d, dist: %lld, dep: %d, par: %d\n", node, dist[node], dep[node], rod);
    for (pii sus:ls[node]) if (sus.X!=rod) dfs(sus.X, node, d+sus.Y);
    R[node]=dtime;
}

void build() {
    for (int i=1; i<LOG; ++i) {
        for (int j=0; j<n; ++j) {
            par[i][j]=par[i-1][par[i-1][j]];
        }
    }
}

int lca(int x, int y) {
    if (dep[x]<dep[y]) swap(x, y);
    for (int i=LOG-1; i>=0; --i) {
        if (dep[x]-(1<<i)>=dep[y]) x=par[i][x];
    }
    if (x==y) return x;
    for (int i=LOG-1; i>=0; --i) {
        if (par[i][x]!=par[i][y]) x=par[i][x], y=par[i][y];
    }
    return par[0][x];
}

void Init(int N, int A[], int B[], int D[]) {
    n=N;
    for (int i=0; i<n-1; ++i) {
        ls[A[i]].pb(mp(B[i], D[i]));
        ls[B[i]].pb(mp(A[i], D[i]));
    }
    dfs(0, 0, 0);
    build();
    for (int i=0; i<n; ++i) da[i]=db[i]=MAX;
}

ll rek() {
    ll ret=MAX;
    int node=v[ind++];
    while (ind<(int)v.size() && L[v[ind]]>=L[node] && L[v[ind]]<R[node]) {
        int sus=v[ind];
        ret=min(ret, rek());
        da[node]=min(da[node], da[sus]+dist[sus]-dist[node]);
        db[node]=min(db[node], db[sus]+dist[sus]-dist[node]);
    }
    ret=min(ret, da[node]+db[node]);
    return ret;
}

ll Query(int S, int x[], int T, int y[]) {
    v.clear();
    for (int i=0; i<S; ++i) {
        v.pb(x[i]);
        da[x[i]]=0;
    }
    for (int i=0; i<T; ++i) {
        v.pb(y[i]);
        db[y[i]]=0;
    }

    sort(v.begin(), v.end(), cmp);
    for (int i=1; i<S+T; ++i) {
        int lc=lca(v[i-1], v[i]);
        v.pb(lc);
    }


    v.pb(0);
    sort(v.begin(), v.end());
    v.erase(unique(v.begin(), v.end()), v.end());
    sort(v.begin(), v.end(), cmp);
//    for (int x:v) printf("%d ", x);
//    printf("\n");

    ind=0;
    ll sol=rek();

    for (int x:v) da[x]=db[x]=MAX;
    return sol;
}

• Subtask 1: 15 / 15
• Subtask 2: 18 / 18
• Subtask 3: 67 / 67