Cod sursa(job #51240)

Utilizator scapryConstantin Berzan scapry Data 10 aprilie 2007 16:58:28
Problema Stramosi Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.34 kb
#include <stdio.h>
#include <string.h>

enum { maxn = 250001 };
int firstbit[maxn];

int n;
int dad[maxn];
bool v[maxn];
int ans[maxn][32]; //powers of 2 (we only need log2(E) < 19)

void precalc_firstbit()
{
	int i, bit = 0;

	for(i = 0; i < maxn; i++)
	{
		if(i >= (1 << bit)) bit++;
		firstbit[i] = bit;
	}
}

int calc(int node, int rank)
{
	if(rank > n) return 0;

	int bit;

	while(true)
	{
		if(node == 0) return 0;
		if(rank == 0) return node;

		bit = firstbit[rank];

		node = ans[node][bit];
		rank &= (~(1 << bit));
	}
}

void precalc(int node)
{
	if(v[node]) return;

	//terminal case - root
	if(dad[node] != 0)
	{
		if(!v[ dad[node] ]) precalc(dad[node]);

		for(int i = 0; i < 31; i++)
		{
			ans[node][i] = calc(dad[node], (1 << i) - 1);
			if(ans[node][i] == 0) break; //all will be zero from now on
		}
	}

	v[node] = true;
}

void wraprintf(FILE *f, int val)
{
	fprintf(f, "%d\n", val);
}

int main()
{
	int queries, i, who, rank;
	FILE *fi = fopen("stramosi.in", "r"),
	     *fo = fopen("stramosi.out", "w");
	if(!fi || !fo) return 1;

	fscanf(fi, "%d%d", &n, &queries);
	for(i = 1; i <= n; i++)
		fscanf(fi, "%d", &dad[i]);

	precalc_firstbit();

	for(i = 1; i <= n; i++)
		precalc(i);

	for(i = 0; i < queries; i++)
	{
		fscanf(fi, "%d %d ", &who, &rank);
		fprintf(fo, "%d\n", calc(who, rank));
	}

	fclose(fi);
	fclose(fo);
       	return 0;
}