数论数学
淘7的悟空
这个作者很懒,什么都没留下…
展开
-
luogu P3599 Koishi Loves Construction
题面题解构造题。对于x=1的情况,显然,当ai=na_i=nai=n的时候,iii只能等于1。然后尝试构造前缀和%n\%n%n为0,1,−1,2,−3...0,1,-1,2,-3...0,1,−1,2,−3...的情况,发现aaa数组为n,1,n−2,3,n−3...n,1,n-2,3,n-3...n,1,n−2,3,n−3...,发现只有nnn为偶数时有解,然后就构造成功了。对于x=2的情况,发现ai=na_i=nai=n的时候iii只能等于nnn,不然在iii之后的所有前缀积就全是0了,尝原创 2020-08-03 12:56:02 · 187 阅读 · 0 评论 -
[CEOI2011]Hotel
题面题面,懒得改LaTeX了题解这题是cy大佬讲的,用了一个叫做凸优化的东西。具体想法是,先想假如没有o的限制,你可以选任意个,那么最优策略就是先让容量最小的房间选订单,选它能够接受的报酬最高的订单,维护报酬最高的订单的这个操作可以用优先队列。那么考虑如果有o的限制,应该怎么做呢?对于同一组数据,f(o)f(o)f(o)表示不同的ooo时答案的大小,那么我们从样例中可以得知f(1)−f(0)>f(2)−f(1)f(1)-f(0)>f(2)-f(1)f(1)−f(0)>f(2)−原创 2020-07-15 16:04:05 · 168 阅读 · 0 评论 -
[六省联考2017]组合数问题
题面给定 n,p,k,rn,p,k,rn,p,k,r ,求∑i=0∞Cnkik+rmodp\sum_{i=0}^{\infin}C_{nk}^{ik+r} mod pi=0∑∞Cnkik+rmodpn≤109n\le 10^9n≤109,r<k≤50r<k\le 50r<k≤50,2≤p≤230−12\le p\le 2^{30}-12≤p≤230−1题解看数据范围发现大概率是矩阵乘。C有一个递推式,Cmn=Cm−1n+Cm−1n−1C_{m}^{n}=C_{m-1}^{n}原创 2020-05-28 10:34:21 · 199 阅读 · 0 评论 -
[SDOI2017]数字表格
题面∏i=1n∏j=1mfgcd(i,j)\prod_{i=1}^{n}\prod_{j=1}^mf_{gcd(i,j)}i=1∏nj=1∏mfgcd(i,j)题解这题神题啊!把所有莫比乌斯反演题里的sigma改成了pi,然后用相同的思路得出答案,这题非常考验对莫比乌斯反演题的熟练程度,同时需要一定的变通能力,属于很有意思的数学题(这里再吹一下山东省选的数论题,好多都是神题)∏i=1...原创 2020-03-22 23:28:40 · 104 阅读 · 0 评论 -
[CQOI2015]选数
题面∑LH∑LH...∑LH[gcd(a1,a2,...,an)=k]\sum_{L}^{H}\sum_{L}^{H}...\sum_{L}^{H}[gcd(a_1,a_2,...,a_n)=k]∑LH∑LH...∑LH[gcd(a1,a2,...,an)=k](n个sigma)题解这题比较创新的地方是选出n个数,但是由于之前做过选4个数的题,看到这道就直接秒了。最后的式子是这...原创 2020-03-22 22:07:55 · 105 阅读 · 0 评论 -
[SDOI2014]数表
题面∑i=1n∑j=1m[gcd(i,j)≤a]∗gcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{m}[gcd(i,j)\le a]*gcd(i,j)i=1∑nj=1∑m[gcd(i,j)≤a]∗gcd(i,j)题解看到式子明显要莫比乌斯反演,于是变成:∑x=1min(n,m)⌊nx⌋⌊mx⌋∑d∣xaμ(xd)d\sum_{x=1}^{min(n,m)}\lflo...原创 2020-03-22 10:40:48 · 200 阅读 · 0 评论 -
[SDOI2015]约数个数和
题面∑i=1n∑j=1md(ij)\sum_{i=1}^n\sum_{j=1}^md(ij)i=1∑nj=1∑md(ij)d(x)d(x)d(x)表示x的约数个数题解其实写这篇题解主要是为了证明下面这个式子:d(ij)=∑x∣i∑y∣j[gcd(i,j)=1]d(ij)=\sum_{x|i}\sum_{y|j}[gcd(i,j)=1]d(ij)=x∣i∑y∣j∑[gcd(i...原创 2020-03-14 23:43:45 · 204 阅读 · 0 评论 -
萌新自己整理的关于OI中的数论题技巧
首先,我认为OI中的数论题和数竞中的数论题的本质区别在于,数竞中的数论题往往需要一个结论性的答案,而OI中的数论题则需要考虑你的式子推导到这一步的复杂度,然后得出一个并非结论性的答案。所以OI中的数论题往往需要互换求和号、莫比乌斯反演等等的技巧,但是数竞中的数论题需要各种分类讨论,以及和OI数论题完全不同的技巧。所以说到这里,我总结出一些套路性的公式和技巧:d(ij)=∑x∣i∑y∣jgcd...原创 2020-03-14 14:47:18 · 521 阅读 · 0 评论 -
[国家集训队]Crash的数字表格
题目背景数学题没有背景题面∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n}\sum_{j=1}^mlcm(i,j)i=1∑nj=1∑mlcm(i,j)题解直接推式子吧:∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^{n}\sum_{j=1}^mlcm(i,j)i=1∑nj=1∑mlcm(i,j)=∑i=1n∑j=1mi∗jgcd(i,j)=\s...原创 2020-03-11 10:01:02 · 111 阅读 · 0 评论 -
[EER1]礼物
题目背景Update:时限扩大到 3 秒。题目描述小 Z 送了你一个数列,具体的,有 a1=1a_1=1a1=1,a2=2a_2=2a2=2,ai=2ai−1+kai−2(3≤i≤n)a_i=2a_{i-1}+ka_{i-2}(3\le i\le n)ai=2ai−1+kai−2(3≤i≤n),其中 n 是数列的长度,k 是她设定的一个正整数参数。小 Z 告诉你一个秘密,这个数...原创 2020-03-09 11:22:25 · 220 阅读 · 2 评论 -
[NOI2010]能量采集
题面:数学题,就不放题面了题解首先得出式子:2∑i=1n∑j=1mgcd(i,j)−n∗m2\sum_{i=1}^{n}\sum_{j=1}^{m}gcd(i,j)-n*m2i=1∑nj=1∑mgcd(i,j)−n∗m然后如果还记得这个式子(有的人管这个叫欧拉反演)的话n=∑d∣nφ(d)n=\sum_{d|n}\varphi (d)n=d∣n∑φ(d)那这题就O(N)切了(...原创 2020-03-03 22:20:39 · 141 阅读 · 0 评论 -
数论数学:斐波那契与黄金分割数
斐波那契与黄金分割数详见生成函数(一),里面有对斐波那契数列通项公式的推导。可以得出斐波那契通项公式为:f(n)=15.[(1+52)n−(1−52)n]f(n)=\frac1{\sqrt5}.[(\frac{1+\sqrt5}{2})^n-(\frac{1-\sqrt5}{2})^n]f(n)=51.[(21+5)n−(21−5)n]由于∣1−52∣<1|\frac{...原创 2020-02-12 16:28:26 · 734 阅读 · 0 评论 -
洛谷P6046 纯粹容器 senpai
题目描述白王制造了 nnn 个容器,并将它们排成了一队,从左到右依次编号为 1∼n1 \sim n1∼n。第 iii 个容器的强度为 aia_iai,保证 aia_iai 互不相同。为了挑选出最纯粹的容器,白王会进行 n−1n-1n−1 轮操作,每轮操作中,他会等概率随机挑选两个 位置相邻 且 未被击倒 的容器,令它们进行决斗,在一次决斗中,强度较小的容器将会被击倒并移出队列。显然最后留下...原创 2020-02-05 22:09:36 · 290 阅读 · 0 评论 -
prufer序列(HNOI2004树的计数以及HNOI2008明明的烦恼)
什么是prufer序列prufer序列是对于每一棵树的生成序列,同时对于每个prufer序列都可以生成一棵树。树生成prufer序列的原则每次找到图上度数为一的点(由于是树,所以必然能找到)中标号最小的那个,删除,并把连向它的点加入序列。...原创 2020-02-05 21:58:19 · 120 阅读 · 0 评论 -
洛谷P6016-[CSGRound3]出游
题目描述学校组织了一次暑期出游活动,报名将在第 TT 天截止。一共有 nnn 位同学,第 iii 位同学有 aia_iai 位朋友。朋友关系是单向的,换句话说,小 Z 有一个朋友是小 Y,并不意味着小 Y 一定也有一个朋友是小 Z。另外,自己也可能是自己的朋友。第 000 天时,每位同学会决定自己是否参加活动。第iii位同学有$ p_i$的概率决定参加,1−pi1-p_i1−pi的概率决...原创 2020-01-31 21:38:04 · 201 阅读 · 0 评论 -
[WC2011]最大XOR和路径
题目描述考虑一个边权为非负整数的无向连通图,节点编号为 111 到 NNN,试求出一条从 111 号节点到 NNN 号节点的路径,使得路径上经过的边的权值的 XOR 和最大。路径可以重复经过某些点或边,当一条边在路径中出现了多次时,其权值在计算 XOR 和时也要被计算相应多的次数,具体见样例。输入格式输入文件 xor.in 的第一行包含两个整数 NNN 和 MMM, 表示该无向图中点的数目...原创 2019-12-17 14:46:16 · 280 阅读 · 0 评论 -
[BJWC2011]元素
题目描述相传,在远古时期,位于西方大陆的 Magic Land 上,人们已经掌握了用魔法矿石炼制法杖的技术。那时人们就认识到,一个法杖的法力取决于使用的矿石。一般地,矿石越多则法力越强,但物极必反:有时,人们为了获取更强的法力而使用了很多矿石,却在炼制过程中发现魔法矿石全部消失了,从而无法炼制出法杖,这个现象被称为“魔法抵消” 。特别地,如果在炼制过程中使用超过一块同一种矿石,那么一定会发生“...原创 2019-12-10 10:56:48 · 95 阅读 · 0 评论 -
[cqoi2013]新Nim游戏
题目描述传统的Nim游戏是这样的:有一些火柴堆,每堆都有若干根火柴(不同堆的火柴数量可以不同)。两个游戏者轮流操作,每次可以选一个火柴堆拿走若干根火柴。可以只拿一根,也可以拿走整堆火柴,但不能同时从超过一堆火柴中拿。拿走最后一根火柴的游戏者胜利。本题的游戏稍微有些不同:在第一个回合中,第一个游戏者可以直接拿走若干个整堆的火柴。可以一堆都不拿,但不可以全部拿走。第二回合也一样,第二个游戏者也有这...原创 2019-12-09 16:31:20 · 149 阅读 · 0 评论 -
杜教筛
狄利克雷卷积与一些积性函数积性函数:e(x)=[x==1]e(x)=[x==1]e(x)=[x==1]I(x)=1I(x)=1I(x)=1id(x)=xid(x)=xid(x)=xφ(x)=∑i=1x[gcd(i,x)==1]\varphi(x)=\sum_{i=1}^{x}[gcd(i,x)==1]φ(x)=∑i=1x[gcd(i,x)==1]μ(x)(定义迷模)\mu(x)(定义迷...原创 2019-11-27 17:52:43 · 109 阅读 · 0 评论 -
伪题解[Cnoi2019]最终幻想
这里貌似要写题解,但是不可能的,正解是打表但是还好这题是个数学题,所以可以记录一下式子。题面异常简单,但令人难以动手:有一个 NNN 维超球, 求使用 KKK 个 N−1N-1N−1 维超平面可以将这个 NNN 维超球划分成多少个 NNN 维块。答案对 998244353998244353998244353 取模。然后我用了半个小时推式子,用了一个小时想写法,最后看题解发现正解是分块打...原创 2019-11-21 20:26:35 · 378 阅读 · 0 评论 -
【SDOI2008】仪仗队
数学题简而言之就是∑i=1n∑j=1n[gcd(i,j)==1]\sum_{i=1}^{n}\sum_{j=1}^{n}[gcd(i,j)==1]i=1∑nj=1∑n[gcd(i,j)==1]很明显这玩意就是:2×∑i=1nϕ(i)+12\times\sum_{i=1}^{n}\phi(i)+12×i=1∑nϕ(i)+1然后题目范围N是4e4,所以这是个线性筛板子。。然后你可能觉得本...原创 2019-11-11 19:54:06 · 313 阅读 · 0 评论 -
[国家集训队2012] 和与积
一道数论好题(突然发现我还没写莫比乌斯反演的博客,先挖个坑吧。。)原创 2019-10-18 21:38:38 · 146 阅读 · 0 评论 -
数论数学:生成函数(一)
嗯,生成函数是个函数。。生成函数即母函数,是组合数学中尤其是计数方面的一个重要理论和工具。生成函数的应用简单来说在于研究未知(通项)数列规律,用这种方法在给出递推式的情况下求出数列的通项,生成函数是推导Fibonacci数列的通项公式方法之一。简而言之,生成函数就是把计数问题(排列组合)转化成代数问题的这么一个函数,神奇的是它可以推例如斐波那契之类的数列的通项公式。|------------...原创 2019-09-29 23:18:00 · 486 阅读 · 0 评论 -
数论数学:欧拉恒等式的证明
今天突然想到我应经写过泰勒展开了!正好又遇到了欧拉恒等式,就顺便在博客上记录一下啦。。题目:试证明:eπi+1=0试证明:e^{\pi i}+1=0试证明:eπi+1=0先将exe^xex这个东西泰勒展开:eπi=∑k=0∞πkikk!=1+πi−π2/2−π3i/6+π4/24……e^{\pi i}=\sum_{k=0}^\infty \frac{\pi^k i^k}{k!}=1...原创 2019-09-23 20:04:50 · 2812 阅读 · 0 评论 -
UVA10529 Dumb Bones -(数学期望,DP)
额,看见很多大佬写了DP,然后这里提供一个记忆化搜索的方法最开始的想法对于一段区间[i,j],往里面插入第j−i+2j-i+2j−i+2块骨牌,最后复杂度是O(n3)O(n^3)O(n3),转移方程:f[i][j]=min(f[i][k−1]×1−pr1−pl−pr+f[k+1][j]×1−pl1−pl−pr+11−pl−pr)(i<k<j)f[i][j]=min...原创 2019-08-20 13:17:46 · 177 阅读 · 1 评论 -
HAOI 2016 放棋子 -(组合数学)
题目描述给你一个N*N的矩阵,每行有一个障碍,数据保证任意两个障碍不在同一行,任意两个障碍不在同一列,要求你在这个矩阵上放N枚棋子(障碍的位置不能放棋子),要求你放N个棋子也满足每行只有一枚棋子,每列只有一枚棋子的限制,求有多少种方案。输入格式第一行一个N,接下来一个N*N的矩阵。N<=200,0表示没有障碍,1表示有障碍,输入格式参考样例输出格式一个整数,即合法的方案数。这题...原创 2019-08-20 13:19:30 · 452 阅读 · 3 评论 -
数论数学:生成函数(二)
继续延续上次生成函数(一)中探讨的解递推式的问题|----------------------------------------------------------------------------------------------------------------------|前置知识卷积二项式定理对于两个数列<Ai><A_i><Ai>...原创 2019-10-01 17:26:09 · 254 阅读 · 0 评论 -
数论数学:所有自然数之和为-1/12的证明
这篇博客也是因为泰勒展开而水出来的一篇求证:1+2+3+4+....=−1121+2+3+4+....=-\frac{1}{12}1+2+3+4+....=−121证明:设S=1+2+3+4+....S=1+2+3+4+....S=1+2+3+4+....则根据泰勒展开:11−x=1+x+x2+x3...\frac{1}{1-x}=1+x+x^2+x^3...1−x1=1+x+x2+x3...原创 2019-10-02 22:05:02 · 18166 阅读 · 0 评论 -
数论数学:生成函数(三)
看之前建议重温数论数学:生成函数(一),因为今天的内容有一定前置知识|----------------------------------------------------------------------------------------------------------------------|前置知识二元一次方程组求解(真的)|----------------------...原创 2019-10-03 22:15:55 · 148 阅读 · 0 评论 -
数论数学:泰勒展开
WHAT IS 泰勒展开!!??泰勒展开是个好东西呀!有什么用捏?他可以把不可能计算的各种非人函数的估值计算出来。例如:f(x)=sin(x)和f(x)=cos(x)f(x)=sin(x)和f(x)=cos(x)f(x)=sin(x)和f(x)=cos(x)然后把他们化成一个:f(x)=∑0∞aixif(x)=\sum_0^\infin a_ix^if(x)=0∑∞aixi...原创 2019-06-12 19:53:50 · 656 阅读 · 0 评论