結果

問題 No.1878 union-find の数え上げ
ユーザー SSRS
提出日時 2022-03-18 21:55:04
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 1,287 bytes
コンパイル時間 2,351 ms
コンパイル使用メモリ 207,416 KB
最終ジャッジ日時 2025-01-28 10:20:46
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 5 WA * 7 RE * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
const long long MOD = 998244353;
int main(){
  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);
  }
  vector<int> bfs;
  vector<int> d(N);
  d[0] = 0;
  queue<int> Q;
  Q.push(0);
  while (!Q.empty()){
    int v = Q.front();
    Q.pop();
    bfs.push_back(v);
    for (int w : c[v]){
      d[w] = d[v] + 1;
      Q.push(w);
    }
  }
  reverse(bfs.begin(), bfs.end());
  vector<vector<long long>> dp(N);
  long long sum = 0;
  for (int v : bfs){
    if (v != 0){
      sum += d[v];
      assert(sum <= 10000000);
      dp[v] = vector<long long>(d[v], 0);
      int cnt = c[v].size();
      vector<long long> P(cnt);
      for (int i = 0; i < cnt; i++){
        P[i] = dp[c[v][i]][d[v]];
      }
      for (int i = d[v] - 1; i >= 0; i--){
        for (int j = 0; j < cnt; j++){
          P[j] += dp[c[v][j]][i];
          P[j] %= MOD;
        }
        dp[v][i] = 1;
        for (int j = 0; j < cnt; j++){
          dp[v][i] *= P[j];
          dp[v][i] %= MOD;
        }
      }
    }
  }
  long long ans = 1;
  for (int v : c[0]){
    ans *= dp[v][0];
    ans %= MOD;
  }
  cout << ans << endl;
}
0