結果
| 問題 |
No.2598 Kadomatsu on Tree
|
| ユーザー |
umezo
|
| 提出日時 | 2024-01-03 00:42:50 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 311 ms / 2,000 ms |
| コード長 | 1,406 bytes |
| コンパイル時間 | 2,526 ms |
| コンパイル使用メモリ | 197,640 KB |
| 最終ジャッジ日時 | 2025-02-18 15:59:21 |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 60 |
ソースコード
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define ALL(v) v.begin(),v.end()
typedef long long ll;
#include<bits/stdc++.h>
using namespace std;
const int MOD=998244353;
const ll c2=499122177;
const int MAX=200200;
vector<int> G[MAX];
ll dp[MAX],ans[MAX],A[MAX];
void dfs1(int v,int p){
dp[v]=1;
for(auto nv:G[v]){
if(nv==p) continue;
dfs1(nv,v);
dp[v]+=dp[nv];
}
}
void dfs2(int v,int p){
ll ls=0,ls2=0,rs=0,rs2=0;
for(auto nv:G[v]){
if(A[nv]<A[v]){
ls=(ls+dp[nv])%MOD;
ls2=(ls2+dp[nv]*dp[nv]%MOD)%MOD;
}
else if(A[nv]>A[v]){
rs=(rs+dp[nv])%MOD;
rs2=(rs2+dp[nv]*dp[nv]%MOD)%MOD;
}
ans[v]=((ls*ls%MOD-ls2+MOD)%MOD*c2%MOD+(rs*rs%MOD-rs2+MOD)%MOD*c2%MOD)%MOD;
}
int N=G[v].size();
vector<ll> L(N),R(N);
rep(i,N) L[i]=R[i]=dp[G[v][i]];
for(int i=1;i<N;i++) L[i]=(L[i]+L[i-1])%MOD;
for(int i=N-2;i>=0;i--) R[i]=(R[i]+R[i+1])%MOD;
rep(i,N){
if(G[v][i]==p) continue;
dp[v]=1;
if(i) dp[v]=(dp[v]+L[i-1])%MOD;
if(i<N-1) dp[v]=(dp[v]+R[i+1])%MOD;
dfs2(G[v][i],v);
}
}
int main(){
ios::sync_with_stdio(false);
std::cin.tie(nullptr);
int n;
cin>>n;
rep(i,n-1){
int u,v;
cin>>u>>v;
u--,v--;
G[u].push_back(v);
G[v].push_back(u);
}
rep(i,n) cin>>A[i];
dfs1(0,-1);
dfs2(0,-1);
ll sum=0;
rep(i,n) sum=(sum+ans[i])%MOD;
cout<<sum<<endl;
return 0;
}
umezo