結果

問題 No.847 Divisors of Power
ユーザー wunderkammer2
提出日時 2020-01-23 13:00:59
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 12 ms / 2,000 ms
コード長 1,559 bytes
コンパイル時間 896 ms
コンパイル使用メモリ 104,488 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 04:27:32
合計ジャッジ時間 1,869 ms
ジャッジサーバーID
(参考情報)
judge1 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<algorithm>
#include<cmath>
#include<cstdio>
#include<functional>
#include<iomanip>
#include<iostream>
#include<map>
#include<numeric>
#include<queue>
#include<set>
#include<string>
#include<utility>
#include<vector>

using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, ll> P;
const ll MOD = 1000000007;
#define rep(i,n) for(int i=0;i<n;i++)
#define repl(i,s,e) for(int i=s;i<e;i++)
#define reple(i,s,e) for(int i=s;i<=e;i++)
#define revrep(i,n) for(int i=n-1;i>=0;i--)
#define all(x) (x).begin(),(x).end()

int N, K, M;

map<int, ll> primeNumberFactor(int n)
{
    map<int, ll> res;

    for (int i = 2; i * i <= n; i++)
    {
        while (n % i == 0)
        {
            ++res[i];
            n /= i;
        }
    }
    if (n != 1) res[n] = 1;
    return res;
}

ll countFactors(vector<P> list, ll num, int idx)
{
    //リストを最後まで巡回してもubを越えなかった
    if (idx >= list.size()) return 1;

    auto x = list[idx];

    ll count = 0;

    reple(i, 0, x.second)
    {
        //numは単調増加するため、一度でもubを超えたら終了
        if (num > M) break;

        //次の素因数へ再帰
        count += countFactors(list, num, idx + 1);

        num *= x.first;
    }

    return count;
}

int main()
{
    cin >> N >> K >> M;

    auto primes = primeNumberFactor(N);

    vector<P> counts;

    for (auto p : primes)
    {
        counts.emplace_back(p.first, p.second * K);
    }

    cout << countFactors(counts, 1, 0) << endl;

    return 0;
}
0