记录编号	357523	评测结果	AAAAAAAAAA
题目名称	[HAOI 2011]问题B	最终得分	100
用户昵称	AntiLeaf	是否通过	通过
代码语言	C++	运行时间	4.514 s
提交时间	2016-12-11 14:13:01	内存使用	0.72 MiB

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=50010;
void mu_table(int);
int solve(int,int);
bool isp[maxn];
int prime[maxn],cnt=0,mu[maxn];
int T,a,b,c,d,k;
int main(){
#define MINE
#ifdef MINE
	freopen("b.in","r",stdin);
	freopen("b.out","w",stdout);
#endif
	mu_table(50000);
	scanf("%d",&T);
	while(T--){
		scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
		printf("%d\n",solve(b,d)-solve(a-1,d)-solve(b,c-1)+solve(a-1,c-1));
	}
#ifndef MINE
	printf("\n-------------------------DONE-------------------------\n");
	for(;;);
#endif
	return 0;
}
void mu_table(int n){
	fill(isp+1,isp+n+1,true);
	mu[1]=1;
	for(int i=2;i<=n;i++){
		if(isp[i]){
			prime[++cnt]=i;
			mu[i]=-1;
		}
		for(int j=1;j<=cnt;j++){
			if(i*prime[j]>n)break;
			isp[i*prime[j]]=false;
			if(i%prime[j])mu[i*prime[j]]=-mu[i];
			else{
				mu[i*prime[j]]=0;
				break;
			}
		}
		mu[i]+=mu[i-1];
	}
}
int solve(int n,int m){
	n/=k;m/=k;
	int ans=0,last=0;
	for(int i=1;i<=n&&i<=m;i=last+1){
		last=min(n/(n/i),m/(m/i));
		ans+=(n/i)*(m/i)*(mu[last]-mu[i-1]);
	}
	return ans;
}
/*
	if(n>m) swap(n,m);
	for(i=1;i<=n;i=last+1){
		last=min(n/(n/i),m/(m/i));
		re+=(n/i)*(m/i)*(sum[last]-sum[i-1]);
	}
	return re;
话说这是为啥......我有必要好好想想......
*/