結果

問題 No.847 Divisors of Power
コンテスト
ユーザー ganmodokix
提出日時 2019-07-05 22:23:21
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 39 ms / 2,000 ms
コード長 3,153 bytes
コンパイル時間 2,775 ms
コンパイル使用メモリ 190,704 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-11-14 21:30:01
合計ジャッジ時間 3,835 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize ("O3")
#pragma GCC target ("tune=native")
#pragma GCC target ("avx")

#include <bits/stdc++.h>

// 汎用マクロ
#define ALL_OF(x) (x).begin(), (x).end()
#define REP(i,n) for (long long i=0, i##_len=(n); i<i##_len; i++)
#define RANGE(i,is,ie) for (long long i=(is), i##_end=(ie); i<=i##_end; i++)
#define DSRNG(i,is,ie) for (long long i=(is), i##_end=(ie); i>=i##_end; i--)
#define UNIQUE(v) { sort((v).begin(), (v).end()); (v).erase(unique((v).begin(), (v).end()), (v).end()); }
template<class T> bool chmax(T &a, const T &b) {if (a < b) {a = b; return true;} return false; }
template<class T> bool chmin(T &a, const T &b) {if (a > b) {a = b; return true;} return false; }
#define INF 0x7FFFFFFF
#define LINF 0x7FFFFFFFFFFFFFFFLL
#define Yes(q) (q ? "Yes" : "No")
#define YES(q) (q ? "YES" : "NO")
#define DUMP(q) cerr << "[DEBUG] " #q ": " << (q) << " at " __FILE__ ":" << __LINE__ << endl
#define DUMPALL(q) { cerr << "[DEBUG] " #q ": ["; REP(dumpall_i, (q).size()) { cerr << q[dumpall_i] << (dumpall_i == (q).size() - 1 ? "" : ", "); } cerr << "] at " __FILE__ ":" << __LINE__ << endl; }
template<class T> T gcd(const T &a, const T &b) { return a < b ? gcd(b, a) : b ? gcd(b, a % b) : a; }
template<class T> T lcm(const T &a, const T &b) { return a / gcd(a, b) * b; }

// gcc拡張マクロ
#define popcount __builtin_popcount
#define popcountll __builtin_popcountll

// エイリアス
using  ll = long long;
using ull = unsigned long long;
using  ld = long double;
using namespace std;

// モジュール
vector<ll> factorize(ll n) {

    if (n < 2) return vector<ll>();

    vector<ll> ans;
    ll n0 = n;
    while (n % 2 == 0) {
        ans.push_back(2);
        n /= 2;
    }
    for (ll i = 3; i * i <= n0; i += 2) {
        while (n % i == 0) {
            ans.push_back(i);
            n /= i;
        }
    }
    if (n > 1) ans.push_back(n);
    return ans;

}

ll reconstruct(const vector<pair<ll, ll>> &p, ll limit = -1) {

    ll r = 1;
    REP(i, limit < 0 ? p.size() : min(limit, (ll)p.size())) {
        for (ll j = p[i].second; j--; ) {
            r *= p[i].first;
        }
    }
    return r;

}

ll solve(vector<pair<ll, ll>> &s, const vector<pair<ll, ll>> &pfac, ll m, ll d) {
    
    if (d >= pfac.size()) {

        ll x = reconstruct(s);
        // REP(i, s.size()) cerr << "[" << s[i].first << " ^ " << s[i].second << "] "; cerr << endl;
        return x <= m;

    } else {

        ll ans = 0;
        REP(i, pfac[d].second + 1) {
            s[d].second = i;
            ll rec = reconstruct(s, d + 1);
            if (m < rec) break;
            ans += solve(s, pfac, m, d + 1);
        }
        return ans;

    }

}

// 処理内容
int main() {
    
    ll n, k, m; cin >> n >> k >> m;
    
    vector<ll> nfac = factorize(n);
    ll facsize = nfac.size();

    vector<pair<ll, ll>> pfac;
    for (ll i = 0; i < facsize; ) {
        ll j = i;
        while (j < facsize && nfac[i] == nfac[j]) j++;
        pfac.emplace_back(nfac[i], (j - i) * k);
        i = j;
    }

    vector<pair<ll, ll>> s = pfac;
    for (auto &si : s) si.second = 0;
    cout << solve(s, pfac, m, 0) << endl;
    
}
0