数论
beginendzrq
主要负责给各位大佬端茶送水...
展开
-
PE 432
果然还是自己太弱了,看的别人的题解QAQ 题目大意:定义 S(n,m)=∑mi=1φ(ni)S(n,m) = \sum_{i=1}^m\varphi(ni),求S(510510,1011)S(510510,10^{11}) 题解: 1)μ(n)!=0\mu(n)!=0 有φ(nk)=∑d|(n,k)φ(nd)φ(k)\varphi(nk) = \sum_{d|(n,k)}\varphi(\f原创 2016-03-30 20:34:20 · 343 阅读 · 0 评论 -
PE 216 【二次同余】
作为一个AFO的选手我又回到了机房( ̄口 ̄) 能继续学OI当然很开心,但是看到机房里的人少了一半,连含爷也走了,外带坑了我妈两万软妹币,结果还是只能狗牌退出,还是好难过。。。各种纠结啊有木有 算了,做道数学题冷静一下(说好的数据结构呢...———–正题—————–题目大意:设 t(n)=2n2−1t(n) = 2n^2 - 1,求t(n)中素数的个数2<=n<=5∗1072 <= n <= 5*原创 2016-04-16 21:29:51 · 604 阅读 · 0 评论 -
跪烂数竞大爷orz
证明一个命题: if(a,b)=1,then ∀x∈N and x>ab−a−b,都能表示成x=pa+qb的形式(p,q∈N) 说真的感觉在哪里见过这个结论但是又想不起来了…… 我们来简单证明一下: 贝祖定理: ∵∃x,y∈Z,ax+by=1 ∴ ∀n∈Z,∃p=nx,q=ny s.t. pa+qb=n 注意到,p可以用减掉一个b的代价使得q加上a。 注意到,p可以转载 2016-08-15 22:58:37 · 618 阅读 · 1 评论 -
PE 201
题目大意:给定一个集合 S={12,22,⋯,1002}S = \{1^2,2^2,\cdots,100^2\},从 S 中选择 50 个数求和,设得到的和的集合为 A,求出 A 中哪些数只出现了一次,并求出这些数的和。。。 这道题窝不会啊,,,只会暴力dp。。。 如果泥萌谁会了记得教我2333(滑稽【答案】115039000原创 2016-08-16 00:36:05 · 478 阅读 · 0 评论 -
PE 211
题目大意:设σ2(n)=∑d|nd2\sigma_2(n) = \sum_{d|n}d^2,求所有σ2(n)\sigma_2(n)为完全平方数的n的和。。。大半夜不睡觉愉快地划水。。。【答案】1922364685#include <iostream>#include <cmath>#include <cstdio>#define LL long longusing namespace std原创 2016-08-19 00:55:53 · 1268 阅读 · 0 评论 -
PE 157
哎,开学前最后划一下水,下次泥萌见到窝就是高三狗了TATATATAT题目大意:求有多少对a<=b满足1a+1b=p10n\frac 1a + \frac 1b = \frac p{10^n},1 <= n <= 9PE的难度果真有毒……设 a = gx,b = gy,g = gcd(a,b) 有 x+ygxy=p10n\frac {x + y}{gxy} = \frac p{10 ^ n}原创 2016-08-25 13:56:03 · 338 阅读 · 0 评论 -
一些证明
费马小定理证明:窝们构造出 p-1 个数:a,2a,3a,…,(p-1)a先证明这 p-1 个数两两不同余:若存在两个数同余则有 (i−j)a≡0 (mod p)(i - j)a \equiv 0 \ (mod \ p) 因为 gcd(a,p)=1\gcd (a,p) = 1,所以有 i−j≡0 (mod p)i - j \equiv 0 \ (mod \ p) 显然不成立所以这 p-1 个数构原创 2016-08-26 13:53:39 · 616 阅读 · 0 评论 -
PE 202 【莫比乌斯反演】
题目大意:给一个正三角形的镜子,光线可以从一个顶角射入和射出,给出反射次数,求光路的方案数先挖坑。。。去背阿房宫赋TAT【答案】1209002624原创 2016-09-16 22:58:58 · 461 阅读 · 0 评论 -
分圆多项式
这是一个坑,而且估计这辈子都不会来填。。。定义 xn−1=0x^n - 1 = 0 的解为n次单位根 学过法法塔知道它在复数域上存在n个根(k∈[0,n)k \in [0,n)) wkn=e2πkn=cos2πkn+isin2πknw_n^k = e^{\frac {2\pi k}n} = \cos \frac {2\pi k}n + i\sin \frac {2\pi k}n 而这n个解对应原创 2016-12-01 17:13:55 · 6100 阅读 · 0 评论 -
HDU 5393 【数论】
最近总被虐…被虐…虐…%%%fsf大爷,太强辣 欢迎神犇打脸(求轻喷题目大意:给出 x0x_0,每次操作将 x 变成 (kx + b) mod p,求最少多少次操作能回到 x0x_0…首先 k = 0,1 要特判 …n 次操作后的值就是 ∑n−1i=0kib+knx\sum_{i=0}^{n-1}k^ib + k^nx 所以就能得到∑i=0n−1kib+knx≡x(modp)kn−1k−1(b+原创 2016-06-16 15:33:13 · 793 阅读 · 0 评论 -
spoj GCDMAT
题目大意:给定 a,b,c,da,b,c,d,求 ∑bi=a∑dj=cgcd(i,j)\sum_{i=a}^b\sum_{j=c}^dgcd(i,j) ∑i=1n∑j=1mgcd(i,j)=∑d=1nd∑i=1⌊nd⌋∑j=1⌊md⌋[gcd(i,j)==1]=∑d=1nd∑i=1⌊nd⌋∑j=1⌊md⌋∑g|i,jμ(g)=∑d=1nd∑g=1⌊nd⌋μ(g)⌊ngd⌋⌊mgd⌋\begin{原创 2016-06-10 21:44:44 · 542 阅读 · 0 评论 -
PE 372 【类欧几里德】【Stern-Brocot Tree】
题目大意:求 ⌊y2x2⌋\lfloor\frac {y^2}{x^2}\rfloor 为奇数的个数,其中 2∗106<x,y<=1092*10^6<x,y <= 10^9设 k=⌊y2x2⌋k = \lfloor\frac {y^2}{x^2}\rfloor,有 x2k<=y2<x2(k+1)x^2k <= y^2 <x^2(k+1) 得 xk√<=y<xk+1−−−−√x\sqrt k <=原创 2016-05-13 19:11:40 · 531 阅读 · 1 评论 -
PE 441
#include #include#include#include#include#include#include#include#include#define D doubleusing namespace std;const int n = 10000000,N = n + 20;D f[N],g[N],t[N];int miu[N],prime[N];bool原创 2016-03-23 12:52:52 · 385 阅读 · 0 评论 -
PE 530
orz rzz 题目大意:设 f(n)=∑d|ngcd(d,nd)f(n) = \sum_{d|n}gcd(d,\frac{n}{d}),F(n)=∑ni=1f(i)F(n) = \sum_{i = 1}^{n}f(i),求F(1015)F(10^{15})的值 题解:F(n)=∑ni=1∑d|igcd(d,id)=∑d∑nd2i=1∑d′|i[gcd(d′,id′)==1]=∑d∑d′dμ(d原创 2016-03-26 22:19:15 · 681 阅读 · 0 评论 -
PE 512
orz rzz 题目大意:设f(n)=∑ni=1ϕ(ni)mod(n+1)f(n) = \sum_{i=1}^n\phi(n^i) mod (n+1),g(n)=∑ni=1f(i)g(n) = \sum_{i=1}^nf(i),求g(5∗108)g(5*10^8) 题解:f(n)=ϕ(n)∑ni=1(−1)i−1f(n) = \phi(n)\sum_{i=1}^n(-1)^{i-1} 当 n原创 2016-03-26 22:35:20 · 427 阅读 · 0 评论 -
PE 443
题目大意:给定g(4)=13g(4) = 13,g(n)=g(n−1)+gcd(n,g(n−1))(n>=4)g(n) = g(n-1) + gcd(n, g(n-1))(n >= 4),求g(1015)g(10^{15})打表发现大多数gcd的值都为1,考虑去找下一个gcd不为1的值 设当前位置f(n)=gf(n) = g 即求gcd(n+1+d,g+d)!=1gcd(n + 1 + d,g原创 2016-04-25 11:52:11 · 588 阅读 · 1 评论 -
PE 521【DP】
题目大意:设ms(i)ms(i)代表i的最小质因数,求∑ni=2ms(i)\sum_{i=2}^nms(i),n=1012mod109n = 10^{12}\;mod\;10^9对于pi>n−√pi>\sqrt n,pi只会出现一次,所以可以先处理<=n−√<=\sqrt n的素数 设 F[i,j] 表示不含p1..ip_{1..i}的因数,范围在[1,j]的数的个数 设 G[i,j] 表示不含原创 2016-04-25 18:50:41 · 715 阅读 · 1 评论 -
PE 439 【莫比乌斯反演】【杜教筛】
垃圾题,坑了我一天QAQ题目大意:求∑Ni=1∑Nj=1σ(ij)\sum_{i=1}^N\sum_{j=1}^N\sigma(ij),N=1011N = 10^{11}∑Ni=1∑Nj=1σ(ij)\sum_{i=1}^N\sum_{j=1}^N\sigma(ij) =∑Ni=1∑Nj=1∑x|i∑y|jxjy[gcd(x,y==1)]\sum_{i=1}^N\sum_{j=1}^N\sum_{原创 2016-04-18 14:45:51 · 2405 阅读 · 2 评论 -
PE 75
题目描述:对于给定周长n,有且只有一种方案使得所围成的直角三角形各边长为整数,求有多少个这样的n(1<=n<=1.5∗1061.5*10^6)和圆上的整点那道题差不多……可以先求出各边长互质的三角形 a2+b2=c2a^2 + b^2 = c^2 a2=(c+b)(c−b)a^2 = (c + b)(c - b) 则有c+bc + b,c−bc - b为完全平方数或完全平方数的两倍 然后枚举原创 2017-06-10 00:00:22 · 455 阅读 · 0 评论