#include<iostream>
#include<cstring>
using namespace std;
const int mod=998244353;
const int N=2e5+1000,M=2e5+100,B=510;
#define int long long
int f[N],finv[N];
inline int qpow(int a,int b){
    int ret=1;
    while(b){
        if(b&1) ret=ret*a%mod;
        a=a*a%mod;b>>=1;
    }return ret;
}
inline int comb(int n,int m){
    if(n<m||m<0) return 0;
    return f[n]*finv[m]%mod*finv[n-m]%mod;
}
int n,m;
int g[510][N],inv[N];
inline void init(){
    for(int m=0;m<=M;m+=B){
        int i=m/B;
        g[i][0]=1;
        g[i][1]=(m-1)*g[i][0]%mod;
        for(int n=2;n<=M;n++){
            g[i][n]=inv[n]*(g[i][n-1]*(n+m-2)%mod+g[i][n-2])%mod;
        }
    }
}
// inline int DP(int n,int m){
    // if(n-2*m<0) return 0;
    // return (f[n-m]*f[m]%mod*dp[m+1][n-2*m]%mod+mod)%mod;
// }
int s[4];
inline int pval(int x,int y){
    if(x<0) return 0;
    if(x%B==0)return g[x/B][y];
    int pre=x/B+1;
    s[0]=g[pre][y],s[1]=g[pre][y+1],s[2]=g[pre][y+2];
    int n=y+2;
    for(int m=pre*B-1;m>=x;m--){
        s[2]=(s[2]-s[1])%mod;
        s[1]=(s[1]-s[0])%mod;
        s[0]=(n*s[2]-(n+m-2)*s[1])%mod;
    }
    return s[0];
}
int x,y;
inline void solve(){
    cin>>n>>m;
    int prod=f[n-m]*f[m]%mod;
    x=m+1,y=n-2*m;
    cout<<(prod*pval(x,y)%mod+mod)%mod<<'\n';
}
int32_t main(){
    // freopen("qwq.in","r",stdin);freopen("qwq.out","w",stdout);
    f[0]=finv[0]=1;
    for(int i=1;i<=M;i++) f[i]=f[i-1]*i%mod;
    finv[M]=qpow(f[M],mod-2);
    for(int i=M-1;i>=1;i--) finv[i]=finv[i+1]*(i+1)%mod;
    for(int i=1;i<=M;i++) inv[i]=f[i-1]*finv[i]%mod;
    inv[M+1]=qpow(M+1,mod-2);
    cerr<<"???"<<endl;
    init();
    int T;ios::sync_with_stdio(false);cin.tie(0),cout.tie(0);
    cin>>T;
    while(T--) solve();
}