首页 > 其他分享 >数学数论全集

数学数论全集

时间:2022-10-03 18:33:16浏览次数:43  
标签:数论 omega varphi mu perp 数学 prod 全集 equiv

扩展欧几里得定理

\[ax+by=(a,b); bx_0+(a \bmod b)y_0=(a,b); x = y_0; y = (a/b) y_0 + b(x_0) \]

证:\({a}x+{b}y=(a,b)=(b,a \bmod b)=bx+(a \bmod b)y=bx_0+(a-(a/b)b)y_0={{b}}x_0+{{a}}y_0-(a/b){{b}}y_0=\mathbf{{a}}(y_0)+\mathbf{{b}}(x_0+(a/b)y_0)\),化简后式为 \(ax+by\) 形式,对应得证。

欧拉定理

\[x^{\varphi(p)} \equiv 1 \pmod p(x \perp p) \]

引理一:对于两个与 \(p\) 互质的数 \(a,b\),\(p \perp ab\)。
证:显然有 \(a,b\) 均无 \(p\) 共同因子,设 \(\mathscr S x\) 为 \(x\) 的所有不同质因子集,则有 \(p \notin \mathscr S a, p \notin \mathscr S b\),而且显然有 \(\mathscr S (ab) = \mathscr S a \cup \mathscr S b\)。所以有 \(p \notin \mathscr S (ab)\),故 \(p \perp ab\)。
引理二:任取 \(k \in [0, p-1] \wedge k \perp p\),令 \(S\) 为在 \([0, p-1]\) 之内所有与 \(p\) 互质的数,令 \(T\) 为 \(S\) 内每个数乘 \(k\) 的集合 \(\bmod p\)(不可重),则有 \(T \subseteq S\)。
令 \(a \in S\) 则 \(ak \in T\)。此时 \(a,k\) 均 \(\perp p\),故 \(ak \in S\),所以 \(ak \in T \Longrightarrow ak \in S\),故得证。
引理三:任取 \(k \in [0, p-1] \wedge k \perp p\),\(S,T\) 在上个引理定义,\(S=T\)。
显然只需要证明 \(|T| = |S|\)。反证法,假设不等大,则必有两个不同元素 \(a,b\),满足 \(ak \equiv bk\),又因为 \(k \perp p\),故 \(k^{-1}\) 存在,于是两边乘 \(k^{-1}\),即可推出矛盾。
原命题证明:引理三套用 \(k = x\),得到 \(S_x = T_x\),对集合内所有元素求积得到:

\[\prod S \equiv \prod S \times x^{|S|} \pmod p \]

,又因为 \(\prod S \perp p\)(反复套用引理一),于是可以对 \(\prod S\) 取逆元,而且根据 \(\varphi\) 的定义,有 \(|S| = \varphi(p)\),得证。
(常用变式:费马小定理,也就是满足 \(p \in \mathbf P\) 的情况下,有 \(x \equiv x^p \pmod p\),证明显然。)
(常用变式:扩展欧拉定理,将 \(x^y\) 分成 \(x^tx^{y-t}\) 的形式,其中 \(t \perp p\),套用欧拉定理。)
扩展欧拉定理的形式:

\[a^b \equiv \begin{cases} a^{b \bmod \varphi(m)}, &a \perp m, \\ a^b, &a \not\perp m \wedge b < \varphi(m) \\ a^{(b \bmod \varphi(m)) + \varphi(m)}, &a \not\perp m \wedge b \ge \varphi(m) \end{cases}\]

,证明烦琐,不做过多证明。

CRT

方程组 \(x \equiv a_i \pmod{b_i}\) 一定有一解

\[x = \sum a_i \epsilon_i \]

,其中

\[\epsilon_i = M_i \times M_{i_{b_i}}^{-1} \]

,其中

\[M_i = \dfrac{\prod b}{b_i} \]


证明:这个方程组的 \(a\) 可以抽象为一个向量 \(\overrightarrow{a}\),且有 \(\overrightarrow{a}\) 的解加 \(\overrightarrow{b}\) 的解就是 \(\overrightarrow{a} + \overrightarrow{b}\) 的解,于是我们只需要找出基向量 \([1,0,0,\dots,0,0],[0,1,0,\dots,0,0],\dots,[0,0,0,\dots,0,1]\) 即可。于是对于第 \(i\) 个基向量 \(\varepsilon_i\)(和 \(\epsilon_i\) 区分开,实际上,我们马上就会看到 \(\varepsilon_i\) 的解就是 \(\epsilon_i\)),有 \(M_i \mid \varepsilon_i\) 的解,以及 \(\varepsilon_i\) 的解 \(\equiv 1 \pmod{b_i}\),于是有 \(\varepsilon_i\) 的解为 \(k \times M_i\),而由于其 \(\equiv 1 \pmod{b_i}\),可以发现 \(k = M_{i_{b_i}}^{-1}\),于是得到 \(\varepsilon_i\) 的解就是 \(\epsilon_i\)。得证。

ExCRT

发现要求 \(M_i \perp b_i\)。如果不满足怎么办呢?
则两两合并:

\[x = a_1 + pb_1 = a_2 + qb_2 \]

,可以通过 exgcd 解出 \(p,q\),套入可以将两个方程合并。

莫反

定理:

\[\mu * 1 = \epsilon \]


定义 \([a_1\quad a_2\quad a_3\quad \cdots] = \prod p_i^{a_i}\)。若 \([a_1\quad a_2\quad a_3\quad \cdots] = r\),则记 \(r_i = a_i\)。定义 \(\omega(n) = \sum [n_i \ge 1]\)。定义 \(t(n) \iff (\max n_i) \ge 2\)。定义 \(S(n) = \{p_i \mid S_i \ge 1\}\)。显然有 \(|S(n)| = \omega(n)\)。
则显然有 \(\omega(n) = \omega(m) \Longrightarrow (\mu * 1)(n) = (\mu * 1)(m)\)。
\(\mu\) 可以解释为 \(\mu(n) = \left\{\begin{matrix}0 & t(n) \\ -1^{\omega(n)} & \text{otherwise}\end{matrix}\right.\),于是对于 \(\neg (t(n) \vee t(m))\),有 \(\omega(n) = \omega(m) \Longrightarrow \mu(n) = \mu(m)\),于是建立 \(S(n) \leftrightarrow S(m)\) 的双射 \(f(x)\),因为对于相同的 \(K_i \in [0,1]\),\(\omega(\prod\limits_{i \in S(n)} f(i)^{K_i})=\omega(\prod\limits_{i \in S(n)} i^{K_i})\),所以可以建立一个双射 \(g(x) : \{d \mid n \wedge \neg t(d)\} \leftrightarrow \{d \mid m \wedge \neg t(d)\}\),满足 \(\mu(g(x)) = \mu(f(x))\),而 \(\forall x \in \{d \mid m \wedge t(d)\} \cup \{d \mid n \wedge t(d)\}, \mu(x) = 0\),所以根据 \((1-1)^n\) 的展开,得证。
(选自我自己写的博客)

小发现

把函数写成

\[f(n) = \prod_i f'(p_i,n_i) \]

的形式,则会有

\[(f*g)(n) = \prod_i \sum_{j=0}^{n_i} f'(p_i, j) \times g'(p_i, n_i - j) \]

,于是可以搞一道原创题:求 \((1*1*\cdots*1*1)(n)\),只不过我已经公开了(笑)
还有一些常见函数的表示:
\(\epsilon'(p, a) = [a = 0]\),\(1'(p, a) = 1\),\(id'(p, a) = p^a\),\(\mu'(p, a) = [a < 2](1 - 2a)\),\(d'(p, a) = (a + 1)\),\(\sigma'(p, a) = p^a+1\),\(\varphi'(p, a) = [a < 2](1 - 2a)p^a\)

关于质数

标签:数论,omega,varphi,mu,perp,数学,prod,全集,equiv
From: https://www.cnblogs.com/lhx-oier/p/16742937.html

相关文章

  • 阿里巴巴全球数学竞赛品牌传播分析(上)
    9月26日,2022阿里巴巴全球数学竞赛获奖名单公布,4座金杯分别由平均年龄25岁,来自美国麻省理工学院、美国布朗大学、北京大学在读数学博士斩获。77位获奖者中00后超五成引热议,因......
  • 数学知识1.3
    一、简述本文章主要介绍欧拉函数以及快速幂的相关算法。二、欧拉函数定义\(1∼N\)中与\(N\)互质的数的个数被称为欧拉函数,记为\(\phi(N)\)。若在算数基本定理中,\(N......
  • 2022.9.11 2022年全国高中数学联赛A卷加试第二题另解
    二.设整数\(n(n>1)\)恰有\(k\)个互不相同的素因子,记\(n\)的所有正约数之和为\(\sigma(n)\),证明\(\sigma(n)|(2n-k)!\)(\(2022\)年全国高中数学联赛加试第二题)解析思路是很......
  • 函数学习
    1.函数分为:库函数和自定义函数查找库函数的网站:http://www.cplusplus.com   或着​​https://legacy.cplusplus.com/reference/​​中文版的C语言函数网站:cpprefenrer......
  • 个人数论专题总结
    中国剩余定理(CRT)证明与应用问题定义给定一组同余方程:\[(S):\begin{cases}x≡a_1(\text{mod}m_1)\\x≡a_2(\text{mod}m_2)\\……\\x≡a_n(\text{mod}m_n)\\\end{cases......
  • 数论笔记
    倍数若\(a,b,k\in\mathbbN\),且\(a\timesk=b\),那么\(b\)是\(a\)的倍数,称\(a\)整除\(b\),记作\(a\midb\)。\([1,n]\)中\(x\)的倍数有\(\left\lfl......
  • 2022高联P2数论
    P2:设整数\(n(n>1)\)恰有k个互不相同的质因子,记n的所有正约数之和为\(\sigma(n)\).求证:\(\sigma(n)|(2n-k)!\).既然已给出了质因子个数和\(\sigma(n)\),那么就可设出\(n\)......
  • CF1514B AND 0, Sum Big-不学数学,你怎么会有进步呢?qaq
    Problem-1514B-Codeforces题意:给定n,k元素为[0,2^k-1],(一共2^k个数),选n个数组成数列,使得1:所有元素的&的和为02:所有元素的和尽量大解:假如n=2,k=2.可知元素只有2位0:0......
  • 《有个数学模型想咨询下吧友的意见》 回复
    《有个数学模型想咨询下吧友的意见》    https://tieba.baidu.com/p/8049516112      风消云散亦无悔:假设某游戏某种材料,记为材料A爆率是20%不是腾......
  • 算法数学笔记-零、常用数表及杂项
    目录零、常用数表及杂项常用数表牛顿迭代牛顿广义二项式定理一些结论范德蒙德卷复数相乘突然发现博客园可以存笔记,这样就可以避免出门没带电脑而又想看笔记的情况了,还方便......