![](https://img-blog.csdnimg.cn/20201014180756930.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
LZRcqbz
这个作者很懒,什么都没留下…
展开
-
【生成函数】五边形数定理与整数划分问题详解
Part 0. 前置知识简单的公式推导生成函数大量数学公式警告Part 1. 什么是五边形数五边形数看一张图吧:XP 系统带的画图真的挺好用!我们不难得出五边形数的数列:P={1,5,12,22,…}P = \{1, 5, 12, 22,\ldots\}P={1,5,12,22,…}再找一下规律,我们就可以得到五边形数的通项公式:Pn=n(3n−1)2P_n = \...原创 2019-12-07 17:50:33 · 2953 阅读 · 2 评论 -
【数论】中国剩余定理与扩展中国剩余定理详解
Part. 0 前置知识扩展欧几里得算法;模运算相关知识。Part. 1 中国剩余定理中国剩余定理是用于求解形如:{x≡a1(mod m1)x≡a2(mod m2)…x≡an(mod mn)\begin{cases}x\equiv a_1(\mod m_1)\\x\equiv a_2(\mod m_2)\\\ldots\\x\equiv a_n(\mod m_n)\...原创 2019-11-14 16:04:39 · 424 阅读 · 0 评论 -
乘法逆元详解【费马小定理+扩展欧几里得算法】
乘法逆元何为乘法逆元?对于两个数a,pa,pa,p若gcd(a,p)=1gcd(a,p)=1\gcd(a,p)=1则一定存在另一个数bbb,使得ab≡1(modp)ab≡1(modp)ab\equiv1(\mod p),并称此时的bbb为aaa关于111模ppp的乘法逆元。我们记此时的bbb为inv(a)inv(a)inv(a)或a−1a−1a^{-1}。举个例子:5×3≡1(mod1...原创 2018-08-05 19:58:41 · 11189 阅读 · 6 评论 -
【AtCoder】【数学】1974 いろはちゃんとマス目 / Iroha and a Grid
AtCoder 1974 いろはちゃんとマス目 / Iroha and a Grid题目 ◇题目传送门◆题目大意 给定H,W,A,BH,W,A,BH,W,A,B,H,WH,WH,W分别表示矩形的长宽,A,BA,BA,B表示位于左下角的A×BA×BA\times B的矩形不能通过,求从左上角到右下角的方案数,答案对109+7109+710^9+7取模。思路...原创 2018-08-06 20:21:51 · 338 阅读 · 0 评论 -
【CodeForces】【数学】616E-Sum of Remainders
CodeForces 616E Sum of Remainders题目 ◇题目传送门◆题目大意 给定n,mn,mn,m求∑mi=1(nmodi)∑i=1m(nmodi)\sum_{i=1}^{m}(n\mod i)。思路暴力???1013101310^{13},早就炸了。。。我们就用数学方法分析一下:Step1.变模为求和众所周知,nmodinmod...原创 2018-08-02 21:32:53 · 278 阅读 · 0 评论 -
【数论】Pollard-Rho 算法总结
Pollard-Rho 算法总结预备知识质因数分解;生日悖论。然后就愉快的开始吧。Part.I 如何对一个大合数进行质因数分解设有一个大合数N=p×q,(p≠q)N=p\times q,(p\not= q)N=p×q,(p=q),那么我们可以用什么方法去把它分成两个数的乘积呢?显然我们可以只找到ppp或者qqq,另一个用NNN去除一下就好了。那么如何去找这个数?显然可以用...原创 2019-09-17 22:35:04 · 542 阅读 · 0 评论 -
大步小步算法(BSGS)详解
BSGS(Baby Step Gaint Step)算法详解BSGS算法是用于求解形如Ax≡B(mod P)A^x\equiv B(\mod P)Ax≡B(modP)的高次同余方程,其中A,B,PA,B,PA,B,P为已知,且要求gcd(A,P)=1\gcd(A,P)=1gcd(A,P)=1是一个质数。前置知识哈希算法分析先证明一个结论:若原方程有解,则x∈[0,P−1)x\in...原创 2019-09-23 22:27:51 · 1539 阅读 · 1 评论 -
扩展大步小步算法详解
扩展BSGS算法详解前置知识BSGS(可以参见我的博客BSGS算法详解)辗转相除法算法分析扩展BSGS主要用于求解高次同余式Ax≡B(mod P)A^x\equiv B(\mod P)Ax≡B(modP),其中gcd(A,P)\gcd(A,P)gcd(A,P)为任意的数。我们先把这个同余式改写成A⋅Ax−1+kP=B,k∈ZA\cdot A^{x-1}+kP = B,k\in ...原创 2019-09-24 21:31:48 · 261 阅读 · 0 评论