欧拉函数就是指:对于一个正整数n,小于n且和n互质的正整数(包括1)的个数,记作φ(n) 。
欧拉函数的通式:φ(n)=n*(1-1/p1)(1-1/p2)(1-1/p3)*(1-1/p4)……(1-1/pn),其中p1, p2……pn为n的所有质因数,n是不为0的整数。φ(1)=1(唯一和1互质的数就是1本身)。
晦涩难懂一:代码中的 ans-=ans/i; ans = ans(1 - 1/n) 这一步就是对应欧拉函数的通式~还是应该比较容易看懂的吧?
晦涩难懂二:while(n%i==0) n/=i; 这一个语句是为了保证完全消除我们刚才得到的那个i因子。确保我们下一个得到的i是n的素因子。
晦涩难懂三:if(n>1)ans-=ans/n; 这个语句是为了保证我们已经除完了n的所有的素因子,有可能还会出现一个我们未除的因子,如果结尾出现n>1 ,说明我们还剩一个素因子木有除。
①N>1,不大于N且和N互素的所有正整数的和是 1/2Neular(N)。 推荐题目:HDOJ 3501
②若(N%a == 0 && (N/a)%a==0) 则有:E(N)=E(N/a) x a;
③若(N%a==0 && (N/a)%a!=0) 则有:E(N)=E(N/a)(a-1);
#include<bits/stdc++.h>
using namespace std;
const int N = 3e6 + 1;
long long int E[N];
int euler(int n)
{
int ans = n;
for(int i = 2; i * i <= n; i++)
{
if(n % i == 0)
{
ans -= ans / i;
while(n % i == 0)
{
n /= i;
}
}
}
if(n > 1) ans -= ans / n;
return ans;
}
//打表求欧拉函数
void get_euler()
{
E[1] = 1;
for(int i = 2; i < N; i++) E[i] = i;
for(int i = 2; i < N: i++)
{
if(E[i] == i)
{
for(int j = i; j < N; j += i)
E[j] -= E[j] / i;
}
}
}
int main()
{
int n;
get_euler();
int m;
while(cin >> n >> m&& n != 0)
{
cout << E[m] - E[n - 1] << endl;
}
return 0;
}