Problem b【HAOI2011】【莫比乌斯反演进阶】

版权声明：这是gigo写的QAQ https://blog.csdn.net/qq_42835815/article/details/85014385

传送门：https://www.luogu.org/problemnew/show/P2522

这道题，，其实就是zap-queries 的进阶版

在其基础上加了个差分，其他只字未变

上代码

// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define in read()
using namespace std;

inline int in{
	int cnt=0,f=1; char ch=0;
	while(!isdigit(ch)){
		ch=getchar();
		if(ch=='-')f=-1;
	}
	while(isdigit(ch)){
		cnt=cnt*10+ch-48;
		ch=getchar();
	}
	return cnt*f;
}
#define N 50003
int mu[N],prim[N],cnt,vis[N];
long long sum[N];
inline void mumu(int n){
	mu[1]=1;
	for(register int i=2;i<=n;i++){
		if(!vis[i]){
			prim[++cnt]=i;
			mu[i]=-1;
		}
		for(register int j=1;j<=cnt&&prim[j]*i<=n;j++){
			vis[prim[j]*i]=1;
			if(i%prim[j]==0)break;
			else mu[prim[j]*i]=-mu[i];
		}
	}
	for(register int i=1;i<=n;i++)sum[i]=sum[i-1]+(long long)mu[i];
}int a,b,c,d,k;
inline long long query(int a,int b){
	if(a>b)swap(a,b);long long ans=0;
	for(register int l=1,r;l<=a;l=r+1){
		r=min(a/(a/l),b/(b/l));
		ans+=(1ll*a/(1ll*l*k))*(1ll*b/(1ll*l*k))*(sum[r]-sum[l-1]); 
	}
	return ans;
}
int main(){
	int t=in;mumu(N);
	while(t--){
		a=in;b=in;c=in;d=in;k=in;
		printf("%lld\n",query(b,d)-query(b,c-1)-query(a-1,d)+query(a-1,c-1));
	}
	return 0;
} 

注意1ll的重要性，，longlong嘛毕竟。

公式推法和上次一样。

另外，，本题如果不开O2会T6个点。。ovo