#include<bits/stdc++.h>
using namespace std;
constexpr int mod=998244353,N=1e5+5;
int n;vector<int> e[N];
int sz[N],cnt[N],val[N],tot[N];bool sb[N];
void dfs(int p,int f){
    sz[p]=1,cnt[p]=(e[p].size()==1);
    for(int x:e[p])if(x!=f) dfs(x,p),sz[p]+=sz[x],cnt[p]+=cnt[x];
    tot[min(sz[p],n-sz[p])]++;
    if(e[p].size()<=2) sb[n-sz[p]]=1;
    if(e[f].size()<=2) sb[sz[p]]=1;
    val[sz[p]]=cnt[p];
}
inline int qpow(int a,int b){
    int ans=1;
    while(b){
        if(b&1) ans=ans*1ll*a%mod;
        a=a*1ll*a%mod,b>>=1;
    }
    return ans;
}
int fac[N],inv[N];
inline int C(int n,int m){if(n<m||m<0)return 0;return fac[n]*1ll*inv[m]%mod*inv[n-m]%mod;}
int main(){
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    cin>>n;
    for(int i=1,u,v;i<n;i++) cin>>u>>v,e[u].emplace_back(v),e[v].emplace_back(u);
    for(int i=fac[0]=1;i<=n;i++) fac[i]=fac[i-1]*1ll*i%mod;
    inv[n]=qpow(fac[n],mod-2);
    for(int i=n;i>=1;i--) inv[i-1]=inv[i]*1ll*i%mod;
    dfs(1,0);
    long long ans=2;
    for(int i=1;i<=n/2;i++)if(tot[i]){
        bool f=(tot[i]>1||tot[i-1]||sb[i]||i*2==n);
        ans+=(f?C(n,i):C(n,i)+mod-C(n-cnt[1],n-i-(val[i]?:cnt[1]-val[n-i])))*(i*2==n?1:2);
    }
    cout<<ans%mod<<'\n';
}