結果

問題 No.2127 Mod, Sum, Sum, Mod
ユーザー LyricalMaestroLyricalMaestro
提出日時 2024-02-25 02:10:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 172 ms / 2,000 ms
コード長 2,731 bytes
コンパイル時間 257 ms
コンパイル使用メモリ 81,700 KB
実行使用メモリ 73,448 KB
最終ジャッジ日時 2024-02-25 02:10:08
合計ジャッジ時間 3,652 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
53,588 KB
testcase_01 AC 39 ms
53,588 KB
testcase_02 AC 145 ms
73,448 KB
testcase_03 AC 41 ms
53,588 KB
testcase_04 AC 40 ms
53,588 KB
testcase_05 AC 42 ms
53,588 KB
testcase_06 AC 40 ms
53,588 KB
testcase_07 AC 42 ms
53,588 KB
testcase_08 AC 42 ms
53,588 KB
testcase_09 AC 42 ms
53,588 KB
testcase_10 AC 41 ms
53,588 KB
testcase_11 AC 41 ms
53,588 KB
testcase_12 AC 42 ms
53,588 KB
testcase_13 AC 42 ms
53,588 KB
testcase_14 AC 40 ms
53,588 KB
testcase_15 AC 40 ms
53,588 KB
testcase_16 AC 40 ms
53,588 KB
testcase_17 AC 41 ms
53,588 KB
testcase_18 AC 41 ms
53,584 KB
testcase_19 AC 42 ms
53,588 KB
testcase_20 AC 40 ms
53,588 KB
testcase_21 AC 38 ms
53,588 KB
testcase_22 AC 172 ms
73,448 KB
testcase_23 AC 140 ms
73,192 KB
testcase_24 AC 41 ms
53,588 KB
testcase_25 AC 40 ms
53,588 KB
testcase_26 AC 97 ms
73,448 KB
testcase_27 AC 41 ms
53,588 KB
testcase_28 AC 140 ms
73,448 KB
testcase_29 AC 42 ms
53,588 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