我们可以发现，只要一个原子在某一秒被触发，那么后面都会一直被触发。因为我们只在 \(m\) 个原子上投放中子，所以可以先把不可能被触发的原子从中删去，并删去它与其它原子相连的边。因为一个原子只要被触发，不管被多少中子激发，都只会发送恒定的一个中子。所以在一定时间后每秒释放的能量是恒定的。所以我们可以用一个数组记录这个原子每秒释放的能量，在它被访问到的时候，也就是每秒释放的能量变化的时候统计这个原子从上一次统计到现在变化之前，释放的总能量。在搜索完后，对每个原子加上现在到未来第 \(k\) 秒所释放的能量。若在第 \(k\) 秒时，还没有到达恒定值，即还没有搜索完，就直接停止搜索，然后输出答案即可。
接下来是代码环节：
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int MOD = 998244353;
int n, k, m, a[500005], w[500005], num[500005], cnt, ans[500005], ad[500005], la[500005];
vector<pair<int, int>> v[500005];
bool mark[500005], flag[500005];
signed main() {ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);cin >> n >> k >> m;for (int i = 1; i <= m; i ++) {cin >> w[i];mark[w[i]] = true;}for (int i = 1; i <= n; i ++) {cin >> a[i];for (int j = 1; j <= a[i]; j ++) {int x;cin >> x;v[i].push_back({x, ++cnt});num[x] ++;}}queue<int> q;for (int i = 1; i <= n; i ++) {if (num[i] == 0 && !mark[i]) {q.push(i);}}while (!q.empty()) {int u = q.front();q.pop();for (auto i : v[u]) {num[i.first] --;if (num[i.first] == 0 && !mark[i.first]) {q.push(i.first);}flag[i.second] = true;}}queue<pair<int, int>> qq;for (int i = 1; i <= m; i ++) {qq.push({w[i], 1});}while (!qq.empty()) {int u = qq.front().first, tim = qq.front().second;if (tim > k) {break;}if (ad[u] == 0) {ad[u] = a[u];} else {ans[u] += (tim - la[u]) * ad[u] % MOD;}ad[u] ++;la[u] = tim;qq.pop();for (auto i : v[u]) {if (!flag[i.second]) {flag[i.second] = true;qq.push({i.first, tim + 1});}}}for (int i = 1; i <= n; i ++) {ans[i] += (k - la[i] + 1) % MOD * ad[i] % MOD;cout << ans[i] << " ";}return 0;
}