【51nod1584】加权约数和（莫比乌斯反演）（线性筛）

传送门

---

题解：

一看这种题就是大力上莫比乌斯反演。

首先$max$并不是能够处理的东西。

设：
$$\begin{gathered} A(n)=\sum _{i=1}^n\sum _{j=1}^{i}i\sigma (ij) \\ B(n)=\sum _{i=1}^{n}i\sigma (i^2) \end{gathered}$$

则答案为$$Ans=2\cdot A(n)-B(n)$$

化简：
$$\begin{array}{rlrl}A(n)& =& & \sum _{i=1}^{n}i\sum _{j=1}^i\sum _{p\mid i}\sum _{q\mid j}[gcd(p,q)=1]\frac{pj}{q}\\ & =& & \sum _{d=1}^n\mu (d)\sum _{d\mid p}p\sum _{p\mid i}i\sum _{d\mid q,q\le i}\sum _{q\mid j,j\le i}\frac{j}{q}\\ & =& & \sum _{d=1}^n\mu (d)\sum _{p=1}^{\lfloor \frac{n}{d}\rfloor }dp\sum _{p\mid i,i\le \lfloor \frac{n}{d}\rfloor }id\sum _{q=1}^i\sum _{q\mid j,j\le i}\frac{j}{q}\\ & =& & \sum _{d=1}^n\mu (d)d^2\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }i\sigma (i)S_{\sigma }(i)\\ & & & \\ & & & \\ B(n)& =& & \sum _{i=1}^{n}i\sigma (i^2)\\ & =& & \sum _{i=1}^{n}i\sum _{p\mid i}\sum _{q\mid i}[gcd(p,q)=1]\frac{ip}{q}\\ & =& & \sum _{d=1}^n\mu (d)\sum _{i=1}^n\sum _{d\mid p\mid i}\sum _{d\mid q\mid i}\frac{ip}{q}\\ & =& & \sum _{d=1}^n\mu (d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }id\sum _{p\mid i}\sum _{q\mid i}\frac{idp}{q}\\ & =& & \sum _{d=1}^n\mu (d)d^2\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }i{\sigma }^2(i)\end{array}$$

其中$S_{\sigma }$表示$\sigma$的前缀和。

则我们有$$\begin{gathered} Ans=2\cdot A(n)-B(n) \\ =\sum _{d=1}^n\mu (d)d^2\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }i\sigma (i)(2\cdot S_{\sigma }(i)-\sigma (i)) \end{gathered}$$

整除分块并不像是可以过的样子。

但是我们可以考虑调和级数预处理。

枚举每一个$i$和$d$，加到$ans[i\ast d]$，我们发现，它对所有$n\ge d\ast i$的询问都产生影响，求一遍前缀和即可。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define gc get_char
#define cs const

namespace IO{
	inline char get_char(){
		static cs int Rlen=1<<18|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	
	template<typename T>
	inline T get(){
		char c;
		while(!isdigit(c=gc()));T num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	}
	inline int getint(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

cs int mod=1e9+7;
inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline void Inc(int &a,int b){(a+=b)>=mod&&(a-=mod);}
inline void Dec(int &a,int b){(a-=b)<0&&(a+=mod);}

cs int P=1e6+7;

bool mark[P];
int pr[P],pcnt;
int sigma[P],sum_sigma[P],mu[P];
int dm[P],dmt[P];

int f[P],g[P],ans[P];

inline void linear_sieves(int n=1e6){
	sigma[1]=mu[1]=1;
	for(int re i=2;i<=n;++i){
		if(!mark[i]){
			pr[++pcnt]=i;
			mu[i]=mod-1;sigma[i]=i+1;
			dm[i]=1;dmt[i]=i+1;
		}
		for(int re j=1;i*pr[j]<=n;++j){
			int k=i*pr[j];
			mark[k]=true;
			if(i%pr[j]){
				dm[k]=i;dmt[k]=pr[j]+1;
				mu[k]=dec(0,mu[i]);
				sigma[k]=mul(sigma[i],pr[j]+1);
			}
			else {
				dm[k]=dm[i];mu[k]=0;
				dmt[k]=add(mul(dmt[i],pr[j]),1);
				sigma[k]=mul(dmt[k],sigma[dm[i]]);
				break;
			}
		}
	}
	for(int re i=1;i<=n;++i){
		f[i]=mul(mu[i],mul(i,i));
		sum_sigma[i]=add(sum_sigma[i-1],sigma[i]);
		g[i]=mul(mul(i,sigma[i]),dec(add(sum_sigma[i],sum_sigma[i]),sigma[i]));
	}
	for(int re i=1;i<=n;++i){
		for(int re j=1,lj=n/i;j<=lj;++j)
		Inc(ans[i*j],mul(f[i],g[j]));
	}
	for(int re i=1;i<=n;++i)Inc(ans[i],ans[i-1]);
}

signed main(){
#ifdef zxyoi
	freopen("sigma.in","r",stdin);
#endif
	linear_sieves();int T=getint();
	for(int re i=1;i<=T;++i)cout<<"Case #"<<i<<": "<<ans[getint()]<<"\n";
	return 0;
}