結果
| 問題 | 
                            No.3118 Increment or Multiply
                             | 
                    
| コンテスト | |
| ユーザー | 
                             | 
                    
| 提出日時 | 2025-04-20 16:59:25 | 
| 言語 | C++17  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 20 ms / 2,000 ms | 
| コード長 | 1,333 bytes | 
| コンパイル時間 | 2,108 ms | 
| コンパイル使用メモリ | 194,836 KB | 
| 実行使用メモリ | 7,848 KB | 
| 最終ジャッジ日時 | 2025-04-20 16:59:29 | 
| 合計ジャッジ時間 | 4,055 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge5 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 1 | 
| other | AC * 35 | 
ソースコード
#include <bits/stdc++.h>
using namespace std;
using int64 = long long;
const int64 MOD = 998244353;
int64 mod_mul(int64 a, int64 b, int64 mod = MOD) {
    return (__int128)a * b % mod;
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(NULL);
    int T;
    cin >> T;
    const int64 INV2 = (MOD + 1) / 2;
    while (T--) {
        int64 N, A;
        cin >> N >> A;
        int64 ans = 0;
        if (A == 1) {
            int64 nmod = N % MOD;
            ans = mod_mul(mod_mul(nmod, (nmod - 1 + MOD) % MOD), INV2);
        } else {
            int64 cur = N;
            while (cur > 0) {
                if (cur < A) {
                    int64 cmod = cur % MOD;
                    int64 term = mod_mul(mod_mul((cmod - 1 + MOD) % MOD, cmod), INV2);
                    ans = (ans + term) % MOD;
                    break;
                }
                int64 m = cur / A;
                int64 r = cur % A;
                int64 term1 = mod_mul(m % MOD, (r + 1) % MOD);
                ans = (ans + term1) % MOD;
                int64 D = cur - m;
                int64 D_mod = D % MOD;
                int64 term2 = mod_mul(mod_mul((D_mod - 1 + MOD) % MOD, D_mod), INV2);
                ans = (ans + term2) % MOD;
                cur = m;
            }
        }
        cout << ans << '\n';
    }
    return 0;
}