点击查看代码
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define endl "\n"
#define ios ios::sync_with_stdio(0),cin.tie(0),cout.tie(0)
const int N = 2e6+10;
const int p = 998244353;
int f[N],inv[N],yz[N],dep[N];
int km(int a,int b)
{
int res=1;
while(b)
{
if(b&1) res=res*a%p;
a=a*a%p;
b>>=1;
}
return res;
}
int n,d[N],dp[N],mx=0;
vector<int> h[N];
bool vis[N];
void add(int a,int b)
{
h[a].push_back(b);
h[b].push_back(a);
}
void bfs()
{
queue<int> q;
q.push(1);
vis[1]=1,d[1]=0;
while(q.size())
{
int u=q.front();
q.pop();
for(auto v:h[u])
{
if(vis[v]) continue;
vis[v]=1;
d[v]=d[u]+1;//算出深度
dep[d[v]]++;//纪录每个深度点的个数
if(h[v].size()==1&&v!=1) yz[d[v]]++;//纪录叶子节点个数
mx=max(mx,d[v]);//纪录下最大深度
q.push(v);
}
}
}
void solve()
{
cin>>n;
for(int i=1;i<n;i++) {
int u,v;
cin>>u>>v;
add(u,v);
}
bfs();
dp[0]=1;
int ans=0;
for(int i=1;i<=mx;i++)//递推求答案过程,注意求逆元时取模
{
int a=dep[i],b=yz[i];
ans+=i%p*dp[i-1]%p*b%p*km(a,p-2)%p;
ans%=p;
dp[i]=dp[i-1]%p*(a-b)%p*km(a,p-2)%p;
dp[i]%=p;
}
cout<<ans<<endl;
}
signed main()
{
ios;
int t=1;
//cin>>t;
while(t--)
{
solve() ;
}
return 0;
}