求逆元

我们来看一点清真的东西：
<关于逆元>
一． p特别小，b也特别小（即$p*b < long int$）
$(a/b)\mod p = a \mod (b*p)/b;$

---

二． p为质数
费马小定理
$a^{p-2}$是$a$的逆元

---

三． p不为质数
$a^{\phi(p)} \mod p = 1$;
$a^{\phi(p)-1}$ 为$p$的逆元
1. 求一次$\phi(p)$
利用公式$\phi(p)=p*(p_1-1/p_1)*(p_2-1/p_2)*……*(p_i-1/p_i);$
其中$p_i$为$p$的质因子。
2. 若$\phi(p)$有多次变化
$O(p)$ 线性筛 $\phi(1-p)$

---

四．线性筛欧拉函数
在求解素数的同时求出欧拉函数。
利用性质：
1.p,q互质时：$\phi(p*q) = \phi(p)*\phi(q);$
2.p为素数：$\phi(p)=p-1;$
3.当p为素数，如果i % p = 0：$\phi(i*p)=\phi(i)*p$

int p,cnt;
bool mark[MAXN];
int prime[MAXN],phi[MAXN];
int main()
{
    scanf("%d",&p);
    for(int i=2; i<=p; i++)
    {
        if(!mark[i])
        {
            prime[cnt++]=i;
            phi[i]=i-1;
        }
        for(int j=0; j<cnt && prime[j]*i<=p; j++)
        {
            mark[i*prime[j]]=1;
            if(!(i%prime[j]))
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
                /*phi[prime[j]]=prime[j]-1;*/
        }
    }
}

---

五．p为质数
$\phi(p^k) = \phi(p)*p^{k-1} = (p-1)*p^{k-1};$

---

六．$a^{\phi(p)} \mod p = 1$
$a^b \mod p = a^{b\mod\phi(p)} \mod p$

---

七．Exgcd - 扩展gcd
$a*x = 1(\mod b)$
$a*x+b*y=1;$
$a*x+b*y=gcd(a,b)=gcd(b,a \mod b);$
$x=y',y=x'-a/b*y';$