prime_sol.cpp