HDU 2588 欧拉函数

题目链接

题意：
给定两个整数$N$,$M$，求有多少个整数$X$($1\le X\le N$)，满足$M\le gcd(X,N)$

思路：

题意很简单，思路却有点难想。

首先$N$是一个已知的整数，而$gcd(X,N)$一定是$N$的一个约数
而$N$的约数并不会很多，故我们可以逐一枚举，并找到每一个不小于$M$的约数，此处设该约数为P
则：
$P|N$ 且 $M\le P$

而对于枚举出的每一个P，原问题得到了简化：
对于一个已知的整数P，求有多少个整数$X$($1\le X\le N$)，满足$gcd(X,N)=P$

想到这里，解法其实离我们已经很近了。
让我们设简化后的问题答案为$Ans$：
则：
$$Ans=\sum _{i=1}^N[gcd(i,N)==P]$$
$$=\sum _{i=1}^{\lfloor \frac{N}{P}\rfloor }[gcd(i\ast P,N)==P]$$
$$=\sum _{i=1}^{\lfloor \frac{N}{P}\rfloor }[gcd(i,\lfloor \frac{N}{P}\rfloor )==1]$$

显然，上面的式子正是欧拉函数的定义式
即$\phi (\lfloor \frac{N}{P}\rfloor )$

对于不太理解欧拉函数的朋友，可以接着看下面的推导：

回到刚才的式子：
$$Ans=\sum _{i=1}^{\lfloor \frac{N}{P}\rfloor }[gcd(i,\lfloor \frac{N}{P}\rfloor )==1]$$

基于容斥的思想，我们引入艾弗森约定：
$$\sum _{d|n}\mu (d)=(n==1)$$
假设：$T=\lfloor \frac{N}{P}\rfloor$

故原式可转化为：
$$Ans=\sum _{i=1}^T[gcd(i,T)==1]$$
$$=\sum _{i=1}^T\sum _{d|gcd(i,T)}\mu (d)$$
$$=\sum _{d|T}\mu (d)\ast \frac{T}{d}$$

由公式：
$$\sum _{d|T}\frac{\mu (d)}{d}=\frac{\phi (T)}{T}$$

最后推出：
$$Ans=\phi (T)$$
此题得解。

---

代码：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;

ll euler(ll x){
    ll res = x;
    for(int i=2 ;i*i<=x ;i++){
        if(x%i == 0){
            res = res/i*(i-1);
            while(x%i==0) x/=i;
        }
    }
    if(x>1) res = res/x*(x-1);
    return res;
}

int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        ll n,m;
        scanf("%I64d%I64d",&n,&m);

        ll ans = 0;
        for(ll i=1 ;i*i<=n ;i++){
            if(n%i == 0){
                if(i >= m){
                    ans += euler(n/i);
                }
                if((n/i)>=m && n/i != i){
                    ans += euler(i);
                }
            }
        }
        printf("%I64d\n",ans);
    }
    return 0;
}