数论
文章平均质量分 77
Flyppy_White
这个作者很懒,什么都没留下…
展开
-
HDU6434 Problem I. Count(欧拉函数)
HDU6434 Problem I. Count(欧拉函数)题目大意求∑i=1n∑j=1i−1[gcd(i−j,i+j)=1](n≤2e7)\sum_{i=1}^n\sum_{j=1}^{i-1}[gcd(i-j,i+j)=1](n\le 2e7)i=1∑nj=1∑i−1[gcd(i−j,i+j)=1](n≤2e7)T(T≤1e5T\le 1e5T≤1e5)组数据.解题思路(误)...原创 2019-12-04 14:25:11 · 328 阅读 · 0 评论 -
Loj 6053. 简单的函数(Min_25筛)
Loj 6053. 简单的函数(Min_25筛)题目大意有一积性函数,当p为质数时有f(pc)=p XOR cf(p^c)=p\ XOR\ cf(pc)=p XOR c求∑i=1nf(i) mod (109+7)\sum_{i=1}^nf(i)\ mod\ (10^9+7)∑i=1nf(i) mod (1...原创 2019-05-28 00:20:38 · 229 阅读 · 0 评论 -
min25筛学习笔记&模板详解
min25筛算法推导&模板详解问题引入引入函数σ0(n)\sigma _0(n)σ0(n)为n的正因数的数量,求S(n,k)=∑i=1nσ0(ik)S(n,k)=\sum_{i=1}^n\sigma_0(i^k)S(n,k)=∑i=1nσ0(ik),n,k≤1010n,k\le 10^{10}n,k≤1010简化分析与引入对于σ0(n)\sigma_0(n)σ0(n)函数,...原创 2019-05-22 23:20:35 · 3140 阅读 · 2 评论 -
P5325 【模板】Min_25筛(min25筛)
P5325 【模板】Min_25筛(min25筛)题目大意定义积性函数f(x),且f(pk)=pk(pk−1)f(p^k)=p^k(p^k-1)f(pk)=pk(pk−1)(pp是一个质数),求∑i=1nf(x)\sum_{i=1}^n f(x)i=1∑nf(x)对10^9+7取模。解题思路min25筛版题https://blog.csdn.net/baiyifeifei/...原创 2019-05-24 12:25:25 · 535 阅读 · 0 评论 -
HDU5528 Count a b(欧拉函数&数的分解)
HDU5528 Count a * b(欧拉函数&数的分解)题目大意设f(n)f(n)f(n)为小于n的a,b使得有a×b mod n≠0a\times b\ mod \ n =\not 0a×b mod n≠0的a,b的组数。现定义g(n)=∑d∣nf(d)g(n)=\sum_{d|n}f(d)g(n)=∑d∣nf(d),现给出n求g(...原创 2019-05-08 00:30:10 · 194 阅读 · 0 评论 -
CF 920F SUM and REPLACE(线性递推&线段树)
CF 920F SUM and REPLACE(线性递推&线段树)题目大意设以函数D(i)D(i)D(i)为i的因数的数量。现在给出一段序列{a}。有两种操作,操作1为将l和r之间的所有的数aia_iai变为D(ai)D(a_i)D(ai),另一种操作求l,r之间aia_iai的和解题思路首先线性地递推出 所有的D(i)D(i)D(i),用一棵线段树维护这一段序列,每次更新一...原创 2019-05-03 00:47:49 · 161 阅读 · 0 评论 -
CF 842C Ilya And The Tree(数论)
CF 842C Ilya And The Tree(数论)题目大意给出一棵树,其上每个节点有一个对应的值,将一个节点的美丽值定义为这个点到根节点的路径上所有的节点的值的gcd,现在可以将一些节点的值变为0,问每个节点最大的美丽值为多少解题思路分为两种情况:是否设定根节点为0将根节点变为0,则直接求所有的节点到根节点的gcd即可不将根节点设定为0,则每个节点的美丽值必然为根节点的因数,设...原创 2019-05-03 00:38:25 · 323 阅读 · 0 评论 -
南昌网络赛tsy's number(莫比乌斯反演&线性递推)
南昌网络赛tsy’s number(莫比乌斯反演&线性递推)题目大意给出n,求∑i=1n∑j=1n∑j=1nφ(i)φ(j2)φ(k3)φ(i)φ(j)φ(k)φ(gcd(i,j,k))\sum_{i=1}^n\sum_{j=1}^n\sum_{j=1}^n{\varphi(i)\varphi(j^2)\varphi(k^3)\over\varphi(i)\varphi(j)\va...原创 2019-04-26 13:25:12 · 599 阅读 · 0 评论 -
CF896D Nephren Runs a Cinema(卡特兰数&不互素数取模下除法)
Nephren Runs a Cinema(卡特兰数&不互素数求逆元)题目大意要求构造长度为n的序列{an}\{a_n\}{an}(令其和序列为SnS_nSn)使得有{∀i Si≥0∀iai=0 or 1 or−1l ≤Sn≤r\begin{cases}\forall_i \ S_i\ge 0\\\forall_i a_i=...原创 2019-04-13 01:07:41 · 589 阅读 · 0 评论 -
CF906D Power Tower(欧拉降幂)
CF906D Power Tower(欧拉降幂)题目大意给出一段有n个数的序列a1,a2,...ana_1,a_2,...a_na1,a2,...an给出q个询问.第i个询问给出li,ril_i,r_ili,ri,每组询问求alal+1...ar mod ma_l^{a_{l+1}^{...^{a^r}}}\ mod\ malal+1...ar ...原创 2019-04-08 22:56:02 · 678 阅读 · 4 评论 -
Congruence Equation(数论)
Congruence Equation(数论)题目大意给出a,b,p,x求解方程n⋅an≡b(mod p)n\cdot a^n\equiv b(mod\ p)n⋅an≡b(mod p)中所有n≤xn\le xn≤x 的解的数量解题思路显然n⋅ann\cdot a^nn⋅an有循环节p⋅(p−1)p\cdot(p-1)p⋅(p−1),因此找出最小的解minn之后(...原创 2019-04-04 10:44:39 · 659 阅读 · 0 评论 -
Uoj #188. 【UR #13】Sanrd(min_25筛)
Uoj #188. 【UR #13】Sanrd(min_25筛)题目大意现在有一函数f(x)f(x)f(x)其函数值为x的次大素因数,如为素数则为0,现求∑i=lrf(i)\sum_{i=l}^rf(i)∑i=lrf(i)解题思路注意到min_25筛的更新中有一步为g(n,m)=∑m<p≤n,pe≤n,e≥1,p∈Pf(pe)([e>1]+g(⌊npe,p⌋...原创 2019-05-28 00:35:11 · 246 阅读 · 0 评论 -
HDU6537(2019湘潭邀请赛F)Neko and functio(min_25筛+容斥原理)
HDU6537(2019湘潭邀请赛F)Neko and functio(min_25筛+容斥原理)题目大意给出一个函数f(n,k)f(n,k)f(n,k)其值为选出k个大于1的数字使得有∑i=1kai=n(ai>1)\sum_{i=1}^ka_i=n(a_i>1)∑i=1kai=n(ai>1)的方案数。先对其求前缀和。需要注意的是,不同的顺序,为不同的方...原创 2019-07-05 00:12:57 · 452 阅读 · 0 评论 -
杭电多校第一场1007Meteor(类欧几里得&莫比乌斯反演&Stern-Brocot树&分数二分 )
杭电多校第一场1007Meteor(类欧几里得&莫比乌斯反演&Stern-Brocot树&分数二分 )题目大意求分子分母均小于n且互素的第K小的数字解题思路AC代码#include<bits/stdc++.h>#define int long longusing namespace std;typedef long long LL;struc...原创 2019-07-23 22:26:16 · 278 阅读 · 0 评论 -
Comet OJ - Contest #10 C 鱼跃龙门(质因数分解&扩展欧几里得)
Comet OJ - Contest #10 C 鱼跃龙门(质因数分解&扩展欧几里得)题目大意给定一个正整数 n,一共有 n 座龙门,跳过第 j (j<n) 座龙门将会到达第 j+1座龙门前,特殊地,跳过第 n座龙门后将会到达第 1 座龙门前。胖头鱼一开始在第一座龙门前,接下来,第 i 个时刻内它会向前跳 i 次,每次跳过 1 座龙门,求最小的正整数 x 满足第 x 个时刻结束...原创 2019-09-09 20:07:23 · 441 阅读 · 1 评论 -
luogu P2480 [SDOI2010]古代猪文(Lucas定理&中国剩余定理)
luogu P2480 [SDOI2010]古代猪文(Lucas定理&中国剩余定理)题目大意iPig打算研究古时某个朝代的猪文文字。根据相关文献记载,那个朝代流传的猪文文字恰好为远古时期的k分之一,其中k是N的一个正约数(可以是1和N)。不过具体是哪k分之一,以及k是多少,由于历史过于久远,已经无从考证了。iPig觉得只要符合文献,每一种能整除N的k都是有可能的。他打算考虑到所有可能...原创 2019-08-20 23:13:41 · 133 阅读 · 0 评论 -
luogu P3327 [SDOI2015]约数个数和(rng58-clj等式)
luogu P3327 [SDOI2015]约数个数和(rng58-clj等式)题目大意设d(x)d(x)d(x)为x的因子的个数,给出N和M求∑i=1n∑j=1md(i,j)\sum_{i=1}^n\sum_{j=1}^md(i,j)i=1∑nj=1∑md(i,j)解题思路根据rng58-clj等式,有∑i=1n∑j=1md(i,j)=∑i=1n∑j=1m[gcd(i,j)=...原创 2019-08-14 00:16:27 · 284 阅读 · 0 评论 -
HDU5608 function(杜教筛)
HDU5608 function(杜教筛)题目大意给出一个函数其满足N2−3N+2=∑d∣Nf(d)N^2-3N+2=\sum_{d|N}f(d)N2−3N+2=d∣N∑f(d)求∑i=1nf(i)\sum_{i=1}^nf(i)∑i=1nf(i)解题思路很明显的杜教筛,它甚至帮你把需要构造的卷积给出来了,对着做就行了。可以发现原函数本质就是莫比乌斯函数卷积上那个多项式,为了方...原创 2019-08-14 00:15:15 · 218 阅读 · 0 评论 -
2019牛客多校第7场 K.Function(min_25筛)
2019牛客多校第7场 K.Function(min_25筛)题目大意定义一个函数csl(p,x)={3e+1x=pe and ∃a,b p=a2+b21x=pe and ̸∃a,b p=a2+b20othercsl(p,x)=\begin{cases}3e+1&x=p^e\ and\ \exist a,b\ p...原创 2019-08-08 23:03:00 · 532 阅读 · 0 评论 -
牛客多校第三场 D Big Integer(欧拉函数&计数)
牛客多校第三场 D Big Integer(欧拉函数&计数)题目大意给出一个函数A(n)=∑i=0n−110iA(n)=\sum_{i=0}^{n-1}10^iA(n)=∑i=0n−110i,现在给出一个素数p,以及范围n,m.求有多少(i,j)(i,j)(i,j)满足1≤i≤n,1≤j≤m,A(ij)≡0(mod p)1\le i\le n,1\le j\le m,A(i...原创 2019-07-25 23:35:23 · 191 阅读 · 0 评论 -
2019牛客多校第5场 generator 2(BSGS)
2019牛客多校第5场 generator 2(BSGS)题目大意已知x0x_0x0,a,b,p,现有递推式xn=(xn−1∗a+b)%px_n=(x_{n-1}*a+b)\%pxn=(xn−1∗a+b)%p,每次询问一个v,问最小的n使得有xn=vx_n=vxn=v。解题思路我们首先可以发现这个数列可以写出通项xn=x0⋅an+ban−1a−1(mod p...原创 2019-08-03 00:45:07 · 299 阅读 · 3 评论 -
HDU6607Easy Math Problem(min_25筛+杜教筛+拉格朗日插值)
HDU6607Easy Math Problem(min_25筛+杜教筛+拉格朗日插值)题目大意求∑i=1n∑j=1ngcd(i,j)klcm(i,j)[gcd(i,j)∈prime]\sum_{i=1}^n\sum_{j=1}^ngcd(i,j)^klcm(i,j)[gcd(i,j)\in prime]i=1∑nj=1∑ngcd(i,j)klcm(i,j)[gcd(i,j)∈pri...原创 2019-07-30 17:00:06 · 787 阅读 · 0 评论 -
一些与gcd有关的前缀和公式
一些与gcd有关的前缀和公式给定n求小于等于n且与n互素的数之和∑d=1nd[gcd(d,n)=1]\sum_{d=1}^nd[gcd(d,n)=1]d=1∑nd[gcd(d,n)=1]这个问题可以通过高斯求∑i=1ni\sum_{i=1}^ni∑i=1ni的方式来求解可以发现,当n≥2n\ge2n≥2时,有gcd(d,n)=1gcd(d,n)=1gcd(d,n)=1则必然有gcd...原创 2019-07-30 00:43:37 · 769 阅读 · 0 评论 -
杭电多校第一场1011Function(莫比乌斯反演)
杭电多校第一场1011Function(莫比乌斯反演)题目大意求∑i=1ngcd(⌊i3⌋,i) mod 998244353\sum_{i=1}^ngcd(\left\lfloor\sqrt[3]{i}\right\rfloor,i)\ mod\ 998244353i=1∑ngcd(⌊3i⌋,i) mod 998244353解题思路以下极大...原创 2019-07-23 23:00:25 · 268 阅读 · 0 评论 -
CF 687B Remainders Game(质因数分解&扩展中国剩余定理应用)
CF 687B Remainders Game(质因数分解&扩展中国剩余定理应用)题目大意现选定一个k与x,x未知,给出n个数c,可否根据x与c之间模数得出x模k?解题思路据扩展中国剩余定理,可以知道我们总可以将两个同余式子{x≡a1(mod m1)x≡a2(mod m2)\begin{cases}x\equiv a_1 (mod\ m_1)\\x\equ...原创 2019-03-13 16:40:55 · 372 阅读 · 0 评论 -
Divisors CodeForces - 1033D(数论)
Divisors CodeForces - 1033D(数论)题目大意给出n个数,每个数有3到5个因数,问n个数的积有多少因数解题思路由于每个数有3-5个因子,在除去了数本身之外就只有这个几种情况三个因子:为完全平方数四个因子:为完全三次方数,或为两个不同的素数的乘积五个因子:为完全四次方数首先对原数进行开二,三,四方,若能开完全则直接分解,否则进行进一步的分解对需要进一步分...原创 2019-03-02 18:44:21 · 464 阅读 · 0 评论 -
上帝与集合的正确用法 HYSBZ - 3884(欧拉降幂)
上帝与集合的正确用法 HYSBZ - 3884(欧拉降幂)题目大意给出一个n求无限个2{2222...%n无限个2\begin{cases}2^{2^{2^{2^{...}}}}\end{cases}\%n无限个2{2222...%n解题思路对式子应用欧拉降幂f(n)=222...%n=222...%ϕ(n)+ϕ(n)%nf(n)=2^{2^{2^{...}}}\%n\\=2...原创 2019-02-21 16:45:50 · 162 阅读 · 0 评论 -
CodeForces - 582C Superior Periodic Subarrays(数论&计数)
CodeForces - 582C Superior Periodic Subarrays(数论&计数)题目大意给出一段纯循环序列a,循环节的长度为n.求有多少对数对(l,s)总是满足以s为循环节,循环节中数字为al,al+1..al+s−1a_l,a_{l+1}..a_{l+s-1}al,al+1..al+s−1的序列设为b使得总是有bi−l+1≥aib_{i-l+1}\ge ...原创 2019-02-15 19:25:58 · 227 阅读 · 0 评论 -
Fox And Jumping CodeForces - 510D(数论&记忆化搜索)
Fox And Jumping CodeForces - 510Dhttp://codeforces.com/contest/510/problem/D题目大意给出n个片段(n≤300n\le 300n≤300)每个片段有一定的长度与价值,购买下一个片段后可以在任意位置x通往x+l或x-l处,l为片段的长度,问最少需要花费多少钱才可以保证可以到达任意一个坐标点解题思路题目大意即为求花费尽...原创 2019-02-07 15:17:03 · 217 阅读 · 0 评论 -
CodeForces - 492E Vanya and Field(逆元)
CodeForces - 492E Vanya and Field(逆元)题目大意给出一个n×nn\times nn×n的坐标纸,其上有m棵苹果树,给出dx,dydx,dydx,dy假设当前位置为x,yx,yx,y则每次行动之后的位置会变成(x+dx)mod&nbsp;n,(y+dy)mod&nbsp;n(x+dx)mod\ n,(y+dy)mod\ n(x+dx)mod&nbsp;n,(y+...原创 2019-02-15 15:00:36 · 119 阅读 · 0 评论 -
CodeForces - 225E Unsolvable(梅森素数)
CodeForces - 225E Unsolvable(梅森素数)CodeForces - 225E 题目大意问不满足条件z=[x2]+y+xyz=[\frac{x}{2}]+y+xyz=[2x]+y+xy的数字z从小到大第n个是什么解题思路对x分开奇偶考虑当x为偶数时,不妨令x=2kz=k+y+2ky2z+1=2k+2y+4ky+12z+1=(2k+1)(2y+1)2k...原创 2019-02-14 20:13:31 · 278 阅读 · 3 评论 -
Fibonacci Number CodeChef - FN (BSGS+cipolla算法)
CodeChef - FN Fibonacci Number(BSGS+cipolla)CodeChef - FN 题目大意给出一个数问其最小与斐波那契数列数列的第几项同余于p解题思路斐波那契数列有着通项Fn=15[(1+52)n−(1−52)n]F_n=\frac{1}{\sqrt5}[(\frac{1+\sqrt5}{2})^n-(\frac{1-\sqrt5}{2})^n]F...原创 2019-02-14 16:34:38 · 296 阅读 · 0 评论 -
A Modular Arithmetic Challenge POJ - 3530(递归求解同余不等式)
A Modular Arithmetic ChallengePOJ - 3530 https://vjudge.net/problem/POJ-3530DescriptionGiven four positive integers M, D, L, R, your task is to find the minimal non-negative integer x which satisfi...原创 2019-02-01 00:33:28 · 531 阅读 · 2 评论 -
The Luckiest number POJ - 3696(欧拉函数)
The Luckiest number POJ - 3696题目描述Chinese people think of ‘8’ as the lucky digit. Bob also likes digit ‘8’. Moreover, Bob has his own lucky number L. Now he wants to construct his luckiest number wh...原创 2019-01-26 22:03:03 · 174 阅读 · 0 评论 -
Taxes CodeForces - 735D(哥德巴赫猜想)
Taxes CodeForces - 735D题目描述Mr. Funt now lives in a country with a very specific tax laws. The total income of mr. Funt during this year is equal to n (n ≥ 2) burles and the amount of tax he has to p...原创 2019-01-26 21:50:42 · 264 阅读 · 0 评论 -
二次剩余Cipolla算法
欧拉准则欧拉准则:设p为奇素数,则ap−12≡(ap)(mod p)a^{\frac{p-1}{2}}\equiv (\frac{a}{p})(mod\ p)a2p−1≡(pa)(mod p)其中ap\frac{a}{p}pa为勒让德符号证明:假设a为二次剩余,如a≡b2(mod p)a\equiv b^2(mod\ p)a≡b2(mod ...原创 2019-01-26 18:50:38 · 1073 阅读 · 2 评论 -
CodeForces - 488E Prefix Product Sequence(逆元&威尔逊定理)
CodeForces - 488E Prefix Product Sequence(逆元&威尔逊定理)题目大意给出一个1到n的置换a[1…n]使得序列a1mod n,a1a2mod n,...,∏i=1nai mod n{a_1 mod\ n,a_1a_2mod\ n,...,\prod_{i=1}^na_i\ mod\ n}a1mod&nbs...原创 2019-02-15 21:40:40 · 253 阅读 · 0 评论 -
Mister B and Astronomers CodeForces - 819D(数论)
Mister B and Astronomers CodeForces - 819D(数论)题目大意有n个观察员,第一个观察员在0秒开始观察星空,随后第i个观察员会在第i-1个观察员之之后aia_iai秒进行观察,第一个观察员也会在第n个观察员观察后a1a_1a1秒观察,有一颗星星其闪烁的周期为T,其第一次闪烁的时间不确定.问每个观察员有多少种可能成为第一个观察到这颗星星的人解题思路假...原创 2019-02-16 22:20:12 · 244 阅读 · 0 评论 -
codeforce 1114E Arithmetic Progression(数论&随机数生成)
codeforce 1114E Arithmetic Progression(数论&amp;amp;随机数生成)题目大意交互题.现有一等差数列,现在需要通过程序进行提问,根据提问所得的答案得到等差数列的首项和公差.现可对评测姬进行最多60次,两种提问.提问方式一为给出一数x,数列中是否存在严格大于这个数的.提问方式二为询问数列第i项为多少解题思路首先通过二分法,找出数列的最大值.然后用完剩下的询问...原创 2019-02-11 21:18:57 · 281 阅读 · 0 评论 -
Divisors CodeForces - 1033D(pollard_rho算法)
Divisors CodeForces - 1033D(pollard_rho算法)题目大意给出n个数,每个数有3到5个因数,问n个数的积有多少因数解题思路pollard_rho算法分解质因数,然后答案就是∏i∣a(cnt[i]+1)\prod_{i|a}(cnt[i]+1)∏i∣a(cnt[i]+1)AC代码#include<bits/stdc++.h> using...原创 2019-03-02 12:39:30 · 313 阅读 · 0 评论