Pagini recente » Cod sursa (job #586242) | Cod sursa (job #1927656) | Cod sursa (job #1555600) | Cod sursa (job #820471) | Cod sursa (job #425910)
Cod sursa(job #425910)
#include<stdio.h>
#include<stdlib.h>
#include<cmath>
FILE *f=fopen("desc.in","r");
FILE *g=fopen("desc.out","w");
long long divs[3001],n,k,nrdivs,i,j,p,x,y,u,m;
int a[3001][3001];
int cmp(const void *a,const void *b)
{
long long x = *((long long*)a);
long long y = *((long long*)b);
if (x < y)
return -1;
if (x == y)
return 0;
return 1;
}
void divide(){
long long d;
for(d=2;d<=sqrt(n);d++)
if(n%d==0){
divs[0]++;
divs[divs[0]]=d;
}
nrdivs=divs[0];
for(d=1;d<=nrdivs;d++)
if(n/divs[d]!=divs[d])
{ divs[0]++;
divs[divs[0]]=n/divs[d];
}
divs[0]++;
divs[divs[0]]=n;
nrdivs=divs[0];
qsort(divs+1, nrdivs, sizeof(long long),cmp);
}
long long caut(long long y){
if(y==1)
return 1;
p=1;
u=nrdivs;
while(p<=u){
m=(p+u)/2;
if(divs[m]>y)
u=m-1;
else if(divs[m]<y)
p=m+1;
else return m;}
return 0;
}
/*
long long caut1(long long y){
if(y==1)
return 0;
p=1;
u=nrdivs;
while(p<=u){
m=(p+u)/2;
if(divs[m]>y)
u=m-1;
else if(divs[m]<y)
p=m+1;
else return m;}
return 0;
} */
void recons(int lvl, int k,long long i){
int j,p;
for(j=i;j<=lvl;j++){
if(divs[lvl]%divs[j]==0){
if(lvl==j)
fprintf(g,"%lld ",divs[lvl]);
else {
p=caut(divs[lvl]/divs[j]);
if(a[p][j]<k)
k-=a[p][j];
else { fprintf(g,"%lld ",divs[j]);
recons(p,k,j);
break;
}
}
}
}
}
int main(){
//citire
fscanf(f,"%lld%lld",&n,&k);
//divsizorii lui n
divide();
//initializam prima linie cu 1
//for(i=1;i<=nrdivs;i++)
// a[0][i]=1;
for(i=1;i<=nrdivs;i++)
for(j=i;j>=1;j--){
a[i][j]=a[i][j+1]; //valoarea anterioara
if(divs[i]%divs[j]==0) //daca se imparte
{ y=divs[i]/divs[j]; //trebuie sa cautam pozitia divsizorului p
if(y==1)
a[i][j]+=1;
else{
x=caut(y);
a[i][j]+=a[x][j];
}
}
}
//reconstituire
fprintf(g,"%d\n",a[nrdivs][1]);
recons(nrdivs,k,1);
return 0;
}