Pagini recente » Borderou de evaluare (job #929952) | Borderou de evaluare (job #2399583) | Borderou de evaluare (job #2045353) | Borderou de evaluare (job #148322) | Cod sursa (job #3304635)
#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;
}