結果

問題 No.2127 Mod, Sum, Sum, Mod
ユーザー LyricalMaestroLyricalMaestro
提出日時 2024-02-25 02:10:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 139 ms / 2,000 ms
コード長 2,731 bytes
コンパイル時間 352 ms
コンパイル使用メモリ 82,040 KB
実行使用メモリ 75,152 KB
最終ジャッジ日時 2024-09-29 10:37:04
合計ジャッジ時間 2,781 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
53,912 KB
testcase_01 AC 35 ms
53,448 KB
testcase_02 AC 123 ms
73,612 KB
testcase_03 AC 38 ms
52,528 KB
testcase_04 AC 36 ms
52,920 KB
testcase_05 AC 36 ms
53,080 KB
testcase_06 AC 39 ms
53,972 KB
testcase_07 AC 35 ms
52,744 KB
testcase_08 AC 39 ms
53,468 KB
testcase_09 AC 38 ms
53,344 KB
testcase_10 AC 36 ms
53,792 KB
testcase_11 AC 37 ms
53,196 KB
testcase_12 AC 39 ms
53,244 KB
testcase_13 AC 38 ms
53,548 KB
testcase_14 AC 38 ms
54,124 KB
testcase_15 AC 37 ms
52,856 KB
testcase_16 AC 38 ms
54,488 KB
testcase_17 AC 36 ms
52,584 KB
testcase_18 AC 35 ms
53,008 KB
testcase_19 AC 37 ms
53,196 KB
testcase_20 AC 36 ms
53,408 KB
testcase_21 AC 40 ms
53,768 KB
testcase_22 AC 139 ms
74,536 KB
testcase_23 AC 138 ms
75,152 KB
testcase_24 AC 37 ms
53,632 KB
testcase_25 AC 35 ms
52,960 KB
testcase_26 AC 84 ms
74,512 KB
testcase_27 AC 37 ms
53,428 KB
testcase_28 AC 127 ms
74,420 KB
testcase_29 AC 36 ms
53,576 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

## https://yukicoder.me/problems/no/2127

import math

MOD = 998244353

def main():
    N, M = map(int, input().split())

    sqrt_n = int(math.sqrt(N))

    inv_pow2 = pow(2, MOD - 2, MOD)
    inv_pow4 = pow(4, MOD - 2, MOD)
    inv_pow12 = pow(12, MOD - 2, MOD)

    # j がsqrt(N)以下のケース
    answer1 = 0
    for j in range(1, min(M, sqrt_n) + 1):
        r = N % j
        q = N // j

        ans1 = ((j - 1) * j) % MOD
        ans1 *= inv_pow2
        ans1 %= MOD
        ans1 *= q
        ans1 %= MOD

        ans2 = ((r + 1) * r) % MOD
        ans2 *= inv_pow2
        ans2 %= MOD
        
        ans = (ans1 + ans2) % MOD
        answer1 += ans
        answer1 %= MOD
    
    # q　がsqrt(N)以下のケース
    answer2 = 0
    if M > sqrt_n:
        def sum1(n):
            ans = (n * (n + 1)) % MOD
            ans *= (n - 1)
            ans %= MOD
            ans *= pow(6, MOD - 2, MOD)
            ans %= MOD
            return ans
        
        def sum2(n, N, q):
            ans1 = (q * q) % MOD
            ans1 *= n
            ans1 %= MOD
            ans1 *= (n + 1)
            ans1 %= MOD
            ans1 *= (((2 * n) % MOD) + 1) % MOD
            ans1 %= MOD
            ans1 *= inv_pow12
            ans1 %= MOD

            ans2 = (2 * N + 1) % MOD
            ans2 *= q
            ans2 %= MOD
            ans2 *= n
            ans2 %= MOD
            ans2 *= (n + 1)
            ans2 %= MOD
            ans2 *= inv_pow4
            ans2 %= MOD

            ans3 = (N * (N + 1)) % MOD
            ans3 *= n
            ans3 %= MOD
            ans3 *= inv_pow2
            ans3 %= MOD

            answer = (ans1 - ans2) % MOD
            answer += ans3
            answer %= MOD

            return answer

        add_start = False
        sqrt_n2 = N // sqrt_n
        for q in range(sqrt_n2 + 2):
            if q == 0:
                lower = N
                if M <= lower:
                    continue
                add_start = True
                upper = M
            else:
                lower = max(N // (q + 1), sqrt_n)
                if not add_start:
                    if M <= lower:
                        continue
                    add_start = True
                    upper = M
                else:
                    upper = max(N // q, sqrt_n)

            ans = (sum1(upper) - sum1(lower)) % MOD
            ans *= q
            ans %= MOD


            ans2 = (sum2(upper, N, q) - sum2(lower, N, q)) % MOD
            ans += ans2
            ans %= MOD
            answer2 += ans
            answer2 %= MOD

    answer = (answer1 + answer2) % MOD
    print(answer)
            

        
    
























if __name__ == "__main__":
    main()
0