oj.uz

Submission #93640

# Submission time Handle Problem Language Result Execution time Memory
93640 2019-01-10T13:56:16 Z samzhang Construction of Highway (JOI18_construction) C++14
100 / 100
 817 ms 23836 KB 
//#pragma comment(linker, "/stack:200000000")
//#pragma GCC optimize("Ofast")
//#pragma GCC optimize(3)
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
//#pragma GCC target("sse3","sse2","sse")
//#pragma GCC target("avx","sse4","sse4.1","sse4.2","ssse3")
//#pragma GCC target("f16c")
//#pragma GCC optimize("inline","fast-math","unroll-loops","no-stack-protector")
//#pragma GCC diagnostic error "-fwhole-program"
//#pragma GCC diagnostic error "-fcse-skip-blocks"
//#pragma GCC diagnostic error "-funsafe-loop-optimizations"
//#pragma GCC diagnostic error "-std=c++14"
#include "bits/stdc++.h"
//#include "ext/pb_ds/tree_policy.hpp"
//#include "ext/pb_ds/assoc_container.hpp"
#define PB push_back
#define PF push_front
#define LB lower_bound
#define UB upper_bound
#define fr(x) freopen(x,"r",stdin)
#define fw(x) freopen(x,"w",stdout)
#define iout(x) printf("%d\n",x)
#define lout(x) printf("%lld\n",x)
#define REP(x,l,u) for(ll x = l;x<u;x++)
#define RREP(x,l,u) for(ll x = l;x>=u;x--)
#define complete_unique(a) a.erase(unique(a.begin(),a.end()),a.end())
#define mst(x,a) memset(x,a,sizeof(x))
#define all(a) a.begin(),a.end()
#define PII pair<int,int>
#define PLL pair<ll,ll>
#define MP make_pair
#define sqr(x) ((x)*(x))
#define lowbit(x) ((x)&(-(x)))
#define lson (ind<<1)
#define rson (ind<<1|1)
#define se second
#define fi first
#define sz(x) ((int)x.size())
#define EX0 exit(0);

typedef  long long ll;
typedef unsigned long long ull;
typedef double db;
typedef long double ld;
using namespace std;
typedef vector<ll> VLL;
typedef vector<int> VI;
const int block_size = 320;
typedef complex<ll> point;
const ll mod = 1e9+7;
const ll inf = 1e9+7;
const ld eps = 1e-9;
const db PI = atan(1)*4;
template<typename T>
inline int sign(const T&a) {
    if(a<0)return -1;
    if(a>0)return 1;
    return 0;
}
string to_string(string s) {
    return '"' + s + '"';
}

string to_string(const char* s) {
    return to_string((string) s);
}

string to_string(bool b) {
    return (b ? "true" : "false");
}

template <typename A, typename B>
string to_string(pair<A, B> p) {
    return "(" + to_string(p.first) + ", " + to_string(p.second) + ")";
}

template <typename A>
string to_string(A v) {
    bool first = true;
    string res = "{";
    for (const auto &x : v) {
        if (!first) {
            res += ", ";
        }
        first = false;
        res += to_string(x);
    }
    res += "}";
    return res;
}

void debug_out() { cerr << endl; }

template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T) {
    cerr << " " << to_string(H);
    debug_out(T...);
}

#ifndef ONLINE_JUDGE
#define dbg(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)
#else
#define dbg(...) {}
#endif

template<typename T,typename S>inline bool upmin(T&a,const S&b){return a>b?a=b,1:0;}
template<typename T,typename S>inline bool upmax(T&a,const S&b){return a<b?a=b,1:0;}

template<typename T> inline void in(T &x) {
    x = 0;
    T f = 1;
    char ch = getchar();
    while (!isdigit(ch)) {
        if (ch == '-') f = -1;
        ch = getchar();
    }
    while (isdigit(ch))  {
        x = x * 10 + ch - '0';
        ch = getchar();
    }
    x *= f;
}

ll twop(int x) {
    return 1LL<<x;
}

template<typename A,typename B > inline void in(A&x,B&y) {
    in(x);
    in(y);
}
template<typename A,typename B,typename C>inline void in(A&x,B&y,C&z) {
    in(x);
    in(y);
    in(z);
}
template<typename A,typename B,typename C,typename D> inline void in(A&x,B&y,C&z,D&d) {
    in(x);
    in(y);
    in(z);
    in(d);
}
ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}
namespace SOLVE {
    
    int color[100010];
    int relabel[100010];
    struct SegTree{
        static const int maxn = 100010;
        
        struct node{
            int l,r;
            int cover;
        };
        
        node no[maxn*4];
        void push_up(int ind){
            
        }
        void push_down(int ind){
            if(no[ind].cover >= 1){
                no[lson].cover = no[rson].cover = no[ind].cover;
                no[ind].cover = 0;
            }
        }
        void build(int l,int r,int ind){
            no[ind].l = l;
            no[ind].r = r;
            if(l == r){
                (no[ind].cover = relabel[l]);
            }else{
                int mid = (l+r)/2;
                build(l,mid,lson);
                build(mid+1,r,rson);
                push_up(ind);
            }
        }
        void update(int l,int r,int ind,int val){
            if(l>no[ind].r || r<no[ind].l)return;
            if(l<=no[ind].l && no[ind].r <= r){
                no[ind].cover = val;
            }else{
                push_down(ind);
                update(l,r,lson,val);
                update(l,r,rson,val);
                push_up(ind);
            }
        }
        void query(int l,int r,int ind,vector<PII>& ans){
            if(l>no[ind].r || r<no[ind].l)return;
            if(no[ind].cover){
                ans.PB(MP(min(r,no[ind].r) - max(l,no[ind].l)+1,no[ind].cover));
            }else{
                push_down(ind);
                query(l,r,lson,ans);
                query(l,r,rson,ans);
                push_up(ind);
            }
        }
    };
    
    
    SegTree tree;

    namespace HLD {
        //0不能被使用
        struct edge {
            int to;
            edge(int x):to(x){}
        };
        
        const int root = 1;
        const int maxn = 100010;
        vector<edge>adj[maxn];
        int dfnToID[maxn],dfn[maxn],head[maxn],fa[maxn],dep[maxn],size[maxn],heavy[maxn],r[maxn],cnt = 1;
        void firstDfs(int cur,int _fa) {
            dep[cur] = dep[_fa]+1;
            size[cur]=1;
            fa[cur] = _fa;
            for(auto e:adj[cur]) {
                if(e.to!=_fa) {
                    firstDfs(e.to,cur);
                    size[cur]+=size[e.to];
                }
            }
            int heavyChild = 0;
            for(auto e:adj[cur]) {
                if(e.to!=_fa) {
                    if(size[e.to]>size[heavyChild]) {
                        heavyChild = e.to;
                    }
                }
            }
            heavy[cur] = heavyChild;
        }
        
        
        void secondDfs(int cur,int _fa) {
            if(cur!=heavy[_fa]) {
                head[cur] = cur;
            } else {
                head[cur] = head[_fa];
            }
            dfn[cur] = cnt++;
            r[cur] = dfn[cur];
            dfnToID[dfn[cur]] = cur;
            if(!heavy[cur])return;
            secondDfs(heavy[cur],cur);
            r[cur] = r[heavy[cur]];
            for(auto e:adj[cur]) {
                if(e.to==_fa||e.to==heavy[cur])continue;
                secondDfs(e.to,cur);
                r[cur] = r[e.to];
            }
        }
        void init() {
            firstDfs(root,0);
            secondDfs(root,0);
        }
        
        int kthFather(int k,int cur) {
            while(k) {
                if(head[cur] == cur) {
                    k--;
                    cur = fa[head[cur]];
                } else {
                    if(dep[cur]-dep[head[cur]]<=k) {
                        k-=dep[cur]-dep[head[cur]];
                        cur = head[cur];
                    } else {
                        return dfnToID[dfn[cur]-k];
                    }
                }
            }
            return cur;
        }
        int LCA(int u,int v) {
            while(head[u]!=head[v]) {
                if(dep[head[u]]>dep[head[v]])swap(u,v);
                v = fa[head[v]];
            }
            if(dep[u]<dep[v])return u;
            return v;
        }
        
        int dis(int u,int v){
            return dep[u]+dep[v]-2*dep[LCA(u, v)];
        }
        void getColors(int u,vector<PII>&ans){
            if(u == 0)return;
            getColors(fa[head[u]], ans);
            tree.query(dfn[head[u]], dfn[u], 1, ans);
        }
        bool cmp(PII a,PII b){
            return a.se>b.se;
        }
        
        ll gao(vector<PII>::iterator s,vector<PII>::iterator t){
            if(t == s+1)return 0;
            auto mid = s + (t-s)/2;
            ll ans = gao(s, mid) + gao(mid, t);
            ll cnt = 0;
            auto cur = s;
            for(auto i = mid;i!=t;i++){
                while(cur!=mid and cur->se > i->se){
                    cnt += cur->fi;
                    ++cur;
                }
                ans += cnt * i->fi;
            }
            inplace_merge(s, mid, t,cmp);
            
            
            
            
            
            
            return ans;
        }
        
        void solve(int u){
            vector<PII>ans;
            getColors(u, ans);
            int last_color = ans.back().se;
            ans.pop_back();
            cout<<gao(begin(ans), end(ans))<<endl;
            while(u!=0){
                tree.update(dfn[head[u]], dfn[u],1,last_color);
                u = fa[head[u]];
            }
        }
        
    }
    vector<PII>pr;
    void main(){
        int n;in(n);
        REP(i,1,n+1)in(color[i]);
        REP(i,1,n){
            int a,b;in(a,b);
            HLD::adj[a].PB(b);
            pr.PB(MP(a,b));
        }
        HLD::init();
        REP(i,1,n+1){
            relabel[HLD::dfn[i]] = color[i];
        }
        tree.build(1, n, 1);
        for(auto i:pr){
            HLD::solve(i.se);
        }
        
        
        
        
        
        
        
        
        
        
    }
}


signed main() {
    
    
    
    
    
    int t;
//    in(t);
    t = 1;
    while(t--){
        SOLVE::main();
    }
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    return 0;
}

Subtask #1
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 2680 KB Output is correct
2 Correct 4 ms 2680 KB Output is correct
3 Correct 4 ms 2680 KB Output is correct
4 Correct 4 ms 2808 KB Output is correct
5 Correct 5 ms 2808 KB Output is correct
6 Correct 5 ms 2808 KB Output is correct
7 Correct 5 ms 2812 KB Output is correct
8 Correct 5 ms 2808 KB Output is correct
9 Correct 5 ms 2808 KB Output is correct
10 Correct 5 ms 2808 KB Output is correct
11 Correct 5 ms 2808 KB Output is correct
12 Correct 5 ms 2936 KB Output is correct
13 Correct 5 ms 2812 KB Output is correct
14 Correct 4 ms 2808 KB Output is correct
15 Correct 5 ms 2808 KB Output is correct
16 Correct 5 ms 2812 KB Output is correct
17 Correct 5 ms 2808 KB Output is correct
18 Correct 5 ms 2808 KB Output is correct
19 Correct 5 ms 2808 KB Output is correct
20 Correct 5 ms 2940 KB Output is correct
21 Correct 5 ms 2808 KB Output is correct
22 Correct 5 ms 2808 KB Output is correct
23 Correct 5 ms 2808 KB Output is correct
24 Correct 5 ms 2808 KB Output is correct
25 Correct 5 ms 2808 KB Output is correct
26 Correct 4 ms 2808 KB Output is correct
27 Correct 4 ms 2808 KB Output is correct
28 Correct 4 ms 2936 KB Output is correct
29 Correct 4 ms 2808 KB Output is correct
30 Correct 5 ms 2808 KB Output is correct
31 Correct 5 ms 2808 KB Output is correct
32 Correct 5 ms 2780 KB Output is correct
33 Correct 5 ms 2808 KB Output is correct
34 Correct 5 ms 2808 KB Output is correct
35 Correct 5 ms 2808 KB Output is correct
36 Correct 5 ms 2808 KB Output is correct
37 Correct 5 ms 2808 KB Output is correct
38 Correct 5 ms 2808 KB Output is correct
39 Correct 5 ms 2808 KB Output is correct
40 Correct 5 ms 2808 KB Output is correct
41 Correct 5 ms 2808 KB Output is correct
42 Correct 6 ms 2808 KB Output is correct
43 Correct 5 ms 2808 KB Output is correct
44 Correct 5 ms 2808 KB Output is correct

Subtask #2
9.0 / 9.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 2680 KB Output is correct
2 Correct 4 ms 2680 KB Output is correct
3 Correct 4 ms 2680 KB Output is correct
4 Correct 4 ms 2808 KB Output is correct
5 Correct 5 ms 2808 KB Output is correct
6 Correct 5 ms 2808 KB Output is correct
7 Correct 5 ms 2812 KB Output is correct
8 Correct 5 ms 2808 KB Output is correct
9 Correct 5 ms 2808 KB Output is correct
10 Correct 5 ms 2808 KB Output is correct
11 Correct 5 ms 2808 KB Output is correct
12 Correct 5 ms 2936 KB Output is correct
13 Correct 5 ms 2812 KB Output is correct
14 Correct 4 ms 2808 KB Output is correct
15 Correct 5 ms 2808 KB Output is correct
16 Correct 5 ms 2812 KB Output is correct
17 Correct 5 ms 2808 KB Output is correct
18 Correct 5 ms 2808 KB Output is correct
19 Correct 5 ms 2808 KB Output is correct
20 Correct 5 ms 2940 KB Output is correct
21 Correct 5 ms 2808 KB Output is correct
22 Correct 5 ms 2808 KB Output is correct
23 Correct 5 ms 2808 KB Output is correct
24 Correct 5 ms 2808 KB Output is correct
25 Correct 5 ms 2808 KB Output is correct
26 Correct 4 ms 2808 KB Output is correct
27 Correct 4 ms 2808 KB Output is correct
28 Correct 4 ms 2936 KB Output is correct
29 Correct 4 ms 2808 KB Output is correct
30 Correct 5 ms 2808 KB Output is correct
31 Correct 5 ms 2808 KB Output is correct
32 Correct 5 ms 2780 KB Output is correct
33 Correct 5 ms 2808 KB Output is correct
34 Correct 5 ms 2808 KB Output is correct
35 Correct 5 ms 2808 KB Output is correct
36 Correct 5 ms 2808 KB Output is correct
37 Correct 5 ms 2808 KB Output is correct
38 Correct 5 ms 2808 KB Output is correct
39 Correct 5 ms 2808 KB Output is correct
40 Correct 5 ms 2808 KB Output is correct
41 Correct 5 ms 2808 KB Output is correct
42 Correct 6 ms 2808 KB Output is correct
43 Correct 5 ms 2808 KB Output is correct
44 Correct 5 ms 2808 KB Output is correct
45 Correct 7 ms 2808 KB Output is correct
46 Correct 17 ms 3164 KB Output is correct
47 Correct 17 ms 3192 KB Output is correct
48 Correct 17 ms 3192 KB Output is correct
49 Correct 13 ms 3576 KB Output is correct
50 Correct 13 ms 3576 KB Output is correct
51 Correct 14 ms 3576 KB Output is correct
52 Correct 13 ms 3320 KB Output is correct
53 Correct 14 ms 3448 KB Output is correct
54 Correct 14 ms 3320 KB Output is correct
55 Correct 14 ms 3320 KB Output is correct
56 Correct 14 ms 3320 KB Output is correct
57 Correct 20 ms 3192 KB Output is correct
58 Correct 20 ms 3192 KB Output is correct
59 Correct 21 ms 3192 KB Output is correct
60 Correct 20 ms 3192 KB Output is correct
61 Correct 15 ms 3320 KB Output is correct
62 Correct 15 ms 3448 KB Output is correct
63 Correct 15 ms 3320 KB Output is correct
64 Correct 17 ms 3064 KB Output is correct
65 Correct 18 ms 3196 KB Output is correct
66 Correct 39 ms 3172 KB Output is correct
67 Correct 17 ms 3192 KB Output is correct
68 Correct 14 ms 3576 KB Output is correct
69 Correct 15 ms 3320 KB Output is correct
70 Correct 14 ms 3320 KB Output is correct
71 Correct 14 ms 3320 KB Output is correct
72 Correct 21 ms 3192 KB Output is correct
73 Correct 21 ms 3192 KB Output is correct
74 Correct 16 ms 3448 KB Output is correct
75 Correct 16 ms 3192 KB Output is correct
76 Correct 16 ms 3192 KB Output is correct
77 Correct 16 ms 3192 KB Output is correct
78 Correct 16 ms 3192 KB Output is correct
79 Correct 15 ms 3192 KB Output is correct
80 Correct 16 ms 3192 KB Output is correct
81 Correct 18 ms 3196 KB Output is correct
82 Correct 17 ms 3192 KB Output is correct
83 Correct 18 ms 3192 KB Output is correct
84 Correct 17 ms 3196 KB Output is correct
85 Correct 17 ms 3240 KB Output is correct
86 Correct 18 ms 3192 KB Output is correct

Subtask #3
84.0 / 84.0

# Verdict Execution time Memory Grader output
1 Correct 4 ms 2680 KB Output is correct
2 Correct 4 ms 2680 KB Output is correct
3 Correct 4 ms 2680 KB Output is correct
4 Correct 4 ms 2808 KB Output is correct
5 Correct 5 ms 2808 KB Output is correct
6 Correct 5 ms 2808 KB Output is correct
7 Correct 5 ms 2812 KB Output is correct
8 Correct 5 ms 2808 KB Output is correct
9 Correct 5 ms 2808 KB Output is correct
10 Correct 5 ms 2808 KB Output is correct
11 Correct 5 ms 2808 KB Output is correct
12 Correct 5 ms 2936 KB Output is correct
13 Correct 5 ms 2812 KB Output is correct
14 Correct 4 ms 2808 KB Output is correct
15 Correct 5 ms 2808 KB Output is correct
16 Correct 5 ms 2812 KB Output is correct
17 Correct 5 ms 2808 KB Output is correct
18 Correct 5 ms 2808 KB Output is correct
19 Correct 5 ms 2808 KB Output is correct
20 Correct 5 ms 2940 KB Output is correct
21 Correct 5 ms 2808 KB Output is correct
22 Correct 5 ms 2808 KB Output is correct
23 Correct 5 ms 2808 KB Output is correct
24 Correct 5 ms 2808 KB Output is correct
25 Correct 5 ms 2808 KB Output is correct
26 Correct 4 ms 2808 KB Output is correct
27 Correct 4 ms 2808 KB Output is correct
28 Correct 4 ms 2936 KB Output is correct
29 Correct 4 ms 2808 KB Output is correct
30 Correct 5 ms 2808 KB Output is correct
31 Correct 5 ms 2808 KB Output is correct
32 Correct 5 ms 2780 KB Output is correct
33 Correct 5 ms 2808 KB Output is correct
34 Correct 5 ms 2808 KB Output is correct
35 Correct 5 ms 2808 KB Output is correct
36 Correct 5 ms 2808 KB Output is correct
37 Correct 5 ms 2808 KB Output is correct
38 Correct 5 ms 2808 KB Output is correct
39 Correct 5 ms 2808 KB Output is correct
40 Correct 5 ms 2808 KB Output is correct
41 Correct 5 ms 2808 KB Output is correct
42 Correct 6 ms 2808 KB Output is correct
43 Correct 5 ms 2808 KB Output is correct
44 Correct 5 ms 2808 KB Output is correct
45 Correct 7 ms 2808 KB Output is correct
46 Correct 17 ms 3164 KB Output is correct
47 Correct 17 ms 3192 KB Output is correct
48 Correct 17 ms 3192 KB Output is correct
49 Correct 13 ms 3576 KB Output is correct
50 Correct 13 ms 3576 KB Output is correct
51 Correct 14 ms 3576 KB Output is correct
52 Correct 13 ms 3320 KB Output is correct
53 Correct 14 ms 3448 KB Output is correct
54 Correct 14 ms 3320 KB Output is correct
55 Correct 14 ms 3320 KB Output is correct
56 Correct 14 ms 3320 KB Output is correct
57 Correct 20 ms 3192 KB Output is correct
58 Correct 20 ms 3192 KB Output is correct
59 Correct 21 ms 3192 KB Output is correct
60 Correct 20 ms 3192 KB Output is correct
61 Correct 15 ms 3320 KB Output is correct
62 Correct 15 ms 3448 KB Output is correct
63 Correct 15 ms 3320 KB Output is correct
64 Correct 17 ms 3064 KB Output is correct
65 Correct 18 ms 3196 KB Output is correct
66 Correct 39 ms 3172 KB Output is correct
67 Correct 17 ms 3192 KB Output is correct
68 Correct 14 ms 3576 KB Output is correct
69 Correct 15 ms 3320 KB Output is correct
70 Correct 14 ms 3320 KB Output is correct
71 Correct 14 ms 3320 KB Output is correct
72 Correct 21 ms 3192 KB Output is correct
73 Correct 21 ms 3192 KB Output is correct
74 Correct 16 ms 3448 KB Output is correct
75 Correct 16 ms 3192 KB Output is correct
76 Correct 16 ms 3192 KB Output is correct
77 Correct 16 ms 3192 KB Output is correct
78 Correct 16 ms 3192 KB Output is correct
79 Correct 15 ms 3192 KB Output is correct
80 Correct 16 ms 3192 KB Output is correct
81 Correct 18 ms 3196 KB Output is correct
82 Correct 17 ms 3192 KB Output is correct
83 Correct 18 ms 3192 KB Output is correct
84 Correct 17 ms 3196 KB Output is correct
85 Correct 17 ms 3240 KB Output is correct
86 Correct 18 ms 3192 KB Output is correct
87 Correct 45 ms 3960 KB Output is correct
88 Correct 144 ms 6264 KB Output is correct
89 Correct 587 ms 14568 KB Output is correct
90 Correct 608 ms 14616 KB Output is correct
91 Correct 561 ms 14700 KB Output is correct
92 Correct 315 ms 23788 KB Output is correct
93 Correct 352 ms 23788 KB Output is correct
94 Correct 335 ms 23836 KB Output is correct
95 Correct 331 ms 18792 KB Output is correct
96 Correct 325 ms 19180 KB Output is correct
97 Correct 320 ms 19172 KB Output is correct
98 Correct 320 ms 19436 KB Output is correct
99 Correct 316 ms 18516 KB Output is correct
100 Correct 710 ms 14740 KB Output is correct
101 Correct 817 ms 14956 KB Output is correct
102 Correct 771 ms 15084 KB Output is correct
103 Correct 804 ms 15072 KB Output is correct
104 Correct 429 ms 18616 KB Output is correct
105 Correct 405 ms 18412 KB Output is correct
106 Correct 407 ms 18412 KB Output is correct
107 Correct 592 ms 13676 KB Output is correct
108 Correct 593 ms 13932 KB Output is correct
109 Correct 588 ms 14184 KB Output is correct
110 Correct 303 ms 23280 KB Output is correct
111 Correct 323 ms 18796 KB Output is correct
112 Correct 337 ms 18668 KB Output is correct
113 Correct 336 ms 17772 KB Output is correct
114 Correct 723 ms 14908 KB Output is correct
115 Correct 766 ms 14440 KB Output is correct
116 Correct 407 ms 17784 KB Output is correct
117 Correct 395 ms 16112 KB Output is correct
118 Correct 356 ms 15468 KB Output is correct
119 Correct 452 ms 15132 KB Output is correct
120 Correct 416 ms 15772 KB Output is correct
121 Correct 374 ms 15284 KB Output is correct
122 Correct 447 ms 14860 KB Output is correct
123 Correct 465 ms 16280 KB Output is correct
124 Correct 442 ms 15608 KB Output is correct
125 Correct 544 ms 15208 KB Output is correct
126 Correct 497 ms 15916 KB Output is correct
127 Correct 475 ms 15336 KB Output is correct
128 Correct 532 ms 14864 KB Output is correct