数论
文章平均质量分 87
cervoliu
这个作者很懒,什么都没留下…
展开
-
关于一个欧拉函数的性质的证明
性质对于任意的 n∈N∗n∈N* ,有:∑d|nφ(d)=n\sum_{d|n}\varphi(d)=n证明方法一设集合 M={1,2,3,⋯,n−1,n}M=\{1,2,3,\cdots,n-1,n\} 我们尝试将集合中的数分类。 每个数都能按照其与nn的最大公因数来分。 不妨设我们当前讨论MM中与nn的最大公因数为dd的数有多少个,d|nd|n。 假设dx∈Mdx∈M并且gcd(dx,原创 2016-02-18 22:03:09 · 2382 阅读 · 0 评论 -
【JZOJ 5163】PS的烦恼
Description话说PS总是有着各种各样的烦恼,这天,他又在为自己失败的感情史烦恼着。这时,他心中的女神,魔法少女小圆从天而降,她对他说,如果你能帮我解决一个问题,我就让你永远没有烦恼。 问题是这样的: 寻找一个最大的k,使得存在一个x使得x^k=y,那么f(y)=k,即y最多可以开k次方根。 小圆的要求是求出从a到b的f值之和(包括a和b)。 100%的数据满足:2<=a<=b<=1原创 2017-06-25 12:55:18 · 422 阅读 · 0 评论 -
杜教筛
实质就是这么一条等式 ∑i=1n∑d|if(d)=∑i=1n∑j=1⌊n/i⌋f(j)=∑i=1nf(i)⌊n/i⌋\sum_{i=1}^n\sum_{d|i}f(d)=\sum_{i=1}^n\sum_{j=1}^{\lfloor n/i \rfloor }f(j)=\sum_{i=1}^nf(i)\lfloor n/i \rfloor 通常左边不容易求,那么我们就将等式卷上一个简单的积性函数原创 2017-07-20 11:14:39 · 400 阅读 · 0 评论 -
【JZOJ 5250】 质数
Description求∑i=1n2f(i)mod 998244353\sum_{i=1}^n2^{f(i)} mod\ 998244353 其中f(i)f(i)表示ii的不同质因子个数 n≤1012n\leq 10^{12}Analysis设g(i)=2f(i)g(i)=2^{f(i)},容易发现g(i)g(i)是积性函数,可以用线性筛法筛出来,不过没什么用 答案的表示有很多种,仔细想想原创 2017-08-11 22:24:55 · 385 阅读 · 1 评论 -
拉格朗日插值法
前言在数值分析中,拉格朗日插值法是一种多项式插值方法。拉格朗日插值法下文,本蒟蒻可能有说的不妥甚至不对的地方, 欢迎大神来打脸。 定义多项式函数f(x)f(x)f(x),则点值运算就是给定xxx,求y=f(x)y=f(x)y=f(x) 插值是点值的逆运算,就是给定nnn次多项式的n+1n+1n+1个点值表达(xi,yi)(xi,yi)(x_i,y_i),要求出f(x)f(x)f...原创 2016-06-20 21:31:43 · 3280 阅读 · 0 评论 -
Miller Rabin素数测试 复习小记
Preface为什么是复习小记?因为忘了个精光QAQMiller Rabin素数测试基于以下两个原理:费马小定理若(a,p)=1(a,p)=1且pp是质数,则ap−1≡1(modp)a^{p-1}≡1\ (mod\ p),反之不一定成立二次探测定理若pp是素数,则x2≡1(modp)x^2≡1\ (mod\ p)的解只有x=1x=1或x=p−1(−1)x=p-1(-1) 也即如果x2≡1原创 2018-01-03 17:26:02 · 637 阅读 · 1 评论 -
自然数幂和 斯特林数
2018 UPD: 其实第二类斯特林数做自然数幂和更简单,这里简单写一下: 由一个基本式子出发 nk=∑i=0k{k i}[n]in^k=\sum_{i=0}^k\big\{^k_{\ i}\big\}[n]_i 考虑对nn求和 Ans=∑i=0nik=∑i=0n∑j=0k{kj}[i]jAns=\sum_{i=0}^ni^k=\sum_{i=0}^n\sum_{j=0}^k\big\{原创 2017-07-12 22:50:32 · 1534 阅读 · 0 评论 -
原根相关
定义设pp是正整数,gg是整数,若gg模pp的阶等于φ(p)φ(p),则称gg为模pp的一个原根。 简单来说,如果gg是pp的原根,那么gg的1…φ(p)\varphi(p)次幂mod pmod\ p的结果一定互不相同。应用NTT时如果模数pp可以写成p=k2c+1p=k2^c+1,可以求出pp的原根gg,并使gp−1ng^{\frac{p-1}{n}}等价于ωn\omega_n原创 2018-01-16 21:36:47 · 382 阅读 · 1 评论 -
Pollard's Rho 快速质因数分解 复习小记
Description为什么又是复习小记?因为又忘了个精光QAQPollard’s Rho分治思想我们实现过程find(n)find(n)表示对nn进行质因数分解。 如果能找到任意一个d|n,d≠1,d≠nd|n,d\neq 1,d\neq n,那么就可以转化成两个子问题find(d)find(d)与find(n/d)find(n/d)。当然如果nn本身就是质数那么肯定是原创 2018-01-05 17:19:57 · 1445 阅读 · 0 评论 -
【Burnside引理】【Pólya计数法】【Pollard's Rho】【JZOJ 5457】 项链
Description现在有m 种颜色的珠子。定义一个长度为n的项链为一个顺次连接n个珠子的环, 将所有旋转和翻转看作是等价的。 比如说, [1,2,3,4]通过旋转等价于[2,3,4,1],[3,4,1,2], [4,1,2,3]; [1, 2,3,4] 通过翻转等价于[1,4,3,2], [3,2,1,4], [2,1,4,3],[4,3,2,1]。 同时, 你还可以进行一种颜色转换操原创 2017-11-07 15:38:10 · 471 阅读 · 0 评论 -
数论/组合数学 杂项
数论相关欧拉定理及扩展欧拉定理Lucas定理当ppp是质数的时候,有 (mn)≡(⌊m/p⌋⌊n/p⌋)(mmodpnmodp)(modp)(mn)≡(⌊m/p⌋⌊n/p⌋)(mmodpnmodp)(modp){m\choose n}≡{{\lfloor m/p\rfloor}\choose{\lfloor n/p\rfloor}}{{m\bmod p}\choose {...原创 2017-08-11 16:53:24 · 507 阅读 · 0 评论 -
浅谈一类积性函数的前缀和
杜教筛的好文 链接 必须顶 有时间要去搞(现在姿势水平太弱QAQ)转载 2016-12-12 19:10:05 · 572 阅读 · 0 评论 -
【SDOI2014】数表
Description求∑i=1n∑j=1mf(gcd(i,j))[f(gcd(i,j))≤a]\sum_{i=1}^n\sum_{j=1}^mf(gcd(i,j))[f(gcd(i,j))\leq a] 其中f(n)=∑d|ndf(n)=\sum_{d|n}d QQ组数据,每组数据给出n,m,an,m,a 1≤n,m≤105,1≤Q≤2∗1041\leq n,m\leq 10^5,1\leq原创 2017-03-10 18:39:29 · 502 阅读 · 0 评论 -
【JZOJ 3851】reward
DescriptionBsny最近公司运作不佳,本年度利润才m元,但员工的奖金还是要发的,公司有n个员工,怎么发奖金这个完全由老板Bsny自己决定。Bsny想要么把这m元全发了,激励一下员工,但具体怎么分配方案有很多。比如m=1, n=2, 那么可以员工1发1元,员工2发0元;也可以员工1发0元,员工2发1元,有两种方案。 但其实,Bsny还是有点吝啬的,他想这m元不一定全部作为奖金,可以部分留给原创 2016-06-24 20:26:53 · 1844 阅读 · 0 评论 -
【组合数取模】【SDOI2013】【BZOJ 3129】【JZOJ 3214】方程
Description给定方程 X1+X2+…+Xn=mX_1+X_2+…+X_n=m 我们对第1…n11\ldots n_1个变量 进行一些限制 : X1≤A1X_1≤A_1 X2≤A2X_2≤A_2 …… Xn1≤An1X_{n1}≤A_{n1} 我们对第n1+1,n1+2…n1+n2n_1+1,n_1+2\ldots n_1+n_2个变量进行一些限制 : Xn1+1≥An1+1原创 2016-06-23 19:29:01 · 1163 阅读 · 0 评论 -
【NOI 2010】【BZOJ 2005】【JZOJ 2225】能量采集
Description 对于100%的数据:1≤n,m≤100,0001 ≤ n, m ≤ 100,000Analysis首先,让我们想一想,线段上会碰到多少个植物其实是什么? 假设坐标为(dn,dm),d=gcd(n,m),x=nd,y=md(dn,dm),d=gcd(n,m),x=\dfrac{n}{d},y=\dfrac{m}{d}。 那么对于所有的xi,yi(xi≤n,yi≤m)x原创 2016-06-17 17:49:13 · 504 阅读 · 0 评论 -
【BZOJ 2818】 gcd(附φ的线性筛法预处理)
Description给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的 数对(x,y)有多少对. 1<=N<=10^7Analysis设gcd(px,py)=p,pgcd(px,py)=p,p为素数。 满足上式的x,yx,y必然互质。 因此 Ans=∑1≤p≤n,p为素数(1+2∗∑x=2⌊np⌋φ(x))Ans=\sum_{1\leq p\leq n,p为素数}(1+2*\su原创 2016-06-15 12:32:26 · 813 阅读 · 0 评论 -
【JZOJ 4496】【GDSOI 2016】第一题 互补约数 (两种解法)
Description求∑i=1n∑d|igcd(d,id)\sum_{i=1}^n\sum_{d|i}gcd(d,\dfrac id) 对于100%的数据,n≤1011n\leq 10^{11}。Analysis方法一这是一个不需要莫比乌斯反演的方法。 ∑i=1n∑d|igcd(d,id)\sum_{i=1}^n\sum_{d|i}gcd(d,\dfrac id) =∑i∗j≤ngcd(i,原创 2016-06-13 19:33:05 · 860 阅读 · 0 评论 -
中国剩余定理
中国剩余定理X≡a1(modX≡a1(mod m1)m1) ............ X≡ak(modX≡ak(mod mk)mk) 的解。 令M=m1m2...mkM=m1m2...mk,若m1,m2...mkm1,m2...mk两两互质,则同余方程组在[0,M)[0,M)有唯一解。该解可以以如下方法构造。 令Mi=M/mi(i∈原创 2016-03-27 12:47:21 · 771 阅读 · 0 评论 -
【JZOJ 4714】 公约数
Description给定一个正整数,在[1,n]的范围内,求出有多少个无序数对(a,b)满足gcd(a,b)=a xor b。 对于100%的数据满足n<=10^7Analysis正解异常简单。 对于任意数对(a,b)(a,b),不妨设a>ba>b。 首先a xor b=c⇒a xor c=ba\ xor\ b=c\Rightarrow a\ xor\ c=b 设a−b=ca-b=c g原创 2016-08-19 17:16:31 · 514 阅读 · 0 评论 -
【JZOJ 4813】【NOIP2016提高组 五校联考2】running
Description小胡同学是个热爱运动的好孩子。 每天晚上,小胡都会去操场上跑步,学校的操场可以看成一个由n 个格子排成的一个环形,格子按照顺时针顺序从0 到n -1 标号。 小胡观察到有m 个同学在跑步,最开始每个同学都在起点(即0 号格子),每个同学都有个步长ai,每跑一步,每个同学都会往顺时针方向前进ai 个格子。由于跑道是环形的,如果 一个同学站在n - 1原创 2016-10-06 22:54:06 · 520 阅读 · 0 评论 -
【2011集训队出题】【BZOJ 2154】【JZOJ 1938】Crash的数字表格
Description莫比乌斯反演 数论求∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^n\sum_{j=1}^mlcm(i,j) n,m≤107n,m\leq 10^7Analysis设n≤mn\leq m 变形 ∑i=1n∑j=1mijgcd(i,j)\s原创 2016-07-05 12:25:23 · 563 阅读 · 1 评论 -
【JZOJ 4919】 神炎皇
Description神炎皇乌利亚很喜欢数对,他想找到神奇的数对。 对于一个整数对(a,b),若满足a+b<=n且a+b是ab的因子,则成为神奇的数对。请问这样的数对共有多少呢? 对于100%的数据n<=10^14。Analysis设a=a′x,b=b′x,n=gcd(a,b),ab=a′b′x2a=a'x,b=b'x,n=gcd(a,b),ab=a'b'x^2 那么gcd(a′,b′)=1,原创 2016-12-12 19:05:08 · 486 阅读 · 0 评论 -
【JZOJ 5746】一道比较强的 自然数幂和 板题
Description给定m,km,km,k,共TTT次询问,每次输入一个nnn,求∑ni=1ik∑i=1nik\sum_{i=1}^ni^k在 modmmodm\bmod m意义下的值 mmm的最大质因子≤3∗105≤3∗105\leq 3*10^5 2≤n,m,k≤1018,1≤T≤3∗1032≤n,m,k≤1018,1≤T≤3∗1032\leq n,m,k\leq 10^{18},1...原创 2018-05-28 22:16:33 · 288 阅读 · 0 评论