結果

問題 No.644 G L C C D M
ユーザー merom686merom686
提出日時 2018-02-03 01:12:27
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 17 ms / 2,000 ms
コード長 1,146 bytes
コンパイル時間 969 ms
コンパイル使用メモリ 80,244 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-31 08:26:38
合計ジャッジ時間 1,902 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

constexpr int P = 1000000007;

vector<int> prime_vector(int n) {
    n = max(n, 6);
    int m = (int)sqrt(n - 1);

    vector<bool> q(n);
    for (int i = 2; i <= m; i++) {
        if (q[i]) continue;
        for (int j = i * i; j < n; j += i) {
            q[j] = true;
        }
    }

    vector<int> p; p.reserve(n / log(n * 0.329));
    for (int i = 2; i < n; i++) {
        if (!q[i]) p.push_back(i);
    }
    return p;
}

int main() {
    int n, m;
    cin >> n >> m;

    int l = n / m;

    auto pv = prime_vector((int)sqrt(l) + 1);

    int64_t r = 0;
    for (int i = 2; i <= l; i++) {
        int k = i, t = i;
        for (int p : pv) {
            if (t < p * p) break;
            if (t % p == 0) {
                while (t % p == 0) t /= p;
                k -= k / p;
            }
        }
        if (t > 1) k -= k / t;
        r += k;
    }

    r *= 2;
    for (int i = 1; i <= n - 2; i++) {
        r *= i;
        r %= P;
    }

    cout << r << endl;

    return 0;
}
0