#include <bits/stdc++.h>

using namespace std;

#define deg( x ) cout<<""#x"="<<x<<endl
#define int long long
#define endl '\n'

const int N=1e5+10;
const int mod=998244353;

vector<int> G[ N ];

int sz[ N ];
int sum[ N ];
int faa[ N ];
int hh[ N ];

void init( int n ){
    for( int i=1; i<=n; i++ ){
        G[ i ].clear();
        hh[ i ]=faa[ i ]=sz[ i ]=sum[ i ]=0;
    }
}

void dfs( int rt , int fa ){
    sz[ rt ]=1;
    faa[ rt ]=fa;
    int re=0;

    for( auto son:G[ rt ] ){
        if( son==fa ) continue;
        dfs( son , rt );

        hh[ rt ]+=sz[ son ]*sz[ son ];
        hh[ rt ]%=mod;
        sz[ rt ]+=sz[ son ];
        sz[ rt ]%=mod;
        re+=sum[ son ];
        re%=mod;
    }

    sum[ rt ]+=re+sz[ rt ]*sz[ rt ];
    sum[ rt ]%=mod;
}

int exgcd( int a , int b , int &x, int &y ) //扩展欧几里得求同余方程解； x,y及为 x*a+y*b=1 的一组可行解；
{
    if( b==0 ){
        x=1,y=0;
        return a;
    }

    int com=exgcd( b , a%b , x , y );

    int _x=x , _y=y ;
    y=_x-(a/b)*_y; x=_y;

    return com;
}

int getinv( int x , int m ) // exgcd获取x在模m情况下的逆元
{
    int ans,y;
    exgcd( x , m , ans , y );
    ans=( ( ans%m )+m )%m;
    return ans;
}


void solve(){
    int n;
    cin>>n;

    init( n );

    for( int i=1,u,v; i<n; i++ ){
        cin>>u>>v;
        G[ u ].emplace_back( v );
        G[ v ].emplace_back( u );
    }

    dfs( 1 , 0 );

    int ans=0;
    for( int i=1; i<=n; i++ ){
        int cnt=0;

        int rt=i;
        for( auto son:G[ rt ] ){
            if( son==faa[ rt ] ) continue;

            cnt=sz[ son ]*sz[ son ]; cnt%=mod;
            cnt*=( n-sz[ son ] ); cnt%=mod;
            cnt*=( n-sz[ son ] ); cnt%=mod;
            ans+=cnt; ans%=mod;

            cnt=0;
            cnt=2*( n-sz[ son ] )*( n-sz[ son ] );
            cnt%=mod;

            int e=( sum[ son ]+mod*mod-sz[ son ]*sz[ son ] );
            e%=mod;
            cnt*=e; cnt%=mod;
            ans+=cnt; ans%=mod;

            cnt=sz[ son ]*sz[ son ]; cnt%=mod;

            e=( sum[ rt ]-sum[ son ]-sz[ rt ]*sz[ rt ]+mod*mod+mod )%mod;
            cnt*=2*e; cnt%=mod;
            ans+=cnt; ans%=mod;
        }
    }

    int jian=0;

    for( int i=1; i<=n; i++ ){
        int rt=i;
        for( auto son:G[ rt ] ){
            if( son==faa[ rt ] ) continue;
            int cnt=sz[ son ]*sz[ son ];
            int e=hh[ rt ]-cnt+mod;
            e%=mod;
            jian+=cnt*e;
            jian%=mod;
        }
    }

    ans=ans-jian+mod;
    ans%=mod;

//    for( int i=1; i<=n; i++ ){
//        cout<<sz[ i ]<<" ";
//    }
//    cout<<endl;
//
//    for( int i=1; i<=n; i++ ){
//        cout<<sum[ i ]<<" ";
//    }
//    cout<<endl;

//    ans*=2; ans%=mod;

//    deg( ans );
    int inv=( ( n )*( n-1 ) )/2;
    inv%=mod;

    inv*=inv;
    inv%=mod;

    inv=getinv( inv , mod );

    ans*=inv;
    ans%=mod;

    cout<<ans<<endl;

//    ans=918384806;
//    ans*=100;
//    ans%=mod;
//    deg( ans );
}

signed main() {
    ios::sync_with_stdio( 0 );
    cin.tie( 0 ); cout.tie( 0 );

    int T=1;
    cin>>T;

    while( T-- ){
        solve();

    }

    return 0;
}