[BZOJ 2301] HAOI 2011 Problem b · 莫比乌斯

这题和Bzoj1101有区别吗。。。囧。。。

题解：狂戳这里

#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;

const int N=1e5+5;
int p[N],tot,check[N+10];
int moblus[N],sum[N];
int a,b,c,d,k,T,ans;

void init(){
	memset(check,1,sizeof check);
	moblus[1]=1;
	for (int i=2;i<=N;i++){
		if (check[i]) {
			p[++tot]=i;
			moblus[i]=-1;
		}
		for (int j=1;j<=tot && p[j]*i<=N;j++){
			check[i*p[j]]=0; 
			if (i%p[j]==0) {
				moblus[i*p[j]]=0;
				break;
			}
			else moblus[i*p[j]]=-moblus[i];
		}
	}
	for (int i=1;i<=N;i++) sum[i]=sum[i-1]+moblus[i];
}

int calc(int n,int m){	//求[1,n][1,m]区间内互质的(x,y)的对数
	int ret=0;
	if (n>m) swap(n,m);
	for (int L=1,R=0;L<=n;L=R+1){
		R=min(n/(n/L),m/(m/L));		// 分段
		ret+=(sum[R]-sum[L-1])*(n/L)*(m/L);
	}
	return ret;
}

int main(){
	init();
	scanf("%d",&T);
	while (T--){
		scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
		ans=calc(b/k,d/k)-calc((a-1)/k,d/k)-calc(b/k,(c-1)/k)+calc((a-1)/k,(c-1)/k);
		printf("%d\n",ans);
	}
	return 0;
}