結果
問題 |
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])