結果

問題 No.847 Divisors of Power
ユーザー ikd
提出日時 2019-07-05 22:10:18
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,074 bytes
コンパイル時間 1,041 ms
コンパイル使用メモリ 85,668 KB
実行使用メモリ 9,728 KB
最終ジャッジ日時 2024-10-06 21:52:14
合計ジャッジ時間 2,017 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 25 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <vector>
#include <set>
#include <tuple>

using namespace std;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
int64_t N, K, M;
vector<pair<int64_t, int64_t>> factors;
set<int64_t> cand;

vector<pair<int64_t, int64_t>> f(int64_t n) {
  vector<pair<int64_t, int64_t>> ret;
  for (int64_t i = 2; i * i <= n; i++) {
    int cnt = 0;
    while (n % i == 0) {
      n /= i;
      cnt += 1;
    }
    if (cnt > 0) {
      ret.emplace_back(i, cnt);
    }
  }
  if (n > 1) {
    ret.emplace_back(n, 1);
  }
  return ret;
}

void dfs(int i, int64_t val) {
  if (val > M or val < 0) return;
  cand.insert(val);
  if (i >= factors.size()) return;
  int64_t v, cnt;
  tie(v, cnt) = factors[i];
  dfs(i + 1, val);
  int64_t p = v;
  rep(j, cnt) {
    if (val * p < 0 or val * p > M) break;
    dfs(i + 1, val * p);
    p *= v;
    if (p < 0 or p > M) break;
  }
}

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

  factors = f(N);
  for (auto &p: factors) {
    p.second *= K;
  }
  dfs(0, 1);
  cout << cand.size() << endl;

  return 0;
 }
0