欧拉函数学习

欧拉函数用希腊字母φ表示,φ(N)表示N的欧拉函数.

欧拉函数的一些性质:

1.对于素数p, φ(p)=p-1，对于对两个素数p,q φ（pq）=pq-1

2.对于一个正整数N的素数幂分解N=P1^q1*P2^q2*...*Pn^qn.

   φ(N)=N*(1-1/P1)*(1-1/P2)*...*(1-1/Pn).

3.除了N=2,φ(N)都是偶数.

4.设N为正整数,∑φ(d)=N (d|N).

那么我们是不是能在筛素数的同时求出所有数的欧拉函数呢.

答案是可以.

φ(n)=n*（1-1/p1)(1-1/p2)....(1-1/pk),其中p1、p2…pk为n的所有素因子。

比如：φ(12)=12*(1-1/2)(1-1/3)=4。
利用这个就比较好求了，可以用类似求素数的筛法。
先筛出N以内的所有素数，再以素数筛每个数的φ值。
比如求10以内所有数的φ值：
设一数组phi[11]，赋初值phi[1]=1,phi[2]=2...phi[10]=10；
然后从2开始循环，把2的倍数的φ值*(1-1/2)，则phi[2]=2*1/2=1,phi[4]=4*1/2=2,phi[6]=6*1/2=3....；
再是3，3的倍数的φ值*(1-1/3)，则phi[3]=3*2/3=2,phi[6]=3*2/3=2，phi[9]=.....；
再5，再7...因为对每个素数都进行如此操作，因此任何一个n都得到了φ(n)=n*（1-1/p1)(1-1/p2)....(1-1/pk)的运算
觉得这个“筛”还是比较好用的，以前求数的所有因子之和也是用的它。
代码如下：

<span style="font-size:18px;"><span style="font-size:18px;">void Init(){       
     euler[1]=1;      
     for(int i=2;i<Max;i++)      
       euler[i]=i;      
     for(int i=2;i<Max;i++)      
        if(euler[i]==i)      
           for(int j=i;j<Max;j+=i)      
              euler[j]=euler[j]/i*(i-1);//先进行除法是为了防止中间数据的溢出       
}  </span></span>   

文出自 ： http://blog.csdn.net/sentimental_dog/article/details/52002608

之前比赛一道题，暴力过了，其实原理是欧拉函数，当时如果我知道了我也分析不出来，现在知道思路自己模拟一下子

POJ - 3090

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<cmath>
using namespace std;
int n, m;
int phi[1005];
int main()
{
    phi[1] = 1;//边打表边筛选
    for(int i = 2; i < 1005; i ++)phi[i] = i;
    for(int i = 2; i < 1005; i ++)
        if(phi[i] == i)
        for(int j = i; j < 1005; j += i)
            phi[j] = phi[j] / i * (i - 1);

    while (scanf("%d", &n) != EOF)
    {
        phi[1] = 1;
        for (int i = 1; i <= n; ++ i)
        {
            scanf("%d", &m);
            int sum = 0;
            for (int j = 1; j <= m; ++ j )
            {
                sum += phi[j];
            }
            printf("%d %d %d\n", i, m, 2 * sum + 1);
        }
    }
    return 0;
}