Pagini recente » Cod sursa (job #1528683) | Cod sursa (job #1171145) | Cod sursa (job #1879867) | Cod sursa (job #1057201) | Cod sursa (job #956833)
Cod sursa(job #956833)
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <queue>
#include <vector>
#define MAXN 505
using namespace std;
queue <int> q;
int P[MAXN];
int res[MAXN*7][MAXN];
int sss[MAXN*7][MAXN];
int p;
int N;
int gg,ss;
double dp[MAXN][MAXN];
vector <int> vec;
bool isPrime(int x)
{
for(int i=0; i < p && P[i] * P[i] <= x;i++)
{
if((x%P[i]) == 0) return 0;
}
return 1;
}
double solve(int k,int s)
{
if(k==p || !s) return 0;
double &ret = dp[k][s];
if(ret != -1) return ret;
int t=P[k];
ret=solve(k+1,s);
for(int i=1;t < s && t != ss;t*=P[k],i++)
{
if(ret < log(t)+solve(k+1,s-t))
{
ret=log(t)+solve(k+1,s-t);
res[k][s]=t;
sss[k][s]=i;
}
}
return ret;
}
void mem()
{
while(1) q.push(2);
}
int main()
{
freopen("nummst.in","r",stdin);
freopen("nummst.out","w",stdout);
scanf(" %d",&N);
for(int i=2;;i++)
{
if(!(N%i))
{
ss=i;
gg=N/i;
break;
}
}
for(int i=2;i<=ss;i++)
if(isPrime(i)) P[p++]=i;
for(int i=0;i<500;i++)
for(int j=0;j<500;j++)
dp[i][j]=-1;
solve(0,ss);
return 0;
for(int i=0;i<p;i++)
{
if(sss[i][ss])
{
printf("%d ",res[i][ss]*gg);
}
ss-=res[i][ss];
}
if(ss > 1000 || ss < 0) mem();
while(ss--) printf("%d ",gg);
printf("\n");
return 0;
}