Submission #1265000

#	Time	Username	Problem	Language	Result	Execution time	Memory
1265000		_ll	Harbingers (CEOI09_harbingers)	C++20	20 / 100	3 ms	3264 KiB

harbingers

#include<bits/stdc++.h>// persistent stack - binary lifting - convex hull trick (Obra de arteeeeee)
using namespace std;
typedef long long ll;
typedef long double ld;
vector<array<ll, 2>> g[10000]; 
array<ll, 2> all[10000];
ld oo = 1e12, re[10000];
struct no {
        array<ld, 4> inf = {1, oo, -1, -1}; // coeficientes e intersecção
        array<ll, 20> bf = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};// binary lifting
        vector<ll> f;     // filhos
};
vector<no> ps;

array<ld, 2> inter(array<ld, 2> a, array<ld, 2> b){
        if(b[0] == a[0]) return {-1, -1}; // pelas ruas serem não negativas isso nunca acontece (o multiplicador sempre diminui)
        ld x = (b[1] - a[1]) / (a[0] - b[0]);
        return {x, b[0] * x + b[1]};
}

ll men(array<ld, 4> a, array<ld, 2> b){ // me diz se é menor
        array<ld, 2> xy = inter({a[0], a[1]}, b); // preciso da intersecção de a com b;
        return(xy[0] <= a[2]);
}

void nps(ll &ip, no &at){ // vou criar um ip auxiliar pois vou mudar as informações
    at.bf[0] = ip;
    for(ll i = 1; i < 20; i++){
        at.bf[i] = ps[at.bf[i - 1]].bf[i - 1];
    }
    ps[ip].f.push_back(ps.size());
    ip = ps.size();
    ps.push_back(at);
}

void ins(ll &at, array<ld, 2> a){
        while(ps[at].bf[0] > 0 && men(ps[ps[at].bf[0]].inf, a)) // olho o pai do ultimo se o ponto de termino for maior 
                at = ps[at].bf[0]; // que o ponto de intersecção sei q o atual continua o do pai
        if(at > 0){
            no aux = ps[at];
            array<ld, 2> xy = inter({ps[at].inf[0], ps[at].inf[1]}, a);
            aux.inf[2] = xy[0], aux.inf[3] = xy[1];
            ps[aux.bf[0]].f.push_back(ps.size());
            at = ps.size();
            ps.push_back(aux);
        }
        no f;
        f.inf = {a[0], a[1], oo, oo};
        nps(at, f); // arruma o binarylifting, coloca ele na pilha e muda o at
} // agora falta o consulta na posição x a partir de um indice atual da pilha

ll vl(ll iat, ll x){
    for(ll i = 20; i >= 0; i--) if(ps[ps[iat].bf[i]].inf[2] >= x)
        iat = ps[iat].bf[i];
    return (ll) (ps[iat].inf[0] * x + ps[iat].inf[1]); // acha o minimo em um x usando a queue com iat no topo
}

void dfs(ll p, ll at, ll del, ll ia){
        ll res = vl(ia, all[at][1]), c = all[at][0] + all[at][1] * del; // olha a beleza disso
        re[at] = res + c;
        ins(ia, {(ld) -del, (ld) res + c});
        for(auto [i, j] : g[at]) if(i != p)
                dfs(at, i, del + j, ia);// tenho q somar isso nas arestas
}

int main(){
        ios_base::sync_with_stdio(0);
        cin.tie(0); cout.tie(0);

        ll n; cin >> n;
        for(ll i = 1, a, b, c; i < n; i++){
                cin >> a >> b >> c; a--, b--;
                g[a].push_back({b, c});
                g[b].push_back({a, c});
        }
        for(ll i = 1; i < n; i++)
                cin >> all[i][0] >> all[i][1];
        no ini, z;
        ps.push_back(ini);
        z.inf = {0, 0, oo, oo};
        ll at = 0;
        nps(at, z);
        dfs(-1, 0, 0, at);
        for(ll i = 1; i < n; i++)
                cout << to_string((ll) re[i]) << " ";
        cout << "\n";
} 

• Subtask 1: 20 / 100