factorial_sol.cpp