数论
WAautomaton
我好弱啊~
展开
-
反演定理及其应用
一、莫比乌斯反演莫比乌斯反演公式: 若 F(n)=∑i|nf(i)F(n)=∑i|nf(i)F(n)=\sum_{i|n}f(i) 则 f(n)=∑i|nμ(i)F(ni)f(n)=∑i|nμ(i)F(ni)f(n)=\sum_{i|n}\mu(i)F(\frac{n}{i}) 为什么呢?请看证明: 由于莫比乌斯函数具有性质∑i|nμ(i)=[n=1]∑i|nμ(i)=[n=1]\...原创 2018-08-04 13:09:07 · 7851 阅读 · 0 评论 -
[Codeforces 1111E] Tree(虚树+二项式反演)
题目链接题目大意给定一棵树,有一些询问。每次询问给出kkk个点和两个数m,rm,rm,r,表示让原树以rrr为根,把这kkk个点分成至多mmm组,每组内不存在一个点是另一个点的祖先。求方案数膜1000000007.n,Q≤105,∑k≤105,m≤min(k,300)n,Q\le 10^5,\sum k\le 10^5,m\le min(k,300)n,Q≤105,∑k≤105,m≤min(...原创 2019-02-20 16:36:00 · 557 阅读 · 0 评论 -
SDOI2018 反回文串(莫比乌斯反演+Pollard-Rho)
题目链接题目大意求所有的串,满足其所有循环同构串中至少有一个串是回文串。n≤1018n\le 10^{18}n≤1018题解第一步我就想偏了orz……我以为要分析这样串的性质……考虑所有回文串,共有k⌈n2⌉k^{\left\lceil\frac{n}{2}\right\rceil}k⌈2n⌉个,我们考虑把每个回文串第一个字符挪到后面,不断进行这样的操作,直到形成了新的回文串为止。假...原创 2019-02-13 20:18:07 · 282 阅读 · 0 评论 -
SDOI2017 硬币游戏(数论)
题目链接题目大意给定一些长度相同的01串,接下来从第一位开始等概率随机01,直到出现某个给定的串为止。问每个串出现的概率是多少。题解神仙题orz……直接判定两两之间的关系比较困难,考虑加入辅助变量xxx表示当前没有任何串出现。比如A=001,B=100那么xxx+011=A赢+(B赢+1)+(B赢+01)即x/8=pA+pB/2+pB/4x/8=p_A+p_B/2+p_B/4x/8...原创 2019-02-12 17:11:55 · 446 阅读 · 0 评论 -
HNOI2017 抛硬币(组合数学+扩展卢卡斯)
题目链接题目大意给定a,ba,ba,b,考虑所有长度为aaa的01串和长度为bbb的01串,统计所有前者中的1严格比后者多的方案数。答案模10k10^k10k。b≤a≤1015,k≤9,a−b≤10,000b\le a\le 10^{15},k\le 9,a-b\le 10,000b≤a≤1015,k≤9,a−b≤10,000题解神仙思路……a=ba=ba=b的情况很简单,只需要考虑...原创 2019-02-21 21:57:43 · 251 阅读 · 0 评论 -
NOI2018 冒泡排序 解题记录
【题目背景】最近,小 S 对冒泡排序产生了浓厚的兴趣。为了问题简单,小 S 只研究对 1到 n的排列的冒泡排序。 下面是对冒泡排序的算法描述。 输入:一个长度为 n 的排列 p[1…n] 输出:p 排序后的结果。 for i = 1 to n do for j = 1 to n - 1 do if(p[i] > p[i + 1]) 交换 p[i] 与 p[i + 1] 的值...原创 2018-07-29 11:14:03 · 1871 阅读 · 3 评论 -
SDOI2018 旧试题(莫比乌斯反演+三元环计数)
题目链接题目大意求:∑i=1A∑j=1B∑k=1Cd(ijk)\sum_{i=1}^A\sum_{j=1}^B\sum_{k=1}^Cd(ijk)i=1∑Aj=1∑Bk=1∑Cd(ijk)其中d(i)d(i)d(i)为iii的因数个数。A,B,C≤105A,B,C\le 10^5A,B,C≤105题解居然真的会有这么美妙的结论qwq……d(ijk)=∑a∣i∑b∣j∑c∣k[gc...原创 2019-02-14 22:55:37 · 349 阅读 · 0 评论 -
类欧几里得算法推导
第一种f(a,b,c,n)=∑i=0nai+bcf(a,b,c,n)=\sum_{i=0}^n\frac{ai+b}{c}f(a,b,c,n)=i=0∑ncai+b情况一:a≥c or b≥ca\ge c~or~b \ge ca≥c or b≥cf(a,b,c,n)=f(a%c,b%c,c,n)+n(n+1)2⋅ac+bc⋅(n+1)f(a,b,c...原创 2019-01-07 21:25:17 · 262 阅读 · 0 评论 -
容斥原理证明及应用
普通的容斥原理例题给定平面上n个多边形,请求出其覆盖的总面积。n≤10, 边数≤50,000n\le 10,~边数\le 50,000n≤10, 边数≤50,000解决方案1——自适应辛普森积分该方法适应于大部分求覆盖面积的问题,但是由于精度问题,不易于实现。解决方案2——按点坐标离散化该方法实现较为复杂,在多边形数量多的时候占优势,但是边数很多的情况下就显得力不...原创 2018-09-30 23:28:43 · 2549 阅读 · 0 评论 -
埃氏筛法+线性筛法+杜教筛+min25筛总结
埃氏筛法这个筛法是最朴素的筛法了,可以在O(nloglogn)O(nloglogn)O(nloglogn)的时间内(基本O(n)O(n)O(n))筛出[1,n]中所有素数。实现非常简单,从2开始遍历,对于每个质数都暴力算出它的所有倍数并筛掉,根据欧拉的调和级数定理,这个时间是O(nlogn)O(nlogn)O(nlogn)级别的,但是只有质数才需要计算倍数,然后不知怎么回事复杂度就变成O(nl...原创 2018-09-12 20:48:28 · 3051 阅读 · 2 评论 -
快速傅里叶变换(FFT)与快速数论变换(NTT)+例题
多项式的系数表示法考虑多项式A(x)=∑i=0naixiA(x)=∑i=0naixiA(x)=\sum\limits_{i=0}^na_ix^i,其中{a0,a1,…,an}{a0,a1,…,an}\{a_0,a_1,\dots,a_n\}被称为多项式A(x)A(x)A(x)的系数向量。每个多项式都有唯一的系数向量,每个系数向量都对应唯一的多项式。多项式的点值表示法我们可以把多项式...原创 2018-08-07 16:09:49 · 7640 阅读 · 1 评论