![](https://img-blog.csdnimg.cn/20201014180756738.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
文章平均质量分 81
200815147
这个作者很懒,什么都没留下…
展开
-
[BZOJ]4403: 序列统计 Lucas定理
Description给定三个正整数N、L和R,统计长度在1到N之间,元素大小都在L到R之间的单调不降序列的数量。输出答案对10^6+3取模的结果。首先考虑一个长度为n的序列,把第i个数+i,问题就转化为求单调上升的序列数量,元素的取值范围也由[l,r]变为[l+1,r+n],所以方案数为C(r+n-l,n),然后把所有的加起来就行了,最后要求的变为C(n+r-l+1,r-l+1)原创 2017-06-30 20:50:07 · 515 阅读 · 0 评论 -
4737: 组合数问题 lucas定理+数位DP
Description 组合数C(n,m)表示的是从n个物品中选出m个物品的方案数。举个例子,从(1,2,3)三个物品中选择两个物品可以有(1,2),(1,3),(2,3)这三种选择方法。根据组合数的定义,我们可以给出计算组合数C(n,m)的一般公式: C(n,m)=n!/m!*(n?m)! 其中n!=1×2×?×n。(额外的,当n=0时,n!=1) 小葱想知道如果给定n,m和k,对于所...原创 2018-03-19 09:29:29 · 812 阅读 · 0 评论 -
2017 计蒜之道 复赛 E 商汤智能机器人 lucas定理
题解:推式子过程可以看官方题解,推出式子后再用lucas拆一下,发现有一个部分是可以递归算的,另外一个部分是公共的,然后就可以了,这种lucas+递归的套路以前也见过,如bzoj4591: [Shoi2015]超能粒子炮·改。代码:#include<bits/stdc++.h>using namespace std;#define LL long long#d...原创 2018-03-19 09:33:40 · 341 阅读 · 0 评论 -
5027: 数学题 扩展欧几里得
Description 给出a,b,c,x1,x2,y1,y2,求满足ax+by+c=0,且x∈[x1,x2],y∈[y1,y2]的整数解有多少对?题解:直接exgcd求出一组解之后乱搞就行了,注意对a,b等于0、正负数的处理就行了。代码:#include<bits/stdc++.h>using namespace std;#define LL long lo...原创 2018-03-14 14:39:37 · 253 阅读 · 0 评论 -
2694: Lcm 莫比乌斯反演
题解:莫比乌斯反演复习。 题目即求: ∑i=1A∑j=1B[μ(gcd(i,j))!=0]lcm(i,j)∑i=1A∑j=1B[μ(gcd(i,j))!=0]lcm(i,j)\sum_{i=1}^A\sum_{j=1}^B[\mu(gcd(i,j))!=0]lcm(i,j) 枚举gcd: ∑d=1min(A,B)[μ(d)==1]d∑i=1⌊Ad⌋∑j=1⌊Bd⌋[gcd(i,j)==...原创 2018-04-10 11:01:51 · 359 阅读 · 1 评论 -
2986: Non-Squarefree Numbers 莫比乌斯函数+二分+容斥
Description 一个正整数K被称为squarefree,如果它没有一个D^2(D>1)这样的约数。题解:这道题和2440几乎一样啊。首先二分答案,问题变成求小于等于x有多少个squarefree,然后用容斥原理来算的话,就是至少有一个平方数因子的-至少有两个平方数因子的+至少有三个平方数因子的……然后对于一个完全平方数x=i2x=i2x=i^2,他的贡献就是μ(i)×⌊...原创 2018-01-09 11:20:44 · 251 阅读 · 0 评论 -
burnside引理
前言今天大概学习了一下burnside引理,下面来小结一下。参考资料南京外国语学校 陈瑜希 的集训队论文 polya计数法的应用基本概念(由于作者懒,直接截图) 置换群顾名思义,就是一个元素都是置换的群。对于置换群,它的二元运算就是置换之间的连接,比如说(13213244)∗(14233241)=(12243341)(12343124)∗(12344321)=...原创 2018-08-10 21:42:53 · 3572 阅读 · 0 评论 -
#2578. 「TJOI2018」教科书般的亵渎 拉格朗日插值
题解:显然这题就是求一个自然数幂和减去一些东西,自然数幂和是可以用拉格朗日插值求的,这个东西其实很简单,自己推也能推出来,主要是第一次写,记录一下。代码其实可以优化成O(Tm3)O(Tm3)O(Tm^3),但是过了就懒得优化了。代码:#include<cstdio>#include<cstring>#include<iostream>#incl...原创 2018-08-14 20:35:36 · 304 阅读 · 2 评论 -
「FJOI2016」建筑师 DP 数学
题解:先考虑暴力。把数字从大到小往序列中插,就可以得到一个DP的方法:f[i][j][k]f[i][j][k]f[i][j][k]表示当前插了前iii大的数,左边能看到jjj个,右边能看到kkk个,转移显然:f[i][j][k]=f[i−1][j][k]×(i−2)+f[i−1][j−1][k]+f[i−1][j][k−1]f[i][j][k]=f[i−1][j][k]×(i−2)+f[i−1...原创 2018-08-31 11:38:13 · 239 阅读 · 0 评论 -
2749: [HAOI2012]外星人 数论
题解:这个题的关键是要找出一个与操作次数相关的量。这个量就是222在取ϕ\phiϕ过程中的222因子个数。因为如果一开始没有222,那么在某一个大于222的质数−1-1−1后就会出现,一开始有222,那么这次就会把这个222变成111。所以算一下过程中222的个数即可。代码:#include<bits/stdc++.h>using namespace std;#define ...原创 2018-09-25 16:50:57 · 192 阅读 · 1 评论 -
[BZOJ]3142: [Hnoi2013]数列 数学
Description小T最近在学着买股票,他得到内部消息:F公司的股票将会疯涨。股票每天的价格已知是正整数,并且由于客观上的原因,最多只能为N。在疯涨的K天中小T观察到:除第一天外每天的股价都比前一天高,且高出的价格(即当天的股价与前一天的股价之差)不会超过M,M为正整数。并且这些参数满足M(K-1)<N。小T忘记了这K天每天的具体股价了,他现在想知道这K天的股价有多少种可能题解:...原创 2018-12-23 19:59:28 · 168 阅读 · 0 评论 -
杜教筛
复习一波。应用在O(n23)O(n^{2\over 3})O(n32)的时间内求某些积性函数的前缀和。怎么搞设我们要求和的积性函数为f(i)f(i)f(i),S(n)=∑i=1nf(i)S(n)=\sum_{i=1}^nf(i)S(n)=∑i=1nf(i)。先给它卷上另外一个积性函数,设为g(i)g(i)g(i)。设T=∑i=1n∑d∣if(d)g(i/d)T=\sum_{i=1}...原创 2019-01-10 13:13:22 · 256 阅读 · 0 评论 -
[BZOJ]5020 & [LOJ]#2289: [THUWC 2017]在美妙的数学王国中畅游 LCT+泰勒展开
Description数字和数学规律主宰着这个世界。机器的运转,生命的消长,宇宙的进程,这些神秘而又美妙的过程无不可以用数学的语言展现出来。这印证了一句古老的名言:“学好数理化,走遍天下都不怕。”学渣小R被大学的数学课程虐得生活不能自理,微积分的成绩曾是他在教室里上的课的最低分。然而他的某位陈姓室友却能轻松地在数学考试中得到满分。为了提升自己的数学课成绩,有一天晚上(在他睡觉的时候...原创 2019-01-11 13:44:13 · 174 阅读 · 0 评论 -
Codeforces 960G. Bandit Blues 分治NTT+第一类斯特林数+DP
题解:考虑最简单的DP,fi,jf_{i,j}fi,j表示i个数的排列,上升序列长度为j的方案数,考虑从大到小放数字,即最后一次放最小的数字,容易得到转移fi,j=fi−1,j−1+(i−1)fi−1,jf_{i,j}=f_{i-1,j-1}+(i-1)f_{i-1,j}fi,j=fi−1,j−1+(i−1)fi−1,j,这其实是第一类斯特林数的递推式。枚举n放在哪里,答案为∑i=1n...原创 2019-01-02 17:26:17 · 152 阅读 · 0 评论 -
5131: [CodePlus2017年12月]可做题2 矩阵乘法 exgcd
题解:现在看挺水的啊……我当初还不会做呢……就找一下规律,矩阵乘法一下,然后用exgcd求一下同余方程的解的个数就好了,但还是有几个地方要注意的:1、方程的负号可以直接忽略。2、判断无解。3、要注意%,比如a1a1a1,上来就要%一下。代码:#include<bits/stdc++.h>using namespace std;#define LL long long#...原创 2018-03-12 11:55:25 · 165 阅读 · 0 评论 -
1129: [POI2008]Per 中国剩余定理
Description 给你一个序列s,你把这个序列的所有不同排列按字典序排列后,求s的排名mod m题解:完全不会啊……代码都几乎是照着栋老师打的……这题思维难度不算太大(但我依然没有想出来),就是一位一位固定,算每一位上比当前排列小的有多少个,然后就可以推出式子(具体的式子可以看别人的blog),然后用中国剩余定理合并就行了,但这个东西我只在学扩展lucas的时候用过,已经忘了,于是...原创 2018-03-11 21:53:01 · 292 阅读 · 0 评论 -
[BZOJ]1485: [HNOI2009]有趣的数列 Catalan数
Description 我们称一个长度为2n的数列是有趣的,当且仅当该数列满足以下三个条件: (1)它是从1到2n共2n个整数的一个排列{ai}; (2)所有的奇数项满足a132n-1,所有的偶数项满足a242n; (3)任意相邻的两项a2i-1与a2i(1≤i≤n)满足奇数项小于偶数项,即:a2i-12i。 现在的任务是:对于给定的n,请求出有多少原创 2017-06-26 12:34:04 · 252 阅读 · 0 评论 -
bzoj 3629: [JLOI2014]聪明的燕姿
求所有约数和等于S的数。强大的dfs。 先要知道d是怎么求的。 d(n)=∏ki(1+p1+p2+p3……+pa)所以……直接上代码void cal(int last,int tot,int num){ if(tot==1){ans[++n]=num;return;} if(tot-1>sqrts&&isprime(tot-1))转载 2017-08-17 16:13:57 · 183 阅读 · 0 评论 -
[BZOJ]1965: [Ahoi2005]SHUFFLE 洗牌 扩展欧几里得
扩展欧几里得原创 2017-09-29 15:08:45 · 237 阅读 · 0 评论 -
[BZOJ]2751: [HAOI2012]容易题(easy) 暴力+快速幂
暴力+快速幂原创 2017-10-01 11:33:20 · 235 阅读 · 0 评论 -
扩展欧几里得(exgcd)
exgcd原创 2017-09-30 08:05:33 · 420 阅读 · 0 评论 -
[BZOJ]3738: [Ontak2013]Kapitał 扩展lucas
扩展lucas原创 2017-12-16 14:40:35 · 517 阅读 · 0 评论 -
[BZOJ]2142: 礼物 扩展lucas
扩展lucas原创 2017-12-16 14:37:59 · 246 阅读 · 0 评论 -
[BZOJ]2226: [Spoj 5971] LCMSum 欧拉函数(或莫比乌斯反演)
欧拉函数(或莫比乌斯反演)原创 2017-12-30 10:18:53 · 288 阅读 · 0 评论 -
数分块
数分块原创 2018-01-01 21:19:01 · 162 阅读 · 0 评论 -
欧拉定理
感觉比费马小定理要强一点。FYC转载 2018-01-02 15:16:30 · 202 阅读 · 0 评论 -
小于等于n且与n互质的数的和
小于等于n且与n互质的数的和原创 2017-12-30 10:06:10 · 1963 阅读 · 0 评论 -
Codeforces919E
前言:惭愧,昨天竟然没有做出来……题解:其实挺简单。首先像这种题目肯定是有个周期的,对于这题,一个周期的长度就是p×(p−1)" role="presentation" style="position: relative;">p×(p−1)p×(p−1)p\times(p-1)。因为n%p" role="presentation" style=原创 2018-02-01 10:39:53 · 545 阅读 · 0 评论 -
20180315模拟题C
n,m,r<=106n,m,r<=106n,m,r原创 2018-03-15 19:06:39 · 276 阅读 · 0 评论 -
[BZOJ]5193: [Usaco2018 Feb]Cow Gymnasts 数学+乱搞
Description厌倦了农场生活,奶牛们变卖了所有尘世的财产,加入了一支巡回马戏团。到现在为止,奶牛们已经能够进行简单的表演了:耍火把、走钢丝、骑独轮车——没什么是拥有灵巧蹄子的奶牛做不到的。但是,马戏团指挥想要为他们的下一次演出创作一个更加激动人心的表演。新的表演所用的舞台由N个围成一圈的平台组成。在每个平台上,1至N头奶牛叠成一摞,一头叠在一头上面。当指挥给出信号的时候,每摞奶牛同时顺时...原创 2019-02-19 14:52:23 · 330 阅读 · 0 评论