結果

問題 No.1276 3枚のカード
ユーザー ayaoniayaoni
提出日時 2020-10-31 01:25:38
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 923 ms / 2,000 ms
コード長 681 bytes
コンパイル時間 195 ms
コンパイル使用メモリ 81,920 KB
実行使用メモリ 77,296 KB
最終ジャッジ日時 2024-07-22 04:12:55
合計ジャッジ時間 30,295 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
52,224 KB
testcase_01 AC 41 ms
52,608 KB
testcase_02 AC 39 ms
52,096 KB
testcase_03 AC 40 ms
52,224 KB
testcase_04 AC 40 ms
52,736 KB
testcase_05 AC 40 ms
51,968 KB
testcase_06 AC 41 ms
52,224 KB
testcase_07 AC 40 ms
52,580 KB
testcase_08 AC 41 ms
52,480 KB
testcase_09 AC 40 ms
52,352 KB
testcase_10 AC 796 ms
76,416 KB
testcase_11 AC 416 ms
76,288 KB
testcase_12 AC 509 ms
76,800 KB
testcase_13 AC 528 ms
76,976 KB
testcase_14 AC 116 ms
76,544 KB
testcase_15 AC 439 ms
76,416 KB
testcase_16 AC 855 ms
76,800 KB
testcase_17 AC 498 ms
76,288 KB
testcase_18 AC 392 ms
76,544 KB
testcase_19 AC 594 ms
76,416 KB
testcase_20 AC 378 ms
76,544 KB
testcase_21 AC 388 ms
76,544 KB
testcase_22 AC 875 ms
76,544 KB
testcase_23 AC 113 ms
76,644 KB
testcase_24 AC 564 ms
76,416 KB
testcase_25 AC 240 ms
76,672 KB
testcase_26 AC 683 ms
76,624 KB
testcase_27 AC 780 ms
76,544 KB
testcase_28 AC 749 ms
76,868 KB
testcase_29 AC 461 ms
76,544 KB
testcase_30 AC 154 ms
76,416 KB
testcase_31 AC 110 ms
76,800 KB
testcase_32 AC 252 ms
76,288 KB
testcase_33 AC 87 ms
76,788 KB
testcase_34 AC 198 ms
76,672 KB
testcase_35 AC 190 ms
76,288 KB
testcase_36 AC 80 ms
74,624 KB
testcase_37 AC 183 ms
76,620 KB
testcase_38 AC 144 ms
76,500 KB
testcase_39 AC 259 ms
76,416 KB
testcase_40 AC 827 ms
76,288 KB
testcase_41 AC 696 ms
76,672 KB
testcase_42 AC 746 ms
76,416 KB
testcase_43 AC 662 ms
76,288 KB
testcase_44 AC 921 ms
76,544 KB
testcase_45 AC 837 ms
76,416 KB
testcase_46 AC 835 ms
76,672 KB
testcase_47 AC 697 ms
76,416 KB
testcase_48 AC 671 ms
76,416 KB
testcase_49 AC 750 ms
76,504 KB
testcase_50 AC 822 ms
76,416 KB
testcase_51 AC 715 ms
77,296 KB
testcase_52 AC 628 ms
76,416 KB
testcase_53 AC 728 ms
76,928 KB
testcase_54 AC 828 ms
76,800 KB
testcase_55 AC 654 ms
76,544 KB
testcase_56 AC 635 ms
76,288 KB
testcase_57 AC 794 ms
76,800 KB
testcase_58 AC 829 ms
76,416 KB
testcase_59 AC 728 ms
76,800 KB
testcase_60 AC 923 ms
76,672 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
def I(): return int(sys.stdin.readline().rstrip())


N = I()
mod = 10**9+7

if N == 3:
    print(0)
    exit()


def f(M):  # (1の約数の個数) + … + (Mの約数の個数)、つまり M//1 + … + M//M
    if M == 0:
        return 0
    res = 0
    for i in range(1,int(M**.5)+1):
        res += M//i
        res %= mod
    res *= 2
    res %= mod
    res -= int(M**.5)**2
    res %= mod
    return res


ans = 0
for x in range(2,int(N**.5)+1):
    y = N//(x+1)
    z = N//x
    ans += (x-1)*(N+1-x)*(z-y)
    ans %= mod
    ans -= (x-1)*(f(z)-f(y))
    ans %= mod

for b in range(1,y+1):
    ans += ((N//b)-1)*(N+1-(N//b)-f(b)+f(b-1))
    ans %= mod

print(ans)
0