我用了一种正确的玄学做法，不知道名称

f[Ch]就是答案，这个方法类似反向DP

’’’

for(int k=n;k>=1;k--){

	g[k]=(n/k)*(m/k);

	f[k]=g[k];

	for(int j=2;j<=(n/k);j++)f[k]-=f[j*k];

	//ans+=f[k]*(2*k-1);

}

’’’