Cod sursa(job #532514)

Utilizator lacalutlacalut lacalut Data 11 februarie 2011 20:53:58
Problema Descompuneri Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.19 kb
#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;
}