結果

問題 No.106 素数が嫌い！2
コンテスト
ユーザー @abcde
提出日時 2019-06-16 18:23:50
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,487 ms / 5,000 ms
コード長 1,632 bytes
コンパイル時間 1,879 ms
コンパイル使用メモリ 174,868 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-24 12:07:29
合計ジャッジ時間 9,321 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 13
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

// Efficient program to print all prime factors of a given number
// https://www.geeksforgeeks.org/print-all-prime-factors-of-a-given-number/
// 与えられた正整数についての素因数分解を計算.
// @param X: 素因数分解を行う整数.
// @return: 素因数分解 の 結果 を 返却.
map<int, int> div(int X){
    
    // 1. X を 2で割り切れなくなるまで割っていく.
    map<int, int> ret;
    while(X % 2 == 0) ret[2]++, X >>= 1;
    
    // 2. X を 3以上の奇数で, 割り切れなくなるまで順次割っていく.
    for(int i = 3; i <= sqrt(X); i += 2){
        while(X % i == 0){
            ret[i]++;
            X /= i;
        }
    }
    
    // 3. X が 2 より 大きな素数であれば, 追加.
    if(X > 2) ret[X]++;
    
    // 4. 出力.
    return ret;
}

int main() {
    
    // 1. 入力情報取得.
    int N, K;
    scanf("%d %d", &N, &K);
    
    // 2. K が 8以上ならば, ゼロを返却.
    // 2 * 3 * 5 * 7 * 11 * 13 * 17 * 19 = 9699690
    if(K >= 8){
        printf("%d\n", 0);
        return 0;
    }
    
    // 4. 2 ～ N を それぞれ素因数分解.
    int ans = 0;
    for(int i = 2; i <= N; i++){
        map<int, int> divisors = div(i);
        // K個以上 と 指定されていることに注意(※ロジック修正).
        // if(divisors.size() == K) ans++;
        if(divisors.size() >= K) ans++;
    }
    
    // 5. 出力.
    // ex.
    // [入力例]
    // 2000000 7
    // -> 58 で OK?.
    // 2000000 8
    // -> 0 で OK?.
    printf("%d\n", ans);
    return 0;
    
}
0