【总结】noip数学汇总

最新推荐文章于 2021-01-28 22:55:34 发布

测试运算符

最新推荐文章于 2021-01-28 22:55:34 发布

阅读量1.4k

点赞数 8

分类专栏：总结 --------------数学-------------- 文章标签：数学 noip

本文链接：https://blog.csdn.net/qq_40515553/article/details/83783780

版权

--------------数学-------------- 同时被 2 个专栏收录

24 篇文章 0 订阅

订阅专栏

总结

17 篇文章 0 订阅

订阅专栏

noip临近，有点怕需要数学知识定理才能做的题

快速幂

用途

$O(\log t)$ 时间内求解 $x^t$

证明&过程

由于 $x^a\cdot x^b=x^{a+b}$ ，而且每个自然数 $t$ 都可以拆分为不超过 $log_2t$ 个二次幂的和，即 $t=\sum_{i=0} k_i2^i$

则 $x^t=\sum_{i=0}x^{k_i2^i}$

其中 $k_i$ 为 $0 / 1$ ， $x^{2^i}$ 可以递推求得

代码

int qpow(int A,int B){
	int res(1);while(B){
		if(B&1)res=1ll*res*A%p;
		A=1ll*A*A%p,B>>=1;
	}return res;
}

一些定理

威尔逊定理：对于任意质数 $p$ ，有 $(p-1)!\equiv p-1\pmod p$
费马小定理：对于任意质数 $p$ ，有 $a^{p-1}\equiv 1\pmod p$

扩展欧几里得

用途

求解形似 $p\cdot a+q\cdot b=\gcd(a,b)$ 方程的解

证明&过程

由于 $\gcd(a,b)=\gcd(b,a\bmod b)$ ，则有

$p\cdot a+q\cdot b=\gcd(a,b)=\gcd(b,a\bmod b)=p\cdot b+q\cdot (a-\frac ab\cdot b)$

$=q\cdot a+(p-\frac ab\cdot q)\cdot b$

递归求解即可

代码

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

线性求逆元

用途

在 $O (n)$ 时间内求解区间 $[1, n]$ 中所有数的逆元

证明&过程

求 $i$ 关于 $p$ 的逆元

一般的对于 $i = 1$ ，有 $i^{-1}\equiv 1\pmod p$

对于 $i > 1$ 的情况，设 $p=k\cdot i+r$

放在模意义下， $k\cdot i+r\equiv 0\pmod p$

同乘 $r^{-1}\cdot i^{-1}$ ，得 $k\cdot r^{-1}+i^{-1}\equiv 0\pmod p$

由上面的定义式，有 $k=\lfloor \cfrac pi\rfloor,r=p\bmod i$

则有 $i^{-1}\equiv -\lfloor \cfrac pi\rfloor \cdot (p\bmod i)^{-1}\pmod p$

后面的 $(p\bmod i)^{-1}$ 可以通过调用前面的数组得到，前面的部分相当于转移系数

代码

inv[1]=1;
for(int i=2;i<=n;++i)
	inv[i]=1ll*(p-p/i)*inv[p%i]%p;

中国剩余定理

用途

求解同余方程组

$\begin{cases} x\equiv a_1\pmod {p_1}\\ x\equiv a_i\pmod {p_i}\\ x\equiv a_n\pmod {p_n} \end{cases}$

证明&过程

将式子拆开为多个形似下面的式子

$\begin{cases} x\equiv 0\pmod {p_1}\\ x\equiv a_i\pmod {p_i}\\ x\equiv 0\pmod {p_n} \end{cases}$

则对于每个方程化为方程： $k\cdot \prod_{j\not =i}p_j+r\cdot p_i=1$ ，最后再乘上 $a_i$

扩欧求解，解出第 $i$ 个方程的解为 $t_i$ ，则将答案合并得 $Ans=\sum t_i \pmod {\prod p_i}$

代码

void CRT(){
	read(n);
	int lcm=1;
	for(int i=1;i<=n;++i){
		read(a[i]),read(p[i]);
	    a[i]%=p[i];lcm*=p[i];
	}
	int x,y,ans=0;
	for(int i=1;i<=n;++i){
		int kl=lcm/p[i],t;
	    exgcd(kl,p[i],x,y);
	    t=1ll*x*kl*a[i]%p;
	    ans=(ans+t)%lcm;
	}
	printf("%d\n",ans);
}

斐波那契数列

递推公式

$F_i=F_{i-1}+F_{i-2}$

通项公式

$F_i=\frac {\sqrt 5}5\cdot \big[(\frac {1+\sqrt 5}2)^n-(\frac {1-\sqrt 5}2)^n\big]$ （在这篇博客里有生成函数的推导）

性质

$\sum_{i=1}^nF_i=F_{n+2}-1$
$\sum_{i=1}^nF_i^2=F_n\cdot F_{n+1}$
$\sum_{i=1}^ni\cdot F_i=n\cdot F_{n+2}-F_{n+3}+2$
$\sum_{i=1}^nF_{2i-1}=F_{2n}$
$\sum_{i=1}^nF_{2i}=F_{2n+1}-1$

卡特兰数

公式

$f_n=\sum_{i=0}^{n-1}f_i\cdot f_{n-1-i}$
$f_n=\frac {f_{n-1}\cdot (4n-2)}{n+1}$
$f_n=\frac {\binom {2n}n}{n+1}$
$f_n=\binom {2n}n-\binom{2n}{n-1}$

模型

括号序列计数
二叉树形态计数
凸多边形三角切分计数

欧拉函数&定理

欧拉函数

定义&性质

$\phi (x)$ 表示小于 $x$ 的正整数中，与 $x$ 互质的数的个数

若 $x=\prod p_i^{k_i}$ ，则 $\phi (x)=x\cdot \prod \frac {p_i-1}{p_i}$

用途

在枚举 $\gcd$ 时有奇效

线性筛

phi[1]=1;
for(rg int i=2;i<N;++i){
	if(!vis[i])phi[pri[++top]=i]=i-1;
    for(rg int j=1;j<=top&&i*pri[j]<N;++j){
    	vis[i*pri[j]]=1;
        if(i%pri[j]==0){phi[i*pri[j]]=phi[i]*pri[j];break;}
        else phi[i*pri[j]]=phi[i]*(pri[j]-1);
    }
}

欧拉定理

定义

若 $a, p$ 互质（不要求为质数），则 $a^{\phi(p)}\equiv 1 \pmod p$

用途

①一种求逆元的方法

②降幂：

由于 $a^{\phi(p)}\equiv 1\pmod p$

$a^t\equiv \frac {a^t}{a^{\phi (p)}}\equiv a^{t\bmod \phi(p)}\pmod p$

错排

定义

将 $n$ 个物品重新排列，求使得每个物品不在原来位置上的方案数，用 $D_n$ 表示

$D_i=(n-1)(D_{i-1}+D_{i-2})$

伪の二维错排

考过一次，求二维平面中物品都不能在原来的行列上，直接计算 $D_n^2$ 即可

组合数学

定义式

$\binom nm=\frac {n!}{m!(n-m)!}$

性质

$\binom nm=\binom n{n-m}$

$\binom nm=\binom{n-1}{m-1}+\binom{n-1}m$

$m\binom nm=n\binom{n-1}{m-1}$

$\sum_{i=0}^r\binom {n+i}{i}=\binom {n+r+1}r$

$\binom nm\cdot \binom md=\binom nd\cdot \binom {n-d}{m-d}$

$\sum_{i=0}^n(-1)^i\binom ni=0$

$\sum_{i=0}^{\frac n2}\binom n{2i}=2^{n-1}$

$\sum_{i=0}^ni\cdot \binom ni=n\cdot 2^{n-1}$

$\sum_{i=0}^ni^2\cdot \binom ni=(n^2+n)\cdot 2^{n-2}$

$\sum_{i=0}^n\binom ni^2=\binom {2n}n$

$\frac 1{(n-i)!(i-t)!}=\binom {n-t}{n-i}\cdot \frac 1{(n-t)!}$

测试运算符

关注

8
点赞
踩
11

收藏

觉得还不错? 一键收藏
0
评论
【总结】noip数学汇总

noip临近，有点怕需要数学知识定理才能做的题快速幂用途O(log⁡t)O(\log t)O(logt)时间内求解xtx^txt证明&amp;amp;amp;amp;过程由于xa⋅xb=xa+bx^a\cdot x^b=x^{a+b}xa⋅xb=xa+b，而且每个自然数ttt都可以拆分为不超过log⁡2t\log_2tlog2t个二次幂的和，即t=∑i=0ki2it=\sum_{i=0} k_i2^it=∑i...
复制链接

扫一扫

专栏目录