数论
QAQwdd
这个作者很懒,什么都没留下…
展开
-
LuoguP3807 【模板】卢卡斯定理 Lucas+费马小定理求逆元
Lucas定理: Cmn=C⌊m%p⌋⌊n%p⌋∗Cm/pn/pCnm=C⌊n%p⌋⌊m%p⌋∗Cn/pm/pC_{n}^{m} = C_{\lfloor n\%p \rfloor}^{\lfloor m\%p \rfloor} *C_{n/p}^{m/p} 具体计算组合数时,由于涉及到除法取模,所以可以用费马小定理求逆元。 由于gcd(a,p)=1gcd(a,p)=1gcd(a, p) =...原创 2018-08-05 15:10:49 · 213 阅读 · 0 评论 -
NOIP2014 解方程 数论+模拟
题目大意:求一个多项式方程在[1,m][1,m][1, m]的整数解。 题解:没想到吧,这题竟然是暴力(溜 好吧也没有那么简单 要用到一个算法:秦九韶算法,就是减少多项式的计算次数 然后暴力枚举[1,m][1,m][1, m]就行了 ? 不行,由于系数太大,我们还要取模。将原数分别模多个质数,如果答案都为0是就可以近似认为是答案了。这样的话在洛谷上可以AC,不过在BZOJ上会TLE。为什...原创 2018-08-22 22:36:06 · 221 阅读 · 0 评论 -
HNOI2009 有趣的数列 卡特兰数
传送门 题意:题目已经说的很清楚了吧 题解: #include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int MAX = 2000001; const int MAXP原创 2018-08-23 16:05:14 · 240 阅读 · 0 评论 -
SDOI2013 随机数生成器 BSGS+数论
题解: 按照XnXnX_n的递推式,得到XnXnX_n的通项公式,移项后得到 an−1≡(Xn+b∗inv(a−1))∗inv(X1+b∗inv(a−1)) (mod p)an−1≡(Xn+b∗inv(a−1))∗inv(X1+b∗inv(a−1)) (mod&原创 2018-08-20 21:13:12 · 227 阅读 · 0 评论 -
BSGS模板
BSGS即baby step, giant step,用于求解形如ax≡b(modp)ax≡b(modp)a^{x} ≡ b \pmod{p}的方程。 具体做法:设x=i∗t−jx=i∗t−jx = i * t - j,则方程变为ai∗t≡b∗aj(modp)ai∗t≡b∗aj(modp)a^{i * t} ≡ b * a ^ {j}\pmod{p}。先把每个b∗aj(modp)b∗aj(mod...原创 2018-08-20 21:30:24 · 292 阅读 · 0 评论 -
CF869A The Artful Expedient 结论题+数论
传送门 题目描述 Tommy和Karen在玩一个游戏。 他们各自准备了一个长度为n的数组,Tommy的数组称作a,Karen的数组称作b。 保证这2n个数互不相同,设这2n个数构成的集合为S。 现在他们想知道,有多少对(i,j)满足a[i] xor b[j]的结果在集合S中。 当这个数为奇数的时候,Tommy获胜,否则Karen获胜。 现在给你n和这两个数组,你需要告诉他们是谁赢了。 输入 第一行...原创 2018-09-23 21:34:23 · 295 阅读 · 0 评论 -
高斯消元模板
高斯消元,实际上就是用加减消元法解n元一次方程组。 百度百科上已经说的很清楚了,先把方程写成矩阵的形式,依次消去各个系数x,y,z...x,y,z...x,y,z...,可以得到一个上三角矩阵。这时已经可以直接读出最后一个方程的解了,再回代到上面的方程就行了。 具体实现: 1.消元 对于xixixi,找到xixixi系数最大的一个方程,以减少精度误差。 然后,把该方程的xixixi系数转化为1,带...原创 2018-11-02 16:17:28 · 152 阅读 · 0 评论