同余(看了也没用)

同余

定义

若整数$a$和整数$b$除以正整数$m$的余数相同则称$a$,$b$ $mod$ $m$ 同余记作$a\equiv b(mod$ $m)$.

同余的一些小性质

同加性：若$a\equiv b(modm)$，则$a+c\equiv b+c(modm)$
同乘性：若$a\equiv b(modm)$，则$a\times c\equiv b\times c(modm)$，若$a\equiv b(modm),c\equiv d(modm)$，则$a\times c\equiv b\times d(modm)$
同幂性：若$a\equiv b(modm)$，则$a^c\equiv b^c(modm)$
$a\times b\%k=(a\%k)\times (b\%k)\%k$

费马小定理

若$p$是质数,则对于任意整数$a$,有$a^p\equiv a(mod$ $p)$.

扩展欧几里得

Bezout定理
对于任意整数$a$,$b$存在一对整数$x$,$y$满足$ax+by=gcd(a,b)$.
证明:
根据欧几里得定理,我们可以知道$gcd(a,b)=gcd(b,a$ $mod$ $b)$,假设有一组数$x,y$,
满足$b\ast x+(a-b\lfloor \frac{a}{b}\rfloor )y$;可以得到$a\ast y+b(x-\lfloor \frac{a}{b}\rfloor \ast y)$
令$x^{\prime }=y,y^{\prime }=(x-\lfloor \frac{a}{b}\rfloor \ast y)$,则可以得到$a{x}^{\prime }+b{y}^{\prime }=gcd(a,b)$.
可以往下推到$b=0$时,显然有$x=1,y=0$成立.
根据数学归纳法,可知定理成立.
证毕。
乘法逆元
可以用扩欧和费马小定理两种方式
代码如下

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=20001000;
ll n,p;
ll f[N],ni[N],ci[N];//jie存阶层，ni存逆元，ci存P-2次方； 
ll fang(ll a,ll b)
{
	ll sum=1;
	while(b)
	{
		a%=p;
		sum%=p;
		if(b&1) sum=(sum%p*(a%p))%p;
		b=b>>1;
		a=(a%p*(a%p))%p;
	}
	return sum%p;
}
int main()
{
	freopen("invele.in","r",stdin);
	freopen("invele.out","w",stdout);
	cin>>n>>p;
	f[0]=1;
	f[1]=1;
	for(ll i=2;i<=n;i++) 
	f[i]=(f[i-1]*i)%p;
	ci[n]=fang(f[n],p-2);
	for(ll i=n-1;i>=1;i--)
	ci[i]=(ci[i+1]*(i+1))%p;
	for(ll i=1;i<=n;i++)
	{
		ni[i]=(ci[i]*f[i-1])%p;
	}
	for(ll i=1;i<=n;i++)
	printf("%lld\n",ni[i]);
	return 0;
}

版子请忽略freopen.
下面是扩欧

typedef long long ll;
void exgcd(ll a,ll b,ll &x,ll &y,ll &d)
{
    if (!b)
    {
        d=a;x=1;y=0;
        return;
    }
    exgcd(b,a%b,x,y,d);
    ll temp=x;
    x=y;
    y=temp-(a/b)*y;
}
int main()
{
    ll a,b,x,y,d;
    cin>>a>>b;
    exgcd(a,b,x,y,d);
    x=(x%b+b)%b;
    cout<<x<<endl;
    return 0;
}

代码里输出的x就是逆元的值了.

中国剩余定理

又称中国余数定理或孙子定理
简单来说就是这样
$\{\begin{array}{c}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ ⋮\\ x\equiv a_n(mod{m}_n)\end{array}$

对于这样的一个方程组如果任意的两个$m$都互素，则一定有解
所以中国剩余定理就是解但$m$互素的情况

我们尝试来构造这样一组解

设$M=m_1\times m_2\times \cdots \times m_n,M_i=\frac{M}{m_i}$

$M_i$只有在模$m_i$时不为零，模其他$m$都是零

由于我们希望在模$m_i$时得到$a_i$所以乘上$a_i{t}_i$

其中$t_i$是$M_i$在模$m_i$时的逆元

所以最小正整数解就是$x={\sum }_{i=1}^n{a}_i{t}_i{M}_i$，通解就是$x+kM$

如果题目说明了的$m$都是素数则可以用费马小定理
代码

int crt()
{
	int ans=0,M=1;
	for(int i=1;i<=n;i++)
	 M*=m[i];
	for(int i=1;i<=n;i++)
	{
		int Mi=M/m[i],invers=ksm(Mi,m[i]-2,m[i]);
		ans=(ans+a[i]*invers*Mi)%M;
	}
	return ans;
}

也可以用扩展欧几里得

inline int crt()
{
	int M=1,ans=0;
	for(int i=1;i<=n;i++) 
	M*=m[i];
	for(int i=1;i<=n;i++)
	{
		LL Mi=M/m[i],x,y;
		exgcd(Mi,m[i],x,y);
		ans=((ans+a[i]*Mi*x)%M+M)%M;
	}
	return ans;
}

证明中国剩余定理一定有解

我们举这个同余方程为例

$\{\begin{array}{c}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\end{array}$

假设我们分别解出

$\{\begin{array}{c}{x}_1=a_1+k_1\times m_1\\ x_2=a_2+k_2\times m_2\end{array}$

因为 $x_1=x_2$

所以$a_1+k_1\times m_1=a_2+k_2\times m_2$

整理得$k_2\times m_2-k_1\times m_1=a_1-a_2$

这就是$ax+by=d$的形式

所以当$gcd(m_1,m_2)|a_1-a_2$ 时一定有解

又因为$m_1$与$m_2$互质

所以$gcd(m_1,m_2)=1$

显然中国剩余定理一定有解
.