Cod sursa(job #956837)

Utilizator ayakutAlperen Yakut ayakut Data 3 iunie 2013 22:29:52
Problema NumMst Scor 24
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.31 kb
#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);
	for(int i=0;i<p;i++)
	{
		if(sss[i][ss])
		{
			printf("%d ",res[i][ss]*gg);
		}
		ss-=res[i][ss];
	}
	if(ss < 0) mem();
	while(ss--) printf("%d ",gg);
	printf("\n");
	return 0;
}