Pagini recente » Cod sursa (job #2318062) | Cod sursa (job #1260353) | Cod sursa (job #142712) | Cod sursa (job #2010888) | Cod sursa (job #1549412)
#include <cstdio>
#include <cstring>
#include <vector>
#include <cctype>
#include <algorithm>
using namespace std;
const int maxN = 100000;
const int maxSize = 1 << 20;
class ReadStream {
private:
int cursor;
char buf[maxSize];
FILE *f;
inline void move() {
++ cursor;
if(cursor == maxSize) {
cursor = 0;
fread(buf, maxSize, 1, f);
}
}
inline int todig(char c) {
return c - '0';
}
public:
ReadStream() {}
ReadStream(const char *fname) {
f = fopen(fname, "r");
cursor = 0;
fread(buf, maxSize, 1, f);
}
inline ReadStream &operator >> (int& N) {
char sign = '+';
while(!isdigit(buf[cursor])) {
sign = buf[cursor];
move();
}
N = 0;
while(isdigit(buf[cursor])) {
N = N * 10 + todig(buf[cursor]);
move();
}
if(sign == '-')
N = -N;
return *this;
}
inline ReadStream &operator >> (float& N) {
char sign = '+';
while(!isdigit(buf[cursor])) {
sign = buf[cursor];
move();
}
N = 0;
while(isdigit(buf[cursor])) {
N = N * 10 + todig(buf[cursor]);
move();
}
if(buf[cursor] == '.') {
float p = 0.1;
move();
while(isdigit(buf[cursor])) {
N = N + p * todig(buf[cursor]);
p *= 0.1;
move();
}
}
if(sign == '-')
N = -N;
return *this;
}
};
int N;
bool viz[1 + maxN];
vector<int> sons[1 + maxN];
int id[1 + maxN], ac;
int euler[maxN << 1], where[1 + maxN], e;
int rmq[20][maxN << 1];
int Log[(maxN << 1) + 1];
void rmqCompute(int v[], int N) {
for(int i = 0; i < N; ++ i)
rmq[0][i] = v[i];
for(int i = 1; (1 << i) <= N; ++ i)
Log[(1 << i)] = 1;
for(int i = 1; i <= N; ++ i)
Log[i] += Log[i - 1];
int p = 2;
for(int step = 1; p <= N; ++ step, p <<= 1) {
int put = p >> 1;
for(int i = 0; i + put < N; ++ i)
rmq[step][i] = min(rmq[step - 1][i], rmq[step - 1][i + put]);
}
}
int rmQuery(int x, int y) {
int l = y - x + 1;
int step = Log[l];
return min(rmq[step][x], rmq[step][y - (1 << step) + 1]);
}
void dfs(int node = 1) {
if(!viz[node]) {
viz[node] = true;
int nodeAc = ac;
id[ac ++] = node;
where[node] = e;
euler[e ++] = nodeAc;
for (int son: sons[node]) {
dfs(son);
euler[e ++] = nodeAc;
}
}
}
int lca(int x, int y) {
x = where[x];
y = where[y];
return id[rmQuery(min(x, y), max(x, y))];
}
int main() {
ReadStream in("lca.in");
freopen("lca.out", "w", stdout);
int M;
in >> N >> M;
for(int i = 2; i <= N; ++ i) {
int fatherI;
in >> fatherI;
sons[fatherI].push_back(i);
}
dfs();
rmqCompute(euler, e);
while(M --) {
int x, y;
in >> x >> y;
printf("%d\n", lca(x, y));
}
return 0;
}