これは、所望のタイトルを必要とする総重量（通信ブロックの同じ点+通信ブロックポイントに同じ寄与しない全ての貢献）/全実施例（森の数）であります
最初の数を求め、パイルから形成された森の木の組成は、記載purfer配列、\（\ N-）標識されたツリーのためのポイント数（\} {^ N-2-N-）\、基準点次いでが考えられるEGF、EGF-ツリーが設定されているため、\（F（X）\） 、森林の次に生成機能\（E ^ {F（X ）} \）
次いで寄与は、同じ通信ポイント番号をブロックしないに等しい同一の通信ポイントを、ブロックしない\（* M ^ 2 \） 、およびなく、同じ通信ブロックと点の総数に等しい点数の\ （ - \）列挙通信ブロックの数を考慮すべき点のブロックサイズとの通信、の次にサイズが\（N- \）通信ブロックによって寄与されるN ^ {N-2（\ } \ binom { }} {2 N- \）通信ブロックを乗じ、その後も何度も表示され、残りの部分は森林であるので、そのライン上の森の前にこの混乱及びEGFボリューム生成機能
最後に、通信ブロックの同じ点の寄与は、正方形幾分困難なプロセスは、約変換考える\ [\開始{整列} \ sum_ {i = 1} ^ {N} \ sum_ {J = + 1} ^ {n}はDISを（i、j）は^ 2＆ = \ sum_ {i = 1} ^ {N} \ sum_ {J = + 1} ^ {N}（\和_ {（x、y）を\で（I、J）} 1 ）^ 2 \\＆= \ sum_ {i = 1} ^ {N} \ sum_ {J = + 1} ^ {N} \和_ {（X1、Y1）\で（I、J）} \ sum_ { （X2、Y2）が\（IN I、J）} 1 \端{整列} \] がパス上の2つの側面を列挙し、次に結合される。我々は、最初の2つに、列挙を列挙順序を交換しましたまた、2つのエッジを含むどのように多くの経路側、次いで考える。双方が同じである場合、それはそれ以外の場合は前者、列挙に三つのブロックの通信に分割され、2つのブロックにホスト通信部と通信します二つの通信ブロックサイズ\（sz_1、sz_2 \） 、および2つの経路は、出発寄与列挙\（sz_1sz_2の\） 、およびこれら2つの辺の交点と二つのブロックの通信、寄与列挙\を（sz_1sz_2 \） 、総寄与そう2 \）^（（sz_1sz_2）\は、次の3つのブロック通信寄与、すなわち、類似している。（（sz_1sz_2sz_3）^ 2 \）\、。上記のように、上部及びEGFにこのボリューム関数を生成する前に森。3つの通信ブロックの場合は、我々は、2つの側面がある列挙ことに注意してください 順序は、そう\（* 2 \） ;その後、時間の計算がします（（i、j）は\ \ ） の貢献と、\（\（J、I） ）の貢献が来ると考えられているが、その貢献のこの部分があることを\（/ 2 \）
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<cmath>
#include<ctime>
#include<queue>
#include<map>
#include<set>
#define LL long long
#define db double
using namespace std;
const int N=1e6+10,M=(1<<20)+10,mod=998244353,inv2=499122177;
LL rd()
{
LL x=0,w=1;char ch=0;
while(ch<'0'||ch>'9'){if(ch=='-') w=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
return x*w;
}
namespace ct2
{
int fpow(int a,int b){int an=1;while(b){if(b&1) an=1ll*an*a%mod;a=1ll*a*a%mod,b>>=1;} return an;}
int ginv(int a){return fpow(a,mod-2);}
int rdr[M],inv[M],p1[M],p2[M],p3[M],p4[M],p5[M];
void ntt(int *a,int n,bool op)
{
int l=0,x,y;
while((1<<l)<n) ++l;
for(int i=0;i<n;++i)
{
rdr[i]=(rdr[i>>1]>>1)|((i&1)<<(l-1));
if(i<rdr[i]) swap(a[i],a[rdr[i]]);
}
for(int i=1;i<n;i<<=1)
{
int ww=fpow(op?3:332748118,(mod-1)/(i<<1));
for(int j=0;j<n;j+=i<<1)
for(int k=0,w=1;k<i;++k,w=1ll*w*ww%mod)
x=a[j+k],y=1ll*a[j+k+i]*w%mod,a[j+k]=(x+y)%mod,a[j+k+i]=(x-y+mod)%mod;
}
if(!op) for(int i=0,w=ginv(n);i<n;++i) a[i]=1ll*a[i]*w%mod;
}
void polyder(int *aa,int *bb,int n)
{
for(int i=0;i<n-1;++i) bb[i]=1ll*aa[i+1]*(i+1)%mod;
bb[n-1]=bb[n]=0;
}
void polying(int *aa,int *bb,int n)
{
for(int i=1;i<n;++i) bb[i]=1ll*aa[i-1]*inv[i]%mod;
bb[0]=0;
}
void polyinv(int *aa,int *bb,int n)
{
if(n==1){bb[0]=ginv(aa[0]);return;}
polyinv(aa,bb,n>>1);
for(int i=0;i<n;++i) p1[i]=aa[i],p2[i]=bb[i];
ntt(p1,n<<1,1),ntt(p2,n<<1,1);
for(int i=0;i<n<<1;++i) p1[i]=1ll*p1[i]*p2[i]%mod*p2[i]%mod;
ntt(p1,n<<1,0);
for(int i=0;i<n;++i) bb[i]=((bb[i]+bb[i])%mod-p1[i]+mod)%mod;
for(int i=0;i<n<<1;++i) p1[i]=p2[i]=0;
}
void polyln(int *aa,int *bb,int n)
{
polyder(aa,p3,n),polyinv(aa,p4,n);
ntt(p3,n<<1,1),ntt(p4,n<<1,1);
for(int i=0;i<n<<1;++i) p3[i]=1ll*p3[i]*p4[i]%mod;
ntt(p3,n<<1,0);
polying(p3,bb,n);
for(int i=0;i<n<<1;++i) p3[i]=p4[i]=0;
}
void polyexp(int *aa,int *bb,int n)
{
if(n==1){bb[0]=1;return;}
polyexp(aa,bb,n>>1);
polyln(bb,p5,n);
for(int i=0;i<n;++i) p1[i]=bb[i],p5[i]=(-p5[i]+aa[i]+mod)%mod;
p5[0]=(p5[0]+1)%mod;
ntt(p1,n<<1,1),ntt(p5,n<<1,1);
for(int i=0;i<n<<1;++i) p1[i]=1ll*p1[i]*p5[i]%mod;
ntt(p1,n<<1,0);
for(int i=0;i<n;++i) bb[i]=p1[i];
for(int i=0;i<n<<1;++i) p1[i]=p5[i]=0;
}
int fac[N],iac[N],aa[M],f[M],dotf[M],g[M],h[M],h2[M],h3[M],hh[M];
int C(int n,int m){return m<0||n<m?0:1ll*fac[n]*iac[m]%mod*iac[n-m]%mod;}
void wk()
{
int nn=200010,len=1<<18;
inv[0]=inv[1]=1;
for(int i=2;i<=nn;++i) inv[i]=(mod-1ll*(mod/i)*inv[mod%i]%mod)%mod;
fac[0]=1;
for(int i=1;i<=nn;++i) fac[i]=1ll*fac[i-1]*i%mod;
iac[nn]=ginv(fac[nn]);
for(int i=nn;i;--i) iac[i-1]=1ll*iac[i]*i%mod;
aa[1]=1;
for(int i=2;i<=nn;++i) aa[i]=1ll*fpow(i,i-2)*iac[i]%mod;
polyexp(aa,f,len);
f[0]=1;
len<<=1;
for(int i=0;i<=nn;++i) dotf[i]=f[i];
ntt(dotf,len,1);
for(int i=2;i<=nn;++i) g[i]=1ll*C(i,2)*fpow(i,i-2)%mod*iac[i]%mod;
ntt(g,len,1);
for(int i=0;i<len;++i) g[i]=1ll*g[i]*dotf[i]%mod;
ntt(g,len,0);
h[1]=1;
for(int i=2;i<=nn;++i) h[i]=1ll*i*i%mod*fpow(i,i-2)%mod*iac[i]%mod;
ntt(h,len,1);
for(int i=0;i<len;++i) h2[i]=1ll*h[i]*h[i]%mod;
ntt(h2,len,0);
for(int i=1;i<=nn;++i) hh[i]=1ll*h2[i]*inv2%mod;
for(int i=nn+1;i<len;++i) h2[i]=0;
ntt(h2,len,1);
for(int i=0;i<len;++i) h3[i]=1ll*h2[i]*h[i]%mod;
ntt(h3,len,0);
for(int i=1;i<=nn;++i) hh[i]=(hh[i]+h3[i])%mod;
ntt(hh,len,1);
for(int i=0;i<len;++i) hh[i]=1ll*hh[i]*dotf[i]%mod;
ntt(hh,len,0);
for(int i=1;i<=nn;++i)
{
f[i]=1ll*f[i]*fac[i]%mod;
g[i]=1ll*g[i]*fac[i]%mod;
hh[i]=1ll*hh[i]*fac[i]%mod;
}
int T=rd();
while(T--)
{
int n=rd(),m=rd();
m=1ll*m*m%mod;
int ans=(1ll*m*(1ll*f[n]*C(n,2)%mod-g[n]+mod)%mod+1ll*hh[n])%mod;
ans=1ll*ans*ginv(f[n])%mod;
printf("%d\n",ans);
}
}
}
int main()
{
ct2::wk();
return 0;
}