数学杂记

记录一些琐碎的数学知识（也可能不是数学）。以下内容重点放在记结论，以便作者复习。慢慢补充吧，东西会比较杂。。

最大公约数

由公约数的基本性质可得：

inline int gcd(int a,int b){return b?gcd(b,a%b):a;}

扩展欧几里得

求解关于$x,y$的二元一次方程：$$ax+by=gcd(a,b)$$
当等式右端不是$gcd(a,b)$的倍数时，方程无解。
考虑求解：$$b{x}_1+(a\%b)y_1=gcd(b,a\%b)=gcd(a,b)$$
假设已经求出了$x_1,y_1$。转化：
$$b{x}_1+(a-\lfloor \frac{a}{b}\rfloor \ast b)y_1=gcd(a,b)$$
$$a\ast y_1+b\ast (x_1-\lfloor \frac{a}{b}\rfloor \ast y_1)=gcd(a,b)$$
于是有：$x=y_1,y=x_1-\lfloor \frac{a}{b}\rfloor \ast y_1$
当递归到最后一层，有：$a=gcd(a,b),b=0$,于是此时$x=1,y=0$。
注意一下代码中$x$和$y$的螺旋走位。

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

中国剩余定理

求解一组线性同余方程，其中满足模数两两互质。
其实就是对于每一组考虑消去其它方程的影响。
$$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ ⋮\\ x\equiv a_n(mod{m}_n)\end{array}$$
令$M=\prod m_i$，$T_i=\frac{M}{m_i}$，考虑求解$r_i$满足：
$$T_i\ast r_i\equiv 1(mod{m}_i)⟺a_i\ast T_i\ast r_i\equiv a_i(mod{m}_i)$$
而同时，这里有$a_i\ast T_i\ast r_i\equiv 0(mod{m}_k),k\ne i,k\in [1,n]$
于是$x$有解：$x_0={\sum }_{i=1}^n{a}_i\ast T_i\ast r_i$，通解有$x=x_0+k\ast M,k\in \mathbb{Z}$。
模板，比较坑，要龟速乘。

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=15;
int k,a[N],b[N],x,y,mul=1,ans=0;
inline void turn(int &x,int mod){x=(x%mod+mod)%mod;}
inline void exgcd(int a,int b,int &x,int &y){
	if(!b){x=1,y=0;return;}
	exgcd(b,a%b,y,x),y-=a/b*x;
}
inline int Mul(int a,int b,int ret=0){
	for(;b;b>>=1,a=(a+a)%mul)
		if(b&1) (ret+=a)%=mul;
	return ret;
}
signed main(){
	scanf("%lld",&k);
	for(int i=1;i<=k;++i) scanf("%lld",a+i);
	for(int i=1;i<=k;++i) scanf("%lld",b+i),mul*=b[i];
	for(int i=1;i<=k;++i){
		int tmp=mul/b[i];
		exgcd(tmp,b[i],x,y),turn(x,b[i]),turn(a[i],b[i]);
		ans+=Mul(Mul(a[i],tmp),x),turn(ans,mul);
	}cout<<ans<<'\n';
}

扩展中国剩余定理

同上，但不满足模数两两互质。此时$T_i\ast r_i\equiv 1(mod{m}_i)$可能会无解，因为没有互质的条件。
考虑两两合并：
$$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\end{array}⟺\{\begin{array}{l}x=a_1+p\ast m_1\\ x=a_2+q\ast m_2\end{array}$$
将右侧两式相减得到一个二元一次方程组的形式：
$$p\ast m_1-q\ast m_2=a_2-a_1①$$
令$x=\frac{p}{\frac{a_2-a_1}{gcd(m_1,m_2)}},y=\frac{-q}{\frac{a_2-a_1}{gcd(m_1,m_2)}}$，考虑求解：
$$x\ast m_1+y\ast m_2=gcd(m_1,m_2)②$$
首先，如果$gcd(m_1,m_2)\nmid a_2-a_1$，显然①式是无整数解的。否则一定有解。
利用$exgcd$把②式求解出来之后，得到$x,y$，然后$x$和$y$分别乘上$\frac{a_2-a_1}{gcd(m_1,m_2)}$得到$p,-q$，随便取$p$或$-q$代入上面大括号右侧式子（$p$和$-q$代出来是一样的）得到一个满足上述同余式的$x_0$。于是有：
$$\{\begin{array}{l}x\equiv a_1\equiv x_0(mod{m}_1)\\ x\equiv a_2\equiv x_0(mod{m}_2)\end{array}$$
接下来证明：
对于两个非零整数$a,b$，且有整数$x,y$，满足$ax+by=gcd(a,b)$（这是一定存在的）：如果$a|n,b|n$，则$lcm(a,b)|n$
记$a^{\prime }=\frac{a}{gcd(a,b)},b^{\prime }=\frac{b}{gcd(a,b)}$，那么有：
$$n=\frac{n(ax+by)}{gcd(a,b)}=n{a}^{\prime }x+n{b}^{\prime }y$$
又有
$$a|n,b|n⟹a{b}^{\prime }|n{b}^{\prime },a^{\prime }b|n{a}^{\prime }⟹a{b}^{\prime }|n{b}^{\prime }y,a^{\prime }b|n{a}^{\prime }x⟹lcm(a,b)|n{a}^{\prime }x+n{b}^{\prime }y$$
即：
$$lcm(a,b)|n$$

上面有：
$$m_1|x-x_0,m_2|x-x_0⟹lcm(m_1,m_2)|x-x_0$$
即：
$$x-x_0\equiv 0(modlcm(m_1,m_2))⟹x\equiv x_0(modlcm(m_1,m_2))$$

于是这样一个一个合并即可。
模板，注意代码实现里，①式是先把$m_1$，$m_2$和$a_2-a_1$的$gcd(m_1,m_2)$除掉。后面把$p$代入$x$的方程时，再把$m_1$的$gcd(m_1,m_2)$乘回去。不然要$WA$，可能是爆$longlong$。还有无解的时候不能直接$break$掉，因为还要读入。。

#include<cstdio>
#define int long long
#define re register
using namespace std;
inline int gcd(int a,int b){return (b)?gcd(b,a%b):a;}
//ax+by=gcd(a,b)
inline void exgcd(int a,int b,int &x,int &y){
	if(!b){x=1,y=0;return;}
	exgcd(b,a%b,y,x),y-=a/b*x;
}
int n,m1,a1,m2,a2,p,q;
signed main(){
	while(~scanf("%lld",&n)){
		scanf("%lld%lld",&m1,&a1);
		bool flag=1;
		for(int re i=2;i<=n;++i){
			scanf("%lld%lld",&m2,&a2);
			int g=gcd(m1,m2),lcm=m1*m2/g,A=a2-a1;
			if(A%g!=0){flag=0;continue;}
			m1/=g,m2/=g,A/=g,exgcd(m1,m2,p,q),p=(p*A%m2+m2)%m2;
			a1=(a1+p*m1*g)%lcm,m1=lcm;
		}
		if(!flag) puts("-1");
		else printf("%lld\n",a1);
	}
}

高斯消元

解决$n$元一次方程组。复杂度$O(n^3)$。
只要记住一个一个地消元即可，红色是用来消元的位置。
红色下方消为0，消完之后往上求解即可。
下图为用第$i$行第$i$个消元。第$j$行的每个系数整体要减去第$i$行的每个系数乘以$\frac{a[j][i]}{a[i][i]}$，然而$a[i][i]=0$怎么办？换成一个第$i$列非$0$的行即可。
还有一点，就是取对应列最大的一行作为主元来消元可以减小精度误差，洛谷的题解里似乎有解释。
如果选出的主元等于0，说明它下面的系数都为0，在有解的情况下，主元可取任意值，它前面的元素随它的变化而变化，有无穷多解。相当于是少了一个方程。
若某行的系数均为$0$而常数项非$0$，则无解。下面代码只判了是否是唯一解。模板传送门

#include<bits/stdc++.h>
using namespace std;
const int N=105;
const double eps=1e-8;
double a[N][N],ans[N];int n;
inline int sgn(double x){
	if(x>eps) return 1;
	if(x<-eps) return -1;
	return 0;
}
inline int Guass(){
	for(int i=1;i<=n;++i){
		int k=i;
		for(int j=i+1;j<=n;++j)
			if(fabs(a[k][i])<fabs(a[j][i])) k=j;
		swap(a[i],a[k]);
		if(!sgn(a[i][i])) return 0;
		for(int j=i+1;j<=n;++j)
			for(int k=i+1;k<=n+1;++k)
				a[j][k]-=a[j][i]/a[i][i]*a[i][k];
	}
	for(int i=n;i;--i){
		for(int j=i+1;j<=n;++j)
			a[i][n+1]-=ans[j]*a[i][j];
		ans[i]=a[i][n+1]/a[i][i];
	}
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;++i)
		for(int j=1;j<=n+1;++j)
			scanf("%lf",&a[i][j]);
	if(!Guass()) return puts("No Solution"),0;
	for(int i=1;i<=n;++i) printf("%.2lf\n",ans[i]);
}

欧拉定理&广义欧拉定理

欧拉定理：若$a,m$为正整数，且$gcd(a,m)=1$，那么有：
$$a^{\phi (m)}\equiv 1(modm)$$
广义欧拉定理：
$$a^b\equiv \{\begin{array}{l}{a}^{b\%\phi (m)},gcd(a,m)=1\\ a^b,gcd(a,m)\ne 1,b<\phi (m)\\ a^{b\%\phi (m)+\phi (m)},gcd(a,m)\ne 1,b⩾\phi (m)\end{array}$$
在快速幂的时候要仔细想一想指数对不对，毕竟指数不能直接用$m$取模。

原根

参考文章
参考《初等数论》$P241$
定义1：$m⩾1,(a,m)=1$，使：$a^d\equiv 1(modm)$成立的最小正整数$d$称为a对模m的指数（或阶），把它记做${\delta }_m(a)$.
定义2：当${\delta }_m(a)=\phi (m)$时，称$a$是模$m$的原根，简称m的原根。

由欧拉定理：$a^{\phi (m)}\equiv 1(modm)$，可知${\delta }_m(a)⩽\phi (m)$
如何求质数的原根？
一般来说原根不会太大，暴力枚举$[2,\lfloor \sqrt[4]{x}\rfloor ]$即可。
如何验证当前枚举的$a$是不是原根？
首先，循环节长度$s$一定满足$s|\phi (m)$。
那么只需要找出$\phi (m)$的所有质因数$p_1,p_2,\cdots ,p_k$。
如果：$\forall p_i,i\in [1,k]:a^{\frac{\phi (m)}{p_i}}\not\equiv 1(modm)$
则$a$是$m$的一个原根。附上模板

#include<bits/stdc++.h>
using namespace std;
int mod,tot=0,factor[31];
inline int mul(int x,int y){return 1ll*x*y%mod;}
inline int quickpow(int a,int b,int ret=1){for(;b;b>>=1,a=mul(a,a))if(b&1)ret=mul(ret,a);return ret;}
inline bool is_pri_root(int x){
	for(int i=1;i<=tot;++i) 
		if(quickpow(x,(mod-1)/factor[i])==1)
			return false;
	return true;
}
inline int find(int mod){
	int phi=mod-1;
	for(int i=2;i*i<=phi;++i)
		while(phi%i==0) factor[++tot]=i,phi/=i;
	if(phi!=1) factor[++tot]=phi;
	for(int i=1;;++i)
		if(is_pri_root(i)) return i;
}
int main(){
	scanf("%d",&mod),printf("%d\n",find(mod));
}

二项式反演

神犇666
$$f_n=\sum _{i=0}^n(-1{)}^i(\genfrac{}{}{0ex}{}{n}{i})g_i⟺g_n=\sum _{i=0}^n(-1{)}^i\left(\genfrac{}{}{0ex}{}{n}{i}\right)f_i$$
这个式子有很强的对称性。一个更常用的结论：
$$f_n=\sum _{i=0}^n\left(\genfrac{}{}{0ex}{}{n}{i}\right)g_i⟺g_n=\sum _{i=0}^n(-1{)}^{n-i}(\genfrac{}{}{0ex}{}{n}{i})f_i$$
可以用来解决$min\_max$容斥。

错排问题

例题
$a$数组为$1$~$n$的排列，$\forall a[i]!=i$的排列方法有多少种？

令$D[n]$为长度为$n$的排列的错排方案数，有递推公式：
$$D[n]=(n-1)\ast (D[n-1]+D[n-2])$$
证明：
考虑增量法求$D[n]$，假设已知$D[1]$~$D[n-1]$：
现在令所有元素有序排列好，即$a[i]=i$。
第$n$个元素可以选择前$n-1$个位置中的任意一个放进去，假设选了某个位置$k$，那么这个位置的元素$k$就被挤出来了。
如果$k$放在了位置$n$，现在就相当于是$n$和$k$交换了位置，剩下$n-2$个元素需要错排。即为$D[n-2]$种方案。
如果$k$不放在位置$n$，相当于是把$a[k]!=k$的限制变为了$a[n]!=k$，由于$n$的位置已经固定在$k$，所以现在问题转化为了$1$~$n-1$的错排问题，方案数为$D[n-1]$。这$n-1$个限制关系依然是一一对应的。

于是有：$$D[n]=(n-1)\ast (D[n-1]+D[n-2])$$

初始状态$D[1]=0,D[2]=1$

二次剩余

以下船部内容参考：
《初等数论》第四章§5：模为素数的二次剩余
假定$p$是奇素数，设$p\nmid a$，二次同余方程的一般形式是：
$$a{x}^2+bx+c\equiv 0(modp)\cdots (1)$$
由于$p$为奇素数，$p\nmid a$，所以$p\nmid 4a$，所以$(1)$式和同余方程：
$$4a(a{x}^2+bx+c)\equiv 0(modp)$$
的解相同。上式可写为：
$$(2ax+b{)}^2\equiv b^2-4ac(modp)$$
令$y=2ax+b$，于是有：$y\equiv 2ax+b(modp)$，且：
$$y^2\equiv b^2-4ac(modp)$$
显然$x$和$y$是一一对应的。于是问题转化为形如：
$$x^2\equiv d(modp)$$
的同余方程，也就是求解$y$。当$p\mid d$，仅有一解：$x\equiv 0(modp)$。
所以，下面恒假定$p\nmid d$。
定义1：设素数$p>2,d\in \mathbb{Z},p\nmid d$，若$x^2\equiv d(modp)$有解，称d是模p的二次剩余，否则称d是模p的二次非剩余
定理1：在模$p$的一个既约剩余系中，恰有$\frac{p-1}{2}$个模$p$的二次剩余，$\frac{p-1}{2}$个模$p$的二次非剩余。若$d$是模$p$的二次剩余，则$x^2\equiv d(modp)$的解数为$2$。
证明略。
定理2（Euler判别法）：设素数$p>2$，$p\nmid d$，那么：
$d$是模$p$的二次剩余的充分必要条件是：
$$d^{\frac{p-1}{2}}\equiv 1(modp)\cdots (3)$$
$d$是模$p$的非二次剩余的充分必要条件是：
$$d^{\frac{p-1}{2}}\equiv -1(modp)\cdots (4)$$

Wilson定理：设$p$是素数，$r_1,\cdots ,r_{p-1}$是模$p$的既约剩余系，有：
$$r_1\cdots r_{p-1}\equiv -1(modp)$$
特别地，有：
$$(p-1)!\equiv -1(modp)$$

下面证明$Euler$判别法：
首先证明$\forall d\in \mathbb{Z},p\nmid d$，$(3)$或$(4)$有且仅有一式成立（平方差）：
$$d^{p-1}\equiv 1(modp)\Rightarrow (d^{\frac{p-1}{2}}-1)(d^{\frac{p-1}{2}}+1)\equiv 0(modp)$$
所以不可能存在$(3)$、$(4)$都不成立的情况。
若两式都成立，则$p\mid (d^{\frac{p-1}{2}}-1),p\mid (d^{\frac{p-1}{2}}+1)$
$\Rightarrow p\mid (d^{\frac{p-1}{2}}+1)-(d^{\frac{p-1}{2}}-1)\Rightarrow p\mid 2$
由于$p$为奇素数，所以上述情况不可能成立，矛盾。
所以$(3)$或$(4)$有且仅有一式成立。
$(3)$的必要性：若$d$是模$p$的二次剩余，则有$x_0$使得：
$$x_0^2\equiv d(modp)\Rightarrow x_0^{p-1}\equiv d^{\frac{p-1}{2}}(modp)\Rightarrow 1\equiv d^{\frac{p-1}{2}}(modp)$$

$(3)$的充分性：若$(3)$成立，则必有$p\nmid d$.考虑一次同余方程：
$$ax\equiv d(modp)$$
令模$p$的既约剩余系为$S$。那么对于$\forall j\in S$，当$a=j$，必有唯一的$x=x_j\in S$使得方程成立。假设$d$不是模$p$的二次剩余：
那么必有$j\ne x_j$，否则若$j=x_j$，有$j\ast x_j\equiv d\equiv x_j^2(modp)$，推出$d$是模$p$的二次剩余，矛盾。
现在把$S$中的数按照$j,x_j$两两配对，可以配出$\frac{p-1}{2}$组：$j\ast x_j\equiv d(modp)$
于是有：$$(p-1)!\equiv d^{\frac{p-1}{2}}(modp)$$
由$Wilson$定理：
$$-1\equiv d^{\frac{p-1}{2}}(modp)$$
与式$(3)$矛盾。所以必有$j_0$使得$j_0=x_{j_0}$，所以$d$是模$p$的二次剩余。这就证明了充分性。由已证部分可以立即推出剩下的结论。

Cipolla算法

求解关于x的同余方程：
$$x^2\equiv n(modp)$$
$p$为奇素数。算法流程：
①首先找到$a$使得$a^2-n$为模$p$的二次非剩余。根据上面的定理1，在模$p$的一个既约剩余系中，恰有$\frac{p-1}{2}$个模$p$的二次剩余，所以如果随机取这个$a$，找到满足条件的$a$的期望次数为$2$。根据$Euler$判别法，${(a^2-n)}^{\frac{p-1}{2}}\equiv -1(modp)$时满足条件。
②令$\omega =\sqrt{a^2-n}$，求出$x_0=(a+\omega {)}^{\frac{p+1}{2}}$，答案有两个：$x1=x_0,x_2=p-x_0$。

几个需要特判的东西：
①$n=0$，则解为$x_1=x_2=0$。
②$p=2$，则返回$n$即可（$n$只有$0,1$两种取值，两个解相同）

其他情况都有两个解。

具体证明参考

洛谷模板现在写代码越写越丑了。。

#include<bits/stdc++.h>
#define re register
#define cs const

namespace IO{
	cs int Rlen=1<<22|1;
	char buf[Rlen],*p1,*p2;
	inline char gc(){return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;}
	template <typename T>
	inline T get(){
		char ch=gc();T x=0;
		while(!isdigit(ch)) ch=gc();
		while(isdigit(ch)) x=((x+(x<<2))<<1)+(ch^48),ch=gc();
		return x;
	}
	inline int gi(){return get<int>();}
}
using IO::gi;

int T,n,mod,NQ;

inline int mul(int x,int y){return 1ll*x*y%mod;}
inline int add(int x,int y){return x+y>=mod?x+y-mod:x+y;}
inline int dec(int x,int y){return x-y<0?x-y+mod:x-y;}
inline int quickpow(int a,int b,int ret=1){for(;b;b>>=1,a=mul(a,a))if(b&1)ret=mul(ret,a);return ret;}

struct plx{
	int x,y;
	plx(int X=0,int Y=0){x=X,y=Y;}
	friend inline plx operator*(cs plx &a,cs plx &b){
		return plx(add(mul(a.x,b.x),mul(mul(a.y,b.y),dec(mul(NQ,NQ),n))),add(mul(a.x,b.y),mul(a.y,b.x)));
	}
	friend inline plx operator^(plx a,int b){
		plx ret(1,0);
		for(;b;b>>=1,a=a*a) if(b&1)
			ret=ret*a;
		return ret;
	}
};

inline void solve(){
	n%=mod;
	if(!n){puts("0");return;}
	if(mod==2){printf("%d\n",n);return;}
	if(quickpow(n,(mod-1)>>1)==mod-1){puts("Hola!");return;}
	while(1){
		NQ=rand()%mod;
		if(quickpow(dec(mul(NQ,NQ),n),(mod-1)>>1)==mod-1)
			break;
	}
	int x0=(plx(NQ,1)^((mod+1)>>1)).x,x1=mod-x0;
	if(x0>x1) std::swap(x0,x1);
	printf("%d %d\n",x0,x1);
}

int main(){
//	freopen("5491.in","r",stdin);
	T=gi();
	while(T--) n=gi(),mod=gi(),solve();
}

积性函数的卷积

$$I\ast \mu =ϵ$$
$$\mu \ast id=\phi$$
$$I\ast id=\sigma$$
$$I\ast I=d$$
$$I\ast \phi =id$$
注意$id(x)=x,I(x)\equiv 1$，不要搞混了T_T

$1$~$n$中与$n$互质的数的和：
$$\sum _{i=1}^n[gcd(n,i)==1]\ast i=\frac{n\ast \phi (n)+[n==1]}{2}$$
当$n=1$或$n=2$时，上式和即为1.
当$n>2$时，$\phi (n)$为偶数。如果$gcd(n,a)=1$，则$gcd(n,n-a)=1$
也就是说这些与$n$互质的数两两配对凑出了$n$，于是就有了上式。