結果
| 問題 |
No.644 G L C C D M
|
| コンテスト | |
| ユーザー |
square1001
|
| 提出日時 | 2017-08-16 21:02:58 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 578 bytes |
| コンパイル時間 | 1,638 ms |
| コンパイル使用メモリ | 166,220 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2024-10-13 13:01:18 |
| 合計ジャッジ時間 | 2,729 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 14 WA * 11 RE * 2 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
int N, M, mobius[100009], mod = 1000000007;
long long solve(long long x) {
fill(mobius, mobius + x + 1, 1);
for(int i = 2; i <= x; i++) {
for(int j = i; j <= x; j += i) mobius[j] *= -1;
if(i * i <= x) for(int j = i * i; j <= x; j += i * i) mobius[j] = 0;
}
long long ret = 0;
for(int i = 1; i <= x; i++) ret += 1LL * mobius[i] * (x / i) * (x / i);
return ret;
}
int main() {
cin >> N >> M;
int ret = (solve(N / M) - 1) % mod;
for(int i = 1; i <= N - 2; i++) ret = 1LL * ret * i % mod;
cout << ret << endl;
return 0;
}
square1001