2018.7.12（初等数论+组合数学）

最新推荐文章于 2021-05-25 17:55:05 发布

1035719430

最新推荐文章于 2021-05-25 17:55:05 发布

阅读量1.5k

点赞数

分类专栏：组合数学数论

本文链接：https://blog.csdn.net/a1035719430/article/details/81021651

版权

数论同时被 2 个专栏收录

12 篇文章 0 订阅

订阅专栏

组合数学

6 篇文章 0 订阅

订阅专栏

初等数论

素数

大家都会，略过
就是 $O（n）$ 的线性筛

素数定理

唯一分解定理

略过

威尔逊定理

$p$ 为素数，则 $(p-1)! ≡ p-1 (mod$ $p)$ 同时，它的逆定理也成立。

费马小定理

$p$ 为素数， $a$ 为正整数， $a$ 与 $p$ 互质，则 $ap-1≡ 1$ ( $mod$ $p$ ) 。

欧拉小定理

$a^{φ(m)} ≡ 1 ($ mod $m)$ ，( $a$ 与 $m$ 互质)

欧几里得算法

$gcd(a,b)=gcd(b,a$ $mod$ $b)$

扩展欧几里得算法

求解形如ax+by=gcd(a,b)方程的一组解。
当b=0时，原方程的解就是x=1，y=0。
假设 $bx’+(a$ $mod$ $b)y’=gcd(a,b)$ 有解。
$bx' + (a-\lfloor \frac{a}{b} \rfloor b)y' = gcd(a,b)$
$ay'+b(x'-\lfloor \frac{a}{b}\rfloor y')=gcd(a,b)$
所以我们也找到了原方程的一组解，按照欧几里得算法那样递归求解即可

裴蜀定理

若 $a,b$ 是整数,且 $gcd(a,b)=d$ ，那么对于任意的整数 $x,y$ , $ax+by$ 都一定是d的倍数，特别地，一定存在整数 $x,y$ ，使 $ax+by=d$ 成立。
推广:n个整数， $a_1,a_2,a_3......a_n$ 为n个整数， $d$ 是它们的最大公约数，那么存在整数 $x_1......x_n$ 使得 $x_1*a_1+x_2*a_2+...x_n*a_n=d$

逆元

$O(n)$ 递推
$O(logn)$ 求单个数的逆元

中国剩余定理

很显然略过

Baby Step Giant Step(BSGS)

用于求解方程 $A^x mod P = B$ （P为质数）
令 $m = \lfloor \sqrt{p} \rfloor$
$x = i * m + j$
原式转化为 $A^{im + j}≡B$ （ $mod$ $p$ ）
即 $A^{im}≡B（A^{-1}）^j$
左边枚举 $i$ 丢进 $hash$ 表
右边枚举 $j$ 去与 $hash$ 表对应

扩展BSGS

显而易见的，BSGS需要求逆元，也就是要求 $gcd(A,P)=1$
也就是说， $gcd(A,P)≠1$ 时，不能用BSGS。
对于方程 $A^B mod P = C$ ，假设 $d=gcd(A,P)$ ，若 $d|C$ 不成立（特殊条件 $B=1$ 暂不考虑），方程无解。
改写原方程
$\frac{A}{d} × A^{B-1}≡\frac{C}{d} ($ $mod$ $\frac{P}{d} )$
后面任然用BSGS

大数分解

大数分解：将n分解为质因数的积的形式（n可能很大，O(√n)算法无法完成）
大数分解给人的第一感觉：玄学
$Pollard Rho$ 算法：
1.找到一个数p，使得 $p|n$ ，将 $n$ 分解为 $p$ 与 $n/p$ ；
2.如果 $p$ 或 $n/p$ 不是素数，继续递归；
3.如果是素数，记录并退出。

大数分解有两个关键：
第一个是如何判定一个数是不是素数。
这里我们用到了 $Miller Rabin$ 算法：
设我们要判定的数是 $p$ ，如果存在一个底数 $a$ ，使得 $a^{p-1}≡1($ $mod$ $p)$ 不成立，那么 $p$ 肯定不是质数。
随机几次，出错率就比较低了。
为了使出错率更低，还有二次探测的方法。
如果 $p$ 是奇素数，则 $x^2≡1($ $mod$ $p)$ 的解为 $x=1$ 或 $x=p-1($ $mod$ $p)$ 。
对于奇素数 $r×2s+1$ （ $r$ 为奇数），一般对 $a^r$ ， $a^{r*2}$ ，…， $a^{(r*2s-1)}$ 进行二次探测。

另一个是如何找 $p$ 使得 $p|n$ ：
1.找到一个数 $p_1$ ；
2.通过某种玄学手段得到与 $p_1$ 对应的一个数 $p_2$ ；
3.判断 $|p_2-p1|$ 是否能整除 $n$ ，即 $gcd(|p_2-p_1|,n)$ 是否大于1；
4.如果大于 $1$ ， $gcd(|p_2-p_1|,n)$ 即为所求，否则 $p_2$ 作为一个新的 $p_1$ ，继续重复上述过程。
实际使用中，一般用如下公式进行推导： $p_2=(p_1^2+c) mod$ $n$ ，其中c为随机常数。这个公式能保证在形成环之前， $p_2-p_1$ 基本不会相等。
如何判断是否形成环：标记i位于1，标记j位于2，标记i随着过程不断++，碰到标记j后，标记j的位置*2，并记录当前标记i处的值，这样能保证当j-i>环的大小时，一定能找到环。碰到环了就直接退出，换一个随机常数再找。

一些例题

压惊水题

给定一个长度为n的数组，求ai|aj这样的数对(i,j)有多少个？
$n,a_i<=10^6$
其实只要枚举每一个数组中的数及其他们的倍数就行了
$O(nlogn)$

BZOJ1477 青蛙的约会

exgcd裸题
$x+ans*m ≡ y+ans*n (mod$ $L)$
$x-y ≡ ans*(n-m) (mod$ $L)$
剩下扩欧好了

#include<bits/stdc++.h>
using namespace std;
long long x,b,a,y,m,n,l,t,d,ans;
int ex_gcd(long long &x,long long &y,long long a,long long b)
{
    if(!b) return x=d/a,y=0;
    ex_gcd(x,y,b,a%b),t=x;
    return x=y,y=t-a/b*y;
}
int main() 
{
    scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l);
    d=(y-x+l)%l;
    a=(m-n+l)%l;
    b=l;
    long long a1=a,b1=b;
     while(b1)
     {
       t=b1;
       b1=a1%b1;
       a1=t;    
     }
    if(d%a1==0)
     {
      ex_gcd(ans,y,a,b);
      while(ans<0) ans=(ans+l)%l;
      printf("%lld",ans%l);
     }
    else 
     printf("Impossible");
    return 0;
}

BZOJ2257 瓶子和饮料

给定n个瓶子的容量，选出其中k个，最大化k个瓶子之间无度量所能表达的最小体积。
$k≤n≤1000，容量≤10^9$
实际上是答案就是 $k$ 个瓶子体积的 $gcd$
枚举所有的质因子判断一下是否出现次数超过 $k$ 即可
$O(n\sqrt{max(V)})$
https://blog.csdn.net/a1035719430/article/details/81055790

BZOJ2818 GCD

~~当然是莫比乌斯反演啦~~
好啦，N这么小怎么可能毒瘤莫比乌斯呢
莫比乌斯可以优化我们这题的做法做到更优
其实我们只要枚举 $gcd = p$ ，然后 $(x,y)$ 可以等价到 $(\frac{x}{gcd(x,y)} , \frac{y}{gcd(x,y)})$
也就是我们只要在 $[1,\frac{N}{p}]$ 之间直接统计答案即可

BZOJ4173数学

这题还真是鬼畜….
~~我们也可以鬼畜，打表出答案~~ $φ(n)*φ(m)*n*m$
我们考虑证明
第二个 $Σ$ 中的( $n$ $mod$ $k + m$ $mod$ $k ≥ k$ )等价于
$\lfloor \frac{n+m}{k} \rfloor - \lfloor \frac{n}{k} \rfloor - \lfloor \frac{m}{k}\rfloor = 1$
将原式转化
这里写图片描述
那么原式答案就为 $φ(n)*φ(m)*n*m$

其实打表是很好的，只要打表后面那个 $Σ$ 就可以愉快AC了

bzoj3122随机数生成器

如果高中数列你学的没问题你可以很轻松的写出通项公式
蓝后你就可以得到形如： $d=(bc+t)*(x1+bc)^{-1}$ 的式子了
BSGS！

组合数学

加法原理&乘法原理

略过

排列组合

公式略过
经典例题
某保密装置须同时使用若干把不同的钥匙才能打开。现有７人，每人持若干钥匙。须４人到场，所备钥匙才能开锁。问①至少有多少把不同的钥匙？②每人至少持几把钥匙？
解析：
①每３人至少缺１把钥匙，且每３人所缺钥匙不同。故至少共有C(7,3)=35把不同的钥匙。
② 任一人对于其他６人中的每３人，都至少有１把钥匙与之相配才能开锁。故每人至少持C(6,3)＝20把不同的钥匙。

可重复的排列

如果S是一个多重集,它有K个不同的类型元素，各元素分别为 $n1,n2,…,nk$ 个，那么，Ｓ的全排列个数为 $\frac{n!} { n1!*n2!*…*nk!}$

可重复的组合

设Ｓ是一个具有k种类型元素的多重集，每种元素均具有无限的重复数．则S的r-组合数为Ｃ(r+k-1,r).
例题简单略过

二项式定理

$(x+y)^n = \sum_{k = 0}^{n} C_{n}^{n-k} x^{n-k}y^k$
对照杨辉三角理解
特殊情况：
$(1+x)^n=\sum_{k=0}^{n} C_{n}^{n-k} x^k$

组合数性质

⑴ $C_{n}^{r}=C_{n}^{n-r}$
显然

⑵ $C_{n}^{r}=C_{n-1}^{r}+C_{n-1}^{r-1}$
显然

⑶ $\sum_{k = 0}^{r} C_{n + k}^{k} = C_{n + r + 1}^{r}$
略微证明：
由 $C_{n}^{0} = C_{n +1}^{0}$
则 $C_{n}^{0} + C_{n + 1}^{1} + C_{n + 2}^{2} + ………+C_{n+r}^{r}$
$=C_{n + 2}^{1} + C_{n+2}^{2} +…………+C_{n+r}^{r}$
$=C_{n+3}^{2} + …… +C_{n+r}^{r}$
$=C_{n + r + 1}^{r}$
证毕

⑷ $C_{n}^{l} × C_{l}^{r} = C_{n}^{r} * C_{n-r}^{l-r}$
显然

⑸ $\sum_{k = 0}^{n} C_{n}^{k} = 2^n$
一句话证明：将 $x=1$ 带入二项式定理

⑹ $\sum_{k=0}^{n} (-1)^k × C_{n}^{k} = 0$
一句话证明：将 $x = -1$ 带入二项式定理

⑺ $C_{m+n}^{r} = \sum_{k = 0}^{r} C_{m}^{k} ×C_{n}^{r-k}$
显然

⑻ $C_{m + n}^{m} = \sum_{k = 0}^{m} C_{m}^{k} ×C_{n}^{k}$
一句话证明：用 $C_{n}^{r} = C_{n}^{n-r}$ 将每一项的左边（含m的组合数）变换一下即转化为前一个式子

⑼ $k×C_{n}^{k} = n × C_{n-1}^{k-1}$
左右展开即可

⑽ $\sum_{k = 0}^{\frac{n}{2}} C_{n}^{2×k}=\sum_{k = 0}^{\frac{n}{2} -1} C_{n}^{2×k+1} = 2^{n-1}$
简要证明：
首先由前面我们可以知道 $\sum_{k = 0}^{n} C_{n}^{k} = 2^n$ 和
$\sum_{k=0}^{n} (-1)^k × C_{n}^{k} = 0$
式子2告诉我们前面两项相等，式子1告诉我们两项都是 $2^{n-1}$

⑾ $\sum_{k = 1}^{n} k×C_{n}^{k} = n×2^{n-1}$
证明：将k的每一项通过 $k×C_{n}^{k} = n × C_{n-1}^{k-1}$ 换一下
变换得 $n ×\sum_{k = 1}^{n} C_{n-1}^{k-1}$
由 $\sum_{k = 0}^{n} C_{n}^{k} = 2^n$ 得
原式 $=n×2^{n-1}$

⑿ $\sum_{k = 1}^{n} k^2 × C_{n}^{k} = n×(n+1)×2^{n-2}$
证明：
将⑼ $k×C_{n}^{k} = n × C_{n-1}^{k-1}$ 左右同乘 $（k-1）$ 再化简
得 $k(k-1)C_{n}^{k} = n(n-1)C_{n-2}^{k-2}$
然后原始就好处理了。
原式 $=\sum_{k=1}^{n} (k^2-k) C_{n}^{k} +\sum_{k=1}^{n} C_{n}^{k}$
$=\sum_{k=1}^{n} n(n-1)C_{n-2}^{k-2} + n × \sum_{k=1}^{n} C_{n-1}^{k-1}$
$=n(n-1) × 2^{n-2} + n × 2^{n-1}$
$=n(n+1)×2^{n-2}$

⒀ $\sum_{k = 0}^{n} (C_{n}^{k})^2 = C_{2n}^{n}$
证明：
对于⑻ $C_{m + n}^{m} = \sum_{k = 0}^{m} C_{m}^{k} ×C_{n}^{k}$
令 $m=n$ 即可得到原式

⒁ $(1+x)^n≡1+x^n (mod$ $n)$
一句话证明：二项式定理拆开，中间的部分都包括因子n

组合数取模

$n≤1000，m≤1000$
直接上杨辉三角
$C(n,m)=C(n-1,m-1)+C(n-1,m)$

$n≤10^5，m≤10^5$ ，p是素数（比n和m大）
预处理阶乘以及阶乘的逆元

$n≤10^9，m≤10^9，p≤10^5$ 且为素数。
$Lucas$ 定理
C(n,m)=C(n%p,m%p)*C(n/p,m/p)
$C_{n}^{m} = C_{n$ $mod$ $p}^{m$ $mod$ $p}$
只需要递归处理 $C_{\frac{n}{p}}^{\frac{m}{p}}$ \一部分即可。
$Lucas$ 定理的证明：
$C_{n}^{m}$ 就是 $(1+x)^n$ 中 $x^m$ 项的系数
设 $n=ap+b$ , $m=cp+d$
有 $(1+x)^{ap+b}=((1+x)^p)^a(1+x)^b=(1+x^p)^a(1+x)^b$
所以 $(1+x)^{ap+b}中x^{cp+d}$ 项的系数就是 $C_{a}^{c} × C_{n}^{d}$

$n≤10^9，m≤10^9，p=p1^{c1}*p2^{c2}…pk^{ck}，pi^{ci}≤10^5$
对于不同的 $pi^{ci}$ ，我们可以用中国剩余定理定理解决
现在主要的问题就是对 $pi^{ci}$ 取模。
将 $C_{n}^{m}$ 表示为 $k×pi^x$ ，就可以解决问题了。

Catalan数

非线性递推关系
$h(n+1)=\sum_{k=2}^{n} h(k)*h(n-k+2)$
$(n-3)h(n)=(n/2)*[h(3)*h(n-1)+h(4)*h(n-2)+......h(n-2)*h(4)+h(n-1)*h(3)]$ .
线性递推关系
$h(n+1)=[(4n-6)/n]*h(n).$
组合式
$h(n+1)=C(2n-2,n-1)/n.$

Fibonacci数

性质：
$F(1)+F(2)+⋯+F(n)=F(n+2)-1$
$F2(1)+F2(2)+⋯+F2(n)=F(n)*F(n+1)$
$F(1)+2*F(2)+⋯+n*F(n)=n*F(n+2)-F(n+3)+2$
$F(1)+F(3)+⋯+F(2n-1)=F(2n)$
$F(2)+F(4)+⋯+F(2n)=F(2n+1)-1$
均可通过数学归纳法证明

一一对应

将一个复杂问题转化为另一个对应的简单问题，使得两者的所有情况一一对应；或者将一个不熟悉、冗余的问题转化为一个熟悉、清晰的问题。这是我们解题的重要思路和思考方法

prufer编码

假定已知的n个顶点标志为1,2,…,n，假定T是其中的一棵树，叶节点中有标号最小者。设为a1，a1的邻接点为b1，从T中消去点a1和边（a1,b1）。再从余下的树T1中寻找标号最小的叶节点，设为a2，a2的邻接点为b2，从从T1中消去点a2和边（a2,b2）。如此步骤n-2次，直到最后剩下一条边为止。于是一棵树T对应一序列{b1,b2,…,bn-2}，这些数是1~n中的数，并且允许重复。

反过来从b1,b2…bn-2可以恢复树T本身

在序列(1)中找出第一个不出现在序列(2)中的数，这个数显然便是a1，同时形成边(a1,b1)，并从(1)中消去a1,从(2)中消去b1。在余下的(1)和(2)中继续以上的步骤n-2次，直到序列(2)为空集为止。这时序列(1)中剩下的两个数x,y，边(x,y)就是树T的最后一条边

Prufer编码指的就是序列 ${b1,b2,...,bn-2}$