欧拉函数+中国剩余定理入门

欧拉函数

互质：

互质是公约数只有1的两个整数，叫做互质整数。公约数只有1的两个自然数，叫做互质自然数，后者是前者的特殊情形

• 5和5不互质，因为5和5的公因数有1、5。
• 1和任何数都成倍数关系，但和任何数都互质。因为1的因数只有1，而互质数的原则是：只要两数的公因数只有1时，就说两数是互质数。因为1只有一个因数所以1既不是质数（素数），也不是合数，无法再找到1和其他数的别的公因数了。1和-1与所有整数互素，而且它们是唯一与0互素的整数。
• 互质数的写法：如c与m互质，则写作$(c,m)=1$

欧拉函数介绍：

• 在数论，对正整数n，欧拉函数是小于或等于n的正整数中与n互质的数的数目（因此φ(1)=1）。
• 此函数以其首名研究者欧拉命名(Euler’s totient function)，它又称为Euler’s totient function、φ函数、欧拉商数等。 例如φ(8)=4，因为1,3,5,7均和8互质。
• 从欧拉函数引伸出来在环论方面的事实和拉格朗日定理构成了欧拉定理的证明。

公式1：
$$\begin{array}{r}\phi \left(x\right)=x\prod _{i=1}^n\left(1-\frac{1}{p_i}\right)\\ \left(其中p_1,p_2\dots \dots p_n为x的所有质因数，x是不为0的整数\right)\end{array}$$
公式2：
$$\begin{array}{r}\phi \left(n\right)=n-1\\ \left(n为质数\right)\end{array}$$

//常规求法
ll eular(ll n){
    ll ans = n;
    for(ll i=2; i<=sqrt(n); i++){
        if(n%i==0){
            ans = ans/i*(i-1);
            while(n%i==0)n/=i;
        }
    }
    if(n>1)ans = ans/n*(n-1);
    return ans;
}

//记忆化求法
void euler(){  
    for(int i=2;i<maxn;i++){  
        if(!Eu[i])
        for(int j=i;j<maxn;j+=i){  
            if(!Eu[j])Eu[j]=j;  
            Eu[j]=Eu[j]/i*(i-1);  
        }  
    }  
}

例题：

GCD（HDU-2588）

题意：

求$x<=n$且 $gcd(x，n)>=m$ 的 $x$的个数

思路：
$$\gcd \left(x,N\right)\ge M\iff \gcd \left(\frac{x}{M},\frac{N}{M}\right)\ge 1$$
$$\begin{array}{r}\gcd \left(\frac{x}{M},\frac{N}{M}\right)\ge 1\iff \gcd \left(\frac{x}{\alpha },\frac{N}{\alpha }\right)=1\\ \left(\frac{\alpha }{M}\ge 1,\alpha 没必要是M的倍数\right)\end{array}$$
先找出所有满足条件的$\alpha$(首先是n的因数,然后比m大)值，之后利用欧拉函数，求出$\frac{x}{\alpha }$的个数，最后累加个数即可

代码：

#include<iostream>
#include<stdio.h>
#include<cmath>
#define ll long long
using namespace std;
ll eular(ll n){
    ll ans = n;
    int g=(int)sqrt(n+0.5);
    for(ll i=2; i<=g; i++){
        if(n%i==0){
            ans = ans/i*(i-1);
            while(n%i==0)n/=i;
        }
    }
    if(n>1)ans = ans/n*(n-1);
    return ans;
}
int main(){
    int t;
    scanf("%d",&t);
    ll n,m;
    while(t--){
        scanf("%lld%lld",&n,&m);
        ll ans=0;
        for(ll i=1;i*i<=n;i++){
            if(n%i==0){//找出满足条件的α(i)值
                if(i>=m) ans+=eular(n/i);
                if(i*i!=n&&n/i>=m) ans+=eular(i);
            }
        }
        printf("%lld\n",ans);
    }
}

中国剩余定理

解如下方程：
$$\left(S\right):\{\begin{array}{l}x\equiv a_1\left(mod{m}_1\right)\\ x\equiv a_2\left(mod{m}_2\right)\\ ⋮\\ x\equiv a_n\left(mod{m}_n\right)\end{array}$$
先取$M={\prod }_{i=1}^n{m}_i$ ，$M_i=\frac{M}{m_i}$

之后令$t_i=M_i^{-1}$为$M_i$模$m_i$的数论倒数($t_i$为$M_i$模$m_i$意义下的逆元)

那么方程组$(S)$的通解形式为
$$\begin{array}{r}x=a_1{t}_1{M}_1+a_2{t}_2{M}_2+\dots +a_n{t}_n{M}_n+kM\\ =kM+\sum _{i=1}^n{a}_i{t}_i{M}_i,k\in Z\end{array}$$

在模$M$意义下，方程组$(S)$只有一个解：
$$x=\left(\sum _{i=1}^n{a}_i{t}_i{M}_i\right)modM$$

例题：

曹冲养猪

代码：

#include<iostream>
#include<stdio.h>
using namespace std;
typedef long long ll;
int a[15],b[15];
ll m=1,ans=0,p,q;
void exgcd1(ll a,ll b,ll &x,ll &y)
{
    if(b==0){x=1,y=0;return ;}
    exgcd1(b,a%b,y,x);
    y-=x*(a/b);
}//无返回值版 
ll exgcd2(ll a,ll b,ll &x,ll &y)
{
    if(!b){x=1, y=0; return a;}
    int gc=exgcd2(b,a%b,y,x);
    y-=a/b*x;
    return gc;
}//递归求ax+by=gcd(a,b),有返回值版 
int main()
{
    int n;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d%d",&a[i],&b[i]),m*=a[i];
    for(int i=1;i<=n;i++)
    {
        ll x=m/a[i];
        exgcd1(x,a[i],p,q);
        p=(p%a[i]+a[i])%a[i];//p的最小正整数解
        ans+=((p*b[i]*x)%m+m)%m;
        ans%=m;
 
    }
    printf("%lld",ans);
    return 0;
}