hdu 2824：The Euler function（欧拉函数）

hdu 2824：The Euler function（欧拉函数）

---

Problem Description

---

The Euler function phi is an important kind of function in number theory, (n) represents the amount of the numbers which are smaller than n and coprime to n, and this function has a lot of beautiful characteristics. Here comes a very easy question: suppose you are given a, b, try to calculate (a)+ (a+1)+…+ (b)

---

Input

---

There are several test cases. Each line has two integers a, b (2<a<b<3000000).

---

Output

---

Output the result of (a)+ (a+1)+…+ (b)

---

Sample Input

---

3 100

---

Sample Output

---

3042

题解：典型的欧拉函数模板题
给定区间a,b要求a到b的欧拉函数对应的值的和

于是翻开了欧拉函数的word文档：
欧拉函数的定义:
在数论中，对于正整数N,少于或等于N ([1,N]),且与N互质的正整数(包括1)的个数，记作φ(n)。
φ函数的值：
φ(x)=x(1-1/p(1))(1-1/p(2))(1-1/p(3))(1-1/p(4))……(1-1/p(n)) 其中p(1),p(2)…p(n)为x的所有质因数;x是正整数; φ(1)=1(唯一和1互质的数，且小于等于1)。注意：每种质因数只有一个。

例如:
φ(10)=10×(1-1/2)×(1-1/5)=4;
φ(30)=30×(1-1/2)×(1-1/3)×(1-1/5)=8;
φ(49)=49×(1-1/7)=42;

于是找到了两个模板：

欧拉函数模板
(1)直接求小于或等于n,且与n互质的个数:

int Euler(int n)
{
    int ret=n;
    for(int i=2; i<=sqrt(n); i++)
        if(n%i==0)
        {
            ret=ret/i*(i-1);                   //先进行除法防止溢出(ret=ret*(1-1/p(i)))
            while(n%i==0)
                n/=i;
        }
    if(n>1)
        ret=ret/n*(n-1);
    return ret;
}

（2）筛选模板:求[1,n]之间每个数的质因数的个数

#define size 1000001
int euler[size];
void Init()
{
     memset(euler,0,sizeof(euler));
          euler[1]=1;
     for(int i=2;i<size;i++)
       if(!euler[i])
       for(int j=i;j<size;j+=i)
       {
              if(!euler[j])
               euler[j]=j;
               euler[j]=euler[j]/i*(i-1);//先进行除法是为了防止中间数据的溢出
         }
}

开始打算用第一个，结果TLE，于是就用的第二个模板（注意long long）

代码一：

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#define ll long long
using namespace std;
const int size=3000000;
ll euler[size+10];
void init()
{
    memset(euler,0,sizeof(euler));
    euler[1]=1;
    for(ll i=2;i<=size;i++)
    {
        if(euler[i]==0)
        {
            for(ll j=i;j<=size;j+=i)
            {
                if(euler[j]==0)
                {
                    euler[j]=j;
                }
                euler[j]=euler[j]/i*(i-1);
                //先进行除法是为了防止中间数据的溢出
            }
        }
    }
    /*for(ll i=3;i<=size;i++)
    {
        euler[i]+=vis[i-1];
    }*/
}
int main()
{
    init();          //初始化，首先求出每个数对应的质因数个数
    ll a,b;
    while(~scanf("%lld%lld",&a,&b))
    {
        ll sum=0;
        for(ll i=a;i<=b;i++)
            sum+=euler[i];
        printf("%lld\n",sum);
    }
    return 0;
}

代码二： 运用了前缀和思想

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#define ll long long
using namespace std;
const int size=3000000;
ll euler[size+10];
void init()
{
    memset(euler,0,sizeof(euler));
    euler[1]=1;
    for(ll i=2;i<=size;i++)
    {
        if(euler[i]==0)
        {
            for(ll j=i;j<=size;j+=i)
            {
                if(euler[j]==0)
                {
                    euler[j]=j;
                }
                euler[j]=euler[j]/i*(i-1);
                //先进行除法是为了防止中间数据的溢出
            }
        }
    }
    for(ll i=3;i<=size;i++)     //前缀和
    {
        euler[i]+=euler[i-1];
    }
}
int main()
{
    init();          //初始化，首先求出每个数对应的质因数个数
    ll a,b;
    while(~scanf("%lld%lld",&a,&b))
    {
        printf("%lld\n",euler[b]-euler[a-1]);
    }
    return 0;
}