数论
文章平均质量分 75
616156
赛高______↑
展开
-
codejam-Round1A-2008-Numbers
题意求(3+√5)n次方的整数部分后三位分析乍一看觉得就是普通的一道快速幂的题,准备直接开搞,才发现因为√5无限不循环小数,这就使得整数部分即使在千位以上仍然会对百位以下造成影响(1000*0.1=1)。 想到这里,不难发现这是一道数学题。 观察这个式子(3+√5),不难想到一种常见套路——共轭 设a=3+√5a=3+√5,b=3−√5b=3-√5SolutionA原创 2017-07-03 10:23:04 · 185 阅读 · 0 评论 -
【数论】AGC003 Anticube
分析:有点套路的数学题。很显然,如果把每个数质因数分解,那么每个质因数的次数对3取模,归为一类,然后对每一类而言,考虑与其矛盾的类,取max即可。重点就在于如何质因数分解。很显然,如果要取模,则先保证质数的三次方在101010^{10}1010以内,那么质数的范围就大大缩小,到了[2,3∗103][2,3*10^3][2,3∗103]左右。这个范围内质数的个数就更少了,大约34034034...原创 2018-09-18 12:26:25 · 125 阅读 · 0 评论 -
【NTT】【数论】【图论】AGC005F Many Easy Problems
分析:这题最恶心的一点就在开头:首先,在树上,联通块大小为边的大小+1,所以可以算边的贡献:对每条边而言,如果它能造成贡献,那么必然在它两端都有选中的点,设选k个点,那么方案数就是C(n,k)−C(suma,k)−C(sumb,k)C(n,k)-C(sum_a,k)-C(sum_b,k)C(n,k)−C(suma,k)−C(sumb,k),suma,sumbsum_a,sum_bsuma...原创 2018-09-20 17:21:54 · 199 阅读 · 0 评论 -
【DP】2018国庆三校联考D1T1
题意:在N个格子之间,放入D-1个隔板(可以重合),要求每两个相邻隔板之间距离不超过M。求方案数。N,M≤2000N,M\leq 2000N,M≤2000D≤1012D\leq 10^{12}D≤1012分析:尽管D的范围非常大,但其实很多隔板之间的距离都为0,所以可以考虑距离不为0的隔板的放置方案(即不能重合的方案)。再利用组合数求出在所有隔板中选择一定数量的方案数。所以可以设DP...原创 2018-10-07 16:49:16 · 146 阅读 · 0 评论 -
【DP】2018国庆三校联考D2T2
题意:在一个二维平面上,求从(0,0)出发,经过T步到达(n,m)的方案数。模数P是数个质数之积分析:很简单的中国剩余定理+Lucas水题。直接枚举横向和纵向多走了多少步即可。#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#d原创 2018-10-07 17:20:04 · 187 阅读 · 0 评论 -
【数论】2018国庆三校联考D4T2
题意:很简单,给出正整数 n ,求出 n! 在十进制表示下的从最末非零位开始的总共 k 位。n≤1018,k≤3n\leq 10^{18},k\leq 3n≤1018,k≤3分析:有一个暴力的想法:直接从1枚举到N,然后把每个数的2的因数和5的因数除去,剩余部分乘起来mod 10kmod\ 10^kmod 10k然后我们剩下的,没有乘进去的2的整次幂肯定不比5的整次...原创 2018-10-07 19:32:33 · 246 阅读 · 0 评论 -
【数论】Codeforces1027G X-mouse in the Campus
分析:显然,i->ix(mod m)连一条边,则最终一定会形成若干个环,并且,环上每个点与m的gcd值必定相同。并且,gcd值相同的环大小也一定相同。所以,如果能算出对于所有数中,与m的gcd为d的个数f(d)f(d)f(d),并算出相应的当gcd为d时的每个环的大小l(d)l(d)l(d),那么答案就是∑f(d)l(d)\sum \frac {f(d)} {l(d)}∑l(d)f(d)...原创 2018-10-08 07:22:43 · 397 阅读 · 0 评论 -
【数论】【矩阵树定理】BZOJ3534重建
分析:首先,很显然我们要求的是:∑e∈TPe∗∑e′∉T(1−Pe′)\sum_{e\in T}P_e*\sum_{e'\notin T}(1-P_{e'})e∈T∑Pe∗e′∈/T∑(1−Pe′)前者可以用矩阵树定理算出来,问题就是转化后半部分。∑e′∉T(1−Pe′)=∑e(1−Pe)∑e∈T(1−Pe)\sum_{e'\...原创 2019-03-01 21:58:18 · 165 阅读 · 0 评论 -
【BSGS】【数论】TopCoder15277 WrongBase
题意:还是那个船长写的题意:分析:首先,已知的是:Yi≡gxiY_i\equiv g^{x_i}Yi≡gxi然后要求hxih^{x_i}hxi,如果算出T(h=gT)T(h=g^T)T(h=gT)那么hxi=gT×xi=YiTh^{x_i}=g^{T\times x_i}=Y_i^Thxi=gT×xi=YiT#include<cstdio>#includ...原创 2019-03-13 16:12:16 · 207 阅读 · 0 评论 -
【总结】伯努利数
分析:伯努利数满足:B0=1,∑i=0i≤nBiCn+1i=0B_0=1,\sum_{i=0}^{i\leq n} B_iC_{n+1}^i=0B0=1,i=0∑i≤nBiCn+1i=0根据这个性质,可以推导出一些结论:∑i=0i≤nBiCn+1i=0\sum_{i=0}^{i\leq n}B_iC_{n+1}^i=0i=0∑i≤nBiCn+1i=0∑i=0i<...原创 2019-03-27 17:07:49 · 694 阅读 · 0 评论 -
【数论】【容斥原理】【EXGCD】COCI 2017/2018 Round #3 Sažetak
分析:非常板的容斥题。。。考场上时间多点应该还是写得出来的。。。转换一下题目,就是求 满足x≡1(mod&nbsp;ai)且x≡0(mod&nbsp;aj)x≡1(mod&nbsp;ai)且x≡0(mod&nbsp;aj)x\equiv 1(mod\ a_i)且x\equiv 0(mod\ a_j)的x的个数(x≤Nx≤Nx\leq N)。由于N非常大,无法判断求解,只能算贡献。...原创 2018-09-06 18:19:44 · 472 阅读 · 0 评论 -
【总结】【数论】原根和指标
前言:原根和指标在数论中有大量的性质和规律,这里无法一一列举,只能简要写写用到的概率大一些的内容。本文是博主在复习时总结用,许多证明都已经略去,不建议初学者参考(速食主义者除外)。 主要包括:原根的定义及基本性质,求法,指标的定义,以及指标的基本应用。原根定义及性质:要说明原根的完整定义,又要扯到阶那一块,这里就直接用一句话大致表示了:若aaa为模m的原根 那么必然满足(a...原创 2018-09-03 14:05:03 · 2591 阅读 · 0 评论 -
【离散对数BSGS】总结
BSGS算法简述对于一个方程: Ax≡B(mod&nbsp;C)Ax≡B(mod&nbsp;C)A^x\equiv B(mod\space C) 求能满足条件的最小的x(GCD(A,C)=1)x(GCD(A,C)=1)x(GCD(A,C)=1) 由于GCD(A,C)=1GCD(A,C)=1GCD(A,C)=1,所以一定满足Aφ(C)≡1(mod&nbsp;C)Aφ(C)≡1(mod&nb...原创 2018-03-23 21:22:34 · 509 阅读 · 0 评论 -
【数论】Lucas定理与拓展Lucas定理
Lucas定理Lucas定理是用于解决快速计算组合数取模的,即计算 Cmn(mod k)Cnm(mod k)C_n^m (mod \space k) 当然,在n与m均不大时(108左右)(108左右)(10^8左右),可以使用O(n)求值的方法,但是当n与m继续扩大,那么传统方法就不适用了。在p的范围较小时且为质数,就可以使用Lucas定理了 Cmn=Cm/pn/...原创 2018-03-31 09:58:21 · 465 阅读 · 0 评论 -
【容斥原理】(AtCoder Regular Contest 093 F) Dark Horse
题意:有2N2N2^N个选手参与一场比赛,比赛规则是:相邻的两个人比赛一次,败者淘汰掉,胜者继续进行,直到只剩一个人为止。 现在给出1号选手会败给哪些选手(实力摸得很清楚啊) 并且已知其他选手之间均满足:两个选手比赛,编号小的一定会胜利。 现在可以安排每个选手初始的位置,要 钦定 1号选手 Chicken Dinner 最后获胜,求能满足条件的初始位置的方案数。分析:首先...原创 2018-03-28 21:20:43 · 366 阅读 · 0 评论 -
Codeforces483Div1 983A Finite or not?
题意:给出N次询问,每次询问十进制下的p/q能否在b进制下表示成有限小数。 N≤105,p,q,b≤1018N≤105,p,q,b≤1018N\leq 10^5,p,q,b\leq 10^{18}分析:在b进制下,该数若能被表示为有限小数,则必须满足: p∗bk≡0(mod q)p∗bk≡0(mod q)p*b^k\equiv 0 (mod \space q...原创 2018-05-26 13:29:53 · 296 阅读 · 0 评论 -
【总结】概率与期望
前言作为NOIP级的知识点,概率与期望算是比较困难的类型了。 但其实也不是无法解决的难题。本文主要通过作者本人的刷题经历,对概率期望类题目进行总结。概率51Nod1639绑鞋带:有n根鞋带混在一起,每根鞋带有两个鞋带头。现在重复n次以下操作:随机抽出两个鞋带头,把它们绑在一起。求最终只形成一个环的概率?依次考虑每一步操作,现在已经选出来了一个头,它必须和非它所在的链的另一...原创 2018-05-31 15:54:31 · 2506 阅读 · 0 评论 -
【组合计数】ARC061F Card Game for Three
题意:有三堆卡牌,牌数分别为N,M,K 每张牌有一个字母(’a’、’b’、’c’)表示下一个拿哪一堆。 现在要求第一堆首先拿完。求方案数。分析:首先,这道题有很多角度可以入手。但最简单的方法是,根据拿的牌的类型计算。 首先,第一堆拿完,则拿的顺序中必然有n个a,且b的数量比m小,c的数量比k小。 但是最后一个必须限定为a,所以这部分拿了的方案数应为Cn−1n−1+iCn...原创 2018-07-19 21:00:55 · 790 阅读 · 1 评论 -
【组合计数】CF1007B Pave the Parallelepiped
题意: 话说出题人在搞笑么。。。全程4个公告。。。T次询问,每次询问给出一个A、B、C,表示一个A×B×CA×B×CA\times B\times C的长方体 求有多少种方案的a×b×ca×b×ca\times b\times c能够凑出这个大的长方体(1≤a≤b≤c)(1≤a≤b≤c)(1\leq a\leq b \leq c) 要求小的立方体必须以同样的方向来凑。...原创 2018-07-14 14:10:35 · 989 阅读 · 2 评论 -
【数论】HDU6390 GuGuFishtion
分析://详情参见problem B 很明显,如果(x,y)=1(x,y)=1(x,y)=1(“()”表示最大公因数) 那么φ(x∗y)=φ(x)∗φ(y)φ(x∗y)=φ(x)∗φ(y)\varphi(x*y)=\varphi(x)*\varphi(y)那么如果(x,y)>1(x,y)>1(x,y)>1 将x、yx、yx、y唯一分解: x=pkx1x1×pkx2x2×...原创 2018-08-14 10:39:00 · 105 阅读 · 0 评论 -
【数论】Codeforces1025G Company Acquisitions
题意:有n个节点,每个节点有两种状态:选中和未选中。每个选中的点后面都跟着若干个(可能是0个)未选中的点。每个未选中的点都一定跟在某个选中的点后面。每次操作随机选择两个被选中的点,随机将其中一个变成未选中,且跟在另一个后面,同时将跟在他后面的节点全部改为选中。求这样操作下去,直到最后只剩一个选中的点的期望步数。分析:非常有趣的一道数论题: 设一个点后面跟着k个被选...原创 2018-08-23 17:03:58 · 473 阅读 · 0 评论 -
【莫比乌斯反演】HDU6428 Problem C. Calculate
题意:求∑i=1i≤A∑j=1j≤B∑k=1k≤Cφ(gcd(i,j2,k3))∑i=1i≤A∑j=1j≤B∑k=1k≤Cφ(gcd(i,j2,k3))\sum_{i=1}^{i\leq A}\sum_{j=1}^{j\leq B}\sum_{k=1}^{k\leq C}\varphi(gcd(i,j^2,k^3))分析:膜拜cch。。。 狄利克雷卷积:(f∗g)(D)=∑d|D...原创 2018-08-23 22:24:59 · 347 阅读 · 0 评论 -
【数论】【图论】【Exgcd】【LCA】Codeforces500G
题意:给出一棵树,有m组询问,每组询问包含2对点(u,v)和(x,y)。现在一个人从u出发,在u和v之间的路径上来回移动,一个人从x出发,在x和y上的路径来回移动,求什么时候两人能刚好出现在一个点上。如果不能输出-1分析非常棒的一道题啊。。。首先,要把这个在图上的问题转化成数论问题:对于每次询问,我们找到以下信息: 两条路径各自的长度,相交部分的长度,相交部分的两个端...原创 2018-09-07 20:29:07 · 250 阅读 · 0 评论 -
【伯努利数】【MTT】51nod 1258序列求和
分析:Sk(n)=1k+1∑i=0i≤kCk+1iBink+1−iS_k(n)=\frac {1} {k+1}\sum_{i=0}^{i\leq k}C_{k+1}^{i}B_in^{k+1-i}Sk(n)=k+11i=0∑i≤kCk+1iBink+1−i#include<cstdio>#include<algorithm>#include<vect...原创 2019-03-27 17:40:50 · 235 阅读 · 0 评论