제출 #474094

#제출 시각아이디문제언어결과실행 시간메모리
474094julian33Cat in a tree (BOI17_catinatree)C++14
0 / 100
9 ms10216 KiB
#include <bits/stdc++.h> using namespace std; #ifdef LOCAL #define deb(...) logger(#__VA_ARGS__, __VA_ARGS__) template<typename ...Args> void logger(string vars, Args&&... values) { cerr<<vars<<" = "; string delim=""; (...,(cerr<<delim<<values,delim=", ")); cerr<<"\n"; } #else #define deb(...) logger(#__VA_ARGS__, __VA_ARGS__) template<typename ...Args> void logger(string vars, Args&&... values) {} #endif #define pb push_back #define sz(x) (int)(x.size()) typedef long long ll; typedef pair<int,int> pii; typedef pair<ll,ll> pll; template<typename T> inline void maxa(T& a,T b){a=max(a,b);} template<typename T> inline void mina(T& a,T b){a=min(a,b);} mt19937 rng(chrono::steady_clock::now().time_since_epoch().count()); const int mxN=2e5+5,mxL=19; //make sure this is right const int mod=1e9+7; vector<int> graph[mxN],path; int ht[mxN],par[mxN]; priority_queue<pii> q; int up[mxN][mxL],ets[mxN],etn[mxN],st[8*mxN],lazy[8*mxN],cur; void push(int v){ if(lazy[v]==1e9) return; mina(st[v<<1],lazy[v]); mina(st[v<<1|1],lazy[v]); mina(lazy[v<<1],lazy[v]); mina(lazy[v<<1|1],lazy[v]); lazy[v]=1e9; } void update(int v,int l,int r,int lq,int rq,int val){ if(lq>rq) return; if(l>=lq && r<=rq){ mina(st[v],val); mina(lazy[v],val); } else{ int m=(l+r)>>1; push(v); update(v<<1,l,m,lq,min(rq,m),val); update(v<<1|1,m+1,r,max(lq,m+1),rq,val); st[v]=min(st[v<<1],st[v<<1|1]); } } int query(int v,int l,int r,int ind){ if(l>ind || r<ind) return 1e9; if(l==r) return st[v]; push(v); int m=(l+r)>>1; return min(query(v<<1,l,m,ind),query(v<<1|1,m+1,r,ind)); } void dfs(int at,int p=-1){ up[at][0]=~p?p:at; for(int i=1;i<mxL;i++) up[at][i]=up[up[at][i-1]][i-1]; ets[at]=++cur; for(int &i:graph[at]){ if(i==p) continue; ht[i]=ht[at]+1; dfs(i,at); } etn[at]=++cur; q.push({ht[at],at}); } int jump(int a,int k){ for(int i=0;i<mxL;i++){ if(k&(1<<i)) a=up[a][i]; } return a; } int lca(int a,int b){ if(ht[a]>ht[b]) swap(a,b); b=jump(b,ht[b]-ht[a]); if(a==b) return a; for(int i=mxL-1;i>=0;i--){ if(up[a][i]!=up[b][i]){ a=up[a][i]; b=up[b][i]; } } return up[a][0]; } int main(){ cin.sync_with_stdio(0); cin.tie(0); cout.tie(0); #ifdef LOCAL freopen("input.txt","r",stdin); freopen("output.txt","w",stdout); #endif int n,d; cin>>n>>d; for(int i=1;i<n;i++){ cin>>par[i]; graph[par[i]].pb(i); graph[i].pb(par[i]); } par[0]=-1; dfs(0); int ans=0; fill(st+1,st+1+8*n,1e9); fill(lazy+1,lazy+1+8*n,1e9); while(!q.empty()){ int at=q.top().second; q.pop(); int lo=query(1,1,2*n,ets[at]); if(lo+ht[at]<d) continue; for(int i=at;~i && ht[at]-ht[i]<=d;i=par[i]) update(1,1,n,ets[i],etn[i],ht[at]-2*ht[i]); ans++; } cout<<ans<<"\n"; }
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...