Cod sursa(job #778579)

Utilizator ionut_blesneagIonut Blesneag ionut_blesneag Data 15 august 2012 01:03:31
Problema Descompuneri Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.06 kb
#include<fstream>
#include<algorithm>
#include<vector>
using namespace std;

ifstream f("desc.in");
ofstream g("desc.out");

long long n;
int i,j,t;
int k,nrdiv;
vector<long long> v;
long long produs;
int m[3005][3005];


int main()
{f>>n>>k;

nrpoz=0;
for(i=2; i*i<=n; i++)
  if(n%i==0)
  {nrdiv++;  v.push_back(i);   
   if(i*i!=n)
      {nrdiv++;  v.push_back(n/i);} 
   }     
nrdiv++;  
v.push_back(n);
v.push_back(0);
sort(v.begin(),v.end());    

for(i=1; i<=nrdiv; i++)
  m[i][i]=1;
    
for(i=2; i<=nrdiv; i++)
  for(j=i-1; j>=1; j--)
    {m[i][j]=m[i][j+1];
     if(v[i]%v[j]==0)
     {t=1;
      while(v[t]<v[i]/v[j])   t++;
      m[i][j]+=m[t][j];
      }
    }
          
/*for(i=1; i<=nrdiv; i++)
 {for(j=1; j<=nrdiv; j++)
   g<<m[i][j]<<" ";
 g<<endl;}  */

g<<m[nrdiv][1]<<endl;

i=nrdiv;
j=1;

while(n>1)
{if(m[i][j]-m[i][j+1]<k)
   {k=k-m[i][j]+m[i][j+1];
    j++;}
 else
   {
    g<<v[j]<<" ";
    n=n/v[j];
    if(n!=1)
       while(v[i]>n)  i--;
   }     
}
 
f.close();
g.close();
return 0;}