#include #include extern void generate_primes(void); extern unsigned int gcd(unsigned int x, unsigned int y); extern void factorize(unsigned int n, std::vector& product); extern mpz_class foc_reference(unsigned int n); extern mpz_class const& foc(unsigned int n); extern int number_of_primes; extern unsigned int primes[]; inline unsigned int lcm(unsigned int x, unsigned int y) { return (x / gcd(x, y)) * y; }