前言

本篇文章中使用的字母\(p\),指\(\text{任意的} p \in \text{素数集合}\)

应用场景

若函数\(f(x)\)满足,

  1. \(f(x)\)是积性函数
  2. \(f(p)\)可以使用多项式表示。
  3. 已知\(f(p)\),要能在常数级的时间内计算\(f(p^x),x \in N^+\)。

Min_25筛可以在\(\Theta(\frac{n^{\frac{3}{4}}}{log_2n})\)的时间复杂度内计算\(f(x)\)的前缀和

或者说\(\Theta(n ^ {1 - \epsilon})\)?个人倾向上面那种。

算法

分类

\[\sum_{i=1}^n f(i)=\sum_{p=1}^nf(p),p \in prime\ + \sum_{q=1}^nf(q) ,q \notin prime
\]

prime指素数集合

辅助函数\(g(x,y)\)的构造与递推

首先线性筛出质数,线性筛就不多说了(如果您线性筛都不会建议先去学习基础算法)

我们从小到大设\(p_i\)为从小到大排列的第i个质数

例如,\(p_1=2,p_2=3,\dots\)

设素数集合为prime,i的最小质因子为\(MPF_i\) (minimal prime fact of i)。

\(||\)为或者,即\(or\); \([a]\)表示a成立是为\(1\),否则为\(0\)

设函数\(g(x,y)\),使得

\[g(x,y)=\sum_{i=1}^x[i\in prime\ ||\ MPF_i>p_y]f(i)
\]

考虑埃拉托斯特尼筛法,每次选出一个质数,筛掉它的倍数。

注:埃拉托斯特尼筛法的百度百科

我们发现\(g(x,y)\)有一个神奇的性质,我们的\(g(x,y)\)就是\([1,x]\)运行y次筛法以后剩余所有数之和加上所有的\(f(p),p<x\)的和。

我们所求的(prime为质数集合)

\[\sum_{i=1}^x[i\in prime]
\]

实际上等同于\(g(x,|prime|)\),|prime|表示\([1,x]\)之内的质数集合的大小。

首先我们要了解\(g(x,y)\)的初值

\(g(x,0)\)表示所有数的和,也就是把所有数带入\(f(x)\)计算出的结果。

那么接下来就是\(g(x,y)\)的递推了。

  1. \(p_y^2 > x\)。 此时埃筛的第\(j\)次没有筛去任何质数(理论上第\(y\)次应该筛掉的最小质数为\(p_y^2\)),所以\(g(x,y)=g(x,y-1)\)
  2. \(p_y^2 \leq x\)。此时埃筛筛去了所有大于\(p_y\)倍的\(p_y\)的倍数,若我们从\(g(x,y-1)\)递推,显然有多计算的部分,该部分就是

\[f(p_y)\times(g(\frac{x}{p_y},y-1)-\sum_{i=1}^{y-1}f(i))
\]

表示成公式就是:

\[g(x,y) = \left\{ \begin{matrix} g(x,y-1) & p_y^2>x \\ \\ g(x,y-1)-f(p_y)\times(g(\frac{x}{p_y},y-1)-\sum_{i=1}^{y-1}f(i)) & p_y^2 \leq x\end{matrix} \right.
\]

前缀函数\(S(x,y)\)的构造与递推

好了现在我们已经有了一个辅助函数\(g(x,y)\),但这玩意貌似一点用都没有....

我们设

\[S(x,y)=\sum_{i=1}^x[MPF_p \geq p_y]f(i)
\]

讲人话就是所有最小质因子大于等于\(p_j\)的\(f\)值之和。

如果我们要求\(\sum f(i),i \in [1,n]\),我们要求的东西是

\[S(n,1)+f(1)
\]

鉴于所有质数对答案的贡献我们已经完成计算,它的贡献是

\[g(n,j)-\sum_{i=1}^{j-1}f(p_i)
\]

(我们要保证最小质因子大于等于\(p_j\)所以要把小于它的质数减去)

考虑合数对答案的贡献,枚举合数的最小质因子和它的出现次数,然后直接计算。

合数对答案的贡献是:

\[\sum_{i=y}^{p_i^2 \leq x} \sum_{j = 1}^{p_i^{j + 1} \leq x}S(\frac{x}{p_i^j},i+1)+f(p_{i}^{j+1})
\]

总结起来就是

\[S(x,y) = g(x,|prime|) - \sum_{i=1}^{y-1}f(p_i) + \sum_{j=y}^{p_j^2 \leq x} \sum_{k = 1}^{p_j^{k + 1} \leq x}S(\frac{x}{p_j^k},j+1)+f(p_{j}^{k+1})
\]

\(prime\)指素数集合,\(||\)表示取当前的最大质因子。

复杂度比较



借用一下WC2019课件中的一张图片

例题

[LOJ6053]简单的函数

代码

#include <cstdio>
#include <cmath>
#define ll long long
#define MAXN 1000005
#define MOD 1000000007 long long sqrtn; ll read(){
ll x = 0; int zf = 1; char ch = ' ';
while (ch != '-' && (ch < '0' || ch > '9')) ch = getchar();
if (ch == '-') zf = -1, ch = getchar();
while (ch >= '0' && ch <= '9') x = x * 10 + ch - '0', ch = getchar(); return x * zf;
} // the number n, cac the sum of f(i) i belongs to [1, n]
ll n;
// primes number of primes in total
ll p[MAXN], pcnt;
// g is the function g, the sum of the f(prime).
ll g[MAXN];
// the number of the prime less the pn_i
ll pn[MAXN];
// the id of every helpful number
ll id[2][MAXN];
// prefix of p
ll prep[MAXN]; // the different value of zcfk, or (n/i)
// zcfk -> https://www.cnblogs.com/linzhengmin/p/11061244.html
ll w[MAXN], wcnt;
// is not a prime -> array in the Euler Sieve
bool np[MAXN]; // Euler Sieve
void getPrime(int n){
np[0] = np[1] = 0;
for (int i = 2; i <= n; ++i){
if (!np[i]) p[++pcnt] = i, prep[pcnt] = (prep[pcnt - 1] + i) % MOD;
for (int j = 1; j <= pcnt && i * p[j] <= n; ++j){
np[p[j] * i] = 1;
if (!(i % p[j])) break;
}
}
} // cac function s
int S(ll x, int y){
if (x <= 1 || p[y] > x) return 0;
int k = (x <= sqrtn) ? id[0][x] : id[1][n / x];
int res = (((ll)g[k] - pn[k] - prep[y - 1] + y - 1) % MOD + MOD) % MOD;
if (y == 1) res += 2;
for (int i = y; i <= pcnt && (ll)p[i] * p[i] <= x; ++i){
ll p1 = p[i], p2 = (ll)p[i] * p[i];
for (int j = 1; p2 <= x; ++j, p1 = p2, p2 *= p[i])
(res += (1ll * S(x / p1, i + 1) * (p[i] ^ j) % MOD + (p[i] ^ (j + 1))) % MOD) %= MOD;
}
return res;
} int main(){
n = read(); sqrtn = sqrt(n); getPrime(sqrtn);
// zcfk -> caculate function g, init array w, id
for (ll i = 1, j; i <= n; i = j + 1){
w[++wcnt] = n / i, j = n / w[wcnt];
(w[wcnt] <= sqrtn) ? id[0][w[wcnt]] = wcnt : id[1][j] = wcnt;
pn[wcnt] = (w[wcnt] - 1) % MOD, g[wcnt] = (w[wcnt] % MOD) * ((w[wcnt] + 1) % MOD) % MOD;
if (g[wcnt] & 1) g[wcnt] = g[wcnt] + MOD;
g[wcnt] >>= 1; g[wcnt]--;
}
for (int j = 1; j <= pcnt; ++j)
for (int i = 1; i <= wcnt && p[j] * p[j] <= w[i]; ++i){
int k = (w[i] / p[j] <= sqrtn) ? id[0][w[i] / p[j]] : id[1][n/(w[i] / p[j])];
// minus the extra part
(g[i] -= (ll)p[j] * (g[k] - prep[j - 1]) % MOD) %= MOD;
// refresh pn
(pn[i] -= pn[k] - j + 1) %= MOD;
}
// answer is S(n,1) plus f(1) = S(n,1) + 1
int ans = S(n, 1) + 1;
printf("%d\n", (ans + MOD) % MOD);
return 0;
}

[算法]Min_25筛的更多相关文章

  1. 关于 min_25 筛的入门以及复杂度证明

    min_25 筛是由 min_25 大佬使用后普遍推广的一种新型算法,这个算法能在 \(O({n^{3\over 4}\over log~ n})\) 的复杂度内解决所有的积性函数前缀和求解问题(个人 ...

  2. LG5325 【模板】Min_25筛

    P5325 [模板]Min_25筛 题目背景 模板题,无背景. 题目描述 定义积性函数$f(x)$,且$f(p^k)=p^k(p^k-1)$($p$是一个质数),求 $$\sum_{i=1}^n f( ...

  3. min_25筛入门

    目录 1.什么是min_25筛 2.前置知识 2.1.数论函数 2.2.埃拉托色尼筛 2.3.欧拉筛 3.min_25筛 3.1.计算质数贡献 3.2.计算总贡献 3.3.实现 4.例题 4.1.[L ...

  4. Note -「Min_25 筛」“你就说这素因子你要不要吧?你要不要?”

      赛上想写,Track Lost 了属于是. \(\mathscr{Intro}\)   Min_25 筛是用于求积性函数前缀和,同时顺带求出一些"有意思"的信息的筛法.   一 ...

  5. 【UOJ448】【集训队作业2018】人类的本质 min_25筛

    题目大意 给你 \(n,m\),求 \[ \sum_{i=1}^n\sum_{x_1,x_2,\ldots,x_m=1}^i\operatorname{lcm}(\gcd(i,x_1),\gcd(i, ...

  6. Min_25 筛 学习笔记

    原文链接https://www.cnblogs.com/zhouzhendong/p/Min-25.html 前置技能 埃氏筛法 整除分块(这里有提到) 本文概要 1. 问题模型 2. Min_25 ...

  7. UOJ188 Sanrd Min_25筛

    传送门 省选之前做数论题会不会有Debuff啊 这道题显然是要求\(1\)到\(x\)中所有数第二大质因子的大小之和,如果不存在第二大质因子就是\(0\) 线性筛似乎可以做,但是\(10^{11}\) ...

  8. 【SPOJ】DIVCNTK min_25筛

    题目大意 给你 \(n,k\),求 \[ S_k(n)=\sum_{i=1}^n\sigma_0(i^k) \] 对 \(2^{64}\) 取模. 题解 一个min_25筛模板题. 令 \(f(n)= ...

  9. 【51NOD1847】奇怪的数学题 min_25筛

    题目描述 记\(sgcd(i,j)\)为\(i,j\)的次大公约数. 给你\(n\),求 \[ \sum_{i=1}^n\sum_{j=1}^n{sgcd(i,j)}^k \] 对\(2^{32}\) ...

随机推荐

  1. LeetCode算法题-Goat Latin Easy(Java实现)

    这是悦乐书的第322次更新,第344篇原创 01 看题和准备 今天介绍的是LeetCode算法题中Easy级别的第192题(顺位题号是824).给出句子S,由空格分隔的单词组成.每个单词仅由小写和大写 ...

  2. 社工 - By浏览器 - Google搜索技巧 - 汇总

    google基本语法 Index of: 使用它可以直接进入网站首页下的所有文件和文件夹中 intext: 将返回所有在网页正文部分包含关键词的网页 intitle: 将返回所有网页标题中包含关键词的 ...

  3. Vue.js系列:生命周期钩子

    开发人员提供了一个Web开发人员可以在Vue.js应用程序的整个生命周期中使用的各种方法的讨论. 生命周期钩子是在Vue对象生命周期的某个阶段执行的已定义方法.从初始化开始到它被破坏时,对象都会遵循不 ...

  4. 【Qt开发】Qt5.7中文显示乱码解决方法两种

    升级到Qt5.X之后,原先解决的Qt显示中文乱码的方法突然不适用了,找了很多方式来解决这个问题 第一种: 在公司代码里看到的方法,先将对应的cpp文件用windows自带的记事本打开,另存为UTF-8 ...

  5. 手把手教你SOAP访问webservice并DOM解析返回的XML数据(转)

    http://blog.csdn.net/u012534831/article/details/51357111 前言: 目前我们项目组还在采用webservice这种http方式,并且某些网站服务提 ...

  6. k线图的分形

    蜡烛图上的分形指标,作为一种特殊的K线组合形态,通过对价格的一系列的高低点的描述,辅助识别出市场潜在的突破和反转点,预判后期走势. 顶分形:相邻的五根K线,若中间那根K线最高价为这五根K线的最高价,则 ...

  7. [DS+Algo] 007 希尔排序及其代码实现

    步骤 将数组列在一个表(一行多列)中,按特定的步长进行插入排序 步长从 length/2 到 1,每次除 2 将数组转换至表是为了更好地理解这算法,算法本身还是使用数组进行排序 算法性能 (根据步长序 ...

  8. 2019 我的世界多了一个 Python

    大一时学过 C语言,大三时用 C 控制单片机,之后就没有别的和编程的交集了. 大约十天前下定决心学 Python,不开玩笑,版本我选 3. 其实我也不是 100% 的零基础,因为一方面,我学过 C:另 ...

  9. Linear Discriminant Analysis

    Suppose that we model each class density as multivariate Gaussian, in practice we do not know the pa ...

  10. javascript xml转json

    1.代码 //加载xml数据 function loadXml(str) { if (str == null) { return null; } var doc = str; try{ doc = c ...