Pagini recente » Cod sursa (job #3266496) | Cod sursa (job #681120) | Cod sursa (job #2833510) | Cod sursa (job #880758) | Cod sursa (job #587618)
Cod sursa(job #587618)
#include<stdio.h>
#include<math.h>
#define NMAX 3333
#define minim(a,b) (a<b ? a : b)
char prim[NMAX];
int v[NMAX];
int n,ram,plus,sol[NMAX],nrs,nrp;
double d[NMAX][NMAX/5];
char pred[NMAX][NMAX/5];
double lnp[NMAX];
int main ()
{
int i,j,k,lim;
freopen("nummst.in","r",stdin);
freopen("nummst.out","w",stdout);
scanf("%d",&n);
if(n%2==0)
{
printf("%d %d\n",n/2,n/2);
return 0;
}
lim=sqrt(n);
plus=1;
ram=n;
for(i=2;i<=lim;i++)
if(n%i==0)
{
plus=n/i;
ram=i;
break;
}
for(i=2;i<=ram;i++)
if(!prim[i])
{
for(j=i*i;j<=ram;j+=i)
prim[j]=1;
v[++nrp]=i;
}
for(i=1;i<=ram;i++)
lnp[i]=log(i);
for(i=1;i<=ram;i++)
for(j=1;j<=nrp;j++)
{
d[i][j]=d[i][j-1];
pred[i][j]=0;
for(k=v[j];k<=i;k*=v[j])
if(d[i][j]<d[i-k][j-1]+lnp[k])
{
d[i][j]=d[i-k][j-1]+lnp[k];
pred[i][j]=k;
}
}
i=ram;j=nrp;
while(j)
{
if(pred[i][j])
sol[++nrs]=pred[i][j];
i-=pred[i][j];
j--;
}
if(i)
sol[++nrs]=i;
for(i=1;i<=nrs;i++)
printf("%d ",sol[i]*plus);
printf("\n");
return 0;
}