---数论
BAJim_H
比孤独更可悲的事情,就是根本不知道自己很孤独,或者分明很孤独,却把自己都骗得相信自己不孤独。
展开
-
[JZOJ4922] 环
Description 小A有一个环,环上有n个正整数。他有特殊的能力,能将环切成k段,每段包含一个或者多个数字。对于一个切分方案,小A将以如下方式计算优美程度: 首先对于每一段,求出他们的数字和。然后对于每段的和,求出他们的最大公约数,即为优美程度。 他想通过合理地使用他的特殊能力,使得切分方案的优美程度最大。 对于20%的数据,n<=16,ai<=10 对于40%的数原创 2016-12-21 20:06:47 · 388 阅读 · 0 评论 -
[JZOJ5539] psy
Description 有很多n位数(可以有前导0),如果一个n位数X对所有的k(1≤kn)k(1≤k都满足 X∗10kMod10n>XX*10^k Mod 10^n > X,这个X我们就认为它脱团了。现在告诉你n,求出有多少个X脱团了。 题目是这样的,设f(n)是n位数里脱团数的数量(脱团数定义如上),现在让你求出f(1)*1^2+f(2)*2^2+…+f(i)*i^2+…+原创 2018-01-22 15:56:21 · 144 阅读 · 0 评论 -
[JZOJ5509] 兔子
Description Solution我们发现最后的答案是前K个的次幂相乘的形式于是可以矩阵乘法求出对应的指数指数可能会爆,并且底数与模数不一定互质 那么用扩展欧拉定理xc≡{xc,xc%φ(m)+φ(m),cmc≥m(modm)\begin{eqnarray}x^c\equiv\begin{cases}x^c, &c<m\cr x^{c\%φ(m)+φ(m)},原创 2018-01-05 21:26:02 · 263 阅读 · 0 评论 -
[JZOJ5508] 距离
Description Solution首先列出最基本的式子 ∑i=1ni((2i+1)k−(2i−1)k)\sum\limits_{i=1}^{n}i((2i+1)^k-(2i-1)^k)前后加减拆开来 =∑i=1ni(2i+1)k−∑i=1ni(2i−1)k=\sum\limits_{i=1}^{n}i(2i+1)^k-\sum\limits_{i=1}^{n}i(2i-原创 2018-01-05 22:35:59 · 215 阅读 · 0 评论 -
[JZOJ5511] 送你一个DAG
Description 给出一个 n 个点 m 条边的 DAG 和参数 k。 定义一条经过 ll 条边的路径的权值为 lkl^k. 对于 i = 1…n, 求出所有 1 到 i 的路径的权值之和, 对 998244353 取模. 对于前 20% 的数据, n ≤ 2000,m ≤ 5000; 对于另 10% 的数据, k = 1; 对于另 20% 的数据, k ≤原创 2018-01-09 12:09:02 · 918 阅读 · 0 评论 -
[JZOJ5527] Silly
Description求∑i=1nφK(gcd(i,n))∗lcm(i,n)(modm)o\sum\limits_{i=1}^{n}φ^K(\gcd(i,n))*lcm(i,n)\pmod moT组询问 N≤1018,K≤109,T≤300,mo≤109+7N\leq 10^{18},K\leq 10^9,T\leq 300,mo\leq 10^9+7且mo为奇质数Solution原创 2018-01-18 12:15:09 · 120 阅读 · 0 评论 -
【学习小记】狄利克雷卷积+杜教筛
Preface这东西分明就是玄学暴力 用来求简单的数论函数的前缀和,像φ,μφ,\mu这类的东西当然,约数和,约数个数之类的也是可以的Text数论函数是指定义域是整数,陪域是复数的函数Dirichlet 卷积定义两个数论函数f,gf,g 它们的狄利克雷卷积表示f∗gf*g,设卷起来得到的新函数是hh h(i)=∑d|if(i)g(id)h(i)=\sum\li原创 2017-07-13 21:37:03 · 3656 阅读 · 4 评论 -
[51Nod 1778] 小Q的集合
Description小Q有一个集合 S ,它的元素个数 |S|=n|S|=n|S|=n 。 对于 S 的任意一个子集合 T ,定义 f(T)=|T|kf(T)=|T|kf(T)=|T|^k ,定义 T 关于 S 的补集为 S−T 。 小Q想知道,如果他等概率地选择一个 S 的子集 T ,那么 f(T)−f(S−T)f(T)−f(S−T)f(T)−f(S−T) 的方差是多少。 由于这个方...原创 2018-03-15 12:32:10 · 414 阅读 · 0 评论 -
[JZOJ5079] 抛硬币
1原创 2018-03-24 15:17:49 · 361 阅读 · 0 评论 -
[BZOJ3817] Sum
Description给出n,r 求∑i=1n(−1)⌊ir√⌋∑i=1n(−1)⌊ir⌋\Large \sum\limits_{i=1}^{n}(-1)^{\lfloor i\sqrt r\rfloor} T组数据T,r<=10000,n<=10^9Solution转化式子 (−1)a=1−2×(amod2)(−1)a=1−2×(amod2)(-1)^a=1...原创 2018-03-14 17:22:47 · 560 阅读 · 2 评论 -
【数论】[51Nod 1236] 序列求和 V3【无实现】
Description求斐波那契数列前N项的K次方和 N<=10^18,K<=100000 答案对1000000009取模Solution考虑斐波那契通项公式 F(i)=15√[(1+5√2)i−(1−5√2)i]F(i)={1 \over \sqrt 5}[({1+\sqrt 5\over 2})^i-({1-\sqrt 5\over 2})^i]设a=1+5√2,b=1−5√2a={1+\原创 2018-04-23 21:29:22 · 262 阅读 · 0 评论 -
【数论】[51Nod 1847] 奇怪的数学题
Description求∑i=1n∑j=1nsgcd(i,j)m∑i=1n∑j=1nsgcd(i,j)m\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{n}sgcd(i,j)^m 其中sgcd(i,j)sgcd(i,j)sgcd(i,j)表示i,j的次大公约数答案对2^32取模 n≤109,m≤50n≤109,m≤50n\leq 10^9,m\leq 5...原创 2018-05-17 17:40:55 · 219 阅读 · 0 评论 -
费马小定理、欧拉定理与扩展欧拉定理(含证明)
这里就以自己做好的PPT图片的形式给出了:原创 2018-06-14 15:42:43 · 22792 阅读 · 10 评论 -
【数论】【容斥】[JZOJ4392] 幂
Description Solution定义i的指数为cicic_i,满足ici≤Aici≤Ai^{c_i}\leq A且ici+1>Aici+1>Ai^{c_i+1}>A如果对于一个iii,存在一个数jjj,使得jp=i,p≤ci,p∈N+jp=i,p≤ci,p∈N+j^p=i,p\leq c_i,p\in N^+,那么我们可以将它们合在一起做,因为只有这样它们...原创 2018-07-18 21:57:09 · 259 阅读 · 0 评论 -
【杂题】【数论】[JZOJ5134] 三元组(无实现)
Description求∑i=1A∑j=1B∑k=1C[(i,j)=1&nbsp;and&nbsp;(j,k)=1&nbsp;and&nbsp;(i,k)=1]∑i=1A∑j=1B∑k=1C[(i,j)=1&nbsp;and&nbsp;(j,k)=1&nbsp;and&nbsp;(i,k)=1]\sum\limits_{i=1}^{A}\sum\limits_{j=1}^{B}原创 2018-08-21 21:41:41 · 220 阅读 · 0 评论 -
【杂题】[51Nod 1238] 最小公倍数之和 V3【数论】【杜教筛】
Description求∑i=1n∑j=1nlcm(i,j)\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{n} lcm(i,j)i=1∑nj=1∑nlcm(i,j)n≤1010n\leq 10^{10}n≤1010Solution将lcmlcmlcm拆开原式化为∑i=1n∑j=1nijgcd(i,j)\sum\limits_{i=1}^{n}\...原创 2018-10-24 17:21:49 · 204 阅读 · 0 评论 -
[JZOJ6022]【GDOI2019模拟2019.2.15】求和【排列组合】【中国剩余定理】
DescriptionSolution∑i=0N∑j=0M[i,j&nbsp;is&nbsp;even](ij)\sum\limits_{i=0}^{N}\sum\limits_{j=0}^{M}[i,j\ is\ even]{i\choose j}i=0∑Nj=0∑M[i,j&nbsp;is&nbsp;even](ji)=∑j=0M[j&nbsp;is&nbsp;even]∑i=...原创 2019-02-15 22:42:35 · 322 阅读 · 0 评论 -
[JZOJ4919] 神炎皇
Description 神炎皇乌利亚很喜欢数对,他想找到神奇的数对。 对于一个整数对(a,b),若满足a+b<=n且a+b是ab的因子,则成为神奇的数对。请问这样的数对共有多少呢? 对于100%的数据n<=100000000000000。Solution先讲一下我的SB思路 我乱搞弄出这么一个式子 若a+b|aba+b|ab,那么a,ba,b一定可以表示为 pn(n+q),p原创 2016-12-12 19:25:04 · 685 阅读 · 0 评论 -
Miller_Rabin素数测试与Pollard_Rho分解质因数
Miller_Rabin测试如果需要快速测试一个数是否是素数,有筛法与试除法此处介绍的是一种基于费马小定理的不确定性算法,当然,这种算法的出错率是极其微小的,尤其当选择的测试数较多时,因此在OI比赛中成为一种实用的算法好的进入正题。我们知道费马小定理,若pp为素数,那么对于a为任意正整数且a不为p的倍数时,有ap−1≡1(modp)a^{p-1}\equiv 1\pmod p然而它的逆定理并不成立,原创 2017-11-08 21:48:20 · 480 阅读 · 0 评论 -
[JZOJ4161] 于神之怒
Description Solution又是一道反演题。设gcd(i,j)=d,f(d)=∑∑[gcd(i,j)=d]gcd(i,j)=d,f(d)=\sum\sum[gcd(i,j)=d]设n<mn<m 然后就是套路。 g(d)=∑d|i,i≤nf(i)=(所有能被d整除的数对数)=⌊ni⌋⌊mi⌋g(d)=\sum\limits_{d|i,i\leq n}f(i)=(所有能被d整除的数原创 2016-12-23 19:54:15 · 389 阅读 · 0 评论 -
[51Nod 1594] Gcd and Phi
Description 求∑i=1n∑j=1nφ(gcd(φ(i),φ(j)))\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{n}φ(gcd(φ(i),φ(j))) n<=2*10^6Solution原创 2016-12-24 08:15:16 · 510 阅读 · 0 评论 -
[51Nod 1223] 分数等式的数量
Description 有这样一个分数等式:1/X + 1/Y = 1/N,(X,Y,N > 0)。给出L,求有多少满足X < Y <= L的等式。 例如:L = 12,满足条件的等式有3个,分别是:1/3 + 1/6 = 1/2, 1/4 + 1/12 = 1/3, 1/6 + 1/12 = 1/4。 L<=10^11Solution化简一下式子。 1/x+1/y=1/N1/x原创 2017-01-04 20:53:32 · 563 阅读 · 2 评论 -
[JZOJ4496] 【GDSOI 2016】互补约数
Description求∑d=1N∑i|dgcd(i,di)\sum\limits_{d=1}^{N}\sum\limits_{i|d} gcd(i,{d\over i})SolutionLF同志说过一句话。 “看到GCD,我们就要坚信它能够反演!”原式显然可以化。 ∑i=1N∑j=1⌊Ni⌋gcd(i,j)\sum\limits_{i=1}^{N}\sum\limits_{j=1}^{\lfl原创 2016-12-30 18:24:11 · 464 阅读 · 0 评论 -
欧拉函数一个神奇性质
Surface本知识产权来自LYD729、HOWARLI以及Alan_CtyText∑d|nφ(d)=n\sum\limits_{d|n}φ(d)=n证明(方法来自HowarLi)当n=1,显然成立 当n=pkn=p^k,pp为质数 ∑d|nφ(d)=1+()\sum\limits_{d|n}φ(d)=1+()原创 2016-12-26 21:08:37 · 425 阅读 · 0 评论 -
[JZOJ3507] 游乐场
Description 今天是个好日子,小A和他的小伙伴们一起去逛游乐园。这时,游乐园中忽然出现了一个伪装的吸血鬼,小A和他的小伙伴们都惊呆了!小伙伴们马上跑向了游乐园的四面八方。当“吸血鬼”回家吃饭的时候,小A才发现他已经和他的小伙伴们走散了。小A是个路痴,所以他只好站在原地等小伙伴们回来。 我们可以将游乐园视为一个N行M列的矩形,最上面一行为第1行,最左边一列为第1列。每个小伙伴手原创 2017-02-10 22:45:26 · 346 阅读 · 0 评论 -
[JZOJ5163] 【NOIP2017模拟6.25】PS的烦恼
Description 话说PS总是有着各种各样的烦恼,这天,他又在为自己失败的感情史烦恼着。这时,他心中的女神,魔法少女小圆从天而降,她对他说,如果你能帮我解决一个问题,我就让你永远没有烦恼。 问题是这样的: 寻找一个最大的k,使得存在一个x使得x^k=y,那么f(y)=k,即y最多可以开k次方根。 小圆的要求是求出从a到b的f值之和(包括a和b)。 100%的数据满足原创 2017-06-25 16:38:28 · 512 阅读 · 0 评论 -
类欧几里得算法(部分)
Preface欧几里得算法,就是辗转相除法。 gcd(i,j)=gcd(j,i%j)定义定义函数 F(a,b,c,n)=∑i=0n⌊ai+bc⌋F(a,b,c,n)=\sum\limits_{i=0}^n \lfloor {ai+b\over c}\rfloor推导一波显然当a≥ca\geq c或者b≥cb\geq c时,F(a,b,c,n)=∑i=0n(⌊(amodc)i+(bmodc)c⌋+原创 2017-08-09 10:05:50 · 1453 阅读 · 2 评论 -
[JZOJ5250]【GDOI2018模拟8.11】质数
Description求∑i=1n2f(i)\sum\limits_{i=1}^{n}2^{f(i)},其中f(i)表示i的不同质因子个数Solution显然,2^f(i)就是每个质因子选或者不选那么2^f(i)等价于∑d|i[gcd(d,i/d)=1]\sum\limits_{d|i}[gcd(d,i/d)=1],两个互质一定不会有相同的设c(q)=∑d|i[gcd(d,i/d)=q]c(q)=\原创 2017-08-12 08:12:13 · 446 阅读 · 1 评论 -
[JZOJ5261]【GDOI2018模拟8.12】求和
Description Solution里面随便变换一下 原式化为 ∑d=1k∑i=1nfd(i)∑x=1⌊ni⌋∑y=1⌊ni⌋[gcd(x,y)=1]\sum\limits_{d=1}^{k}\sum\limits_{i=1}^{n}f_d(i)\sum\limits_{x=1}^{\lfloor {n\over i}\rfloor}\sum\limits_{y=1}^{\lfloor {原创 2017-08-15 21:19:54 · 332 阅读 · 0 评论 -
[51Nod 2026] Gcd and Lcm
Description设f(x)=∑i|xμ(i)if(x)=\sum\limits_{i|x}\mu(i)i 求 ∑i=1n∑j=1nf[gcd(i,j)]f[lcm(i,j)]\sum\limits_{i=1}^{n}\sum\limits_{j=1}^n f[\gcd(i,j)]f[lcm(i,j)]n≤1e9n\leq 1e9Solution其实可以发现,f是个积性函数。。。把gcd(i原创 2017-10-21 21:28:26 · 150 阅读 · 0 评论 -
[JZOJ5457]【NOIP2017提高A组冲刺11.6】项链
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-08 16:47:36 · 473 阅读 · 0 评论 -
[小结] 数论题的式子的转化
Text设g(i)为i的约数个数 它的求和∑i=1ng(i)\sum\limits_{i=1}^{n}g(i)等价于∑i=1n⌊ni⌋\sum\limits_{i=1}^{n}\lfloor {n\over i}\rfloor 即单独考虑每个约数的贡献经典的式子∑i=1n∑j=1n[gcd(i,j)=1]\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{n}[g原创 2017-08-14 08:11:44 · 419 阅读 · 1 评论 -
[JZOJ4367]【GDKOI2016】小学生数学题(口胡)
Description Solution既然保证最后答案有逆元,那么将原式拆成有逆元和没有逆元两部分 即∑a=1⌊np⌋1a∗p+∑i=0⌊np⌋−1∑j=1p−11i∗p+j(modpk)\sum\limits_{a=1}^{\lfloor {n\over p}\rfloor}{1\over a*p}+\sum\limits_{i=0}^{\lfloor {n\over p}\rfloor-1原创 2017-12-19 17:37:07 · 417 阅读 · 0 评论 -
[JZOJ4695]【GDOI2017模拟8.14】佐助的难题
Description 求在1到n!范围内,与m!互质的数的数量,由于答案太大,只需计算答案对R取模之后的答案即可,保证R是一个质数 对于100%的数据n<=10000000,T<=10000 保证R为质数,m<=n,n < RSolution对于求原创 2017-12-14 16:10:28 · 273 阅读 · 0 评论 -
[JZOJ5395]【NOIP2017提高A组模拟10.6】Count
Description 答案对998244353取模Solution有公式 ∑i=1ni[gcd(i,n)=1]=φ(n)∗n/2\sum\limits_{i=1}^{n} i[\gcd(i,n)=1]=φ(n)*n/2证明很简单 因为若gcd(i,j)=1,i<j\gcd(i,j)=1,i<j,一定有gcd(i,j−i)=1gcd(i,j-i)=1所以上面的式子显而易见。那么f(i)=i原创 2017-10-06 21:39:23 · 358 阅读 · 0 评论 -
[JZOJ5214]【HEOI、SXOI2017】相逢是问候(口胡)
Description Solution指数是不能取模的,这很尴尬但是有欧拉定理 若gcd(c,p)=1gcd(c,p)=1,那么cx≡cxmodφ(p)modpc^x\equiv c^{x \mod φ(p)} \mod p因为根据欧拉定理cφ(p)≡1modpc^{φ(p)}\equiv1 \mod p但是此处c,p并不一定互质有欧拉定理的扩展(为什么还有这种操作。。。。) cx≡c(xm原创 2017-07-09 21:48:53 · 615 阅读 · 1 评论 -
【数论模板】二次剩余Cipolla算法,离散对数BSGS 算法
CipollaLL ksm(LL k,LL n){ LL s=1; for(;n;n>>=1,k=k*k%mo) if(n&1) s=s*k%mo; return s;}namespace number{ LL D; struct Z { LL x,y; Z(LL _x=0,LL _y=0){x=_x,y=_y;} }; Z operator...原创 2019-05-09 11:28:10 · 290 阅读 · 0 评论