杜教筛
_Ark
这个作者很懒,什么都没留下…
展开
-
Min25筛学习 + 【51nod1847】奇怪的数学题(Min_25筛+杜教筛)
Min25筛小结——alpc_qleonardo的博客讲的非常清楚,不过其中大S(n,j)S(n,j)S(n,j)表示的应该是从111累加到nnn的F(i)F(i)F(i),iii要么是质数,要么最小质因子大于等于jjj。这样才满足那个递推式。然后来看一道巧妙的例题。题意给出 N,KN,KN,K ,请计算下面这个式子:∑i=1N∑j=1Nsgcd(i,j)k\sum_{i=1}^N\...原创 2020-01-14 11:03:40 · 371 阅读 · 0 评论 -
BZOJ 3944: Sum (杜教筛)
惊险…卡着时限过了…注意有的地方可能爆int(因为有231−12^{31}-1231−1)有几个地方可能会蜜汁RE,具体看代码注释CODE9000ms卡过2333#include<bits/stdc++.h>using namespace std;typedef long long LL;const int MAXN = 1000005;int prime[MAXN/...原创 2019-04-04 22:30:29 · 122 阅读 · 0 评论 -
[bzoj 4176] Lucas的数论 (杜教筛 + 莫比乌斯反演)
题面设d(x)d(x)d(x)为xxx的约数个数,给定NNN,求 ∑i=1N∑j=1Nd(ij)\sum^{N}_{i=1}\sum^{N}_{j=1} d(ij)i=1∑Nj=1∑Nd(ij)N&amp;lt;=109N&amp;lt;=10^9N&lt;=109题目分析有这样一个结论d(ij)=∑x∣i∑y∣j[(x,y)==1]d(ij)=\sum_{x|i}\sum_{y|...原创 2018-12-10 22:45:03 · 244 阅读 · 0 评论 -
[51 Nod 1584] 加权约数和
题意求∑i=1N∑j=1Nmax(i,j)⋅σ1(ij)\large \sum_{i=1}^N\sum_{j=1}^Nmax(i,j)\cdot\sigma_1(ij)i=1∑Nj=1∑Nmax(i,j)⋅σ1(ij)其中1≤N≤1061\le N\le10^61≤N≤1061≤T≤5⋅1041\le T\le5\cdot10^41≤T≤5⋅104σ1(n)\sigma_1(n)σ...原创 2018-12-16 23:08:29 · 262 阅读 · 0 评论 -
狄利克雷卷积 一些常识
e=μ∗I\large e=\mu*Ie=μ∗Iφ∗I=id→φ=μ∗id\large \varphi*I=id\\\to\varphi=\mu*idφ∗I=id→φ=μ∗idσ0=I∗I→I=μ∗σ0\large \sigma_0=I*I\\\to I=\mu*\sigma_0σ0=I∗I→I=μ∗σ0∑d∣nμ(d)=[n=1]\large \sum_{d|n}\mu(d)=[n=...原创 2018-12-16 13:52:34 · 440 阅读 · 0 评论 -
[51Nod 1222] - 最小公倍数计数 (..怎么说 枚举题?)
题面求∑k=ab∑i=1k∑j=1i[lcm(i,j)==k]\large\sum_{k=a}^b\sum_{i=1}^k\sum_{j=1}^i[lcm(i,j)==k]k=a∑bi=1∑kj=1∑i[lcm(i,j)==k]1&amp;amp;amp;lt;=a&amp;amp;amp;lt;=b&amp;amp;amp;lt;=10111&amp;amp;amp;lt;=a&amp;amp;amp;lt;=b原创 2018-12-09 12:04:45 · 284 阅读 · 1 评论 -
[51Nod 1238] 最小公倍数之和 (恶心杜教筛)
题目描述求∑i=1N∑j=1Nlcm(i,j)\sum_{i=1}^N\sum_{j=1}^Nlcm(i,j)i=1∑Nj=1∑Nlcm(i,j)2&amp;amp;amp;amp;amp;amp;amp;lt;=N&amp;amp;amp;amp;amp;amp;amp;lt;=10102&amp;amp;amp;amp;amp;amp;amp;lt;=N&amp;amp;amp;amp;amp;amp;amp;lt;=10^{1原创 2018-12-06 23:07:10 · 316 阅读 · 0 评论 -
[51Nod 1220] - 约数之和 (杜教筛)
题面令d(n)d(n)d(n)表示nnn的约数之和求∑i=1n∑j=1nd(ij)\large\sum_{i=1}^n\sum_{j=1}^nd(ij)i=1∑nj=1∑nd(ij)题目分析先给结论d(ij)=∑x∣i∑y∣jxj/y[(x,y)==1]\large d(ij)=\sum_{x|i}\sum_{y|j}xj/y[(x,y)==1]d(ij)=x∣i∑y∣j∑xj/...原创 2018-12-13 23:00:27 · 421 阅读 · 0 评论 -
杜教筛时间复杂度证明
假设我们要求S(n)S(n)S(n)且有S(n)=∑i=2nS(⌊ni⌋)S(n)=\sum_{i=2}^nS(⌊\frac ni⌋)S(n)=∑i=2nS(⌊in⌋)于是我们要求在整除分块优化下的时间复杂度,只需要求出n个S(⌊ni⌋)\sqrt n个S(⌊\frac ni⌋)n个S(⌊in⌋)的值就行了假设计算出S(n)S(n)S(n)的复杂度为T(n)T(n)T(n),则有T(...原创 2018-12-03 23:06:41 · 2800 阅读 · 2 评论 -
[51Nod 1244] - 莫比乌斯函数之和 & [51Nod 1239] - 欧拉函数之和 (杜教筛板题)
[51Nod 1244] - 莫比乌斯函数之和求∑i=1Nμ(i)\sum_{i=1}^Nμ(i)∑i=1Nμ(i)开推∑d∣nμ(d)=[n==1]\sum_{d|n}\mu(d)=[n==1]d∣n∑μ(d)=[n==1]移项μ(d)=[n==1]−∑d∣n,d&amp;lt;nμ(d)∴S(N)=∑i=1Nμ(i)=∑i=1N([i==1]−∑d∣i,d&amp;lt;iμ(d...原创 2018-12-03 22:23:52 · 223 阅读 · 0 评论 -
[51Nod 1237] 最大公约数之和 (杜教筛+莫比乌斯反演)
题目描述求∑i=1n∑j=1n(i,j)&amp;nbsp;mod&amp;nbsp;(1e9+7)n&amp;amp;lt;=1010\sum_{i=1}^n\sum_{j=1}^n(i,j)~mod~(1e9+7)\\n&amp;amp;lt;=10^{10}i=1∑nj=1∑n(i,j)&amp;nbsp;mod&amp;nbsp;(1e9+7)n&amp;lt;=1010题目分原创 2018-12-05 23:06:17 · 408 阅读 · 0 评论 -
[SPOJ] DIVCNT2 - Counting Divisors (square) (平方的约数个数前缀和 容斥 卡常)
题目vjudge URL:Counting Divisors (square) Let σ0(n)\sigma_0(n)σ0(n) be the number of positive divisors of nnn.For example, σ0(1)=1\sigma_0(1) = 1σ0(1)=1, σ0(2)=2\sigma_0(2) = 2σ0(2)=2 and σ0(6)=4\...原创 2018-12-08 18:23:31 · 508 阅读 · 1 评论 -
[HDU 5608]Function(莫比乌斯反演 + 杜教筛)
题目描述有N2−3N+2=∑d∣Nf(d)N^2-3N+2=\sum_{d|N} f(d)N2−3N+2=∑d∣Nf(d)求∑i=1Nf(i)\sum_{i=1}^{N} f(i)∑i=1Nf(i) mod 109+7~mod~10^9+7 mod 109+71&lt;=T&lt;=5001&lt;=N&lt;=1...原创 2018-12-05 17:14:34 · 204 阅读 · 0 评论 -
[51Nod 1227] 平均最小公倍数 (杜教筛)
题目描述求∑i=ab∑j=1ilcm(i,j)i\large\sum_{i=a}^b\sum_{j=1}^i\frac{lcm(i,j)}ii=a∑bj=1∑iilcm(i,j)1&amp;lt;=a&amp;lt;=b&amp;lt;=1091 &amp;lt;= a &amp;lt;= b &amp;lt;= 10^91&lt;=a&lt;=b&lt;原创 2018-12-08 15:11:55 · 302 阅读 · 0 评论