数论
ixRic
ixRic
展开
-
[51Nod 1584] 加权约数和(约数和函数性质 + 莫比乌斯反演) | 错题本
文章目录题目分析代码题目[51Nod 1584] 加权约数和分析先把 max{i,j}\max\{i, j\}max{i,j} 去掉:原式=2∑i=1N∑j=1iiσ(ij)−∑i=1Niσ(i2)原式=2\sum_{i = 1}^{N}\sum_{j = 1}^i i\sigma(ij)-\sum_{i = 1}^N i\sigma(i^2)原式=2i=1∑Nj=1∑iiσ(ij)−i=1∑Niσ(i2)推式子前先介绍约数和函数有一个重要性质(和约数个数一样可以转化为 gcd\gc原创 2020-09-17 10:50:44 · 225 阅读 · 0 评论 -
[SDOI 2015] 约数个数和(约数个数函数性质 + 莫比乌斯反演) | 错题本
文章目录题目分析代码题目[SDOI 2015] 约数个数和分析d(xy)=∑d1∣x∑d2∣y[gcd(d1,d2)=1]d(xy)=\sum_{d_1|x}\sum_{d_2|y}[\gcd(d_1,d_2) = 1]d(xy)=d1∣x∑d2∣y∑[gcd(d1,d2)=1] 我们证明每个因子与 gcd(d1,d2)=1\gcd(d_1,d_2) = 1gcd(d1,d2)=1 的双射关系:设 x=∏i=1npiai,y=∏i=1npibi (ai+bi≠0)x原创 2020-09-08 11:08:46 · 239 阅读 · 0 评论 -
杜教筛及其时间复杂度分析
文章目录杜教筛方法举例莫比乌斯函数欧拉函数时间复杂度杜教筛杜教筛用于求一类积性函数的前缀和,时间复杂度可以做到 O(n23)O(n^{\frac{2}{3}})O(n32)。方法设我们要求的是积性函数函数 f(x)f(x)f(x) 的前缀和 S(n)=∑i=1nf(i)S(n) = \sum_{i = 1}^n f(i)S(n)=∑i=1nf(i)。我们取一个积性函数 ggg(ggg 具体是什么由 fff 决定,但这不影响推式子)则 ∑i=1n(f∗g)(i)=∑i=1n∑d∣ig(d)f(i原创 2020-09-07 10:44:15 · 986 阅读 · 3 评论 -
[LOCAL] 过河(数学) | 错题本
文章目录题目分析代码题目一只青蛙想从河的一侧跳到另一侧。河中有一些石子,青蛙必须要经过这些石子才能够顺利过河。由于青蛙一次跳过的距离是总是正整数,我们可以把河中青蛙可能到达的点看成数轴上的一串整点:0,1,⋯ ,L0,1,\cdots,L0,1,⋯,L(其中 LLL 是河的宽度)。当青蛙的坐标大于 LLL 或者小于 000 时,我们就认为青蛙跳出了这条河。一旦青蛙跳出这条河,青蛙就不会选择回到河中。注意青蛙不一定最后要跳到河的另一侧。现在,青蛙将会从 000 坐标开始,每次可以随意选择向正方向跳跃原创 2020-09-01 10:07:02 · 230 阅读 · 0 评论 -
LOJ530 「LibreOJ β Round #5」最小倍数(二分)
文章目录题目分析代码题目「LibreOJ β Round #5」最小倍数分析令n!=p1a1p2a2⋯pkakn! = {p_1}^{a_1}{p_2}^{a_2} \cdots {p_k}^{a_k}n!=p1a1p2a2⋯pkak,那么ai=∑j=1∞⌊npij⌋a_i = \sum\limits_{j = 1}^{\infty} \left\lfloor \dfrac{n...原创 2020-03-28 11:32:25 · 299 阅读 · 0 评论 -
C++ 线性筛逆元
O(n)求1到n所有数模p的逆元。原创 2019-08-04 16:36:06 · 1488 阅读 · 2 评论 -
【51nod】欧拉函数之和(数论,杜教筛)
文章目录题目分析题目1239 欧拉函数之和分析欧拉函数φ(n)\varphi(n)φ(n)表示小于nnn的与nnn互质的数的个数。令答案f(n)=∑i=1nφ(i)f(n)=\sum\limits_{i=1}^{n}\varphi(i)f(n)=i=1∑nφ(i)。可以证明n=∑d∣nφ(d)n=\sum\limits_{d|n}\varphi(d)n=d∣n∑φ(d)。下面是一...原创 2018-12-07 13:28:31 · 1877 阅读 · 5 评论 -
【AtCoder】ABC104 We Love ABC
题目题目大意分析只有ABC的情况加上?后的情况题目beta传送门题目大意对于一个字符串TTT,规定它的ABC number是这样的一个三元组的个数:(i,j,k)(i,j,k)(i,j,k)(1≤i<j<k≤|T|1≤i<j<k≤|T|1\leq iTi=Ti=T_i= A,Tj=Tj=T_j= B且Tk=Tk=T_k= C。...原创 2018-08-10 23:54:42 · 1340 阅读 · 0 评论 -
【CodeForces】988E·Divisibility by 25
题目题目大意分析代码题目传送门题目大意给你一个整数nnn(1≤n≤10181≤n≤10181\leq n\leq 10^{18},不含前导零),你能够将任意两个相邻位置上的数互换位置,每交换一次算一步,每一次交换过后的数不能有前导零。现在要求以最少的步数使这个数变为252525的倍数,输出这个步数。分析比赛的时候疯狂地分类讨论,仿佛灵魂都被净化...原创 2018-07-25 10:03:38 · 798 阅读 · 0 评论 -
【组合数学】【逆元】AtCoder - 1974·いろはちゃんとマス目 / Iroha and a Grid
题目题目大意分析题目传送门题目大意一个H×WH×WH\times W的矩形,左下角的A×BA×BA\times B的矩形区域不能走,问左上角走到右下角的方案数模109+7109+710^9+7。分析如果无视障碍,从(x,y)(x,y)(x,y)走到(m,n)(m,n)(m,n)的方案数是Cn−ym−x+n−yCm−x+n−yn−yC_{m-x+n-y...原创 2018-07-24 11:19:05 · 1027 阅读 · 0 评论 -
【数论】无平方因子的数
问题 B(2642): 无平方因子的数时间限制: 1 Sec 内存限制: 128 MB提交: 213 解决: 55[提交][状态][我的提交]题目描述给出正整数n和m,区间[n, m]内的“无平方因子”的数有多少个?整数p无平方因子当且仅当不存在 k > 1,使得p是k2 的倍数。输入第1行:2个整数n和m (1 输出原创 2017-02-07 16:33:45 · 2615 阅读 · 0 评论 -
【数论】计算组合数
问题 B(2713): [POJ2249]计算组合数时间限制: 1 Sec 内存限制: 128 MB题目描述给定正整数n, k,计算C(n, k)。答案保证在2^31以内。输入多组数据,每组数据仅一行,即2个整数n和k (n>=1) and k (0以2个0结束输入输出对每个数据,输出对应的答案 样例输入Copy (如果复制到控制台原创 2017-02-16 13:42:06 · 857 阅读 · 0 评论 -
【数论】【欧拉函数】【快速幂】洗牌机
3333: 洗牌机时间限制: 2 Sec 内存限制: 512 MB题目描述有2n张牌,放在2n个从1到2n的有序位置上。洗牌机每次可以把第i张牌洗到p(i)的位置上。P(i)的定义如下: 问经过最少多少轮洗牌,才会使所有牌回到原来的位置。输入输入格式:有多组数据,每组数据一个整数n(n输出输出格式:对于每组数据,输出一个整数,表原创 2017-03-02 08:44:49 · 552 阅读 · 0 评论 -
【数论】【NOIP2009】Hankson的趣味题
时间限制: 3 Sec 内存限制: 64 MB题目描述Hanks博士是BT (Bio-Tech,生物技术) 领域的知名专家,他的儿子名叫Hankson。现在,刚刚放学回家的Hankson正在思考一个有趣的问题。 今天在课堂上,老师讲解了如何求两个正整数c1和c2的最大公约数和最小公倍数。现在Hankson认为自己已经熟练地掌握了这些知识,他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问原创 2017-07-09 12:36:58 · 680 阅读 · 0 评论