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

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

题解：狂戳这里
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;

const int N=1e5+5;
int p
,tot,check[N+10];
int moblus
,sum
;
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;
}