Pagini recente » Borderou de evaluare (job #133086) | Cod sursa (job #404125) | Cod sursa (job #860013) | Cod sursa (job #2390701) | Cod sursa (job #778579)
Cod sursa(job #778579)
#include<fstream>
#include<algorithm>
#include<vector>
using namespace std;
ifstream f("desc.in");
ofstream g("desc.out");
long long n;
int i,j,t;
int k,nrdiv;
vector<long long> v;
long long produs;
int m[3005][3005];
int main()
{f>>n>>k;
nrpoz=0;
for(i=2; i*i<=n; i++)
if(n%i==0)
{nrdiv++; v.push_back(i);
if(i*i!=n)
{nrdiv++; v.push_back(n/i);}
}
nrdiv++;
v.push_back(n);
v.push_back(0);
sort(v.begin(),v.end());
for(i=1; i<=nrdiv; i++)
m[i][i]=1;
for(i=2; i<=nrdiv; i++)
for(j=i-1; j>=1; j--)
{m[i][j]=m[i][j+1];
if(v[i]%v[j]==0)
{t=1;
while(v[t]<v[i]/v[j]) t++;
m[i][j]+=m[t][j];
}
}
/*for(i=1; i<=nrdiv; i++)
{for(j=1; j<=nrdiv; j++)
g<<m[i][j]<<" ";
g<<endl;} */
g<<m[nrdiv][1]<<endl;
i=nrdiv;
j=1;
while(n>1)
{if(m[i][j]-m[i][j+1]<k)
{k=k-m[i][j]+m[i][j+1];
j++;}
else
{
g<<v[j]<<" ";
n=n/v[j];
if(n!=1)
while(v[i]>n) i--;
}
}
f.close();
g.close();
return 0;}