結果

問題 No.2127 Mod, Sum, Sum, Mod
コンテスト
ユーザー LyricalMaestro
提出日時 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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17 TLE * 1 -- * 9
権限があれば一括ダウンロードができます

ソースコード

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