记录编号 483531 评测结果 AAAAAAAAAA
题目名称 [HAOI 2011]问题B 最终得分 100
用户昵称 GravatarShirry 是否通过 通过
代码语言 C++ 运行时间 10.292 s
提交时间 2018-01-17 17:13:20 内存使用 0.91 MiB
显示代码纯文本
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=50010;
int n,a,b,c,d,k,cnt,pri[maxn],sum[maxn],mobius[maxn];
bool vis[maxn];
void GetMobius(){
	mobius[1]=1;
	for(int i=2;i<=maxn;i++){
		if(!vis[i])mobius[i]=-1,pri[++cnt]=i;
		for(int j=1;j<=cnt&&i*pri[j]<=maxn;j++){
			vis[i*pri[j]]=1;
			if(i%pri[j])mobius[i*pri[j]]=-mobius[i];
			else{mobius[i*pri[j]]=0;break;}
		}
	}
	for(int i=1;i<=maxn;i++)sum[i]=sum[i-1]+mobius[i];
}
long long work(long long x,long long y){
	x/=k,y/=k;
	if(x>y)swap(x,y);
	if(x==0)return 0;
	long long tot=0,tmp;
	for(long long i=1;i<=x;i=tmp){
		tmp=min(x/(x/i),y/(y/i));
		tot+=(x/i)*(y/i)*(sum[tmp]-sum[i-1]);
		tmp++;
	}
	return tot;
}
int main(){
	freopen("b.in","r",stdin);
	freopen("b.out","w",stdout);
	GetMobius();
	scanf("%d",&n);
	while(n--){
		scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
		long long ans=work(b,d)+work(a-1,c-1)-work(a-1,d)-work(b,c-1);
		printf("%lld\n",ans);
	}
	return 0;
}