这题统计答案你可能用了$floor()$，but会有精度问题，记得

const double eps = 1e-7;
......
ans += phi[i] * (floor(n / i) + eps) * (floor(n / i) + eps);