类似于欧拉函数因子的phi值之和等于n,这里直接求即可
只需要知道x^3+3x^2+x在1到n的和为n(n+1)(n+1)(n+4)/4 然后结合逆元的知识即可完成
#include<cstdio>
typedef unsigned long long ll;
const int MAXN = 2000001;
const ll ha = 1333333;
const int mod = 1000000007;
const int inv_4 = 250000002;
int g[MAXN],p[MAXN],cut,pri[MAXN],...