题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1695
题意:因为题目已经说明固定a=1,c=1,所以就是x在[1,b],y在[1,d],求有多少对(x,y)满足gcd(x,y)=k;(gcd指的是求俩数的最小公倍数)。
求解方式:一共有两种:一个是用容斥原理+欧拉函数,一个是莫比乌斯反演,在这我们用的是容器定理+欧拉函数,莫比乌斯反演我会另开一篇文章来讲。
先来说一下欧拉函数。欧拉函数用phi(n)表示(应该是φ(n),但φ的读音就是phi,在这就用phi表示φ),phi(n)表示小于n的数中有phi(n)个数与n互质。
欧拉函数通式:φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn),
再来说一下容斥原理。公式:
也可表示为:
设S为有限集,
,则
看起来乱乱的公式其实就是小学学的数学知识,仔细看看就理解了。
基础已说完,说一下题目解题思路:
在这先让b为b,d的最小值,若b>d就把b,d的值换一下(这并不影响结果)。
因为是求gcd(x,y)=k,所以干脆把区间范围除以k,相当于求[1,b/k],[1,d/k]中互质个数。(若不明白仔细想想)。
第一步,先求x,y都在[1,a]区间时他们互质的对数,看一下欧拉函数的定义,实质上就是求欧拉函数的前a项和(所以我们要预处理每个数的欧拉函数。。。)。
第二步,再求y在[b+1,d],x在[1,b]中时,有多少对(x,y)互质,那么我们来枚举y,首先将y进行合数分解,得到fatcnt个不同素数,然后套用容斥定理就ok,(请注意为什么最后的返回值是n-ans)。
AC代码:
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cmath>
#include<vector>
#include<map>
#include<queue>
using namespace std;
long long a,b,c,d,k;
/*素数筛选和合数分解*/
const int maxn=100100;
int prime[maxn+1];
void getprime(){
memset(prime,0,sizeof(prime));
for(int i=2;i<maxn;i++){
if(!prime[i]) prime[++prime[0]]=i;
for(int j=1;j<=prime[0]&&prime[j]<=maxn/i;j++){
prime[prime[j]*i]=1;
if(i%prime[j]==0) break;
}
}
}
/*合数分解*/
long long factor[105][2];
int fatcnt;
int getfactor(long long x){
if(x==0)
return 0;
fatcnt=0;
long long tmp=x;
for(int i=1;prime[i]<=tmp/prime[i];i++){
factor[fatcnt][1]=0;
if(tmp%prime[i]==0){
factor[fatcnt][0]=prime[i];
while(tmp%prime[i]==0){
factor[fatcnt][1]++;
tmp/=prime[i];
}
fatcnt++;
}
}
if(tmp!=1){
factor[fatcnt][0]=tmp;
factor[fatcnt++][1]=1;
}
return fatcnt;
}
/*欧拉函数筛*/
int phi[100100];
void phitable_ZY(int n){
memset(phi,0,sizeof(phi));
phi[1]=1;
for(int i=2;i<=n;i++) if(!phi[i])
for(int j=i;j<=n;j+=i){
if(!phi[j]) phi[j]=j;
phi[j]=phi[j]/i*(i-1);
}
}
long long cr(long long i){
getfactor(i);
long long ans=0;
for(int i=1;i<(1<<fatcnt);i++){
long long t=1;
int cut=0;
for(int j=0;j<fatcnt;j++){
if(i&(1<<j))
t*=factor[j][0],cut++;
}
if(cut&1)
ans+=b/t;
else
ans-=b/t;
}
return (long long)b-ans;
}
int main(){
getprime();
phitable_ZY(100010);
int t,cas=0;
scanf("%I64d",&t);
while(t--){
scanf("%I64d%I64d%I64d%I64d%I64d",&a,&b,&c,&d,&k);
printf("Case %d: ",++cas);
if(k==0||k>b||k>d){
printf("0\n");
continue;
}
b/=k;
d/=k;
long long ans=0;
if(b>d)
swap(b,d);
for(int i=1;i<=b;i++)
ans+=phi[i];
for(int i=b+1;i<=d;i++){
ans+=cr(i);
}
printf("%I64d\n",ans);
}
return 0;
}