結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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])
0