数论
文章平均质量分 65
对与数论相关题目进行研究
酷酷的Herio
这个作者很懒,什么都没留下…
展开
-
数论选做-
数论选做P1286 两数之和考虑从小到大确定答案,枚举a1a_1a1。计算出a2a_2a2 ,删掉a1+a2a_1+a_2a1+a2。此时最小的就是a1+a3a_1+a_3a1+a3 ,然后通过a1a_1a1 得出a3a_3a3,再删掉与a3a_3a3 相加的数。时间复杂度:O(mxain2logn)O(mx_{a_i}n^2logn)O(mxain2logn)依次类推。参考代码如下#include<bits/stdc++.h>#define mp ma原创 2021-08-05 20:44:15 · 446 阅读 · 0 评论 -
拉格朗日插值法
拉格朗日插值法重心拉格朗日插值法gu拉插板子cal 离散点拉插O(n2)O(n^2)O(n2)inpo 连续点拉插O(n)O(n)O(n)//Lagrange Interpolation #define il inlineil ll ksm(ll a,ll n,ll m=mod){ll s=1;while(n){if(n&1) s=s*a%m;a=a*a%m;n>>=1;}return s;}struct LR{ ll x[N],y[N];int n; ll fa原创 2021-08-03 14:58:16 · 994 阅读 · 0 评论 -
P6156 简单题(推式子&莫反&整除分块)
P6156 简单题(推式子&莫反&整除分块)https://www.luogu.com.cn/problem/solution/P6156// Problem: P6156 简单题// Contest: Luogu// URL: https://www.luogu.com.cn/problem/P6156// Memory Limit: 250 MB// Time Limit: 1000 ms// Date: 2021-07-07 21:28:20// ------原创 2021-07-07 22:30:00 · 278 阅读 · 0 评论 -
P2522 [HAOI2011]Problem b(莫反&整除分块)
P2522 [HAOI2011]Problem b(莫反&整除分块)∑i=ab∑j=cd[gcd(i,j)=k]\large \sum\limits_{i=a}^b \sum\limits_{j=c}^d [gcd(i,j)=k]i=a∑bj=c∑d[gcd(i,j)=k]就是f(n,m,k)=∑i=1n∑j=cm[gcd(i,j)=k]f(n,m,k)=\large \sum\limits_{i=1}^n \sum\limits_{j=c}^m [gcd(i,j)=k]f(n,m,k)=i原创 2021-07-07 20:25:47 · 210 阅读 · 0 评论 -
P3327 [SDOI2015]约数个数和(莫反&整除分块)
P3327 [SDOI2015]约数个数和(莫反&整除分块)d(ij)=∑x∣i∑y∣j[gcd(x,y)=1]d(ij)=\sum\limits_{x|i}\sum\limits_{y|j}[gcd(x,y)=1]d(ij)=x∣i∑y∣j∑[gcd(x,y)=1]该式子证明部分如下:然后就是简单套路了。∑n∑m∑x∣i∑y∣j∑d∣x,d∣yμ(d)\large \sum\limits^n\sum\limits^m\sum\limits_{x|i}\sum\limits_{y原创 2021-07-07 20:00:34 · 251 阅读 · 2 评论 -
狄利克雷卷积总结
狄利克雷卷积总结积性函数常见的积性函数:φ,μ,σ,d\large \varphi,\mu,\sigma,dφ,μ,σ,d常见的完全积性函数:ϵ,I,id\large \epsilon,I,idϵ,I,id函数名数学表达式欧拉函数φ(n)\varphi(n)φ(n)莫比乌斯函数μ(n)\mu(n)μ(n)元函数ϵ(n)=[n=1]\epsilon(n)=[n=1]ϵ(n)=[n=1]恒等函数I(n)=1I(n)=1I(n)=1约数个数函数d(n)原创 2021-07-07 18:59:04 · 609 阅读 · 0 评论 -
欧拉函数的性质证明(补充)
欧拉函数的性质证明(补充)本文只证明φ(n)=n∏i=1k(1−1pi)\varphi(n)=n\prod\limits_{i=1}^k(1-\dfrac{1}{p_i})φ(n)=ni=1∏k(1−pi1)证明需要的东西:积性函数的性质算术基本定理φ(pk)=pk−pk−1=pk−1(p−1)\varphi(p^k)=p^k-p^{k-1}=p^{k-1}(p-1)φ(pk)=pk−pk−1=pk−1(p−1)关于第三个前置知识的证明:证明φ(n)=∏i=1kφ(piai)\原创 2021-07-07 11:00:02 · 249 阅读 · 0 评论 -
P2158 [SDOI2008]仪仗队(欧拉函数&莫反)
P2158 [SDOI2008]仪仗队(欧拉函数&莫反)考虑不同的斜率个数,则求∑i=1n∑i=1n[gcd(i,j)=1]\sum\limits_{i=1}^n\sum\limits_{i=1}^n [gcd(i,j)=1]i=1∑ni=1∑n[gcd(i,j)=1]转换一下:∑i=1n∑i=1n[gcd(i,j)=1]\sum\limits_{i=1}^n\sum\limits_{i=1}^n [gcd(i,j)=1]i=1∑ni=1∑n[gcd(i,j)=1]φ(i)=∑j=1i原创 2021-07-07 10:48:17 · 336 阅读 · 0 评论 -
2018CCPC-K - Mr. Panda and Kakin(欧拉定理&快速乘&CRT&EXGCD)
2018CCPC-K - Mr. Panda and Kakin(欧拉定理&快速乘&CRT&EXGCD)ab(modn)=ca^b\pmod{n}=cab(modn)=c已知n=pq,c,b=230+3n=pq,c,b=2^{30}+3n=pq,c,b=230+3利用RSA解密,找到n\sqrt{n}n附近的p,qp,qp,q。利用欧拉定理可得:ab(modn)=ab(modφ(n))(modn)=ab(mod(p−1)(q−1))(modn)a^b\pmod{n}=a^{原创 2021-07-05 20:34:15 · 494 阅读 · 2 评论 -
CF616E Sum of Remainders(分块&取模)
CF616E Sum of Remainders(分块&取模)余数求和加强了一下范围,注意要边算边取模。因为n,m≤1013n,m\le 10^{13}n,m≤1013,直接乘法会爆掉llllll,所以先对运算数取模,然后再对结果取模。不然会WA。// Problem: CF616E Sum of Remainders// Contest: Luogu// URL: https://www.luogu.com.cn/problem/CF616E// Memory Limit: 250 M原创 2021-07-02 11:23:18 · 374 阅读 · 0 评论 -
97. 约数之和(分治)
97. 约数之和(分治)求aba^bab的约数和 模ppp的值。σ(ab)(modp)\large \sigma(a^b)\pmod{p}σ(ab)(modp)首先利用算数基本定理:ab=(p1k1p2k2…pmkm)b=p1k1bp2k2b…pmkmb\large a^b=(p_1^{k_1}p_2^{k_2}\dots p_m^{k_m})^b=p_1^{k_1b}p_2^{k_2b}\dots p_m^{k_mb}ab=(p1k1p2k2…pmkm)b=p1k1bp2k2b原创 2021-07-02 10:41:06 · 229 阅读 · 0 评论 -
线性筛约数个数和(证明)
线性筛约数个数和(证明)之前已经证明过了线性筛约数和,现在证这个就灰常简单了。传送门非常类似,也是用一个数组来维护d(i)d(i)d(i)。由算术基本定理:n=p1k1p2k2…pmkmn=p_1^{k_1}p_2^{k_2}\dots p_{m}^{k_m}n=p1k1p2k2…pmkm这里g(n)=k1+1g(n)=k_1+1g(n)=k1+1,即nnn的最小素因子个数+1+1+1。时间复杂度:O(n)O(n)O(n)ll f[N],g[N];int cnt,vis[N]原创 2021-07-01 22:15:38 · 269 阅读 · 0 评论 -
P3708 koishi的数学题(因数和)
P3708 koishi的数学题题目传送门值得学习的点因子和σ(n)=∑d∣nd\sigma(n)=\sum\limits_{d|n}dσ(n)=d∣n∑dx(mody)=x−⌊xy⌋×yx\pmod y = x - \lfloor\dfrac{x}{y}\rfloor \times yx(mody)=x−⌊yx⌋×yf(x)=∑i=1nx(modi)f(x)=\sum\limits_{i=1}^n x\pmod{i}f(x)=i=1∑nx(modi),求f(1),f(2)…,f(n)f原创 2021-07-01 21:40:11 · 412 阅读 · 0 评论 -
LC.810. 黑板异或游戏(异或&数学)
LC.810. 黑板异或游戏(异或&数学)题意若当前玩家的局面为异或和为0的数组则胜。分情况1.开始所有数异或和为0,则先手赢。1.开始所有数异或和不为0. 考虑一个问题,不论先手还是后手,他们选择数时数组的奇偶性都是一样的。当nnn为偶数,先手每次都是从偶数个数选择一个数然后变成奇数个交给后手。假设nnn为偶数,先手输了。则必定存在一个局面先手无论选择那个数,剩下奇数个数异或和为0.即:a1⊕a2⊕a3…an=S≠0a_1\oplus a_2\oplus a_3\do原创 2021-05-22 10:18:47 · 288 阅读 · 2 评论 -
Gcd小练习2(LuoGu)
Gcd小练习2(LuoGu)P1414 又是毕业季II预处理每个因数出现的次数,因为随iii增加,答案是非递增的。所以维护一个mxmxmx即可。int mx=0;scanf("%d",&n);for(int i=1;i<=n;i++){ int x;scanf("%d",&x);mx=max(x,mx); for(int i=1;i*i<=x;i++) if(x%i==0){ a[i]++; if(i*i!=x) a[x/i]++; }原创 2021-05-16 11:08:47 · 260 阅读 · 1 评论 -
P1306 斐波那契公约数(数论&矩阵加速)
P1306 斐波那契公约数(数论&矩阵加速)思路一个重要性质:gcd(f[n],f[m])=f[gcd(n,m)]gcd(f[n],f[m])=f[gcd(n,m)]gcd(f[n],f[m])=f[gcd(n,m)],然后就矩阵快速幂加速即可。证明:不失一般性,不妨设n<m,f[n]=a,f[n+1]=bn<m,f[n]=a,f[n+1]=bn<m,f[n]=a,f[n+1]=b。f[n+2]=a+b,f[n+3]=2a+b…,f[m]=f[m−n−1]a+f[m−n]原创 2021-03-20 15:18:01 · 353 阅读 · 1 评论 -
D - The Number of Pairs(数论)
D - The Number of Pairs(数论)题意 给定c,d,xc,d,xc,d,x,求满足c×lcm(a,b)−d×gcd(a,b)=xc\times lcm(a,b)-d\times gcd(a,b)=xc×lcm(a,b)−d×gcd(a,b)=x的(a,b)(a,b)(a,b)对数。思路为了方便表达,令lcm(a,b)=l,gcd(a,b)=glcm(a,b)=l,gcd(a,b)=glcm(a,b)=l,gcd(a,b)=g。由lg=ablg=ablg=ab化简得:ab原创 2021-03-19 14:32:18 · 525 阅读 · 2 评论 -
P3802 小魔女帕琪(期望&概率论)
P3802 小魔女帕琪(期望&概率论)思路 依题意知,每个长度为777的区间贡献次数是独立的,即不受其他影响。所以我们只需考虑一个长度为777的区间的期望次数,最后乘上(N−6)(N-6)(N−6),因为有(N−6)(N-6)(N−6)个长为777的区间。对于一个长度为777的区间,总方案数是7!7!7!而释放魔法的概率是∏i=17aiN+1−i\prod\limits_{i=1}^7\dfrac{a_i}{N+1-i}i=1∏7N+1−iai,这里即使顺序不同,概率也是相同的,原创 2021-03-16 23:36:08 · 281 阅读 · 0 评论 -
P3937 Changing(Sierpinski三角形&规律)
P3937 Changing(Sierpinski三角形&规律)思路找规律或者生成函数(暂时不会)令ft,kf_{t,k}ft,k表示ttt时刻第kkk盏灯的情况。有递推式:ft,k=ft−1,k⊕ft−1,k+1f_{t,k}=f_{t-1,k}\oplus f_{t-1,k+1}ft,k=ft−1,k⊕ft−1,k+1=ft−1,k+ft−1,k+1=f_{t-1,k}+f_{t-1,k+1}=ft−1,k+ft−1,k+1=ft−2,k+2ft−2,k+1+ft−2,原创 2021-03-16 14:49:30 · 240 阅读 · 0 评论 -
C. Going Home(鸽巢原理&暴力)
C. Going Home(鸽巢原理&暴力)题意 给定nnn个数,求满足ax+ay=az+aw,(x,y,z,wa_x+a_y=a_z+a_w,(x,y,z,wax+ay=az+aw,(x,y,z,w不相同)的一组。思路鸽巢原理,因为ax+ay=az+aw∈[1,5×106]a_x+a_y=a_z+a_w\in[1,5\times 10^6]ax+ay=az+aw∈[1,5×106]而某一个和出现四次必有解,所以最多枚举2×1072\times 10^72×107次。原创 2021-03-14 19:39:33 · 495 阅读 · 2 评论 -
康托展开的学习总结
康托展开的学习总结介绍 康托展开是用于求解有关全排列的字典序排名问题的有效算法,康托展开就是全排列到排名的一个映射,而逆康托展开就是排名到全排列的一个映射,两者是双射关系,一一对应。公式 ans=(∑i=1nA[i]×(n−i)!)+1ans=(\sum\limits_{i=1}^n A[i]\times (n-i)!)+1ans=(i=1∑nA[i]×(n−i)!)+1A[i]A[i]A[i]是j∈[i,n]j\in[i,n]j∈[i,n]中满足a[j]<a[i]a[j]&l原创 2021-03-12 15:49:18 · 279 阅读 · 0 评论 -
平衡三进制学习总结
平衡三进制学习总结讲解OI-WIKI挺细的一个知识点,还挺简单的。简单来说就是一个用标准三进制表示的数进行一个转换。ep:64=(02101)3ep: 64=(02101)_3ep:64=(02101)3原本用(0,1,2)(0,1,2)(0,1,2)来表示的数,我们尝试用(−1,0,1)(-1,0,1)(−1,0,1)表示。为方便,记−1=Z-1=Z−1=Z所以有:64=(1Z101)364=(1Z101)_364=(1Z101)3转换方法:当该位为0,10,10,1满足直接跳过。原创 2021-03-02 17:51:45 · 2204 阅读 · 0 评论 -
孙子定理及扩展孙子定理
孙子定理及扩展孙子定理中国剩余定理 (Chinese Remainder Theorem)CRT是解决下列同余方程问题的,mi,mjm_i,m_jmi,mj两两互质。该定理只需要用到两个取余知识点:1.a%b=c→(a+kb)%b=ca\%b=c\rightarrow (a+kb)\%b=ca%b=c→(a+kb)%b=c2.a%b=c→(ka)%b=kc (kc<b)a\%b=c\rightarrow (ka)\%b=kc\ (kc<b)a%b=c→(ka)%b=kc原创 2021-02-28 22:00:45 · 438 阅读 · 0 评论 -
B.系数 (Lucas定理)
B.系数 (Lucas定理)题意给定n,k,p=3n,k,p=3n,k,p=3,求模ppp意义下f(x)=(x2+x+1)nf(x)=(x^2+x+1)^nf(x)=(x2+x+1)n的第kkk项系数:akxka_kx^kakxk的aka_kak。思路关键:模333意义下+x=−2x+x=-2x+x=−2x所以题意转化为:f(x)=(x2−2x+1)n=(x−1)2nf(x)=(x^2-2x+1)^n=(x-1)^{2n}f(x)=(x2−2x+1)n=(x−1)2n由二项式定理可知:第原创 2021-02-24 20:26:34 · 447 阅读 · 0 评论 -
灯泡开关问题
灯泡开关问题995. K 连续位的最小翻转次数每个位置作为翻转的起始位置最多一次,因为翻转两次等于没有翻转,且翻转与顺序没有关系,所以可以从左往右翻转,因为每次固定翻kkk个,可以利用差分实现区间加来判断该位置是否翻转过。319. 灯泡开关每个位置会被 因数给翻转 一次,所以考虑因数 的奇偶性即可,完全平方数是奇数,其他都是偶数,因为两两成对。672. 灯泡开关 Ⅱ实际上只取决于前三个灯,然后mmm最多枚举444,因为大于444都可以归为444。然后特判。1375. 灯泡开关 III原创 2021-02-18 12:44:50 · 454 阅读 · 1 评论 -
T162491 [Round 1]整除(同余方程&exgcd)
T162491 [Round 1]整除(同余方程&exgcd)solve:solve:solve:{x=10m+n (1)y=am+bn (2)\begin{cases}x=10m+n\ \ (1)\\y=am+bn\ \ (2)\end{cases}{x=10m+n (1)y=am+bn (2),考虑得到x,yx,yx,y的关系。(1)(1)(1)式同乘bbb,消掉bnbnbn。→y=bx+m(a−1原创 2021-02-17 15:27:03 · 1109 阅读 · 2 评论 -
F - GCD or MIN(思维&gcd)
F - GCD or MIN(思维&gcd)题意求nnn个数通过gcd(x,y)gcd(x,y)gcd(x,y)或者min(x,y)min(x,y)min(x,y)的操作每次合并一个数,最终得到不同的值个数。思路考虑最后的值是什么,显然必须小于等于min{ai}min\{a_i\}min{ai},因为最小的数被两个操作都不会使答案递增,显然最小的数是一种答案,因为一直取minminmin即可,接下来只需考虑答案是否可能通过gcdgcdgcd最后得到,可以注意到答案必须是某个aia_ia原创 2021-02-06 23:48:56 · 581 阅读 · 0 评论 -
Brute-force Algorithm(矩阵快速幂&欧拉降幂)
Brute-force Algorithm(矩阵快速幂&欧拉降幂)题意 a,b,ab,ab2…a,b,ab,ab^2\dotsa,b,ab,ab2… 序列的a,ba,ba,b指数是斐波那契形式,求在模ppp意义下的第nnn项。思路考虑先用矩阵快速幂预处理出a,ba,ba,b的指数,因为指数非常大,考虑欧拉降幂,只需要矩阵乘法中间加个特判即可。时间复杂度:O(p+Tlogn)O(p+Tlogn)O(p+Tlogn)代码#include<bits/stdc++.h>原创 2021-02-02 16:14:17 · 575 阅读 · 0 评论 -
P4139 上帝与集合的正确用法(欧拉降幂)
P4139 上帝与集合的正确用法(欧拉降幂)题意求2的无穷层幂塔函数在模ppp意义下的值。思路显然最多递归至logloglog层,所以用欧拉降幂即可。时间复杂度:O(Tlog2p)O(Tlog^2p)O(Tlog2p)#include<bits/stdc++.h>using namespace std;typedef long long ll;const int N=1e3+5,M=2e4+5,inf=0x3f3f3f3f,mod=1e9+7;#define mst(a,原创 2021-02-02 15:17:51 · 262 阅读 · 0 评论 -
牛客练习赛22 E.简单数据结构1(BIT&欧拉降幂)
牛客练习赛22 E.简单数据结构1(BIT&欧拉降幂)题意 给定长为nnn序列,有区间加和区间求幂塔值两个操作。思路 用权值BITBITBIT维护差分数组,然后欧拉降幂,这样区间操作就是单点更新lognlognlogn,区间求幂塔值就是log2nlog^2nlog2n。总复杂度:O(mlog2n)O(mlog^2n)O(mlog2n)需要注意的是,区间更新后的值,要在欧拉取模意义下成立。v=mod(v,m)v=mod(v,m)v=mod(v,m)。代码#include&l原创 2021-02-02 14:50:28 · 287 阅读 · 0 评论 -
幂塔个位数的计算(欧拉降幂)
幂塔个位数的计算(欧拉降幂)题意求给定幂塔函数的个位数思路个位数是一个循环,考虑用欧拉降幂,当nnn很大时,因为是%10\% 10%10,所以很快就会降到000,只需设置一个较小的nnn,而个数的结果,取决于aaa的十位和个位。因为个位数最大的循环是444。代码#include<bits/stdc++.h>using namespace std;typedef long long ll;const int N=1e6+5,M=2e4+5,inf=0x3f3f3f3f;#d原创 2021-02-02 12:41:52 · 901 阅读 · 0 评论 -
D. Power Tower (欧拉降幂&幂塔函数)
D. Power Tower (欧拉降幂&幂塔函数)题意 求给定区间的数的幂塔值。思路 扩展欧拉定理降幂的裸题,主要是扩欧+重新定义取模,欧拉函数可以n\sqrt nn+记忆化,主要是递归的思想,当递归到phi(x)=1phi(x)=1phi(x)=1直接返回,因为任何数模111为0,则(sum)%phi(x)+phi(x)=1(sum)\% phi(x)+phi(x)=1(sum)%phi(x)+phi(x)=1。因此时间复杂度:O(qlog2n)O(qlog^2n)O(qlo原创 2021-02-02 12:09:27 · 1015 阅读 · 0 评论 -
E. Border(数论&裴蜀定理)
E. Border(数论&裴蜀定理)题意 给定nnn个纸币对应的面值,求这些纸币任意组合(且每个纸币可以多个)在kkk进制下的个位的值有多少种。思路 考虑这些纸币的答案集合,由裴蜀定理可知,答案集合为:∑fiai=U(x,[gcd(ai)∣x])\sum f_i a_i=U(x,[gcd(a_i)| x])∑fiai=U(x,[gcd(ai)∣x])而这些值在kkk进制下的个位即x%kx\%kx%k,取模转化为减法即:gcd(a1,a2…,an,k)gcd(a_1,a_2\原创 2021-01-30 16:41:20 · 323 阅读 · 0 评论 -
C. Tile Painting(数论)
C. Tile Painting(数论)题意 给定n≤1012n\le 10^{12}n≤1012,表示nnn个格子,对于格子i,ji,ji,j如果n%∣i−j∣=0n\%|i-j|=0n%∣i−j∣=0且∣i−j∣>1|i-j|>1∣i−j∣>1,则i,ji,ji,j格子要涂同色,问最多能涂多少种颜色。思路 n≤1012n\le 10^{12}n≤1012,考虑n\sqrt{n}n解法,可以想到与素数有关。n=1,ans=1n=1,ans=1n=1,ans=1。如原创 2021-01-30 16:39:17 · 399 阅读 · 0 评论 -
D. Fox And Jumping(裴蜀定理&DP)
D. Fox And Jumping(裴蜀定理&DP)题意 给定nnn个数,每个数a[i]a[i]a[i]有对应的cost[i]cost[i]cost[i],求选出若干个数使得gcd=1gcd=1gcd=1的最小花费。思路 题意是由裴蜀定理转化得到,求最小花费使用dpdpdp,令dp[i]dp[i]dp[i]表示gcd=igcd=igcd=i的最小花费。有转移方程:对于第iii个数。dp[gcd(a[i],j)]=min(dp[gcd(a[i],j)],dp[j]+c[i])dp[原创 2021-01-29 16:42:38 · 530 阅读 · 1 评论 -
D. Nezzar and Board (裴蜀定理&gcd)
D. Nezzar and Board (裴蜀定理&gcd)题意 给定nnn个不同的数,每次操作可以选择两个数x,yx,yx,y (x,yx,yx,y可以相同)新增一个数2x−y2x-y2x−y。给定kkk,问是否能在若干次操作后得到kkk。思路 最后答案的形式肯定是:2x−y=x+(x−y)=k2x-y=x+(x-y)=k2x−y=x+(x−y)=k,即一个数加上(它与另一个数之差)的和组成,答案的形式是ai+∑j,k(aj−ak)a_i+\sum\limits_{j,k原创 2021-01-29 15:13:11 · 889 阅读 · 2 评论 -
5648. 生成乘积数组的方案数(质因数分解&组合数学)
5648. 生成乘积数组的方案数(质因数分解&组合数学)思路:质因数分解+组合数学。将kkk进行质因数分解后,考虑对于某一个质因数的个数cntcntcnt,如何分配给个nnn个位置,可以为空。这等价于把cntcntcnt个球放入nnn个篮子,篮子可以为空的方案。我们可以先将nnn个篮子都加入一个球,问题等价于cnt+ncnt+ncnt+n个球放入nnn个篮子,篮子不能为空的方案数,这可以用隔板法解决。ans=C(cnt+n−1,n−1)=C(cnt+n−1,cnt)ans=C(cnt+n-原创 2021-01-24 13:19:30 · 424 阅读 · 0 评论 -
牛客练习赛76 F.phi and phi(莫比乌斯反演&差分)
牛客练习赛76 F.phi and phi(莫比乌斯反演&差分)题解讲的很清楚了,预处理+差分+前缀和即可解决。#include<bits/stdc++.h>using namespace std;typedef long long ll;const int N=1e6+5,M=2e4+5,inf=0x3f3f3f3f,mod=1e9+7;#define mst(a,b) memset(a,b,sizeof a)#define PII pair<int,int>原创 2021-01-16 18:35:22 · 381 阅读 · 0 评论 -
有趣的无理数
有趣的无理数Conway 常数基本介绍传送门这里需要先引入外观序列(Look-and-say)序列Look-and-say数列是指以下特点的整数序列: 1, 11, 21, 1211, 111221, 312211, 13112221, 1113213211, …… 它以数字1开始,序列的第n项是对第n-1项的描述。比如第5项是111221,描述就是3个1,2个2,1个1, 所以下一项(第六项)就是312211。Conway常数的由来关于这个方程和Conway常数的由来,可以参考下面这篇原创 2021-01-14 22:29:33 · 1955 阅读 · 0 评论 -
F.消减整数(gcd)
F.消减整数(gcd)考虑:最后剩余的数为mmm,此时下一个要减的数为kkk,m<km<km<k,所以+n+n+n下一次会剩余2m,3m…2m,3m\dots2m,3m…知道xm≥kxm\ge kxm≥k,然后减去kkk,又小于kkk,所以最终目的就是求m,km,km,k的lcmlcmlcm。xm=k→x=lcm(m,k)m=kgcd(m,k)xm=k\rightarrow x=\dfrac{lcm(m,k)}{m}=\dfrac{k}{gcd(m,k)}xm=k→x=mlcm(m,原创 2021-01-11 13:54:07 · 347 阅读 · 2 评论