結果

問題 No.2127 Mod, Sum, Sum, Mod
ユーザー LyricalMaestroLyricalMaestro
提出日時 2024-02-25 02:08:40
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 2,767 bytes
コンパイル時間 323 ms
コンパイル使用メモリ 82,480 KB
実行使用メモリ 79,680 KB
最終ジャッジ日時 2024-09-29 10:37:01
合計ジャッジ時間 5,506 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 33 ms
59,912 KB
testcase_01 AC 34 ms
53,864 KB
testcase_02 AC 123 ms
74,032 KB
testcase_03 AC 39 ms
53,580 KB
testcase_04 AC 39 ms
52,656 KB
testcase_05 AC 38 ms
54,088 KB
testcase_06 AC 34 ms
53,124 KB
testcase_07 AC 32 ms
53,524 KB
testcase_08 AC 34 ms
52,896 KB
testcase_09 AC 36 ms
54,300 KB
testcase_10 AC 36 ms
54,032 KB
testcase_11 AC 34 ms
53,068 KB
testcase_12 AC 34 ms
54,004 KB
testcase_13 AC 35 ms
54,340 KB
testcase_14 AC 147 ms
72,876 KB
testcase_15 AC 69 ms
72,864 KB
testcase_16 AC 35 ms
53,452 KB
testcase_17 AC 34 ms
52,884 KB
testcase_18 AC 77 ms
74,304 KB
testcase_19 AC 35 ms
53,336 KB
testcase_20 TLE -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
権限があれば一括ダウンロードができます

ソースコード

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_pow6 = pow(6, 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 = M // 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