Cod sursa(job #2788285)

Utilizator PredescuSebastianIonPredescu Sebastian Ion PredescuSebastianIon Data 25 octombrie 2021 14:35:04
Problema NumMst Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.25 kb
#include <fstream>
#include <cmath>

using namespace std;
ifstream f("nummst.in");
ofstream g("nummst.out");
const int nmax=452,pmax=3192;
int n,m,d,prime[nmax],poz[nmax][pmax],sol[pmax];
double dp[2][pmax];
int main()
{
    f>>n;
    d=2;
    while(n%d>0)d++;
    for(int i=2; i<d; i++)
    {
        int j=2;
        while(j<=i/2 && i%j>0)
        {
            j++;
        }
        if(j==i/2+1)prime[++m]=i;
    }
    for(int i=1; i<=m; i++)
    {
        for(int j=1; j<=d; j++)
        {
            dp[i%2][j]=dp[(i-1)%2][j];
        }
        for(int j=prime[i]; j<=d; j*=prime[i])
        {
            double tmp=log(j);
            for(int k=j; k<=d; k++)
            {
                if(dp[i%2][k]<dp[(i-1)%2][k-j]+tmp)
                {
                    dp[i%2][k]=dp[(i-1)%2][k-j]+tmp;
                    poz[i][k]=j;
                }
            }
        }
    }
    int nr=d;
    for(int i=m; i>=1; i--)
    {
        if(poz[i][nr])
        {
            sol[++sol[0]]=poz[i][nr];
            nr-=poz[i][nr];
        }
    }
    for(int i=1; i<=nr; i++)
    {
        sol[++sol[0]]=1;
    }
    d=n/d;
    for(int i=1; i<=sol[0]; i++)
    {
        g<<sol[i]*d<<" ";
    }
    return 0;
}