Cod sursa(job #3304635)

Utilizator 27augainHan Duy Anh 27augain Data 25 iulie 2025 16:52:39
Problema Descompuneri Scor 30
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 2.13 kb
#include <bits/stdc++.h>
using namespace std;

#define MASK(i) (1LL << (i))
#define BIT(x, i) (((x) >> (i)) & 1LL)

using ll = long long;

const int mxn = 9 + 4e3;
const int inf = 0x3f3f3f3f;
const ll lnf = 0x3f3f3f3f3f3f3f3f;

namespace std {
  template <int D, typename T> struct Vec : public vector<Vec<D - 1, T>> {
    static_assert(D >= 1, "Dimension must be positive");
    template <typename... Args>
    Vec(int n = 0, Args... args) : vector<Vec<D - 1, T>>(n, Vec<D - 1, T>(args...)) {}
  };
  
  template <typename T> struct Vec<1, T> : public vector<T> {
    Vec(int n = 0, T val = T()) : std::vector<T>(n, val) {}
  };
}

Vec<1, int> d;
unordered_map<ll, int> idx;
ll mp[mxn][mxn];

ll solve(int id, int mn) {
  if (d[id] == 1) {
    return 1;
  }
  if (mp[id][mn]) {
    return mp[id][mn];
  }
  ll res = 0;
  for (int j = mn; j < (int) (d.size()); ++j) {
    if (d[id] % d[j] == 0) {
      res += solve(idx[d[id] / d[j]], j);
    }
  }
  return mp[id][mn] = res;
}

void get(int id, int mn, int k) {
  if (d[id] == 1) {
    return;
  }
  for (int j = mn; j < (int) (d.size()); ++j) {
    if (d[id] % d[j] == 0) {
      ll cnt = solve(idx[d[id] / d[j]], j);
      if (cnt >= k) {
        cout << d[j] << " ";
        get(idx[d[id] / d[j]], j, k);
        return;
      } else {
        k -= cnt;
      }
    }
  }
}

void _27augain(void) {
  ll n;
  int k;
  cin >> n >> k;
  for (ll i = 1; i * i <= n; ++i) {
    if (n % i == 0) {
      d.push_back(i);
      if (i * i ^ n) {
        d.push_back(n / i);
      }
    }
  }
  sort(begin(d), end(d));
  for (int i = 0; i < (int) (d.size()); ++i) {
    idx[d[i]] = i;
  }
  cout << solve((int) (d.size()) - 1, 1) << "\n";
  get((int) (d.size()) - 1, 1, k);
}

int32_t main(void) {
#define task "desc"
  cin.tie(0)->sync_with_stdio(0);
  for (string iext : {"in", "inp"}) {
    if (fopen((task "." + iext).c_str(), "r")) {
      freopen((task "." + iext).c_str(), "r", stdin);
      freopen(task ".out", "w", stdout);
    }
  }
  int testcase = 1;
  // cin >> testcase;
  while (testcase--) {
    _27augain();
  }
  return 0;
}