結果
問題 | No.644 G L C C D M |
ユーザー |
|
提出日時 | 2018-02-02 22:12:59 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 5 ms / 2,000 ms |
コード長 | 1,241 bytes |
コンパイル時間 | 1,129 ms |
コンパイル使用メモリ | 106,064 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-31 07:39:31 |
合計ジャッジ時間 | 2,171 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 27 |
ソースコード
#define _USE_MATH_DEFINES#include <cstdio>#include <iostream>#include <sstream>#include <fstream>#include <iomanip>#include <algorithm>#include <cmath>#include <complex>#include <string>#include <vector>#include <list>#include <queue>#include <stack>#include <set>#include <map>#include <bitset>#include <numeric>#include <limits>#include <climits>#include <cfloat>#include <functional>#include <iterator>using namespace std;const int MOD = 1000000007;// オイラー関数φ(x):x以下の自然数で、xと互いに素なものの個数void EulerTotient(int n, vector<int>& phi){vector<bool> prime(n+1, true);phi.assign(n+1, 0);for(int i=1; i<=n; ++i)phi[i] = i;for(int i=2; i<=n; ++i){if(!prime[i])continue;for(int j=1; i*j<=n; ++j){prime[i*j] = false;phi[i*j] -= phi[i*j] / i;}}}int main(){int n, m;cin >> n >> m;vector<int> phi;EulerTotient(n/m, phi);long long ans = 0;for(int i=2; i<=n/m; ++i){ans += phi[i] * 2;ans %= MOD;}for(int i=1; i<=n-2; ++i){ans *= i;ans %= MOD;}cout << ans << endl;return 0;}