学习原根 by OI-wiki
根据 OI-wiki 的讲解,加以自己的理解和简化。偏重于算法竞赛而不是数学竞赛。
前置知识:
费马小定理:\(a^{p-1} \equiv 1(\mod p)\),\(p\) 为质数。
欧拉定理:\(a^{\varphi(m)}\equiv 1(\mod m)\),m 为任意正整数。
拉格朗日定理:\(p\) 为质数,\(n\) 整系数次多项式在模 \(p\) 意义下至多有 \(n\) 个不同解。(即多项式 \(f(x)\),\(f(x)\equiv 0(\mod p)\) 的 \(x\) 的取值至多有 \(n\) 个。)
定义一个东西,称作阶。
阶:满足 \(a^n \equiv 1(\mod m)\) 的最小正整数 \(n\),即为 \(a\) 模 \(m\) 的阶,记作 \(\delta_m(a)\)。由欧拉定理可知,任意 \(m\) 都存在 \(\varphi(m)\) 使得满足这个方程,所以 \(n\) 一定存在且 \(n \leq \varphi(m)\)。
性质 \(1\):\(a^1,a^2,...,a^{\delta_m(a)}\) 模 \(m\) 两两不同余。
证明:若 \(a^i \equiv a^j(\mod m),i < j \leq \delta_m{a}\),则 \(a^{j-i}\equiv 1(\mod m)\)(等式两边同时除以 \(a^j\))。
显然 \(j-i<\delta_m(a)\),根据阶的定义,与阶的最小性矛盾。
性质 \(2\):若 \(a^n \equiv 1(\mod m)\),则 \(\delta_m(a) | n\)。
证明:
考虑把 \(n\) 分解为 \(\delta_m(a)\) 的倍数+\(n\) 模 \(\delta_m(a)\) 得到的余数的形式,然后反证。
还有一堆性质,但我要学的是原根,故跳过。
原根:严谨定义:设 \(m\in \N^*,a\in \Z\)。若 \(\gcd(a,m)=1\),且 \(\delta_m(a)=\varphi(m)\),则称 \(a\) 为模 \(m\) 意义下的原根。一般的用途,对于质数 \(p\),其原根为 \(g\)(可以证明一定存在),\(g^i \mod p,0 \leq i \leq p\) 的值互不相同。这一点在 NTT 里是关键所在,如果以后有机会我会写多项式相关的博客。
如何求原根?
一个结论:任意质数都有原根。
第二个结论:若 \(m\) 有原根,其最小原根不多于 \(m^{0.25}\) 级别。
可以考虑暴力枚举每一个数直到找到原根。判别式为原根的定义式。
判断数 \(x\) 是否为原根,先判断 \(x^{\varphi(m)}\equiv 1(\mod m)\) 是否成立,如果成立,然后根据阶的性质二,枚举 \(\varphi(m)\) 的因数是否成立,若因数全为否,则 \(x\) 为原根。
更简洁地,枚举 \(\varphi(m)\) 除以其中一个质因数,一共质因数个数多个数,这些数是否成立等价于所有因数是否成立。
\(\varphi\) 在 \(m\) 取值多但数值小时可以用线性筛,如果 \(m\) 只有一个,但是很大,考虑用根号算法算出 \(\varphi\) 的值。
以后贴一份求阶的线性筛代码。
学习原根 by OI-wiki的更多相关文章
- 【学习笔记】OI玄学道—代码坑点
[学习笔记]\(OI\) 玄学道-代码坑点 [目录] [逻辑运算符的短路运算] [\(cmath\)里的贝塞尔函数] 一:[逻辑运算符的短路运算] [运算规则] && 和 || 属于逻 ...
- 【学习笔记】OI模板整理
CSP2019前夕整理一下模板,顺便供之后使用 0. 非算法内容 0.1. 读入优化 描述: 使用getchar()实现的读入优化. 代码: inline int read() { int x=0; ...
- SPOJLCS Longest Common Substring
题意 A string is finite sequence of characters over a non-empty finite set Σ. In this problem, Σ is th ...
- SPOJ LCS Longest Common Substring 和 LG3804 【模板】后缀自动机
Longest Common Substring 给两个串A和B,求这两个串的最长公共子串. no more than 250000 分析 参照OI wiki. 给定两个字符串 S 和 T ,求出最长 ...
- 「懒惰的美德」我用 python 写了个自动生成给文档生成索引的脚本
我用 python 写了一个自动生成索引的脚本 简介:为了刷算法题,建了一个 GitHub仓库:PiperLiu / ACMOI_Journey,记录自己的刷题轨迹,并总结一下方法.心得.想到一个需求 ...
- OI学习之路上的宝藏网站/App分享
OI学习之路上的宝藏网站/App分享 想要变强吗少年?这里有各种我平时收集的网站/App,它们可以帮助你更好地学习算法或者找到解题思路.废话不多说,快来打开新世界的大门罢~ 知识学习 觉得各种知识晦涩 ...
- 【算法学习笔记】Meissel-Lehmer 算法 (亚线性时间找出素数个数)
「Meissel-Lehmer 算法」是一种能在亚线性时间复杂度内求出 \(1\sim n\) 内质数个数的一种算法. 在看素数相关论文时发现了这个算法,论文链接:Here. 算法的细节来自 OI w ...
- 【算法学习笔记】动态规划与数据结构的结合,在树上做DP
前置芝士:Here 本文是基于 OI wiki 上的文章加以修改完成,感谢社区的转载支持和其他方面的支持 树形 DP,即在树上进行的 DP.由于树固有的递归性质,树形 DP 一般都是递归进行的. 基础 ...
- 动态规划 Dynamic Programming 学习笔记
文章以 CC-BY-SA 方式共享,此说明高于本站内其他说明. 本文尚未完工,但内容足够丰富,故提前发布. 内容包含大量 \(\LaTeX\) 公式,渲染可能需要一些时间,请耐心等待渲染(约 5s). ...
- CCPC比赛与算法学习的个人分享
大赛简介 中国大学生程序设计竞赛(China Collegiate Programming Contest,简称CCPC)是工业和信息化部教育与考试中心主办的 "强国杯"技术技能大 ...
随机推荐
- 小笨自采集在线电脑壁纸 v2.0-支持移动端自适应,支持https
这周闲着没事,翻着网站,记得去年的发表的一篇文章小笨分享一款高清电脑壁纸API,主要是将孟坤壁纸改造支持https,还加了一个搜索功能.以前的壁纸是这样的: 但是,不支持移动端浏览,于是小笨趁着这周周 ...
- JVM——锁
对象头[每个对象都具有对象头] Mark:对象头的标记(32位),描述对象的 hash.锁信息.垃圾回收标记.年龄:内容包括:①.指向锁记录的指针:②.指向 monitor 的指针:③.GC 标记:④ ...
- Stanford CS 144, Lab 0: networking warmup 实验
Stanford CS 144, Lab 0: networking warmup Finish Stanford CS144 lab0 and pass the test. 2023/03/29 - ...
- 集合-ArrayList 源码分析
1.概述 ArrayList 是一种变长的集合类,基于定长数组实现.ArrayList 允许空值和重复元素,当往 ArrayList 中添加的元素数量大于其底层数组容量时,其会通过扩容机制重新生成一个 ...
- 细节拉满,80 张图带你一步一步推演 slab 内存池的设计与实现
1. 前文回顾 在之前的几篇内存管理系列文章中,笔者带大家从宏观角度完整地梳理了一遍 Linux 内存分配的整个链路,本文的主题依然是内存分配,这一次我们会从微观的角度来探秘一下 Linux 内核中用 ...
- Java 框架面试题-Spring Boot自定义配置与自动配置共存
Spring Boot 是一个快速开发框架,可以简化 Spring 应用程序的开发,其中自定义配置是其中一个非常重要的特性. 在 Spring Boot 中,自定义配置允许开发者以自己的方式来配置应用 ...
- ServletContext 详解(转载)
转载:https://www.cnblogs.com/zjdxr-up/p/7761813.html ServletContext,是一个全局的储存信息的空间,服务器开始,其就存在,服务器关闭,其才释 ...
- mysql迁移:mysqldump导出表结构及数据
问题描述:有需要mysql某几张表的需求,某个数据库某几张表,导出先检查相应的数据库和表是否存在 数据泵用法:默认导出的是表结构以及表中的数据 mysqldump -uroot -p -S /data ...
- 从零开始学Vue(二~三)—— Vue 实例 / 模板语法(插值、指令)
概述 vue.js作为现在笔记热门的JS框架,使用比较简单易上手,也成为很多公司首选的JS框架. 但是对于初学者可能学起来有些麻烦,所以推出<从零开始学Vue>系列博客,本系列计划推出19 ...
- C++ 测试框架 GoogleTest 初学者入门篇 丙
theme: channing-cyan *以下内容为本人的学习笔记,如需要转载,请声明原文链接微信公众号「ENG八戒」https://mp.weixin.qq.com/s/RIztusI3uKRno ...