結果
問題 | No.644 G L C C D M |
ユーザー |
![]() |
提出日時 | 2018-02-02 22:02:47 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 97 ms / 2,000 ms |
コード長 | 965 bytes |
コンパイル時間 | 1,041 ms |
コンパイル使用メモリ | 83,140 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-31 07:34:51 |
合計ジャッジ時間 | 2,194 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 27 |
ソースコード
#include <iostream> #include <algorithm> #include <string> #include <vector> #include <cmath> #include <queue> #include <functional> #include <map> using namespace std; long long n, m, mod = 1000000007; bool primes[1000009]; long long euler_phi(long long p) { long long r = p; for (int i = 2; i * i <= p; i++) { if (p % i == 0) r = 1LL * r * (i - 1) / i; while (p % i == 0) p /= i; } if (p != 1) r = 1LL * r * (p - 1) / p; return r; } long long prime(long long p) { if (p == 0) return 0; long long R = 0; for (int i = 1; i <= p; i++) { R += euler_phi(i); } R--; R *= 2; return R; } int main() { for (int i = 2; i <= 1000008; i++) primes[i] = true; for (int i = 2; i*i <= 1000008; i++) { for (int j = i*i; j <= 1000008; j += i) primes[j] = false; } cin >> n >> m; long long G = prime(n / m); long long ret = 1; for (int i = 1; i <= n - 2; i++) { ret *= i; ret %= mod; } ret *= (G % mod); ret %= mod; cout << ret << endl; return 0; }