結果

問題 No.1098 LCAs
ユーザー SSRS
提出日時 2020-05-22 20:24:43
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,485 bytes
コンパイル時間 1,884 ms
コンパイル使用メモリ 179,024 KB
実行使用メモリ 77,616 KB
最終ジャッジ日時 2024-11-22 16:20:04
合計ジャッジ時間 13,209 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 2
other AC * 2 WA * 8 RE * 17 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
int LOG = 19;
struct lowest_common_ancestor{
	int N;
	vector<int> d;
	vector<vector<int>> p;
	lowest_common_ancestor(vector<int> &P, vector<vector<int>> &C){
		N = P.size();
		d = vector<int>(N, 0);
		queue<int> Q;
		Q.push(0);
		while (!Q.empty()){
			int v = Q.front();
			Q.pop();
			for (int w : C[v]){
				d[w] = d[v] + 1;
				Q.push(w);
			}
		}
		p = vector<vector<int>>(N, vector<int>(LOG, -1));
		for (int i = 0; i < N; i++){
			p[i][0] = P[i];
		}
		for (int i = 1; i < LOG; i++){
			for (int j = 0; j < N; j++){
				if (p[j][i - 1] != -1){
					p[j][i] = p[p[j][i - 1]][i - 1];
				}
			}
		}
	}
	int query(int u, int v){
		if (d[u] > d[v]){
			swap(u, v);
		}
		for (int k = 0; k < LOG; k++){
			if ((d[v] - d[u]) >> k & 1){
				v = p[v][k];
			}
		}
		if (u == v){
		    return u;
		}
		for (int k = LOG - 1; k >= 0; k--){
			if (p[u][k] != p[v][k]){
				u = p[u][k];
				v = p[v][k];
				assert(u != -1);
				assert(v != -1);
			}
		}
		return p[u][0];
	}
};
int main(){
  //想定TLE O(N^2logN)
  int N;
  cin >> N;
  vector<int> p(N, -1);
  for (int i = 1; i < N; i++){
    cin >> p[i];
    p[i]--;
  }
  vector<vector<int>> c(N);
  for (int i = 1; i < N; i++){
    c[p[i]].push_back(i);
  }
  lowest_common_ancestor T(p, c);
  vector<long long> ans(N, 0);
  for (int i = 0; i < N; i++){
    for (int j = 0; j < N; j++){
      ans[T.query(i, j)]++;
    }
  }
  for (int i = 0; i < N; i++){
    cout << ans[i] << endl;
  }
}
0