Submission #369367

# Submission time Handle Problem Language Result Execution time Memory
369367 2021-02-21T12:40:20 Z AmineWeslati Regions (IOI09_regions) C++14
65 / 100
8000 ms 39020 KB
//Never stop trying
#include "bits/stdc++.h"
using namespace std;
#define boost ios_base::sync_with_stdio(false); cin.tie(0); cout.tie(0)

typedef long long ll;
typedef string str;
typedef long double ld;
typedef pair<int, int> pi;
#define fi first
#define se second
typedef vector<int> vi;
typedef vector<pi> vpi;
#define pb push_back
#define eb emplace_back
#define sz(x) (int)x.size()
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)
#define endl "\n"
#define FOR(i,a,b) for (int i = (a); i < (b); ++i)
#define ROF(i,a,b) for (int i = (b)-1; i >= (a); --i)

const int MOD = 1e9 + 7; //998244353
const ll INF = 1e18;
const int MX = 2e5 + 10;
const int nx[4] = {0, 0, 1, -1}, ny[4] = {1, -1, 0, 0}; //right left down up

template<class T> using V = vector<T>;
template<class T> bool ckmin(T& a, const T& b) { return a > b ? a = b, 1 : 0; }
template<class T> bool ckmax(T& a, const T& b) { return a < b ? a = b, 1 : 0; }
ll cdiv(ll a, ll b) { return a / b + ((a ^ b) > 0 && a % b); } // divide a by b rounded up
//constexpr int log2(int x) { return 31 - __builtin_clz(x); } // floor(log2(x))

mt19937 rng(chrono::system_clock::now().time_since_epoch().count());
//mt19937_64 rng(chrono::system_clock::now().time_since_epoch().count());
ll random(ll a, ll b){
    return a + rng() % (b - a + 1);
}
#ifndef LOCAL  
#define cerr if(false) cerr
#endif
#define dbg(x) cerr << #x << " : " << x << endl; 
#define dbgs(x,y) cerr << #x << " : " << x << " / " << #y << " : " << y << endl;
#define dbgv(v) cerr << #v << " : " << "[ "; for(auto it : v) cerr << it << ' '; cerr << ']' << endl;
#define here() cerr << "here" << endl;

void IO() {
#ifdef LOCAL
    freopen("input.txt", "r", stdin); 
    freopen("output.txt", "w", stdout);
#endif
}
/////////////////////////ONLY CLEAN CODES ALLOWED/////////////////////////

const int C=800;
int N,R,Q; 
vi a(MX),adj[MX],cnt(MX,0),vec[MX],vec2[MX],in(MX),out(MX);
//unordered_map<int,int>mp[MX];
vi mp[MX];

//mp[x][y]
void dfs(int u, int c){
	for(auto v: adj[u]){
		mp[c][a[v]]++;
		dfs(v,c);
	}
}
void precompute1(){ //(N/C)*N
	FOR(i,1,R+1) if(cnt[i]>C){
		mp[i].resize(R+1);
		for(auto j: vec[i]) dfs(j,i);
	}
}

//in[u],out[u]
//vec2[c]
int cntr=0;
void euler(int u=1, int p=1){
	in[u]=++cntr;
	for(auto v: adj[u]) euler(v,u);
	out[u]=cntr;
}
void precompute2(){ //O(NlgN)
	euler();
	FOR(i,1,N+1) vec2[a[i]].pb(in[i]);
	FOR(i,1,R+1) sort(all(vec2[i]));
}


int solve(int l, int r, int c){
	if(l>r) return 0;
	return upper_bound(all(vec2[c]),r)-lower_bound(all(vec2[c]),l);
}

int main() {
    boost; IO();

    cin>>N>>R>>Q;
    FOR(i,1,N+1){
    	if(i>1){
    		int j; cin>>j;
    		adj[j].pb(i);
    	}
    	cin>>a[i];
    	vec[a[i]].pb(i);
    	cnt[a[i]]++;
    }

    precompute1();
    precompute2();

    //Q*C*lg(N)
    while(Q--){ 
    	int x,y; cin>>x>>y;
    	if(cnt[x]>C){
    		cout << mp[x][y] << endl;
    		cout.flush();
    	}
    	else{
    		int ans=0;
    		for(auto u: vec[x]){
    			int l=in[u],r=out[u];
    			ans+=solve(l+1,r,y);
    		}
    		cout << ans << endl;
    		cout.flush();
    	}
    }

    return 0;
}
//Change your approach 
# Verdict Execution time Memory Grader output
1 Correct 15 ms 22252 KB Output is correct
2 Correct 16 ms 22252 KB Output is correct
3 Correct 18 ms 22252 KB Output is correct
4 Correct 19 ms 22252 KB Output is correct
5 Correct 22 ms 22252 KB Output is correct
6 Correct 34 ms 22380 KB Output is correct
7 Correct 40 ms 22400 KB Output is correct
8 Correct 42 ms 22380 KB Output is correct
9 Correct 67 ms 22892 KB Output is correct
10 Correct 71 ms 22764 KB Output is correct
11 Correct 165 ms 22892 KB Output is correct
12 Correct 185 ms 23404 KB Output is correct
13 Correct 200 ms 23148 KB Output is correct
14 Correct 268 ms 23596 KB Output is correct
15 Correct 288 ms 26092 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 6780 ms 26088 KB Output is correct
2 Correct 3165 ms 24896 KB Output is correct
3 Execution timed out 8041 ms 27116 KB Time limit exceeded
4 Correct 313 ms 23660 KB Output is correct
5 Correct 455 ms 25184 KB Output is correct
6 Correct 1611 ms 24684 KB Output is correct
7 Correct 1900 ms 25480 KB Output is correct
8 Correct 1514 ms 30116 KB Output is correct
9 Correct 2923 ms 29928 KB Output is correct
10 Correct 5401 ms 34732 KB Output is correct
11 Correct 4332 ms 29300 KB Output is correct
12 Execution timed out 8013 ms 30816 KB Time limit exceeded
13 Execution timed out 8050 ms 29676 KB Time limit exceeded
14 Execution timed out 8039 ms 29356 KB Time limit exceeded
15 Execution timed out 8064 ms 33516 KB Time limit exceeded
16 Execution timed out 8103 ms 39020 KB Time limit exceeded
17 Execution timed out 8058 ms 38016 KB Time limit exceeded