結果

問題 No.847 Divisors of Power
ユーザー phspls
提出日時 2020-04-04 23:42:34
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 15 ms / 2,000 ms
コード長 1,119 bytes
コンパイル時間 1,700 ms
コンパイル使用メモリ 183,248 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-03 07:49:27
合計ジャッジ時間 2,564 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define llong long long

llong multi(llong base, int exp) {
    if(exp == 0) return 1;
    if(exp == 1) return base;
    return multi(base, exp/2) * multi(base, exp/2) * multi(base, exp%2);
}

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

    map<int, llong> primes;
    for(int i=2; i<=(int)sqrt(n); i++) {
        while(n % i == 0) {
            n /= i;
            primes[i] += 1;
        }
        if(n == 1) break;
    }
    if(n > 1) primes[n] = 1;

    for(pair<int, llong> p: primes) {
        primes[p.first] *= min(k, 30);
    }

    vector<int> result;
    result.push_back(1);
    for(pair<int, int> p: primes) {
        int limit = result.size();
        for(int i=1; i<=p.second; i++) {
            llong val = multi(p.first, i);
            if(val > m) break;
            rep(j, limit) {
                if(result[j] * val > m) break;
                result.push_back(result[j] * val);
            }
        }
        sort(result.begin(), result.end());
    }
    cout << result.size() << endl;
}
0