数论
文章平均质量分 55
Fat_tu
这个作者很懒,什么都没留下…
展开
-
洛谷3711:仓鼠的数学题(NTT+伯努利数)
题面 题意:给出a数组,求 ∑nk=0Sk(x)ak\sum_{k=0}^nS_k(x)a_k 所表示多项式的每一项系数。额,直接将伯努利数带进S里,得 =∑k=0nakk+1∑g=0kCgk+1Bgxk+1−g=\sum_{k=0}^n\frac{a_k}{k+1}\sum_{g=0}^k C_{k+1}^g B_g x^{k+1-g}=∑k=0nakk!∑g=0nBgg!xk+1−原创 2018-01-20 10:36:32 · 480 阅读 · 0 评论 -
HDU6063:RXD and math(数论)
题面 题意: ∑i=1nkμ2(i)∗nki−−−√\sum_{i=1}^{n^k}\mu^2(i)*\sqrt{\frac{n^k}{i}}1~n的每个数都可以写成a∗b2a*b^2,都对应1个a。a无平方因子。 循环就相当于枚举a,剩下的是算出有多少个b。 没问题,妥妥的。#include #include #include #include #include #原创 2018-01-19 11:11:13 · 241 阅读 · 0 评论 -
51nod1228:序列求和(伯努利数)
题面题意:求幂和∑ni=0ik\sum_{i=0}^ni^k 模1e9+7 n≤1e18,k≤2e3,5000组数据曾记得L指导koi后教我们杜教筛,普及了下幂和怎么求 一本正经地推试子 原来那就是伯努利数 还说实在不会就写高斯消元由于我还很菜,还在学,就随便吹一些。当某个数列Bi{B_i}指数生成函数为xex−1\frac{x}{e^x-1}, 我们定义它叫伯努原创 2018-01-19 21:43:14 · 399 阅读 · 0 评论 -
BZOJ3601:一个人的数论(莫比乌斯反演+伯努利数)
题面题意:给出d和n(n以分解质因数给出),问所有与n互质的数的d次幂之和,即∑xxd[gcd(x,n)==1]∑xxd[gcd(x,n)==1]\sum_xx^d[gcd(x,n)==1]套路推♂倒=∑i|nμ(i)∗id∗∑x=1nixd=∑i|nμ(i)∗id∗∑x=1nixd=\sum_{i|n}\mu(i)*i^d*\sum_{x=1}^{\frac{n}{i}}x^d fa♂...原创 2018-02-26 21:03:20 · 546 阅读 · 1 评论 -
BZOJ3561:DZY Loves Math VI(莫比乌斯反演+线性筛)
题面题意:∑i=1n∑j=1mlcm(i,j)gcd(i,j)∑i=1n∑j=1mlcm(i,j)gcd(i,j)\sum_{i=1}^n \sum_{j=1}^m lcm(i,j)^{gcd(i,j)} 以下内容十分丑陋 就像背包连着胆假就是平底锅一样 这个O(n)O(n)O(n)的外层枚举连着三个幂和,它这就是个nlog(n)nlog(n)nlog(n) 这就是我用线性筛打出的...原创 2018-02-27 15:02:53 · 251 阅读 · 0 评论 -
BZOJ3309 :DZY Loves Math(莫比乌斯反演)
题面f(i)f(i)f(i)表示将iii分解质因数后的最高次幂。 求∑i=1a∑j=1bf(gcd(i,j))∑i=1a∑j=1bf(gcd(i,j))\sum_{i=1}^a \sum_{j=1}^bf(gcd(i,j)) 老年人数论题,好像没什么好说的一轮推导,有ans=∑i=1a⌊ai⌋⌊bi⌋∑d|if(d)μ(id)ans=∑i=1a⌊ai⌋⌊bi⌋∑d|if(d)μ(id)...原创 2018-02-27 19:57:42 · 205 阅读 · 0 评论 -
BZOJ1408:Robot(数论)
题面 题意:以分解质因数形式给出m 考虑m的所有无平方因子的约数d 对有偶数个质因子的d和有奇数个质因子的d 分别求它们φφ\varphi的和 (不会说了,反正很简单)φ(p)=p−1φ(p)=p−1\varphi(p)=p-1 一算下去就好了#include <iostream>#include <fstream>#include <algo...原创 2018-03-07 21:36:52 · 161 阅读 · 0 评论 -
BZOJ4652:循环之美(数学+杜教筛)
题面 这题天秀 题意:求有多少个互质的x,y,x≤n,y≤m 使得xyxy\frac{x}{y}的小数在k进制下从第一位开始无限循环我看到这题,突然发现和我小学5年级奥数学的内容很像(我小学好厉害)当时讲了怎么把一个循环小数化成分数 从第一位开始循环的小数 循环节长度就是分母中9的个数 分子就是循环节 (证明就是等比数列求和,然后求极限,好像当时也讲了) 比如0.12345...原创 2018-03-12 20:34:22 · 283 阅读 · 0 评论 -
BZOJ3328:PYXFIB(数论+矩阵快速幂)
我觉得这是一道非常细的数学题。题面 题意:模P下,求∑k|i,inCin∗Fi\sum_{k|i,i F为斐波那契数,n≤1e18,p≤1e9,k≤2e4,k|(p−1)n≤1e18,p≤1e9,k≤2e4,k|(p-1)。其实数学题也是有套路的(根据本题现学) 看到菲波那契就应该是矩阵 看到整除P-1就应该是原根有关 看到组合数就应该是二项式先搞个算斐波那契的矩阵AA={1原创 2018-01-19 10:37:58 · 271 阅读 · 0 评论 -
洛谷2312:解方程(初等数论)
听别人说是noip很有趣的一题。 题面我一开始看成N,M都很大,以为是用NTT搞快速插值之类的,就自然想到取模。模1004535809为0,大概多项式就是0了。再仔细看数据范围,发现可以暴力搞,大概选几个四位(质)数作为模数,枚举1~10000,分别在同余的情况下算出结果。然后枚举1~m,看是否同余都为0即可。#include #include #include #in原创 2017-10-23 09:04:25 · 364 阅读 · 0 评论 -
51nod 1258:序列求和 V4(伯努利数+多项式求逆+NTT+CRT)
题面 题意:求Sk(n)S_k(n) n≤1e18,k≤5e4。用多项式逆元求出伯努利数,代公式就可以了。 上维基查了下伯努利数,B+B^+与B−B^-的差别只是B1B_1的正负而已.但我依然不知道什么时候用哪个,两个都试一下,发现这题用B+B^+。所以我做这题是为了学CRT的。假设题意要求模P,三模数为p1,p2,p3。我起初还天真地以为可以算出分别模p1,p2,p3的原创 2018-01-22 16:59:15 · 435 阅读 · 0 评论 -
洛谷3790:文艺数学题(莫比乌斯反演+矩阵树定理)
题面 题意:一个nn个点,mm条带权边的图。问所有(生成树边权GCD)的和。 n≤60,m≤3000n≤60,m≤3000某个大佬曾说过,看到GCD就知道是个反演。 所以这题就是个反演。设f(x)f(x)为gcd为xx的生成树个数 g(x)g(x)为gcd为xx倍数的生成树个数 我们可以通过枚举边权是xx倍数的边,通过矩阵树定理n3n^3地算出g(x)g(x)套路反演 an原创 2018-01-24 15:34:13 · 340 阅读 · 0 评论 -
BZOJ3560:DZY Loves Math V(数论,线性筛)
也许我应该下放一下我的懒惰标记了。 题面题意:n≤1e5,a≤1e7,求 由于积性,对每个质数p求贡献。φ(pe)=pe∗p−1p\varphi(p^e)=p^e*\frac{p-1}{p},对于某个a,我们掏空它里面p的幂,得到贡献 ∑ki=0φ(pi)=p−1p∗(−1+∑ki=0pi)+1\sum_{i=0}^k\varphi(p^i)=\frac{p-1}{p}*(-1+原创 2018-01-17 20:04:07 · 246 阅读 · 0 评论 -
BZOJ4407:于神之怒加强版(线性筛)
题面 题意:对于给定的k,求∑i=1n∑j=1mgcd(i,j)k\sum_{i=1}^n\sum_{j=1}^mgcd(i,j)^k n,m,k≤5e6,2000组数据。根据路人试子的推导(倒)方法,我们枚举gcd,得 dkd^k完全积性,可以先处理质数的结果,再算出全部。 且由于这个完全积性,它卷个μ\mu也可以用同样的方法算出来,然后对于每个询问分块就可以了。#inclu原创 2018-01-01 22:19:21 · 301 阅读 · 0 评论 -
BZOJ3512:DZY Loves Math IV(杜教筛)
很巧妙的一道题 题面 题意:∑i=1n∑j=1mφ(ij)\sum_{i=1}^n\sum_{j=1}^m\varphi(ij) n≤1e5,m≤1e9显然杜教筛是这个问题的子问题,而我又没学过除了杜教筛以外的算法 所以这题就是杜教筛了。我在语文课上突然想起了markdown写公式的恐怖。 这个试子复杂度上界为nm−−√n\sqrt{m},玄学记忆化就可以过。#inc原创 2017-12-29 15:05:38 · 566 阅读 · 0 评论 -
51nod1584:加权约数和(数论+线性筛)
离退役不远了,我依然很菜很菜,所以以后博客不再面向对象了,只为自己看懂。 题面题意:求∑i=1n∑j=1nmax(i,j)σ1(ij)\sum_{i=1}^n\sum_{j=1}^nmax(i,j)\sigma_1(ij) 其中σ1(x)\sigma_1(x)为x的约数和。分类讨论2∑i=1n∑j=1iiσ1(ij)−∑i=1niσ1(i2)2\sum_{i=1}^n\sum_{j=1原创 2017-12-25 21:42:24 · 548 阅读 · 0 评论 -
洛谷3678:简单的数学题(画柿子+杜教筛)
最近看了《魔法少女伊莉雅~雪下的誓言》, 这部番除了百合与萌以外还有其他东西, 就和当初看《fate》的我一样。题面 题意:求∑i=1n∑j=1nij(i,j)\sum_{i=1}^n\sum_{j=1}^nij(i,j)n≤1010n≤10^{10}这就是数论题,自己数学课上半个小时推出来,看Kscla的博客一个下午才懂。根据我两页草稿纸的结果,原式为∑i=1ni3∑j=1原创 2017-12-21 21:23:05 · 442 阅读 · 0 评论 -
BZOJ4514:数字配对(费用流+数学)
《Clannad》告诉我,即使有的事再努力也没有用,也要去拼命争取。 无论结果如何,我也不会后悔那次在坡道上与你的相遇。 最后半年的oi之路,就从这个白色相簿的季节开始吧。题面 题意:n个数,若a[i]/a[j]=p,p为质数,则a[i]和a[j]可以配对,产生c[i]*c[j]的价值,a[i]可以配对b[i]次,问在价值>0的前提下,最多配对多少次。不难,题目条件很符合费用流的格式原创 2017-12-07 14:57:43 · 286 阅读 · 0 评论 -
51nod1355:斐波那契的最小公倍数(数论)
题面 题意给出n个a,问LCM{ f(a) },f为斐波那契数。知乎靠谱的题解记住这两个路人性质就好 ①容斥求LCM lcm{S}=∏T⊆S,T≠∅gcd{T}(−1)|T|+1lcm{S}=∏T⊆S,T≠∅gcd{T}(−1)|T|+1\text{lcm}\{S\}=\displaystyle\prod_{T\subseteq S,T\ne \,\emptyset} \gcd\{T...原创 2018-03-13 09:56:56 · 250 阅读 · 0 评论