結果
| 問題 |
No.644 G L C C D M
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-03-30 22:45:09 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 65 ms / 2,000 ms |
| コード長 | 429 bytes |
| コンパイル時間 | 289 ms |
| コンパイル使用メモリ | 82,184 KB |
| 実行使用メモリ | 68,352 KB |
| 最終ジャッジ日時 | 2024-11-15 14:23:29 |
| 合計ジャッジ時間 | 2,432 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 27 |
ソースコード
N,M = map(int,input().split())
P = 10 ** 9 + 7
import sys
if N < 2 * M:
print(0)
exit()
C = N + 5
fact = [1] * C
for i in range(2,C):
fact[i] = fact[i-1] * i % P
q = N // M
dat = [0] * (q + 1)
for i in range(1,q+1):
m = i * M
j = N // m
dat[i] = j * (j - 1) * fact[N-2] % P
for i in range(q,0,-1):
for j in range(2 * i,q + 1,i):
dat[i] -= dat[j]
dat[i] %= P
print(dat[1])