結果

問題 No.644 G L C C D M
ユーザー merom686merom686
提出日時 2018-02-03 01:12:27
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 16 ms / 2,000 ms
コード長 1,146 bytes
コンパイル時間 1,521 ms
コンパイル使用メモリ 79,008 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-30 02:35:13
合計ジャッジ時間 1,960 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 1 ms
4,380 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 2 ms
4,380 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 1 ms
4,376 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 3 ms
4,380 KB
testcase_14 AC 2 ms
4,376 KB
testcase_15 AC 1 ms
4,376 KB
testcase_16 AC 2 ms
4,380 KB
testcase_17 AC 1 ms
4,376 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 2 ms
4,380 KB
testcase_20 AC 1 ms
4,376 KB
testcase_21 AC 2 ms
4,376 KB
testcase_22 AC 16 ms
4,380 KB
testcase_23 AC 9 ms
4,376 KB
testcase_24 AC 6 ms
4,380 KB
testcase_25 AC 2 ms
4,380 KB
testcase_26 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

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