Submission #1301584

#TimeUsernameProblemLanguageResultExecution timeMemory
1301584EkinOnalRegions (IOI09_regions)C++20
95 / 100
7220 ms48476 KiB
//#pragma GCC optimize("O3,unroll-loops,Ofast") //#pragma GCC target("avx2,bmi,bmi2,popcnt,lzcnt") #include <bits/stdc++.h> // #include <ext/pb_ds/assoc_container.hpp> // #include <ext/pb_ds/tree_policy.hpp> using namespace std; // using namespace __gnu_pbds; #define MAX 200005 #define pb push_back // #define mp make_pair // #define int long long // #define double long double #define f first #define s second #define vi vector<int> #define pii pair<int,int> #define si set<int> #define vpii vector<pair<int,int>> const int mod = 1e9+7; const int INF = 1e9; #define SQRT 560 // myMap.begin()->first : key // myMap.begin()->second : value int epow(int a,int b){int ans=1;while(b){if(b&1) ans*=a;a*=a;b>>=1;ans%=mod;a%=mod;}return ans%mod;} int gcd(int a,int b) {if(a<b)swap(a,b);while(b){int tmp=b;b=a%b;a=tmp;}return a;} int mul(int a,int b){return ((a%mod)*(b%mod))%mod;} int sum(int a,int b){return ((a%mod)+(b%mod))%mod;} //typedef tree<int,null_type,less<int>,rb_tree_tag,tree_order_statistics_node_update>ordered_set; // typedef // tree<int,null_type,less_equal<int>,rb_tree_tag,tree_order_statistics_node_update>ordered_multiset; int tin[MAX]; vpii euler[25005]; int tout[MAX],timer=1,dep[MAX],type[MAX]; vi adj[MAX]; void dfs(int node,int par){ tin[node]=timer++; dep[node]=dep[par]+1; euler[type[node]].pb({tin[node],node}); for(auto u : adj[node]){ if(u!=par) dfs(u,node); } tout[node]=timer-1; } set<int> st[MAX]; void dfs2(int node,int par){ for(auto u : adj[node]){ if(u==par) continue; dfs2(u,node); if(st[u].size()>st[node].size()) swap(st[node],st[u]); for(auto x : st[u]) st[node].insert(x); } } void solve(){ int n,r,q; cin>>n>>r>>q; cin>>type[1]; for(int i=2;i<=n;i++){ int p,tur; cin>>p>>tur; adj[p].pb(i); adj[i].pb(p); type[i]=tur; } dfs(1,0); dfs2(1,0); // for(auto u : calc) cout<<u[0]<<" "<<u[1]<<" "<<u[2]<<endl; vpii ans[25005];; for(int i=1;i<=25000;i++){ if(euler[i].size()<=SQRT) continue; for(int j=1;j<=25000;j++){ int tot=0; for(auto x : euler[i]){ int node=x.s; int l = lower_bound(euler[j].begin(),euler[j].end(),make_pair(tin[node],0))-euler[j].begin(); int r = upper_bound(euler[j].begin(),euler[j].end(),make_pair(tout[node],INF))-euler[j].begin(); tot+=r-l; } if(tot) ans[i].pb({j,tot}); } } while(q--){ int a,b; cin>>a>>b; if(euler[a].size()>SQRT){ auto itr = lower_bound(ans[a].begin(),ans[a].end(),make_pair(b,-1)); if(itr==ans[a].end() || (*itr).f!=b) cout<<"0\n"; else cout<<(*itr).s<<endl; } else{ int ans=0; for(auto u : euler[a]){ int node=u.s; int l = lower_bound(euler[b].begin(),euler[b].end(),make_pair(tin[node],0))-euler[b].begin(); int r = upper_bound(euler[b].begin(),euler[b].end(),make_pair(tout[node],INF))-euler[b].begin(); ans+=r-l; } cout<<ans<<endl; } } } int32_t main(/*int32_t argc, char* argv[]*/){ std::ios_base::sync_with_stdio(0); std::cin.tie(0); // freopen("promote.in", "r", stdin); // freopen("promote.out", "w", stdout); int t=1; // cin >> t; while (t--) solve(); return 0; }
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...