51nod 1575 Gcd and Lcm【Min_25筛】

题目描述：

给出$n\le 10^9$，求：
$$\sum _{i=1}^n\sum _{j=1}^i\sum _{k=1}^i[(i,j),(i,k)]\mathrm{m}\mathrm{o}\mathrm{d}{2}^{32}$$

题目分析：

很好，又是一道筛题，我们开始推式子：

推式子内容详见tangjz的答案，以及佐理慧的blog，省略1万字。。

然后我们自闭了。

然后有人告诉我们${\sum }_{j=1}^i{\sum }_{k=1}^i[(i,j),(i,k)]$是一个关于$i$的积性函数。
$$F(P^k)=(2k+1)(P^{2k}-P^{2k-1})+P^{k-1}$$
很好，直接上Min_25筛。
。。。。。。

Code：

#include<bits/stdc++.h>
#define maxn 100005
using namespace std;
typedef unsigned int ui;
const ui inv3 = 2863311531u;
int T,n,m,sn,a[maxn];
ui s[3][maxn],f[maxn];
inline int ID(int x){return x<=sn?x:m-n/x+1;}
int main()
{
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n),sn=sqrt(n),m=0;
		for(int i=1;i<=n;i++) a[++m]=i=n/(n/i),s[0][m]=i-1,s[1][m]=1ll*i*(i+1)/2-1,s[2][m]=1ll*i*(i+1)/2*(2*i+1)*inv3-1;
		//1到sn必然在a中，即a[i<=sn]=i
		for(int i=2;i<=sn;i++) if(s[0][i]^s[0][i-1])//已经筛去了<i的质数，若s[0][i]!=s[0][i-1],说明i是质数
			for(int j=m;a[j]>=i*i;j--)
				s[0][j]-=s[0][ID(a[j]/i)]-s[0][i-1],
				s[1][j]-=i*(s[1][ID(a[j]/i)]-s[1][i-1]),
				s[2][j]-=i*i*(s[2][ID(a[j]/i)]-s[2][i-1]);
		for(int i=1;i<=m;i++) f[i]=3*s[2][i]-3*s[1][i]+s[0][i];//f(p)=3*p^2-3*p+1
		for(int i=sn;i>1;i--) if(s[0][i]^s[0][i-1])//Min_25筛非递归写法
			for(int j=m;a[j]>=i*i;j--)
				for(ui k=1,p0=1,p1=i,sp=3*s[2][i]-3*s[1][i]+s[0][i];1ll*p1*i<=a[j];k++,p0=p1,p1*=i)
					f[j]+=((2*k+1)*p1*p0*(i-1)+p0)*(f[ID(a[j]/p1)]-sp)+(2*k+3)*p1*p1*i*(i-1)+p1;
		printf("%u\n",f[m]+1);
	}
}