![](https://img-blog.csdnimg.cn/20201014180756780.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学
文章平均质量分 90
STcyclone
这个作者很懒,什么都没留下…
展开
-
【数论算法】扩展欧几里得算法详解+编程求解不定方程
扩展欧几里得算法详解扩展欧几里得算法详解及模板题(洛谷P5656)解题思路题目链接题目描述数据范围扩展欧几里得算法代码模板时间复杂度算法正确性的证明归纳基础归纳步算法时间复杂度的证明b=0b=0b=0 时b≠0b \neq 0b=0 且 a<ba < ba<b 时b≠0b \neq 0b=0 且 a≥ba \ge ba≥b 时引理: ∀a,b∈Z+\forall a,b \in \mathbb{Z}^{+}∀a,b∈Z+,若 a≥ba \ge ba≥b,则有 a mod b≤⌊a2原创 2022-03-15 16:09:32 · 1522 阅读 · 0 评论 -
洛谷p3768 简单的数学题
题目大意求∑i=1n∑j=1nijgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{n}ijgcd(i,j)对质数p取膜的结果n<=10^10 5*10^8<=p<=1.1*10^9分析看到这种两层求和带gcd的式子 可以考虑直接用莫比乌斯反演求解。 但注意到i j的上界都是n,所以用ϕ\phi推导公式也可以解决这题。(感觉可能会简便一些)∑i=1n∑j=1nijgcd(i,原创 2017-06-23 20:06:27 · 1435 阅读 · 0 评论 -
乘法逆元的几种求法总结
乘法逆元对于缩系中的元素,每个数a均有唯一的与之对应的乘法逆元x,使得ax≡1(mod n) 一个数有逆元的充分必要条件是gcd(a,n)=1,此时逆元唯一存在 逆元的含义:模n意义下,1个数a如果有逆元x,那么除以a相当于乘以x。下面给出求逆元的几种方法1 循环找解法给定模m和需要求逆的数x,直接暴力枚举1~m-1 检查是否有x*i=1(mod m)这种算法可以应用与写暴力、对拍、模数较小,原创 2016-07-31 23:58:21 · 18762 阅读 · 0 评论 -
狄利克雷卷积及莫比乌斯反演定理
在数论函数集上,狄利克雷卷积定义为如下的二元运算 (f*g)(n)=d|nΣf(d)g(n/d)狄利克雷卷积运算满足交换律,结合律,分配律在数论函数集上,存在单位元函数ε(n)满足ε(1)=1 ε(n)=0(n>1) 使得(f*ε)=(ε*f)=f而常值函数1(n)则定义为1(n)=1一个数论函数f(n)与1(n)的狄利克雷卷积并不是f(n)自身,而是一个常见的式子 (f*1)(n)=d|n原创 2016-07-31 23:29:57 · 2769 阅读 · 1 评论 -
线性筛法求素数表
scanf("%d",&n); fo(i,2,n){ if (!isp[i]){ pri[++pt]=i; //printf("%d ",i); } fo(j,1,pt){ t=(long long)i*pri[j]; if (t>n)原创 2016-07-31 22:59:45 · 654 阅读 · 0 评论 -
线性筛法求phi函数表
void phi_table(){ phi[1]=1; fo(i,2,mn){//从2开始计算 if (!phi[i]){//等价于是质数 pri[++pt]=i; phi[i]=i-1; } fo(j,1,pt){ //枚举添加的小质数 te=i*pri原创 2016-07-31 22:26:55 · 2994 阅读 · 0 评论 -
VIJOS 1516 N连环
题意九连环应该都知道吧! (1)第 1 环可以自由上下 (2)而上/下第 n 环时(n>1),则必须满足: (a)第 n-1 个环在架上 (b)前 n-2 个环全部在架下 为了让大家多学点知识,特改此题输入2个N连环的状态 保证环数小于等于100 第一行是初状态 第二行是末状态 输出一个数,需要的步数n(保证n不超过qword)分析九连环的数学模型实际上是格雷码我们用一串01串表示N原创 2016-06-27 22:55:14 · 664 阅读 · 0 评论 -
FFT/NTT做题方法与调试技巧(+提高码题效率的一些想法)
FFT/NTT是卷积运算常见而实用的优化但是FFT/NTT的处理过程并不像暴力运算(差不多是多项式乘法)那样能够直观地反映卷积结果的实时变化。因此在使用FFT时将会或多或少地加大调试的难度。如果调试程序时直接跟踪变量,每步手算结果比对,不仅会耽误大量时间,而且效果可能并不理想。直接肉眼查错效率可能也不太高。但也正由于FFT的特点,实际上也有一些很方便而实用的调试方法,能够明显提升调试的效率,减小调试原创 2016-05-07 18:59:28 · 2990 阅读 · 3 评论 -
VIJOS 1999 小h的妹子树三
题意简述给定节点个数n和树高h(树根高为0),求出满足节点数位n,树高为h的二叉树的个数 答案对p取模,T组数据 对于100%的数据,满足1<=n<=1200,1<=h<=1200,1<=T<=2000000,1<=P<=1007。分析这题的数据范围提示我们要将所有询问预处理出来 而二叉树计数问题与卷积有着比较密切的联系,模数较小也提示我们这题正解很有可能就是FFT设有n个节点,高度为h的二叉原创 2016-05-18 00:04:52 · 3791 阅读 · 0 评论 -
CF 615E Hexagons
题目大意给定一个由六边形密铺组成的坐标系中的点的坐标规律。求从原点出发以六边形轨迹行走n步所走到的坐标(0 ≤ n ≤ 10^18)题目链接: http://codeforces.com/contest/615/problem/E分析很明显的找规律题了,我是直接求边长后分情况讨论在哪一条边上,代码写的很丑很长。(同时应该比较容易理解吧)思考:如果我们已知点的坐标,该如何求得走了多少步呢。这或原创 2016-01-14 13:45:08 · 933 阅读 · 0 评论 -
CF615D Multipliers
题目大意一个数n由m个质因子pi相乘而得(n=p1p2…pm),求出n的所有因子的乘积模1e9+7的值 1 ≤ m ≤ 200 000 2 ≤ pi ≤ 200 000分析代码原创 2016-01-11 22:38:55 · 683 阅读 · 0 评论 -
VIJOS1420Valentine’s Seat
题目大意可以给排列成n+1行m+1列矩阵的椅子上色。其中,第一行所有的椅子都已上好粉红色,第一列从第2个开始的椅子都已上好天蓝色。对于其他的椅子,应保证,它的颜色和它左边的或上面的椅子颜色相同。求染色的方案数模19900801的值。分析通过画图或者打表,我们可以发现,最终矩阵肯定被分成2个色块,粉红色块位于矩阵右上方,天蓝色块位于矩阵左下方。每一种染色方案都一一对应于一条从n+1行m+1列的格点地图原创 2016-01-11 22:24:28 · 452 阅读 · 0 评论 -
平方和立方和公式的推导及其拓展
描述原创 2015-10-03 12:39:34 · 3766 阅读 · 2 评论 -
VIJOS1413Valentine’s Present
题目描述今天是情人节,小杉已经想好了要给喜欢的人送一份特殊的情人节礼物。 礼物是n个颜色各异的箱子,每个箱子里装一个蛋糕,蛋糕是可以上色的。 一个可爱的上色方案应该满足如下条件: 1. 任意一个蛋糕上的颜色应与一个箱子相同(可以是装它的那个箱子的颜色)。 2. 任意开启一个箱子,按里面蛋糕的颜色打开对应的箱子,这两个箱子(也可以是同一个)里的蛋糕颜色相同。 小杉现在想知道总共有多少种可爱的原创 2015-12-20 13:51:29 · 769 阅读 · 0 评论