#include<bits/stdc++.h>
using namespace std;
const int N=2e5+5,mod=998244353;
inline int add(int a,int b){return a+b>=mod?a+b-mod:a+b;}
inline int mul(int a,int b){return 1ll*a*b%mod;}
inline int qpow(int a,int b){
    int c=1;
    for(;b;b>>=1,a=mul(a,a))if(b&1)c=mul(c,a);
    return c;
}
int fac[N],ifac[N];
int n,m;
int f[5005][5005],ans[5005][5005];
void dp(){
    f[0][0]=f[1][1]=1;
    for(int i=2;i<=5000;i++){
        for(int j=0;j<=i;j++){
            if(j)f[i][j]=add(f[i-1][j-1],mul(i-1,add(f[i-1][j-1],f[i-2][j-1])));
            f[i][j]=add(f[i][j],mul(i-1,add(f[i-1][j],f[i-2][j])));
        }
    }
    for(int i=0;i<=5000;i++){
        for(int j=0;j<<1<=i;j++){
            ans[i][j]=mul(f[i-j][j],fac[j]);
        }
    }
}
int solve(int n,int m){
    if(m<<1>n)return 0;
    int ans=0;
    for(int i=0;i<=n-m-m;i++){
        int val=1ll*fac[n-m]*fac[n-m-i]%mod*ifac[i]%mod*ifac[n-m-m-i]%mod;
        if(i&1)ans=add(ans,mod-val);
        else ans=add(ans,val);
    }
    return ans;
}
int main(){
    fac[0]=1;
    for(int i=1;i<=200000;i++)fac[i]=mul(fac[i-1],i);
    ifac[200000]=qpow(fac[200000],mod-2);
    for(int i=200000;i>=1;i--)ifac[i-1]=mul(ifac[i],i);
    int T;scanf("%d",&T);
    // while(T--)solve();
    dp();
    while(T--){
        int n,m;scanf("%d%d",&n,&m);
        if(n<=5000&&m<=5000)printf("%d\n",f[n][m]);
        else printf("%d\n",solve(n,m));
    }
    return 0;
}