题意:
给出两个区间,求这个两个区间这样的一对数:(a,b) gcd(a,b)=k 的个数,(a,b)和(b,a)等价。
题解:
要求区间(1,l)(1,r)一对数满足最大公倍数等于k,那么等价于求(1,l/k)(1,r/k)一对数满足互质。继续分析发现可以将问题分成两个部分解决:我们假设l比r更小
1、(1,l)和(1,l)这个区间的某个数i,其中一部分解就是这些i对应在(1,i)区间互质的数的个数和。
2、(1,l)和(l,r)求第二个区间的每个数i在第一个区间中互质的数的个数。
第一部分的解计算方法可以用欧拉函数求解。第二部分的解用容斥解决。
对于第二部分解求法:
对于区间(l,r)某个数i,要求互质的数的个数,我们可以先考虑逆问题,我们先求不互质的数的个数,那么不互质的数肯定满足是这个数i质因子的倍数,那么根据容斥原理去枚举子因子,用dfs求解。然后累加所有的i对应互质数的个数。
注意:k可以等于0
#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<string>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
void cmax(int& a,int b){ if(b>a)a=b; }
void cmin(int& a,int b){ if(b<a)a=b; }
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=1000007;
const int maxn=100005;
ll euler[maxn];
int p[maxn][15];
int num[maxn];
void Init(){
euler[1]=1;
for(int i=2;i<maxn;i++){
if(!euler[i]){
for(int j=i;j<maxn;j+=i){
if(!euler[j]) euler[j]=j;
euler[j]=euler[j]*(i-1)/i;
p[j][num[j]++]=i;
}
}
euler[i]+=euler[i-1];
}
}
ll dfs(int cnt,int x,int a){
ll res=0;
for(int i=cnt;i<num[a];i++){
res+=x/p[a][i]-dfs(i+1,x/p[a][i],a);
}
return res;
}
int main(){
//freopen("E:\\read.txt","r",stdin);
int a,b,k,T,l,r;
ll ans;
Init();
scanf("%d",&T);
for(int cas=1;cas<=T;cas++){
scanf("%d%d%d%d%d",&a,&a,&b,&b,&k);
if(k==0){ printf("Case %d: 0\n",cas);continue; }
l=min(a,b);
r=max(a,b);
l/=k;
r/=k;
ans=euler[l];
for(int i=l+1;i<=r;i++){
ans+=l-dfs(0,l,i);
}
printf("Case %d: %I64d\n",cas,ans);
}
return 0;
}