2301: [HAOI2011]Problem b

花了一个晚上学莫比乌斯，觉得自己的脑子被烧坏了。特别感谢连击哥的谆谆教诲，否则我这么蒟蒻根本学不下去。

接下来简单介绍一下算法。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
#define rep(i,j,k) for(i=j;i<=k;++i)
#define per(i,j,k) for(i=j;i>=k;--i)
#define ll long long
#define pli pair<ll,int>
#define mkp make_pair
#define X first
#define Y second
#define N 50005
ll mu[N],pri[N],cnt,Smu[N];bool vis[N];
void getmu(){
	ll i,j;Smu[1]=mu[1]=1;
	rep(i,2,N-1){
		if(!vis[i])mu[pri[++cnt]=i]=-1;
		for(j=1;j<=cnt&&i*pri[j]<N;j++){
			vis[i*pri[j]]=1;
			if(i%pri[j]==0)break;
			mu[i*pri[j]]=-mu[i];
		}
		Smu[i]=Smu[i-1]+mu[i];
	}
}
ll cal(ll n,ll m,ll k){
	if(!n||!m)return 0;
	n/=k;m/=k;
	ll i,j,sum=0,rg=min(n,m);for(i=1;i<=rg;i=j+1){
		j=min(n/(n/i),m/(m/i));
		sum+=(Smu[j]-Smu[i-1])*(n/i)*(m/i);
	}
	return sum;
}
int main(){
	ll _,a,b,c,d,e;
	scanf("%lld",&_);
	getmu();
	while(_--){
		scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&e);
		printf("%lld\n",cal(b,d,e)+cal(a-1,c-1,e)-cal(b,c-1,e)-cal(a-1,d,e));
	}
	return 0;
}

代码略丑，请见谅。