Найти наименьшее количество спичек для построения модели из N кубиков (С++ -> Pascal)

//---------------------------------------------------------------------------
 
#include <stdio.h>
#include <math.h>
int NCube(int n)
                {
                        int c1=(int)(pow(n,1.0/3)+0.5);
                        int c2=(int)(sqrt(n/c1)+0.5);
                        int c3=n/(c1*c2);
                        int S=3*c1*c2*c3+2*(c1*c2+c1*c3+c2*c3)+(c1+c2+c3);
                        int h=n-c1*c2*c3;
                        if(h!=0) {
                                int d1=(int)(sqrt(h)+0.5);
                                if(d1*d1>h) d1--;
                                int d2=h/d1;
                                S+=3*d1*d2+2*(d1+d2)+1;
                                h-=d1*d2;
                                if(h!=0) S+=3*h+2;
                        }
                        return S;
                }
//---------------------------------------------------------------------------
 
int main()
 {
        long N;
        FILE *f;
 
        f = fopen("C:\\input.txt","r");
         fscanf(f, "%d", &N);
        fclose(f);
 
        f = fopen("C:\\output.txt","w");
         fprintf(f, "%d\n", NCube(N));
        fclose(f);
 
        return 0;
 }
//---------------------------------------------------------------------------