[CF543D]Road Improvement

#include <cstdio>
#include <vector>
#include <cctype>
namespace R {
	int x, ch;
	inline int read() {
		ch = getchar();
		while (isspace(ch)) ch = getchar();
		for (x = ch & 15, ch = getchar(); isdigit(ch); ch = getchar()) x = x * 10 + (ch & 15);
		return x;
	}
}
using R::read;

#define maxn 200010
const int mod = 1e9 + 7;
int head[maxn], cnt;
struct Edge {
	int to, nxt;
} e[maxn << 1];
inline void add(int a, int b) {
	e[++cnt] = (Edge) {b, head[a]}; head[a] = cnt;
}

int n;
int f[maxn], g[maxn], ans[maxn], l[maxn], r[maxn];
void dfs(int u, int fa = 0) {
	f[u] = 1;
	for (int i = head[u]; i; i = e[i].nxt) {
		int v = e[i].to;
		if (v != fa) {
			dfs(v, u);
			f[u] = static_cast<long long> (1 + f[v]) * f[u] % mod;
		}
	}
}
void dfs1(int u, int fa = 0) {
	std::vector<int> S;
	for (int i = head[u]; i; i = e[i].nxt) {
		int v = e[i].to;
		if (v != fa) S.push_back(v);
	}
	if (S.size()) {
		l[*S.begin()] = g[u];
		for (std::vector<int>::iterator it = S.begin() + 1; it != S.end(); it++) {
			l[*it] = static_cast<long long> (l[*(it - 1)]) * (f[*(it - 1)] + 1) % mod;
		}
		r[*(S.end() - 1)] = 1;
		if (S.begin() + 1 != S.end()) {
			for (std::vector<int>::iterator it = S.end() - 2; true; it--) {
				r[*it] = static_cast<long long> (r[*(it + 1)]) * (f[*(it + 1)] + 1) % mod;
				if (it == S.begin()) break;
			}
		}
	}
	for (int i = head[u]; i; i = e[i].nxt) {
		int v = e[i].to;
		if (v != fa) {
			g[v] = (static_cast<long long> (l[v]) * r[v] + 1) % mod; 
			ans[v] = static_cast<long long> (g[v]) * f[v] % mod;
			dfs1(v, u);
		}
	}
}

int main() {
	n = read();
	for (int i = 1, x; i < n; i++) {
		x = read();
		add(i + 1, x);
		add(x, i + 1);
	}
	dfs(1);
	ans[1] = f[1]; g[1] = 1;
	dfs1(1);
	for (int i = 1; i <= n; i++) {
		printf("%d", ans[i]);
		putchar(i == n ? ‘\n‘ : ‘ ‘);
	}
	return 0;
}