-----------数论-----------
MaxMercer
成都七中学生,OIer
展开
-
bzoj2818 gcd 欧拉函数(前缀和)
莱昂哈德 欧拉 Description给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的 数对(x,y)有多少对.Input一个整数NOutput如题Sample Input4 Sample Output4这道题是比较水的一道数论题…… 如果说两人数的gcd为质数,那么同时除以这个质数,gcd就为1,gcd为1的话我们就能明白这两数一定互质,那么我们就枚举n以内的质数,对每一个质原创 2017-06-10 10:31:12 · 941 阅读 · 0 评论 -
FFT学习笔记
一点感慨 今天Doggu讲了, 开始觉得这玩意好强啊… 后来自己看的时候发现蛮zz…(wys表示 赞同) 所以很多自己觉得遥不可及, 一辈子都不可能学会的东西只要慢慢理解, 理解好了再把他化为自己的语言, 归到自己的世界里 – – 就觉得不是很难了, 甚至很简单. 当然开始并不是, 从”好难啊”到”这么简单困大爷我这么久,之前脑子抽了zz”. 这其实是不断简化思维吸收理解的一个过程.正文 hh虽原创 2017-11-30 21:40:40 · 407 阅读 · 5 评论 -
[BZOJ]2219 数论之神 离散对数 + 原根
2219: 数论之神Time Limit: 3 Sec Memory Limit: 259 MB Submit: 1052 Solved: 322 [Submit][Status][Discuss] Description在ACM_DIY群中,有一位叫做“傻崽”的同学由于在数论方面造诣很高,被称为数轮之神!对于任何数论问题,他都能瞬间秒杀!一天他在群里面问了一个神题: 对于给定的3个非负整原创 2017-12-01 15:58:51 · 835 阅读 · 0 评论 -
[BZOJ]3944 Sum 杜教筛
3944: SumTime Limit: 10 Sec Memory Limit: 128 MBSubmit: 4415 Solved: 1170[Submit][Status][Discuss]DescriptionInput一共T+1行第1行为数据组数T(T第2~T+1行每行一个非负整数N,代表一组询问Output一共T行,每原创 2017-12-01 20:08:30 · 457 阅读 · 0 评论 -
[BZOJ]2179 FFT快速傅立叶 FFT模板
2179: FFT快速傅立叶Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 3709 Solved: 1970[Submit][Status][Discuss]Description给出两个n位10进制整数x和y,你需要计算x*y。Input第一行一个正整数n。 第二行描述一个位数为n的正整数x。 第三行描述一个原创 2017-12-01 22:00:41 · 709 阅读 · 0 评论 -
[BZOJ]3301 [USACO2011 Feb] Cow Line 康托展开&逆康托展开
3301: [USACO2011 Feb] Cow LineTime Limit: 10 Sec Memory Limit: 128 MBSubmit: 215 Solved: 108[Submit][Status][Discuss]DescriptionThe N (1 yet another one of their crazy games with Farme原创 2017-12-14 22:15:38 · 351 阅读 · 0 评论 -
[BZOJ]2742 [HEOI2012] Akai的数学作业
附上蒟蒻的sb式推导过程, 如有不严谨请指出, 谢谢~ a0+a1∗pq+....an∗pnqn=0a_0 + a_1 * \frac {p}{q} + .... a_n * \frac {p^n}{q^n} = 0 a0=−a1∗pq...−an∗pnqna_0 = -a_1 * \frac {p}{q} ... - a_n * \frac {p^n}{q^n} 左右同时乘以qn原创 2017-12-11 09:07:59 · 326 阅读 · 0 评论 -
[BZOJ]1257 [CQOI2007] 余数之和sum 分块优化
1257: [CQOI2007]余数之和sumTime Limit: 5 Sec Memory Limit: 162 MBSubmit: 5175 Solved: 2404[Submit][Status][Discuss]Description给出正整数n和k,计算j(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值,其原创 2018-01-15 19:57:21 · 319 阅读 · 0 评论 -
[BZOJ]3160 万径人踪灭 Manacher + FFT
3160: 万径人踪灭Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 1753 Solved: 977[Submit][Status][Discuss]DescriptionInputOutputSample InputSample OutputHINTSource2013湖北互测week1[Submit][Status][Discuss]HOM原创 2018-01-17 21:48:58 · 429 阅读 · 0 评论 -
[HDU]4372 Count the Buildings 第一类斯特林数
Count the BuildingsTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 2460 Accepted Submission(s): 804Problem DescriptionThere a原创 2018-01-19 21:47:23 · 415 阅读 · 0 评论 -
[BZOJ]4555 [TJOI2016&HEOI2016] 求和 第二类斯特林数 + NTT
4555: [Tjoi2016&Heoi2016]求和Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 525 Solved: 418[Submit][Status][Discuss]Description在2016年,佳媛姐姐刚刚学习了第二类斯特林数,非常开心。现在他想计算这样一个函数的值:S(i, j)表示第二原创 2018-01-19 21:57:26 · 661 阅读 · 0 评论 -
[BZOJ]3994 [SDOI2015]约数个数和 莫比乌斯 + 分块
3994: [SDOI2015]约数个数和Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 1096 Solved: 757[Submit][Status][Discuss]Description 设d(x)为x的约数个数,给定N、M,求 Input输入文件包含多组测试数据。第一行,一个整数T原创 2017-10-05 20:47:19 · 307 阅读 · 0 评论 -
扩展欧拉定理证明
改成全角就能行首空两格了… 发现这玩意儿没有欧拉定理那么好证啊… 欧拉定理消去律做一下推推即可. 从丁学长那里学会了之后写一下证明. 论有一个善良的学长的重要性!! 终于会用makdown写数学公式了~扩展欧拉定理aq ≡ aq mod φ(m) + φ(m)(mod m) a^q \equiv a^{q mod \varphi (m) + \varphi(m)} (mod m原创 2017-11-28 21:19:34 · 1976 阅读 · 0 评论 -
bzoj3122 随机数生成器 BSGS+费马小定理求逆元
这道题应该算是数论里面比较难的一道题,我顺便借此学了BSGS(Baby Steps Giant Steps),又称力拔盖世算法(手动滑稽).刚开始看到这道题只能感到深深的无助,完全不会做…对于一道数论题如果不会做的话,大家还是可以尝试很多特判骗分的,这道题就有很多特判.从小角度切入,深入了解问题,往往是很重要的. 由于不知道怎么打出一堆数学符号,所以很难与大家分享我的题解。我在网上找了另一个人的博原创 2017-06-10 16:17:55 · 466 阅读 · 0 评论 -
bzoj2440 完全平方数 二分答案+莫比乌斯函数+容斥原理
Description小 X 自幼就很喜欢数。但奇怪的是,他十分讨厌完全平方数。他觉得这些 数看起来很令人难受。由此,他也讨厌所有是完全平方数的正整数倍的数。然而 这丝毫不影响他对其他数的热爱。 这天是小X的生日,小 W 想送一个数给他作为生日礼物。当然他不能送一 个小X讨厌的数。他列出了所有小X不讨厌的数,然后选取了第 K个数送给了 小X。小X很开心地收下了。 然而现在小 W 却记原创 2017-07-02 09:56:22 · 412 阅读 · 0 评论 -
bzoj3505 组合数
http://www.lydsy.com/JudgeOnline/problem.php?id=3505这道题斜边一开始想成o(n^2m^2)的了……实际上是一道很水的题.注意1000*1000的范围,组合数最好递推./************************************************************** Problem: 3505 User:原创 2017-07-05 20:03:38 · 381 阅读 · 0 评论 -
bzoj1042 硬币购物 容斥原理
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1042 Description 硬币购物一共有4种硬币。面值分别为c1,c2,c3,c4。某人去商店买东西,去了tot次。每次带di枚ci硬币,买s i的价值的东西。请问每次有多少种付款方法。 Input 第一行 c1,c2,c3,c4,tot 下面tot行 d1,d2,d3,d4,原创 2017-06-10 09:23:59 · 530 阅读 · 0 评论 -
齐肯多夫定理--斐波那契数列
简述齐肯多夫定理: 任何正整数都可以表示为若干个不连续的斐波那契数之和. 证明一下用F来表示斐波那契数列.数学归纳法:1.针对小部分的情况 : 1 = F(2), 2 = F(3), 3 = F(4), 这个命题是成立的, 但是我们需要证明对于任何数都是成立的.2.针对任意正整数m (1)若m为斐波那契数, 命题显然成立. (2)若m不为斐波原创 2017-08-31 17:37:52 · 2737 阅读 · 0 评论 -
[BZOJ]2142 礼物 扩展Lucas
2142: 礼物Time Limit: 10 Sec Memory Limit: 259 MBSubmit: 1788 Solved: 748[Submit][Status][Discuss]Description一年一度的圣诞节快要来到了。每年的圣诞节小E都会收到许多礼物,当然他也会送出许多礼物。不同的人物在小E心目中的重要性不同,在小E心中分量越重的人,收到的礼原创 2017-10-02 21:49:14 · 429 阅读 · 0 评论 -
NOIP模拟赛 数论专题 扩展欧几里得 + 组合数 + 容斥原理
Problem 1. payInput file: pay.in Output file: pay.out Time limit: 1 second Memory limit: 256 MB Mr.Hu 开了个饭店,来了两位客人:Alice 和Bob,他们吃完饭要结账时,发现他们需要支付c 元钱,但 是Alice 只有面值为a 的钱,Bob 只有面值为b 的钱(他们每个人的钱的和都大于c,原创 2017-10-19 21:30:24 · 732 阅读 · 0 评论 -
[BZOJ]1101 莫比乌斯+分块优化
1101: [POI2007]ZapTime Limit: 10 Sec Memory Limit: 162 MBSubmit: 2731 Solved: 1163[Submit][Status][Discuss]Description FGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x,yInput原创 2017-10-05 19:51:36 · 369 阅读 · 0 评论 -
NOIP模拟赛 数论 积性函数 + BSGS + Lucas模进制数位Dp
Problem 1. facsum题解我们发现后面那个sigma里都是积性函数, 积性函数相乘也是积性函数. 我们可以线性筛筛除后面那个就可以. 考虑后面那个和式设为积性函数g, p为质数的话, g(p) = 2 - p. g(p ^ 2) = 3 - 2p. g(p ^ 3) = 4 - 3p, 以此类推. 然后用线性筛来筛, 若i 与 pr[j]互质的话, g(i * p[j]) = g[i]原创 2017-10-20 18:43:17 · 630 阅读 · 0 评论 -
数论[模板]
如有错误请指出, 万分感谢!#include<stdio.h>#include<cmath>#include<algorithm> typedef long long lnt;const int mod = 1e9 + 7;const int maxn = 1e6 + 5;bool mark[maxn];lnt fac[maxn];//gcdlnt gcd(lnt a, lnt b原创 2017-11-09 21:00:39 · 279 阅读 · 0 评论 -
[51NOD]1288 序列求和 V4 拉格朗日插值法
就是让你求自然数幂之和嘛... 我们发现指数很小但是n很大. 本来这个东西用第一类斯特林数来搞简直极品, 因为模数都可以不是质数... 但是这道题k是万级别的, O(k ^ 2)过不了, 但是这样伯努利数不是也过不了吗? 不过好像可以NTT搞事情. 不过我们用拉格朗日插值法可以做到O(k), 因为这个和他实际上是一个k+1项多项式, 那我们插1 ~ k +2个值进去, 得到了拉格朗日插值式,原创 2018-01-19 22:08:06 · 520 阅读 · 0 评论