求sum(gcd(i,N)),1 <= i <= N
//1.gcd(i,m * n) 当m与n互质时,= gcd(i,n) * gcd(i,m) 因此他为积性函数
//2.sum(gcd(i,N)) 积性函数的和也为积性函数
//3.f(N) = sum(gcd(i,N)) = sum(p * phi(N / p)) p为N的因子
//4.f(p^r) = r * (p^r - p^r-1) + p^r
//质因子分解N,用积性函数求解
#include <iostream>
#include <cstdio>