Min_25筛

有(没)什么用?

求解积性函数 \(F\) 的前缀和

\[\sum_{i=1}^{n}F(i)\]

做法

首先假设 \(F(i)=i^k\)
设 \(P_i\) 为从小到大的第 \(j\) 个质数
设 \(g(x,j)\)

\(g(x,j)=\sum_{i=1}^{x}[i\)为质数或最小质因子$ > P_j]F(i)$

求解 \(g\)

若 \(P_j^2 > x\),则不存在 \(x\) 以内的合数的质因子大于 \(P_j\)

那么 \(g(x,j)=g(x,j-1)\)

否则,\(P_j^2 \le x\),考虑从 \(g(x,j-1)\) 推过来

显然 \(g(x,j-1)\) 中多了最小质因子为 \(P_j\) 的那些合数的贡献,设为 \(R\)

设这些合数为 \(tP_j\),贡献即 \(t^kP_j^k\)

要满足这些合数的最小质因子为 \(P_j\),\(t\) 要满足最小质因子 \(\ge P_j\)

提出 \(P_j^k\) ,所以 \(R=(g(\frac{x}{P_j},j-1)-\)小于 \(P_j\) 的质数的贡献\()P_j^k\)

也就是 \(R=P_j^k(g(\frac{x}{P_j},j-1)-g(P_{j-1}, j-1))\)

(因为 \(P_{j-1}\) 以内为质数或最小质因子$ > P_{j-1} $ 的只有质数)

总结一下就是

\[
g(x, j) = \left\{
\begin{array}\ g(x,j-1), P_j^2 > x\ g(x,j-1)-P_j^k(g(\frac{x}{P_j},j-1)-g(P_{j-1}, j-1)), P_j^2 \le x\ \end{array}\ \right.
\]

求解前缀和

设 \(h(x)=\sum_{i=2}^{x}[i\) 为质数 \(]F(i)\)

假设 \(a_x\) 为 \(P_{a_x}^2 \le x\) 的最大的数

那么 \(h(x) = g(x, a_{\sqrt x})\)

再设 \(S(n,j)=\sum_{i=1}^{n}[i\)的最小质因子$ \ge P_j]F(i)$

分成两个部分计算

  1. \(i\) 为质数,贡献即为 \(h(n)-h(P_{j-1})\)
  2. \(i\) 为合数:

枚举最小质因子 \(P_k\) 及其的指数 \(e\) (这里的 \(k\) 和 \(F(i)=i^k\) 不是一个)

贡献为

\[\sum_{k\ge j}\sum_{e}^{P_k^{e+1}}(F(P_k^e)S(\frac{n}{p_k^e},k+1)+F(P_k^{e+1}))\]

首先积性函数的性质有前面的一部分

\[\sum_{k\ge j}\sum_{e=1}^{P_k^{e+1} \le n}F(P_k^e)S(\frac{n}{p_k^e},k+1)\]

而这样就没有算到 \(F(P_k^e)\) 的贡献,加回来即可

原文地址:https://www.cnblogs.com/cjoieryl/p/9403579.html

时间: 2024-07-29 21:32:55

Min_25筛的相关文章

SP34096 DIVCNTK - Counting Divisors (general) min_25筛

\(\color{#0066ff}{ 题目描述 }\) \(σ_0(i)\) 表示\(i\) 的约数个数 求\(S_k(n)=\sum_{i=1}^n\sigma_0(i^k)\mod 2^{64}\) \(\color{#0066ff}{输入格式}\) 第一行一个T为数据组数 接下来每组数据一个n,一个k \(\color{#0066ff}{输出格式}\) 每个询问输出一行 \(\color{#0066ff}{输入样例}\) 5 1 3 2 3 3 3 10 3 100 3 \(\color{

P4213 【模板】杜教筛(Sum) min_25筛

\(\color{#0066ff}{ 题目描述 }\) 给定一个正整数\(N(N\le2^{31}-1)\) 求 \(ans_1=\sum_{i=1}^n\varphi(i)\) \(ans_2=\sum_{i=1}^n \mu(i)\) \(\color{#0066ff}{输入格式}\) 一共T+1行 第1行为数据组数T(T<=10) 第2~T+1行每行一个非负整数N,代表一组询问 \(\color{#0066ff}{输出格式}\) 一共T行,每行两个用空格分隔的数ans1,ans2 \(\c

51Nod1222 最小公倍数计数 数论 Min_25 筛

原文链接https://www.cnblogs.com/zhouzhendong/p/51Nod1222.html 题意 给定 $a,b$, 求 $$\sum_{n=a}^b \sum_{i=1}^n \sum_{j=1}^i [{\rm lcm } (i,j) = n]$$ $$a,b\leq 10^{11}$$ $${\rm Time \ Limit } = 6s$$ 题解 本题做法很多. Min_25 筛 先差分一下,转化成求前缀和. 先把原题的统计无序数对转化成统计有序数对,最终 $an

UOJ188. 【UR #13】Sanrd [min_25筛]

传送门 思路 也可以算是一个板题了吧qwq 考虑min_25筛最后递归(也就是DP)的过程,要枚举当前最小的质因子是多少. 那么可以分类讨论,考虑现在这个质因子是否就是次大质因子. 如果不是,那么就是\(S(n/p,k+1)\):如果是,那么剩下的必定是一个更大的质数,那么就需要知道一段区间内有多少个质数. 质数个数显然可以min_25筛给搞出来. 于是就做完了. 代码 #include<bits/stdc++.h> clock_t t=clock(); namespace my_std{ u

LG5325 【模板】Min_25筛

P5325 [模板]Min_25筛 题目背景 模板题,无背景. 题目描述 定义积性函数$f(x)$,且$f(p^k)=p^k(p^k-1)$($p$是一个质数),求 $$\sum_{i=1}^n f(x)$$ 对$10^9+7$取模. 输入输出格式 输入格式: 一行一个整数$n$. 输出格式: 一个整数表示答案. 输入输出样例 输入样例#1: 复制 10 输出样例#1: 复制 263 输入样例#2: 复制 1000000000 输出样例#2: 复制 710164413 说明 $f(1)=1,f(

ZOJ The Sum of Unitary Totient (Min_25筛)

题意: 给你一个函数f(n)=(p1a1-1)(p2a2-1)...(prar-1) ,n=p1a1p2a2...prar  求\[\sum\limits_{i = 1}^n {f(i)} \],\[n <  = 1e9\] 思路: \[f(p) = p - 1,f({p^k}) = {p^k} - 1\],直接用min_25筛就可以了,按道理来讲我感觉min_25筛是稳过的,但我交了几十次才过,要不是在网上看到有人用min_25筛过了我还以为我想错了 #include<bits/stdc++

[51nod1965]奇怪的式子(Min_25筛)

题面 传送门 题解 好毒啊--\(de\)了一个晚上的\(bug\)-- 题目所求为 \[\prod_{i=1}^n\sigma_0(i)^{\mu(i)+i}(mod\ 10^{12}+39)\] 那么我们把式子拆开来,变成 \[(\prod_{i=1}^n\sigma_0(i)^{\mu(i)})(\prod_{i=1}^n\sigma_0(i)^i)\] 左边 现在就是要求\(\prod_{i=1}^n\sigma_0(i)^{\mu(i)}\) 我们发现只有在所有质因子都不超过\(1\)

min_25筛题目总结

最重要的一点,g不要看成是函数,而是埃式筛第j轮后的剩下的数的F之和. 1.求区间素数个数.n≤1E11 1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long int ll; 4 const ll maxn=1E6+5; 5 ll n,prime[maxn],size,sqr,back[maxn],m,g[maxn],id1[maxn],id2[maxn]; 6 bool vis[maxn]; 7 void i

奇怪的数学题(51nod1847)——min_25筛+杜教筛+第二类斯特林数

题面 51nod1847 解析   定义$f(i)$为$i$的次大因数,特别地$f(1)=0$  那么我们就是去求$$\sum_{i=1}^{n}\sum_{j=1}^{n}f^{m}(gcd(i, j))$$ 这种东西的套路就是枚举$gcd$然后用欧拉函数处理, \begin{align*}\sum_{i=1}^{n}\sum_{j=1}^{n}f^{m}(gcd(i, j)) &= \sum_{i=1}^{n}\sum_{j=1}^{\left \lfloor \frac{n}{i}\rig