欧拉定理相关及扩展欧几里得

威尔逊定理、费马定理、欧拉函数、欧拉定理、逆元、exgcd

威尔逊定理：$(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$，当且仅当$p$为质数。
其中： $(p-1)!$表示$p-1$的阶乘,即$1\ast 2\ast 3\ast \cdots \ast (p-1)$。
证明：
先证充分性： 即“$p$为质数”$\to$ $(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$。
设$0<i<p$,因为$p$为质数，可知$(i,p)=1$。根据裴蜀定理中的引理二，必存在一个整数$j(0<j<p)$,使得$i\ast j\%p=1$,称$j$是$i$的逆元，显然$i$的逆元具有唯一性。
而$i$的逆元是有可能等于$i$的。在哪些情况下，$i$的逆元等于$i$呢？
令$i\ast i\%p=1$，我们来求$i$ 的值。
由$$i\ast i\%p=1$$
可得：$$i\ast i-k\ast p=1$$
移项可得：
$$i\ast i-1=k\ast p$$
即
$$(i+1)\ast (i-1)=k\ast p$$
因为$p$为质数，所以：要么$(i+1)$是$p$的倍数，要么$(i-1)$是$p$的倍数。因为$0<i<p$,所以必有：
$$\{\begin{array}{l}i+1=p\\ i-1=0\end{array}$$
即$$\{\begin{array}{l}\begin{array}{rl}{i}_1& =p-1\\ i_2& =1\end{array}\end{array}$$
可知：只有当$i$为$1$或$p-1$时，$i$的逆元是它自身。
换句话说，在区间(1,p-1)中，每一个数都有一个不等于自己的唯一的逆元，并且逆元也在这个区间中。即这个区间中的数，可以划分成若干逆元对。
所以，$(p-1)!\%p=1\ast 2\ast 3\ast \cdots \ast (p-1)\%p=1\ast (p-1)\%p=p-1$.
而p-1和-1是模p等价的。所以充分性得证。
**再证必要性：**即 $(p-1)!\equiv -1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ $\to$ “$p$是质数”。
即证“p不是质数” $\to$ $(p-1)!\not\equiv p-1(\mathrm{m}\mathrm{o}\mathrm{d}p)$.
如果p不是质数，则[2,p-1]中必有p的因子，设其为$i$，则$i$和$p/i$都是$p$的因子。
1.如果$i\ne p/i$,则$1\ast 2\ast 3\ast \cdots \ast (p-1)$一定是p的倍数，所以它模$p$一定为0.
2.如果$i=p/i$,则则$1\ast 2\ast 3\ast \cdots \ast (p-1)$必为$i$的倍数，所以它模$p$必为$i$的倍数。而因为$p$是$i$的倍数，且$i>1$,所以(p-1)不可能是$i$的倍数。所以$(p-1)!\%p\ne (p-1)$
得证。

费马定理: 如果$p$为质数，且$a\%p\ne 0$,则有$a^{p-1}\%p=1$
证明：
$$(a\ast 1)\ast (a\ast 2)\ast (a\ast 3)\ast \cdots \ast (a\ast (p-1))\%p=a^{p-1}\ast (p-1)!\%p\text{\hspace{1em}(1)}$$
又因为：
a与p互质，根据裴蜀定理中的推论，$(a\ast 1)\%p,(a\ast 2)\%p,\dots ,(a\ast (p-1))\%p$的结果均互不相等，且在区间[1,p-1]之间，所以：
$$(a\ast 1)\ast (a\ast 2)\ast (a\ast 3)\ast \cdots \ast (a\ast (p-1))\%p=(p-1)!\%p\text{\hspace{1em}(2)}$$
由(1),(2)可知：
$$a^{p-1}\ast (p-1)!\%p=(p-1)!\%p\text{\hspace{1em}(3)}$$
因为$p$是质数，所以$(p-1)!$与$p$互质，所以等式(3)两边同时除以$(p-1)!$,得到：
$$a^{p-1}\%p=1$$
所以，费马定理得证。

欧拉函数

对一个正整数n，它的欧拉函数记为$\varphi (n)$，表示不超过n且与n互质的正整数的个数。
$\varphi (1)=1$

• 如果$n$为质数，则$\varphi (n)=n-1$.

• 如果$n=a^p$，其中$a$为质数，则$\varphi (n)=a^p-a^{p-1}=a^{p-1}\ast (a-1)=a^p\ast (1-1/a)$

• 令$n=a_1^{p_1}{a}_2^{p_2}\dots a_k^{p_k}$，根据积性函数的性质，

  $\begin{array}{rl}\varphi (n)& =\varphi (a_1^{p_1})\ast \varphi (a_2^{p_2})\ast \cdots \ast \varphi (a_k^{p_k})\\ & =a_1^{p_1}\ast (1-1/a_1)\ast a_2^{p_2}\ast (1-1/a_2)\ast \cdots \ast a_k^{p_k}\ast (1-1/a_k)\\ & =n\ast (1-1/a_1)\ast (1-1/a_2)\ast \cdots \ast (1-1/a_k)\end{array}$。

例1：给出一个整数n，求小于n且与n互质的数的个数。
n<=10000000
分析：本题即为求n的欧拉函数。直接通过上述公式来求。

#include<bits/stdc++.h>
using namespace std;
int n;
int getphi(int x)
{
    int res=x;
    for(int i=2;i*i<=x;i++)
    {
        if(x%i==0)
        {
            res=1ll*res*(i-1)/i; //此处要用1ll。否则有可能超出整型范围。
            while(x%i==0)
            x=x/i;
        }
    }
    if(x>1) //最后要注意判断x是否大于1。因为x中可能包含一个大于其平方根的质因数。
    res=1ll*res*(x-1)/x;
    return res;
}
int main()
{
    scanf("%d",&n);
    printf("%d\n",getphi(n));
    return 0;
}

例2.给出一个整数n，求区间[1,n]中所有数的欧拉函数值。
n<=10000000
用例1的方法求欧拉函数，时间复杂度为$O(\sqrt{n})$.如果本题也用这种方法，肯定超时。
我们可以采用类似于质数的筛法来求欧拉函数。
算法时间复杂度为$O(nlogn)$

#include<bits/stdc++.h>
using namespace std;
#define MAXN 1000005
int phi[MAXN];
int n,a,b;
void getphi(int x)
{
    for(int i=2;i<=x;i++)
    {
        if(phi[i]==0) //i一定是质数。
        {
            phi[i]=i-1;
            for(int k=2*i;k<=x;k+=i)
            {
                if(phi[k]==0)phi[k]=k;
                phi[k]=phi[k]/i*(i-1); //每个数都被它的质因子筛一下。
            }
        }
    }
}
int main()
{
    scanf("%d%d",&a,&b);
    if(a>b)swap(a,b);
    getphi(b);
    for(int i=a;i<=b;i++)
    printf("%d\n",phi[i]);
    return 0;
}

欧拉定理

如果$a$与$p$互质，则$a^{\varphi (p)}\%p=1$
其中$\varphi (p)$表示$p$的欧拉函数。
证明：
设$p$的简化剩余系为 { p 1 , p 2 , … , p k } \{p_1,p_2,\dots,p_k\} {p1​,p2​,…,pk​},因为$a$与$p$互质，所以 { a ∗ p 1 % p , a ∗ p 2 % p , … , a ∗ p k % p } \{a*p_1\%p,a*p_2\%p,\dots,a*p_k\%p\} {a∗p1​%p,a∗p2​%p,…,a∗pk​%p}也构成了$p$的简化剩余系。
所以:$$(a\ast p_1)\ast (a\ast p_2)\ast \cdots \ast (a\ast p_k)\equiv (p_1\ast p_2\ast \cdots \ast p_k)(\mathrm{m}\mathrm{o}\mathrm{d}p)$$
所以$$a^{\varphi }(p)\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

逆元

如果整数a，b满足$a\ast b\%p=1$,则称$a,b$在模p意义下互为逆元。
只有当a与p互质时，在模p意义下a才有逆元，$a$的逆元记作$a^{-1}$.
逆元的性质：若b与p互质，则$a/b\%p=a\ast b^{-1}\%p$
由此可见，通过逆元，可以将除法取余变为乘法取余。
注意此处的“除法”，不是整除，可以理解为分式。
在模意义下，分式取余是一个整数,比如$2/3\%5$，它是有意义的。因为在模意义下，每一个数表示的是一个同余类。在模5的意义下，2表示的是所有除以5余数为2的整数。所以2表示的是 { 2 , 7 , 12 , … , } \{2,7,12,\dots,\} {2,7,12,…,}，于是$2/3\%5$和$12/3\%5$是一样的，结果为4.
这种分式取模，我们既可以采用同余等价类来解决，也可以采用逆元来求。采用逆元是更一般的方法。
3在模5意义下的逆元为7，因为$3\ast 7\%5=1$。
于是有：$2/3\%5=2\ast 3^{-1}\%5=2\ast 7\%5=4$.

区间逆元

求区间[1,n]的逆元，如果采用上述求逆元的方法做，时间复杂度为$O(nlogn)$
有一种递推的方法可以求区间的逆元，时间复杂度为$O(n)$。
设模数为m,$f[n]$表示$n$的逆元，其中1~(m-1)的逆元已经求出，则$f[n]=(-f[m\%n]\ast (m/n)\%m+m)\%m$.
这个原理是什么呢？
1.设$f[i]$表示$i$在模m下的逆元，则有 $f[i]=-f[m-i]$。
证明：因为$$f[i]\ast i\%m=1$$,
所以：
$$f[i]\ast (-i)\%m=-1$$
即:$$-f[i]\ast (m-i)\%m=1$$
所以：$$f[m-i]=-f[i]$$
2. $f[i]=k\ast f[k\ast i]$
证明：
因为
$$f[i\ast k]\ast (k\ast i)\%m=1$$
即$$(f[k\ast i]\ast k)\ast i=1$$
即$$f[i]=k\ast f[k\ast i]$$

扩展欧几里得

扩展欧几里得是用来求不定方程如：$ax+by=c$。其中，$a,b,c$为已知整数，其中$gcd(a,b)\mid c$。这样，不定方程才有解。
扩展欧几里得是用递归的思想来求解的。
我们先求方程$ax+by=gcd(a,b)。$求出来该方程的解，则原方程的解乘上系数$(c/gcd(a,b))$就可以了。
因为：$$ax+by=gcd(a,b)=gcd(b,a\%b)=b{x}^{\prime }+(a\%b)y^{\prime }$$
$$ax+by=b{x}^{\prime }+(a\%b)y^{\prime }=b{x}^{\prime }+(a-a/b\ast b)y^{\prime }$$
将a,b看做变量，移项并合并同类项,可得：
$$\{\begin{array}{l}x=y^{\prime }\\ y=x^{\prime }-y^{\prime }\ast (a/b)\end{array}$$
只需要求出$x^{\prime },y^{\prime }$，则可以求出$x,y$。
而求$x^{\prime },y^{\prime }$，可以继续递归。这样递归下去，$gcd(a,b)$中的参数b最终会变为0，此时gcd(a,0)=a.
于是有:$ax+by=gcd(a,0)=a$，可以求出$x=1,y=0$.这是最底层的$x,y$,然后一层层返回，就可以求出原始的$x$和$y$了。

void exgcd(int a,int b,int &d,int &x,int &y)
{
	if(!b)d=a,x=1,y=0;
	else
	{
		exgcd(b,a%b,d,y,x);
		y-=x*(a/b);
	}
}