Number Theory
Feynman1999
理论都是灰色的,而生活之树常青^_^
展开
-
HDU6588 经典gcd求和
需要求∑i=1ngcd(⌊i3⌋,i)mod  998244353n=1021\sum_{i=1}^n\gcd(\lfloor{\sqrt[3]{i}}\rfloor,i)\mod 998244353 \quad n = 10^{21}i=1∑ngcd(⌊3i⌋,i)mod998244353n=1021思路令m=⌊n...原创 2019-08-22 21:03:21 · 610 阅读 · 0 评论 -
数论概论读书笔记 34.丢番图逼近与佩尔方程
丢番图逼近与佩尔方程求解PellPellPell方程,我们的想法是取两个使x2−Dy2x2−Dy2x^2-Dy^2具有相同数值的数对并“将它们相除”举个例子,取D=13D=13D=13,如果我们知道数对(x1,y1)=(11,3)(x1,y1)=(11,3)(x_1,y_1)=(11,3)和(x2,y2)=(14159,3927)(x2,y2)=(14159,3927)(x_2,y_2)=...原创 2018-08-28 23:08:23 · 1028 阅读 · 0 评论 -
数论概论读书笔记 33.丢番图逼近
丢番图逼近如何求出佩尔方程x2−Dy2=1x2−Dy2=1x^2-Dy^2=1的一个正整数解x,yx,yx,y?因式分解后得 (x−yD−−√)(x+yD−−√)=1(x−yD)(x+yD)=1(x-y\sqrt{D})(x+y\sqrt{D})=1 现在有一个问题,x−yD−−√x−yDx-y\sqrt{D}可以小到什么程度?如果能找到整数x,yx,yx,y使其非常小,则可以期望...原创 2018-08-28 23:07:26 · 2245 阅读 · 0 评论 -
数论概论读书笔记 32.佩尔方程
佩尔方程PellPellPell方程是指具有形式x2−Dy2=1x2−Dy2=1x^2-Dy^2=1的方程,其中DDD是一个固定的正整数且不是完全平方数假设可求得PellPellPell方程的一个解x1,y1x1,y1x_1,y_1,则可以利用上一章中对D=2D=2D=2所描述的同样方法来产生一个新的解。将已知解因式分解为 1=x21−Dy21=(x1+y1D−−√)(x1−y1D−−√...原创 2018-08-28 23:06:54 · 432 阅读 · 0 评论 -
数论概论读书笔记 31.再论三角平方数
再论三角平方数定理 (三角平方数定理) 方程x2−2y2=1x2−2y2=1x^2-2y^2=1的每个正整数解都可通过将3+22–√3+223+2\sqrt{2}自乘得到,即解(xk,yk)(xk,yk)(x_k,y_k)可以通过展开下式得到。xk+yk2–√=(3+22–√)k,k=1,2,3,...xk+yk2=(3+22)k,k=1,2,3,...x_k+y_k\s...原创 2018-08-28 23:05:51 · 577 阅读 · 0 评论 -
数论概论读书笔记 30.方程x^4+y^4=z^4
方程x4+y4=z4x4+y4=z4x^4+y^4=z^4指数为4的费马大定理主要用费马降价法来证明,具体看书原创 2018-08-28 23:05:02 · 977 阅读 · 0 评论 -
数论概论读书笔记 29.原根与指标(指标也被称为离散对数(#^.^#))
原根与指标(指标也被称为离散对数(#^.^#))原根我们清楚了。啥是指标呢?对于模13,2是它的原根,则2x mod 132x mod 132^x\ mod \ 13会取遍[1,12][1,12][1,12],x∈[1,12]x∈[1,12]x∈[1,12]比如24=16≡3 mod( 13)24=16≡3 mo...原创 2018-08-28 23:02:24 · 1098 阅读 · 0 评论 -
数论概论读书笔记 28.幂模p与原根
幂模p与原根如果aaa和ppp互素,费马小定理告诉我们,ap−1≡1 (mod p)ap−1≡1 (mod p)a^{p-1}\equiv 1\ (mod \ p)辣么这个指数p−1p−1p-1是唯一的使得结果为1的嘛我们选一些aaa和ppp看一下对于a=3,p=7a=3,p=7a=3,p=7 :指数只有为6时才取到1。 模p余1的a...原创 2018-08-28 23:01:18 · 2113 阅读 · 0 评论 -
数论概论读书笔记 27.欧拉函数与因数和
欧拉函数与因数和前面我们在说明完全数的时候,引入了因子和函数σ(x)σ(x)\sigma(x) ,通常σ(x)σ(x)\sigma(x)是大于xxx的。现在我们做个猜想,不对因子直接求和,而先求欧拉函数再求和。猜想 设d1,d2,...,drd1,d2,...,drd_1,d_2,...,d_r是整除nnn的数,其中包括1与n,则 φ(d1)+φ(d2)+φ(d3)+...+φ(...原创 2018-08-28 23:00:21 · 859 阅读 · 0 评论 -
HDU 3542 (费马降阶法)
probelmGiven an odd prime p. You are to tell whether p can be written into form a^2+b^2 (0InputThere will be up to 2*10^5 test cases. An odd prime p (3<=p<2*10^9) in each line. End w...原创 2018-08-21 10:59:54 · 513 阅读 · 0 评论 -
数论概论读书笔记 41.连分数与佩尔方程
连分数与佩尔方程通过反复地将小数部分翻到分母上并将整数部分分离,我们可以对任何一个数形成连分数:连分数在视觉上滑向右下方,将它们写成分数却要花费很多的笔墨和空间,由于所有的分子都是1,故我们要做的就是列出分母,将连分数简记为[a0,a1,a2,a3,a4,...][a0,a1,a2,a3,a4,...][a_0,a_1,a_2,a_3,a_4,...] 辣么,当我展开到an...原创 2018-08-30 23:35:07 · 898 阅读 · 0 评论 -
数论概论读书笔记 40.O,多美的一个函数
OOO,多美的一个函数这个是计算机领域中的核心,就不介绍了原创 2018-08-30 23:34:29 · 227 阅读 · 0 评论 -
数论概论读书笔记 39.斐波那契与线性递归序列
斐波那契与线性递归序列比内公式 斐波那契序列FnFnF_n用递归公式描述如下: F1=F2=1,Fn=Fn−1+Fn−2,n=3,4,5,...F1=F2=1,Fn=Fn−1+Fn−2,n=3,4,5,...F_1=F_2=1,\quad F_n=F_{n-1}+F_{n-2}\quad ,n=3,4,5,... 则斐波那契序列的第nnn项可用公式 Fn=15–√{(1+5–√2)...原创 2018-08-30 23:33:46 · 368 阅读 · 0 评论 -
数论概论读书笔记 38.二项式系数与帕斯卡三角形
二项式系数与帕斯卡三角形二项式展开的系数我们应该已经很熟悉了 (A+B)n=∑k=0nCknAn−kBk(A+B)n=∑k=0nCnkAn−kBk(A+B)^n=\sum_{k=0}^{n}C_{n}^{k}A^{n-k}B^k 帕斯卡三角形(杨辉三角)如果让二项式系数模ppp简化(其中ppp为素数),会发生什么?模ppp二项式定理 设ppp为素数二项式系数CkpC...原创 2018-08-30 23:33:14 · 1013 阅读 · 0 评论 -
数论概论读书笔记 6.线性方程定理
线性方程与最大公因数形如ax+by的最小正整数等于gcd(a,b) 欧几里得回带法可以简单证明一定有解 下面要证明解最小,因为gcd(a,b)<=min(a,b)gcd(a,b)<=min(a,b)gcd(a,b)(x1,y1)(x1,y1)(x_{1},y_{1}),它可由欧几里得回带法得到。则方程的每一个解可由(x1+k∗bg,y1−k∗ag)(x1+k∗bg,y...原创 2018-08-25 20:13:03 · 487 阅读 · 0 评论 -
数论概论读书笔记 7.因数分解与算术基本定理
因数分解与算术基本定理引理7.1 令p是素数,假设p整除乘积ab,则p整除a或p整除b(或者p既整除a也整除b) 证明略定理7.2(素数整除性质). 假设素数p整除乘积a1a2...ara1a2...ara_1a_2...a_r ,则p整除a1,a2,...,ara1,a2,...,ara_1 \quad ,a_2 \quad , ...,\quad a_r 中至少一个因数。使...原创 2018-08-26 00:10:52 · 648 阅读 · 0 评论 -
数论概论读书笔记 8.同余式
同余式整除性是很好的性质,这在勾股数组、最大公因数与素数分解中得到了体现。而同余式提供了一种描述整除性质的简便方式。 如果m|(a−b)m|(a−b)m|(a-b),我们就说a与b模m同余并记之为a≡b(mod m)a≡b(mod m)a\equiv b(mod\ m) 特别地,a≡a%m(mod m)a≡a%m(mod m)a\equ...原创 2018-08-26 00:11:58 · 905 阅读 · 0 评论 -
51nod1237(EES解法,省空间)
这里用的是最近比较流行的EES方法解决的不了解的可以看下这篇文章:https://blog.csdn.net/Feynman1999/article/details/82874491这种方法只要能找出f(pe)f(p^e)f(pe)的表达式,并且当e+1e+1e+1时可以O(1)O(1)O(1)维护,就可以用EES筛法去做对于本题而言,考虑函数S(n)=∑i=1ngcd(i,n)S(n)...原创 2018-09-29 16:55:32 · 407 阅读 · 0 评论 -
Loj 6053(EES筛法)
problem思路满足EES筛法的要求考虑f(p)f(p)f(p)即p1p^1p1是啥,当p&amp;gt;2p&amp;gt;2p&gt;2时,f(p)=p−1f(p)=p-1f(p)=p−1,而f(2)=3f(2)=3f(2)=3所以预处理出p1,p0p^1,p^0p1,p0的前缀和即可并注意,当i&amp;gt;=2i&amp;gt;=2i&gt;=2时,这个前缀和少计算了2(即...原创 2018-09-29 12:41:37 · 270 阅读 · 0 评论 -
Min_25筛(EES筛)法
Extended Eratosthenes Sieve 参考链接给出一个积性函数(一些非积性函数也可以搞一搞)fff,且f(p)f(p)f(p)为关于ppp的多项式。求S(n)=∑i=1nf(i)S(n)=\sum_{i=1}^nf(i)S(n)=∑i=1nf(i)∀&amp;amp;amp;amp;nbsp;2≤i≤n\forall \ 2\le i\le n∀&amp;amp;amp;amp;nbsp;2≤i≤n,我们可以将iii分为两类最...原创 2018-09-27 22:28:38 · 969 阅读 · 0 评论 -
18四川省赛G题 分段 数论函数
problem计算 ans=∑i=1n∑j=1i(n&nbsp;mod(i×j))ans=∑i=1n∑j=1i(n&nbsp;mod(i×j))ans =\sum^n_{i=1}\sum^i_{j=1} (n\ mod (i \times j)) 1≤n≤10111≤n≤10111 ≤ n ≤ 10^{11}思路解法一首先写作∑ni=1∑ij=1(n−⌊nij⌋∗ij...原创 2018-09-07 17:21:08 · 268 阅读 · 0 评论 -
51nod 1060 (因子最多数)
因子个数最多的数实际上,这样的数称为反素数:即给定NNN,求[1,N][1,N][1,N]之间最大的反素数(即拥有因子数目最多的数)性质: No.1 一个反素数的质因子必然是从2开始连续的质数。 No.2 p=2^t1 x 3^t2 x 5^t3 x 7^t4…..必然t1>=t2>=t3>=….暴力dfs 时间复杂度大概几个log...原创 2018-09-03 21:49:54 · 288 阅读 · 0 评论 -
POj 3358 二进制小数循环节
二进制循环小数循环节problemLet {x} = 0.a*1*a*2*a*3… be the binary representation of the fractional part of a rational number *z. Suppose that {x} is periodic then, we can write{x} = 0.a*1*a*2…*ar(a**r+1...原创 2018-09-03 20:41:16 · 622 阅读 · 0 评论 -
POJ 3696 最幸运的数字
ProblemChinese 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 which is the minimum among al...原创 2018-09-03 19:42:23 · 665 阅读 · 0 评论 -
CF757E(积性函数)
problem思路满足pq=npq=npq=n且gcd(p,q)=1gcd(p,q)=1gcd(p,q)=1的二元组个数是多少呢?显然是2w(n)2w(n)2^{w(n)}个,其中w(n)w(n)w(n)是不同的素因子的数目那么fr+1(n)=∑uv=nfr(u)+fr(v)2fr+1(n)=∑uv=nfr(u)+fr(v)2f_{r+1}(n)=\sum_{uv=n}\f...原创 2018-09-05 21:13:29 · 834 阅读 · 0 评论 -
HDU5528 积性函数
HDU 5528Marry likes to count the number of ways to choose two non-negative integers aa and bbless than mm to make a×ba×b mod m≠0m≠0. Let’s denote f(m)f(m) as the number of ways to choose two non-n...原创 2018-09-05 11:58:54 · 339 阅读 · 0 评论 -
HYSBZ - 2818
problem给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的 数对(x,y)有多少对. N1e7思路枚举素数(大概几e5个) 对于每个素数暴力算(先预处理筛出欧拉值) 时间复杂度 case∗O(素数个数)case∗O(素数个数)case*O(素数个数)代码示例#include<bits/stdc++.h>using nam...原创 2018-09-08 09:10:38 · 342 阅读 · 0 评论 -
bzoj 2705
Longge的数学成绩非常好,并且他非常乐于挑战高难度的数学问题。现在问题来了:给定一个整数N,你需要求出∑gcd(i, N)(1<=i <=N)。Input一个整数,为N。Output一个整数,为所求的答案。Sample Input6Sample Output15Hint【数据范围】对于60%的数据,0思路枚举哪些数...原创 2018-09-04 20:02:01 · 218 阅读 · 0 评论 -
数论概论读书笔记 10.同余式、幂与欧拉公式
同余式、幂与欧拉公式费马小定理很漂亮 ap−1≡1(mod&nbsp;p)ap−1≡1(mod&nbsp;p)a^{p-1}\equiv 1(mod\ p) 但限制P是素数且p∤ap∤ap\nmid a如果p是合数,即使a、p互质,结论也不正确了因此,研究是否有依赖模m的指数使得 a???≡1(mod&nbsp;m)a???≡1(mod&nbsp;m)a^{???}\equiv...原创 2018-08-26 00:14:37 · 873 阅读 · 0 评论 -
数论概论读书笔记 9.同余式、幂与费马小定理
同余式、幂与费马小定理前面我们讨论了关于同余式、同余方程的一些性质小结一下,互质这个条件我觉得很重要①对于ac≡bc(mod m)ac≡bc(mod m)ac\equiv bc(mod \ m) 如果gcd(c,m)=1gcd(c,m)=1gcd(c,m)=1 则可以消去c 得到a≡b(mod m)a≡b(mod m)a\equiv b(m...原创 2018-08-26 00:13:39 · 615 阅读 · 0 评论 -
数论概论读书笔记 37.无理数与超越数
无理数与超越数稍微偏了点理论,,这里不详细记录了主要是:2–√2\sqrt{2}的无理性定理刘维尔不等式关于ββ\beta的好的逼近的引理ββ\beta的超越性定理...原创 2018-08-30 23:32:19 · 864 阅读 · 0 评论 -
数论概论读书笔记 36.高斯整数与唯一因子分解
高斯整数与唯一因子分解稍微偏了点理论,数域的扩展,这里不详细记录了主要是:高斯整数的唯一分解高斯整数带余除法高斯整数公因数性质高斯素数整除性质勒让德两平方数之和定理D1−D3D1−D3D_1-D_3差定理...原创 2018-08-30 23:31:25 · 1257 阅读 · 0 评论 -
数论概论读书笔记 35.数论与虚数
数论与虚数稍微偏了点理论,数域的扩展,这里不详细记录了主要是:代数学基本定理高斯单位定理范数的积性高斯素数定理高斯整除引理原创 2018-08-30 23:21:01 · 348 阅读 · 0 评论 -
数论概论读书笔记 15.梅森素数与完全数
梅森素数与完全数完全数是等于其真因数之和的数。比如6和28均为完全数。欧几里得《几何原本》第9卷的命题36给出了求某些完全数的方法:15.1 欧几里得完全数公式: 如果2p−12p−12^p-1是素数,则2p−1(2p−1)2p−1(2p−1)2^{p-1}(2^p-1)是完全数证明:质因子分解一下,由等比公式易证。自然要问的是,由15.1给出的求解完全数公式的方法,是否能...原创 2018-08-26 09:36:45 · 1975 阅读 · 0 评论 -
数论概论读书笔记 14.梅森素数
梅森素数本节我们研究形如an−1(n≥2)an−1(n≥2)a^n-1(n\ge 2)的素数。例如31就是这样的数,因为31=25−131=25−131=2^5-1容易知道若a是奇数,则an−1an−1a^n-1是偶数,即不可能是素数。且a−1a−1a-1始终是an−1an−1a^n-1的因子。为啥呢?因为使用几何级求和公式,(展开即可证明该公式) xn−1=(x−1)(xn−1+x...原创 2018-08-26 09:35:57 · 454 阅读 · 0 评论 -
数论概论读书笔记 13.素数的计数
素数的计数素数有无穷多个且算术级数的素数 由狄利克雷定理知也有无数个现在的问题是对于一段区间,其间到底有多少个素数呢?或者说大致比例是多少呢?下表可以直观的感觉一下素数的分布定理13.1(素数定理) 当x很大时,小于x的素数的个数近似等于x/ln(x)x/ln(x)x/ln(x) 即 limx→∞π(x)x/ln(x)=1limx→∞π(x)x/ln(x)=1\l...原创 2018-08-26 09:35:14 · 433 阅读 · 0 评论 -
数论概论读书笔记 12.素数
素数定理12.1 (无穷多素数定理). 存在无穷多个素数证明 欧几里得证明法(见课本p55)定理12.2(模4余3的素数定理). 存在无穷多个模4余3的素数证明 变化的欧几里得证法 见书p57面(思路同上)但是这种思想,对于模4余1的素数不起作用!这并不是无聊的问题,理解论据限制与理解为什么论据有效几乎同样重要。定理12.3(算术级数的素数狄利克雷定理) 设a与...原创 2018-08-26 09:33:30 · 597 阅读 · 0 评论 -
数论概论读书笔记 11.欧拉函数与中国剩余定理
欧拉函数与中国剩余定理欧拉公式 aφ(m)≡1(modm)aφ(m)≡1(modm)a^{\varphi(m)} \equiv 1(mod m) 很美但是,我们使用这个公式时不可避免要计算φ(m)φ(m)\varphi(m) 一种简单的方法是枚举1~m-1来检查是否互质,但是当m几十位上百位时,这种方法就不可行了我们试着去寻找一些规律首先,当m是一个素数时,φ(p)=p...原创 2018-08-26 09:32:41 · 809 阅读 · 0 评论 -
POJ - 3048 区间数的最大质因子
To improve the organization of his farm, Farmer John labels each of his N (1 <= N <= 5,000) cows with a distinct serial number in the range 1..20,000. Unfortunately, he is unaware that the cows ...原创 2018-03-13 16:31:01 · 679 阅读 · 0 评论 -
POJ 2689(区间素数筛)
The branch of mathematics called number theory is about properties of numbers. One of the areas that has captured the interest of number theoreticians for thousands of years is the question of primali...原创 2018-03-12 21:18:02 · 1579 阅读 · 0 评论