数论
文章平均质量分 90
alpha202
这个作者很懒,什么都没留下…
展开
-
莫比乌斯反演
莫比乌斯函数是定义在N+\mathbb N^+N+上的函数μμμ考虑将nnn进行分解质因数:n=∏i=1maipin=∏_{i=1}^ma_i^{p_i}n=i=1∏maipi得到该函数的通解:μ(n)={1n=10∃i∈[1,m] pi≥2(−1)m∀i∈[1,m] pi<2μ(n)=\begin{cases}1&n=1\\0&am...原创 2019-02-13 10:37:09 · 96 阅读 · 0 评论 -
多项式算法6:多项式快速幂
一开始本来想封装多项式乘法,按照普通快速幂的思路去做的Ak=(A⌊k2⌋)2A[(2,n)=1]A^k=\left(A^{\left\lfloor\frac{k}{2}\right\rfloor}\right)^2A^{\left[(2,n)=1\right]}Ak=(A⌊2k⌋)2A[(2,n)=1]但logk\log klogk还是太大,只好考虑将kkk转化为多项式的系数了。首先B...原创 2019-06-02 11:50:52 · 1490 阅读 · 0 评论 -
多项式算法5:多项式指数函数
先从牛顿迭代讲起。已知多项式函数G(z)G(z)G(z),求多项式函数F(x)F(x)F(x)满足G(F(x))≡0(modxn)G(F(x))\equiv0 \pmod{x^n}G(F(x))≡0(modxn)考虑用迭代求解,假设我们已经求得F0(x)F_0(x)F0(x)满足G(F0(x))≡0(modx⌈n2⌉)G(F_0(x))\equiv0\pmod{x^{\left\lcei...原创 2019-06-02 11:47:39 · 2583 阅读 · 0 评论 -
多项式算法4:多项式对数函数
已知多项式F(x)F(x)F(x),求解G(x)G(x)G(x)使得G(x)=lnF(x)G(x)=\ln F(x)G(x)=lnF(x)模板题,我们知道lnF(x)=∫F′(x)F−1(x) dx\ln F(x)=∫F'(x)F^{-1}(x)\ \text{d}xlnF(x)=∫F′(x)F−1(x) dx由指数函数的求导公式(xα)′=αx...原创 2019-06-02 11:42:44 · 627 阅读 · 0 评论 -
多项式算法3:多项式除法
已知次数界分别为n,mn,mn,m的多项式F(x)F(x)F(x)与G(x)G(x)G(x),求解次数界为n−mn-mn−m的两个多项式Q(x)Q(x)Q(x)与R(x)R(x)R(x)使得F(x)=G(x)Q(x)+R(x)F(x)=G(x)Q(x)+R(x)F(x)=G(x)Q(x)+R(x)我们令Ar(x)=xDeg(A)A(x−1)A_r(x)=x^{\text{Deg}(A)}A(...原创 2019-06-02 11:38:25 · 4693 阅读 · 0 评论 -
多项式算法2:多项式求逆
已知多项式F(x)F(x)F(x),求满足条件的G(x)G(x)G(x)使得F(x)G(x)≡1(modxn)F(x)G(x)\equiv1\pmod{x^n}F(x)G(x)≡1(modxn)考虑迭代求解令F(x)G0(x)≡1 (modx⌈n2⌉)F(x)G_0(x)\equiv1\ \pmod{x^{\left\lceil\frac{n}{2}\right\rceil}}F...原创 2019-06-02 11:33:15 · 1501 阅读 · 0 评论 -
多项式算法1:多项式乘法
多项式简介对于数域F\mathbb FF,若有∀i∈{1,2,3,⋯ ,n}\forall i\in\{1,2,3,\cdots,n \}∀i∈{1,2,3,⋯,n},则f(x)=a0+a1x+a2x2+⋯+anx=∑i=1naixif(x)=a_0+a_1x+a_2x^2+\cdots+a_nx=\sum_{i=1}^na_ix^if(x)=a0+a1x+a2...原创 2019-06-02 11:27:26 · 5131 阅读 · 1 评论 -
题解 P1829 [国家集训队]Crash的数字表格 / JZPTAB
题目背景提示:原 P1829 半数集问题 已经迁移至 P1028 数的计算题目描述今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple)。对于两个正整数a和b,LCM(a, b)表示能同时整除a和b的最小正整数。例如,LCM(6, 8) = 24。回到家后,Crash还在想着课上学的东西,为了研究最小公倍数,他画了一张NM的表格。每个格子里写了一个数...原创 2019-04-30 21:55:50 · 285 阅读 · 0 评论 -
题解 P4466 [国家集训队]和与积
题目描述给出NNN,统计满足下面条件的数对(a,b)(a,b)(a,b)的个数:1≤a<b≤N1\le a<b\le N1≤a<b≤Na+ba+ba+b整除a×ba×ba×b输入输出格式输入格式一行一个数NNN输出格式一行一个数表示答案思路讲解比较好的莫比乌斯反演题了。ans=∑a=1N∑b=a+1N[(a+b)∣ab]\text{ans...原创 2019-04-30 21:51:27 · 425 阅读 · 0 评论 -
题解 P3911 最小公倍数之和
题目描述对于A1,A2,...,ANA_1,A_2,...,A_NA1,A2,...,AN,求:∑i=1N∑j=1Nlcm(Ai,Aj)\sum_{i=1}^N\sum_{j=1}^Nlcm(A_i,A_j)i=1∑Nj=1∑Nlcm(Ai,Aj)的值lcm(a,b)lcm(a,b)lcm(a,b) 表示a 和b 的最小公倍数输入输出格式输入格式第1行,1个整数N。第...原创 2019-04-30 21:47:58 · 314 阅读 · 0 评论 -
题解 P2260 [清华集训2012]模积和
题目背景数学题,无背景。题目描述求∑i=1n∑j=1m(nmod  i)(mmod  j),i≠j\sum_{i=1}^n\sum_{j=1}^m(n\mod i)(m\mod j),i=\not ji=1∑nj=1∑m(nmodi)(mmodj),i≠jmod 19...原创 2019-04-30 21:42:05 · 219 阅读 · 0 评论 -
题解 TG2013 火柴排队
题目描述:涵涵有两盒火柴,每盒装有n根火柴,每根火柴都有一个高度。 现在将每盒中的火柴各自排成一列, 同一列火柴的高度互不相同, 两列火柴之间的距离定义为: ∑(ai−bi)2\sum(a_i-b_i)^2∑(ai−bi)2其中aia_iai表示第一列火柴中第iii个火柴的高度,bib_ibi表示第二列火柴中第iii个火柴的高度。每列火柴中相邻两根火柴的位置都可以交换,请你通过交换使...原创 2019-04-14 11:51:08 · 148 阅读 · 0 评论 -
题解 JZOJ 1162.【普及组模拟赛】最大杂置(set)
题目描述:令SSS为nnn个元素的集合,则S有2n−12^{n-1}2n−1个子集(除去空集)。现在要你从这2n−12^{n-1}2n−1个子集中选出最多的子集,使这些子集能构成一个杂置。杂置是指任意两个集合没有包含或被包含的关系。例如对于有3个元素的集合{a,b,c}\{a,b,c\}{a,b,c}。{a,b},{a,c},{b,c}\{a,b\},\{a,c\},\{b,c\}{a,b}...原创 2019-04-06 18:55:05 · 209 阅读 · 1 评论 -
伯努利数
参考文献资料:小蒟蒻yyb的博客百度百科-伯努利数伯努利数(Bernoulli)——学习笔记康复计划#3 简单常用的几种计算自然数幂和的方法Bernoulli number has been found for about 200 years.Its name is named after Bernoulli, its discoverer.Recursive definition:...原创 2019-02-13 19:06:55 · 714 阅读 · 0 评论 -
多项式算法7:多项式三角函数
首先证明一下欧拉公式eiθ=cosθ+isinθe^{i\theta}=\cos\theta+i\sin\thetaeiθ=cosθ+isinθ由麦克劳林展开式,我们有ex=∑n=0∞xnn!e^{x}=\sum_{n=0}^{\infty}\frac{x^n}{n!}ex=n=0∑∞n!xn指数乘上虚数单位iiieix=∑n=0∞(ix)nn!=∑n=0∞(i)nxnn!e^{...原创 2019-06-02 11:55:45 · 2312 阅读 · 0 评论