#include using namespace std; typedef long long ll; const ll MOD = 998244353; // modint struct mint { ll v; mint(): v(0) {} mint(ll _v) { v = (_v % MOD + MOD) % MOD; } mint& operator+=(const mint& o) { v += o.v; if (v >= MOD) v -= MOD; return *this; } mint& operator-=(const mint& o) { v -= o.v; if (v < 0) v += MOD; return *this; } mint& operator*=(const mint& o) { v = (ll)((__int128)v * o.v % MOD); return *this; } friend mint operator+(mint a, const mint& b) { return a += b; } friend mint operator-(mint a, const mint& b) { return a -= b; } friend mint operator*(mint a, const mint& b) { return a *= b; } }; int main() { ios::sync_with_stdio(false); cin.tie(NULL); int T; cin >> T; while (T--) { ll N, A; cin >> N >> A; mint ans = 0; if (A == 1) { // f(K) = N - K, sum = N*(N-1)/2 mint mN = N; mint inv2 = (MOD + 1) / 2; ans = mN * mint(N - 1) * inv2; } else { // base-A digits of N (least significant first) vector digits; ll tmp = N; while (tmp > 0) { digits.push_back(tmp % A); tmp /= A; } int D = digits.size(); // powers of A vector pw(D + 1); pw[0] = 1; for (int i = 1; i <= D; ++i) { pw[i] = pw[i - 1] * A; } // prefix sums of digits: S[t] = sum_{j=0..t-1} digits[j] vector S(D + 1); S[0] = 0; for (int i = 1; i <= D; ++i) { S[i] = S[i - 1] + digits[i - 1]; } mint inv2 = (MOD + 1) / 2; // iterate over possible division counts t for (int t = 0; t <= D; ++t) { if (pw[t] > N) break; ll R = N / pw[t]; ll L = (t + 1 <= D ? N / pw[t + 1] : 0); ll cnt = R - L; if (cnt <= 0) continue; mint mcnt = cnt; // sum f = cnt*(t + S[t]) + cnt*(cnt-1)/2 mint base = mcnt * mint(t + S[t]); mint tail = mint(cnt) * mint(cnt - 1) * inv2; ans += base + tail; } } cout << ans.v; if (T) cout << '\n'; } return 0; }