BZOJ P2301 [HAOI2011] Problem b【莫比乌斯】

$$\begin{gathered} f(n,m)=\sum _{i=1}^n\sum _{j=1}^m(gcd(i,j)==k) \\ =\sum _{i=1}^{\frac{n}{k}}\sum _{j=1}^{\frac{m}{k}}(gcd(i,j)==1) \\ =\sum _{i=1}^{\frac{n}{k}}\sum _{j=1}^{\frac{m}{k}}\sum _{d|i,d|j}\mu (d) \\ =\sum _{d|i,d|j}\mu (d)\sum _{i=1}^{\frac{n}{k}}\sum _{j=1}^{\frac{m}{k}} \\ =\sum _{d=1}^{min(\frac{n}{k},\frac{m}{k})}\mu (d)\sum _{i=1}^{\frac{n}{dk}}\sum _{j=1}^{\frac{m}{dk}} \\ =\sum _{d=1}^{min(\frac{n}{k},\frac{m}{k})}\mu (d)\lfloor \frac{n}{dk}\rfloor \ast \lfloor \frac{m}{dk}\rfloor \\ Ans=f(b,d)+f(a-1,c-1)-f(a-1,b)-f(b,c-1) \end{gathered}$$
由于直接计算$f$时间复杂度仍然高，但是结合$\mu (d)$的性质我们发现其实可以对$\mu (d)$做一个前缀和优化一下。

void Mobius(){
	......
}
LL GetAns(LL N,LL M,LL K){
	LL I,Ans=0,Last;
	N/=K,M/=K;
	if(N>M){
		swap(N,M);
	}
	for(I=1;I<=N;I=Last+1){
		Last=min(N/(N/I),M/(M/I));
		Ans+=(Miu[J]-Miu[I-1])*(N/I)*(M/I);
	}
	return Ans;
}
int main(){
	LL I,J,K;
	N=Read();Mobius();
	for(I=1;I<=N;I++){
		A=Read(),B=Read(),C=Read(),D=Read(),K=Read();
		Ans=GetAns(B,D,K)-GetAns(A-1,D,K)-GetAns(B,C-1,K)+GetAns(A-1,C-1,K);
		......
	}	
	return 0;
}