Разбор алгоритма программы

#include "stdafx.h"
#include <iostream>
 
 
 
using namespace std;
int a[1000],i,n,j,ost=0,num=1,c;
 
int main(){
    cin>>n;
    memset(a,0,sizeof(a));
    if(n==0 || n==1) {cout<<1; return 0;}
    a[0]=1;
    for(i=2;i<=n;i++){
        for(j=0;j<num;j++){
            c=(a[j]*i+ost);
            a[j]=c%1000;
            ost=c/1000;
 
        }
        a[num]=ost%1000;
        ost/=1000;
        a[num+1]=ost%1000;
        ost/=1000;
        a[num+2]=ost;
        ost=0;
        if(a[num+1]!=0)num+=2;
        else if(a[num]!=0)num++;
    
    }
    cout<<a[num-1];
    for(i=num-2;i>=0;i--) {
        if(a[i]>99) cout<<a[i];
        else {
            if (a[i]>9) {cout<<0<<a[i];}
            else cout<<'0'<<'0'<<a[i];
        }
    }
return 0;
}

#include <iostream>
 
unsigned long long fact(int n)
{
      if(n < 0) 
        return 0; 
      if (n == 0 || n == 1) 
        return 1; 
      else 
        return n * fact(n - 1); 
}
 
int main()
{
  std::cout << fact(20) << std::endl;
  
  system("PAUSE");
  return 0;
}