#include<iostream>
#include<vector>
#include<algorithm>

using namespace std;
using ll = long long;

#include<atcoder/modint>
using mint = atcoder::modint998244353;
const int B = 2e6;
mint fac[B+1],ifac[B+1];
#include<map>
int main(){
    cin.tie(nullptr);
    ios::sync_with_stdio(false);

    int q;
    cin>>q;
    fac[0] = 1;
    for(int i = 1;i<=B;i++) fac[i] = fac[i-1] * i;
    ifac[B] = fac[B].inv();
    for(int i = B-1;i>=0;i--) ifac[i] = ifac[i+1] * (i+1);
    const int want = 1e5 + 10;
    vector<int> p(want+1,1);
    for(int i = 2;i<=want;i++){
        if(p[i]!=1) continue;
        for(int j = i;j<=want;j+=i) p[j] = i;
    }
    ll all = 0;
    for(int i = 0;i<q;i++){
        ll a,b;
        cin>>a>>b;
        while(a!=1){
            all++;
            a /= p[a];
        }
        mint ans = 0;
        if(all>=b){
            ll use = all - b;
            ans = fac[use+b-1] * ifac[use] * ifac[b-1];
        }        
        cout<<ans.val()<<endl;
    }
}