oj.uz

Submission #369368

# Submission time Handle Problem Language Result Execution time Memory
369368 2021-02-21T12:50:28 Z AmineWeslati Regions (IOI09_regions) C++14
85 / 100
 4344 ms 32896 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=400,MX2=25000;
int N,R,Q; 
vi a(MX),adj[MX],cnt(MX2,0),vec[MX2],vec2[MX2],in(MX),out(MX),mp[MX2];

//mp[x][y]
void dfs(int u, int c, int nb){
	for(auto v: adj[u]){
		mp[c][a[v]]+=nb;
		dfs(v,c,nb+(a[v]==c));
	}
}
void precompute1(){ 
	FOR(i,1,R+1) if(cnt[i]>C){
		mp[i].assign(R+1,0);
		dfs(1,i,(a[1]==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(){ 
	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(); //(N/C)*N
    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

Subtask #1
15.0 / 15.0

# Verdict Execution time Memory Grader output
1 Correct 7 ms 9196 KB Output is correct
2 Correct 8 ms 9196 KB Output is correct
3 Correct 10 ms 9196 KB Output is correct
4 Correct 11 ms 9324 KB Output is correct
5 Correct 15 ms 9452 KB Output is correct
6 Correct 25 ms 9324 KB Output is correct
7 Correct 36 ms 9324 KB Output is correct
8 Correct 48 ms 9452 KB Output is correct
9 Correct 56 ms 9836 KB Output is correct
10 Correct 64 ms 9708 KB Output is correct
11 Correct 153 ms 10008 KB Output is correct
12 Correct 199 ms 10476 KB Output is correct
13 Correct 183 ms 10092 KB Output is correct
14 Correct 350 ms 10476 KB Output is correct
15 Correct 316 ms 13036 KB Output is correct

Subtask #2
70.0 / 85.0

# Verdict Execution time Memory Grader output
1 Correct 2080 ms 13324 KB Output is correct
2 Correct 2055 ms 11912 KB Output is correct
3 Correct 3558 ms 15280 KB Output is correct
4 Correct 324 ms 10604 KB Output is correct
5 Correct 490 ms 12172 KB Output is correct
6 Correct 669 ms 13360 KB Output is correct
7 Correct 1759 ms 13644 KB Output is correct
8 Correct 1420 ms 22104 KB Output is correct
9 Correct 2617 ms 17024 KB Output is correct
10 Correct 3822 ms 32896 KB Output is correct
11 Incorrect 4344 ms 16336 KB Output isn't correct
12 Correct 1773 ms 17900 KB Output is correct
13 Correct 2250 ms 18284 KB Output is correct
14 Correct 2789 ms 19532 KB Output is correct
15 Incorrect 4119 ms 23364 KB Output isn't correct
16 Incorrect 3743 ms 30576 KB Output isn't correct
17 Correct 3552 ms 30552 KB Output is correct