結果

問題 No.847 Divisors of Power
ユーザー risujirohrisujiroh
提出日時 2019-07-05 21:45:45
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 9 ms / 2,000 ms
コード長 1,023 bytes
コンパイル時間 1,774 ms
コンパイル使用メモリ 175,956 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-06 21:23:02
合計ジャッジ時間 2,351 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using lint = long long;
template<class T = int> using V = vector<T>;
template<class T = int> using VV = V< V<T> >;

template<class Z> map<Z, int> factorize(Z n) {
  map<Z, int> res;
  for (Z i = 2; i * i <= n; ++i) while (n % i == 0) ++res[i], n /= i;
  if (n != 1) ++res[n];
  return res;
}

lint powll(lint a, int n) {
  assert(n >= 0);
  lint res = 1;
  for (; n > 0; a *= a, n >>= 1) if (n & 1) res *= a;
  return res;
}

int main() {
  cin.tie(nullptr); ios::sync_with_stdio(false);
  int n, k, m; cin >> n >> k >> m;
  k = min(k, 30);
  auto mp = factorize(n);
  V<> p;
  for (const auto& e : mp) p.push_back(e.first);
  lint t = 1;
  int res = 0;
  auto dfs = [&](const auto& dfs, int j) -> void {
    if (j == p.size()) {
      ++res;
      return;
    }
    for (int i = 0; i <= mp[p[j]] * k; ++i) {
      if (t * powll(p[j], i) > m) break;
      t *= powll(p[j], i);
      dfs(dfs, j + 1);
      t /= powll(p[j], i);
    }
  };
  dfs(dfs, 0);
  cout << res << '\n';
}
0