【BZOJ3309】DZY Loves Math(线性筛)】的更多相关文章

一通正常的莫比乌斯反演后,我们只需要求出g(n)=Σf(d)*μ(n/d)的前缀和就好了. 考虑怎么求g(n).当然是打表啊.设n=∏piai,n/d=∏pibi .显然若存在bi>1则这个d没有贡献.考虑bi为0和1两种情况.如果只看ai最小的质因子的选取情况,会发现大部分情况下其是0还是1,对f的取值是没有影响的,但会使μ取反,于是就抵消为0.而特殊情况即为所有ai均相同,此时若所有bi都取1会使f减少.与一般情况比较可以得到此时g(n)=(-1)质因子个数+1. 然后就可以线性筛了.记录一…
$\sum\limits_{T=1}^{n}\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor\sum\limits_{d|T}f(d)\mu(\frac{T}{d})$ 求出$g(n)=\sum_{d|T}f(d)\mu(\frac{n}{d})$的前缀和,分块加速. 考虑怎么快速求g.观察什么时候d能对答案产生贡献,显然当且仅当:对于n的每个质因子,d包含这个质因子的次幂数至多比n包含这个质因子的次幂数少1,否则n/d就会包含平方因子. 接下…
bzoj2154||洛谷P1829 https://www.lydsy.com/JudgeOnline/problem.php?id=2154 https://www.luogu.org/problemnew/show/P1829 不妨设n<=m 就是求$ans=\sum_{k=1}^m{\frac{1}{k}}\sum_{i=1}^n\sum_{j=1}^m{ij[(i,j)=k]}$ 把1/k后面的那一部分提出来,设为f(k), 然后莫比乌斯反演得到f(k)较简易的计算式,代回ans,并化简…
莫比乌斯反演得 $ans=\sum g[i]\frac{a}{i}\frac{b}{i}$ 其中$g[i]=\sum_{j|i}f[j]\mu(\frac{i}{j})$ 由f和miu的性质可得 设$n=p[1]^{a[1]}p[2]^{a[2]}...p[k]^{a[k]}$ 若存在$a[i]$不等于$a[j]$,则$g[n]=0$ 否则$g[n]=(-1)^{k+1}$ 线性筛$O(n)$预处理,然后每次询问$O(\sqrt{n})$分块计算 #include<cstdio> typede…
Description 对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0. 给定正整数a,b,求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b).T<=1e4; a,b<=1e7. Solution 一开始没仔细看数据范围然后打了一个每个询问O(n)的,当然T了 (盗一张图) 一开始我按照第二行的做的,里层外层循环都和ab有关,每一层都要sqrt(n)…
题目 对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0. 给定正整数a,b,求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b). 输入格式 第一行一个数T,表示询问数. 接下来T行,每行两个数a,b,表示一个询问. 输出格式 对于每一个询问,输出一行一个非负整数作为回答. 输入样例 4 7558588 9653114 6514903 4451211 742…
link 好久没写数学题了,再这样下去吃枣药丸啊. 找一套应该还比较有意思的数学题来做. [bzoj3309]DZY Loves Math 简单推一下. \[\sum_{i=1}^n\sum_{j=1}^mf(\gcd(i,j))\\=\sum_{d=1}^nf(d)\sum_{i=1}^{n/d}\mu(i)\frac n{id}\frac m{id}\\=\sum_{T=1}^n\frac nT\frac mT\sum_{d|T}f(d)\mu(\frac Td)\] 令\(h(T)=\su…
[BZOJ3309]DZY Loves Math Description 对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0.给定正整数a,b,求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b). Input 第一行一个数T,表示询问数.接下来T行,每行两个数a,b,表示一个询问. Output 对于每一个询问,输出一行一个非负整数作为回答. Sample…
[BZOJ3309]DZY Loves Math(莫比乌斯反演) 题面 求 \[\sum_{i=1}^a\sum_{j=1}^bf(gcd(a,b))\] 其中,\(f(x)\)表示\(x\)分解质因数之后,最高的幂次 题解 完全不会莫比乌斯反演了. 先来推式子 \[\sum_{d=1}^a\sum_{i=1}^{a/d}\sum_{j=1}^{b/d}[gcd(i,j)=1]f(d)\] \[\sum_{d=1}^af(d)\sum_{i=1}^{a/d}\sum_{j=1}^{b/d}[gc…
[BZOJ3309]DZY Loves Math Description 对于正整数\(n\),定义\(f(n)\)为\(n\)所含质因子的最大幂指数.例如\(f(1960)=f(2^3×5^1×7^2)=3\),\(f(10007)=1\),\(f(1)=0\). 给定正整数\(a,b\),求\(\sum\limits_{a_i=1}\sum\limits_{b_j=1}f(\gcd(i,j))\). Input 第一行一个数\(T\),表示询问数. 接下来\(T\)行,每行两个数\(a,b\…