結果
問題 | No.1276 3枚のカード |
ユーザー |
![]() |
提出日時 | 2020-10-30 23:46:07 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 940 ms / 2,000 ms |
コード長 | 984 bytes |
コンパイル時間 | 412 ms |
コンパイル使用メモリ | 82,204 KB |
実行使用メモリ | 76,928 KB |
最終ジャッジ日時 | 2024-07-22 03:28:06 |
合計ジャッジ時間 | 30,952 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 61 |
ソースコード
import syssys.setrecursionlimit(10**7)def I(): return int(sys.stdin.readline().rstrip())def MI(): return map(int,sys.stdin.readline().rstrip().split())def LI(): return list(map(int,sys.stdin.readline().rstrip().split()))def LI2(): return list(map(int,sys.stdin.readline().rstrip()))def S(): return sys.stdin.readline().rstrip()def LS(): return list(sys.stdin.readline().rstrip().split())def LS2(): return list(sys.stdin.readline().rstrip())N = I()mod = 10**9+7if N == 3:print(0)exit()def f(M):if M == 0:return 0res = 0for i in range(1,int(M**.5)+1):res += M//ires %= modres *= 2res %= modres -= int(M**.5)**2res %= modreturn resans = 0for x in range(2,int(N**.5)+1):y = N//(x+1)z = N//xans += (x-1)*(N+1-x)*(z-y)ans %= modans -= (x-1)*(f(z)-f(y))ans %= modfor b in range(1,y+1):ans += ((N//b)-1)*(N+1-(N//b)-f(b)+f(b-1))ans %= modprint(ans)