首页 > 其他分享 >数论

数论

时间:2024-05-02 15:44:24浏览次数:21  
标签:phi 数论 质数 合数 times 逆元 ll

数论

常见筛法

对于任意正整数 \(A\),存在唯一集合 {\((p_1,q_1),(p_2,q_2),\dots,(p_n,q_n)\)} 满足 \(A=\prod^n_{i=1} {p_i}^{q_i}\)

\(\min(a,b)+\max(a,b)=a+b\)

\(\gcd(a,b)\times\operatorname{lcm}(a,b)=ab\)

埃氏筛

在 \(O(n \log\log n)\) 时间内预处理出 \([1,n]\) 内所有的质数。

\([2,n]\) 进行扫描,未标记则加入质数中,并把倍数标记为合数。

  • P7960

埃氏筛思想,将十进制含有 \(7\) 的数看为“类质数”。

  • P1835

仍然考虑埃氏筛,扫描每个 \(\le \sqrt{R}\),并对其在 \([L,R]\) 内的倍数标记为合数。

称为区间筛

线性筛

\(O(n)\) 计算。

令 \(low(n)\) 表示 \(n\) d的最小质因子,对 \([2,n]\) 扫描,对于 \(i\),我们枚举所有 \(\le low(i)\) 的质数 \(j\),并将 \(ij\) 标记为合数。

发现 \(n\) 只会在

可以用来求质因子个数。

\(\phi(n)\) 表示 \([1,n]\) 中与 \(n\) 互质的数。

\(n=\sum_{d|n} \phi(\frac{n}{d})\)

如果 \(\gcd(p,q)=1\),则 \(\phi(pq)=\phi(p)\phi(q)\),反证法。

模意义

光速乘

ll times(ll a,ll b,ll c){
ull t=(long double)a*b/c+0.5;
ll ans=(ull)a*b-t*c;
if(ans<0) ans+=c;
return ans;
}

逆元

模意义下的除法

\(a\) 的逆元是 \(a\) 在模意义下的倒数,逆元关系相互的。

威尔逊定理
对于质数 \(p\),有 \(1\times 2\times \dots \times(p-1)\equiv -1(\mod p)\)

费马小定理
\(p\) 为质数:

\(a^{p-1} \equiv 1(\mod p)\)

\(a^{-1} \equiv a^{p-2}\),用快速幂求

\(p\) 为合数:

欧拉定理

若 \(p \ge 2\) 且 \(\gcd(a,p) =1\) ,则 \(a^(\phi(p)) \equiv 1(\mod p)\)

线性预处理逆元

  1. \(1\) 的逆元为 \(1\)。
  2. 让 \(p\) 对 \(i\) 作带余除法,得到 \(p=qi+r(r<i)\)
  3. 将等式带入

inv[i]=fac[i-1]*ifac[i]

裴属定理
对于所有的 \(ai+bj\),最小的正整数为 \(\gcd(a,b)\)

标签:phi,数论,质数,合数,times,逆元,ll
From: https://www.cnblogs.com/CheZiHe929/p/18170253

相关文章

  • 网课-数论学习笔记
    埃氏筛P7960[NOIP2021]报数P1835素数密度:区间筛。预处理\(\sqrt{R}\)内的质数,然后用埃氏筛筛[L,R]的质数。线性筛-EOF-欧拉函数P10031「CfzRound3」XorwithGcd光速乘用于解决$$llTimes(lla,llb,llc){ ullt=(longdouble)a*b/c+0.5; llans=(u......
  • 数论学习笔记 (4):扩展欧几里得算法
    概述扩展欧几里得算法(\(exgcd\))可以用来求形如\(ax+by=c\)的不定方程的通解。铺垫-\(\small\texttt{ax+by=gcd(a,b)}\)的解\(exgcd\)的思想是在用辗转相除法递归\(gcd(a,b)\)的回溯时求出对应方程\(ax+by=gcd(a,b)\)的解。考虑方程\(ax+by=gcd(a,b)\)。看回辗......
  • 数论模拟(1) 小朋友们,我们今天来找规律
    \(60\)分钟,干出来\(30\)至\(40\)分(满分\(50\)),最后一步没写出来还是有点rz.题目:求最小的整数\(n\),使得对至少两个不同的奇素数\(p\),有\[\sum_{k=1}^{n}(-1)^{v_p(k!)}<0.\]解:根据\(v_p\)函数的性质,可以对所有正整数进行规律性地分块,每块中的\(v_p\)值都是相同的:......
  • CF EDU164-E-数论分块
    link:https://codeforces.com/contest/1954/problem/E有一排怪物,第\(i\)只有\(a_i\)的血,每次攻击可以选择在\(i\)处放一个技能,技能会一直向左/右以相同的\(k\)点伤害扩散,直至遇到边界或已经死亡的怪物。问最少需要放几次技能?需要对所有\(k\)回答答案。\(n,a_i\leq10^......
  • 数论学习笔记 (3):因数与倍数
    \(\texttt{godmoo}\text{}\texttt{の}\text{}\texttt{数论学习笔记之}\text{}\boxed{因数与倍数}\)定义因数/约数,倍数:若\(d\midn\),则\(d\)是\(n\)的因数,\(n\)是\(d\)的倍数。公因数/公约数,公倍数:公共的因数/约数、倍数。最大公因(约)数:\(GreatestCommonDi......
  • 数论习题(2) Legendre公式+高斯函数
    本人独自证明,可能存在一定疏漏.题目:\[m!n!(m+n)!\mid(2m)!(2n)!.\]证明:对于每个素数\(p\),考察式子两边的\(p\)进赋值,即证\[v_p((2m)!(2n)!)\geqv_p(m!n!(m+n)!).\]根据\(p\)进赋值的基本性质与Legendre公式,有\[\begin{align*}v_p((2m)!(2n)!)&=v_p((2m)!)......
  • [学习笔记] 丢番图方程 & 同余 & 逆元 - 数论
    首先,他们几个都是兄弟,有着极大的相似性。另外,他们的各自的思想都能够很好的服务于另外几个,有助于加深理解。线性丢番图方程丢番图不是个图啊!他是个man……——百度百科现在主要说的是二元线性丢番图方程:通用形式为\(ax+by=c\)。其中常数全都为整数。很像不定方程对吧?现在在......
  • 一个人的数论 题解
    Solution令指数为\(k\)正常反演得到\[\sum_{d\midn}\mu(d)d^k\sum_{i=1}^{\fracnd}i^k\]设\(f(x)=\sum_{i=1}^xi^k\),它是一个关于\(x\)的\(k+1\)次多项式求这个多项式可以插值\(\mathcalO(n^2)\)(推荐)高斯消元(待定系数法)\(\mathcalO(n^3)\)直接伯努利数\(\ma......
  • 【数论】最大公因数和最小公倍数(GCD和LCM)
    最大公因数(GCD)两个数的最大公因数很好做,使用内置的库函数即可,注意x和y的类型要相同。llgcd=__gcd(x,y);如果要求多个数的最大公因数,那么初始化为0(因为根据定义,0和任何数x的gcd都是x,所以0是gcd操作的幺元),然后分别进行gcd即可。llgcd=0;for(inti=1;i<=n;++i)......
  • 数论
    求逆元费马小定理扩展欧几里得线性求逆元中国剩余定理形如\[x\equiva_1\pmod{n_1}\]\[x\equiva_2\pmod{n_2}\]\[x\equiva_3\pmod{n_3}\]\[x\equiva_4\pmod{n_4}\]OI-WIKI点击查看代码LLCRT(intk,LLa[],LLr[]){LLn=1,ans=0;......