Pagini recente » Cod sursa (job #1862151) | Cod sursa (job #2671694) | Cod sursa (job #2186998) | Cod sursa (job #187059) | Cod sursa (job #532514)
Cod sursa(job #532514)
#include <algorithm>
#include <stdio.h>
#include <vector>
#define ll long long
#define MAX 2048
#define pb push_back
#define mp make_pair
using namespace std;
ll n, k, rez;
vector <ll> vctDiv;
ll vctSol[MAX][MAX];
int main()
{
freopen("desc.in", "r", stdin);
freopen("desc.out", "w", stdout);
scanf("%lld %lld", &n, &k);
vctDiv.pb(n);
for (int i = 2; i * i <= n; i++)
if (n % i == 0)
{
vctDiv.pb(i);
if (i * i != n)
vctDiv.pb(n / i);
}
int m = vctDiv.size();
sort(vctDiv.begin(), vctDiv.end());
for (int i = 0; i < m; i++)
{
int r = 0;
vctSol[i][i] = 1;
for (int j = i; j >= 0; j--)
{
vctSol[i][j] = vctSol[i][j + 1];
if (vctDiv[i] % vctDiv[j] == 0)
{
for (; vctDiv[r] * vctDiv[j] < vctDiv[i]; r++);
vctSol[i][j] += vctSol[r][j];
}
}
}/*
printf("%lld\n", vctSol[m - 1][0]);
int l = m - 1;
for (int i = 0; n > 1 && i < m; )
if (vctSol[l][i] - vctSol[l][i + 1] >= k)
{
printf("%lld ", vctDiv[i]);
n /= vctDiv[i];
if (n > 1)
for (; vctDiv[l] > n; l--);
}
else k -= vctSol[l][i] - vctSol[l][i + 1], i++;*/
fclose(stdin);
fclose(stdout);
return 0;
}