#include<iostream>
#include<vector>
#include<algorithm>
#include<atcoder/modint>
using namespace std;
using namespace atcoder;
using ll=long long;
using mint=modint998244353;
#define all(v) v.begin(),v.end()
#define rall(v) v.rbegin(),v.rend()
template<class T> bool chmax(T &a, T b){if (a < b){a = b;return true;} else return false;}
template<class T> bool chmin(T &a, T b){if (a > b){a = b;return true;} else return false;}

const int MAX=2000010;
const ll MOD=998244353;
mint fac[MAX+1],finv[MAX+1],inv[MAX+1];

void set_fac(){
    fac[0]=fac[1]=1;
    finv[0]=finv[1]=1;
    inv[1]=1;
    for (int i=2;i<=MAX;i++){
        fac[i]=fac[i-1]*i;
        inv[i]=-inv[MOD%i]*(MOD/i);
        finv[i]=finv[i-1]*inv[i];
    }
}

mint cmb(int n,int r){
    if (r<0||n<r)return 0;
    return fac[n]*finv[r]*finv[n-r];
}

int main(){
    ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    int T;
    cin>>T;
    while(T--){
        ll N,M;
        cin>>N>>M;
        vector dp(M+1,vector<mint>(3));
        dp[0][0]=1;
        for(int i=0;i<M;i++){
            dp[i+1][0]+=dp[i][0]+dp[i][1]+dp[i][2];
            dp[i+1][1]+=dp[i][0]*N+dp[i][1]*(N-1)+dp[i][2]*(N-2);
            if(N>=4){
                dp[i+1][2]+=dp[i][0]*(N*(N-1)/2-N);
                dp[i+1][2]+=dp[i][1]*(N*(N-1)/2-N-(N-3));
                dp[i+1][2]+=dp[i][2]*(N*(N-1)/2-N-2*(N-3)+1);
            }
        }
        mint ans=dp[M][0]+dp[M][1]+dp[M][2];
        cout<<ans.val()<<"\n";
    }
}