min25筛学习笔记

min25筛学习笔记

min25筛用于求

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

其中 \(f\) 是一个积性函数且能快速求 \(f(p^e)\) 。除此之外我们还要找到一个完全积性函数 \(f'\) 可以快速求前缀和且 \(f(p)=f'(p)\)

part0

先线性筛出 \(\le sqrt(n)\) 的质数,设 \(p\) 代表质数, \(p_j\) 是第 \(j\) 个质数。

part1

首先需要求出一个dp数组:

\[g(n,j) = \sum_{i\in p}^n f'(i)+\sum_{i\not \in p}^n [\text{minp}_i>p_j]f'(i) \]

即质数和最小质因子大于\(p_j\)的合数处的 \(f'_i\) 之和。

考虑递推,从 \(j-1\) 转移到 \(j\) 时减去了最小质因子为 \(p_j\) 的合数的 \(f_i'\),由于 \(f’\) 是完全积性函数,我们可以提出一个 \(f'(p_j)\) ,则

\[g(n,j)=g(n,j-1)-f'(p_j)\left(g\left(\frac{n}{p_j},{j-1}\right)-g\left(p_{j-1},j-1\right)\right) \]

然后 \(\frac{n}{x}\) 只有\(\sqrt n\) 种,可以离散化,\(g\left(p_{j-1},j-1\right)\)就预处理一下质数的 \(f_p'\) 前缀和。

part2

接着计算答案,我们考虑把质数和合数的答案分开,设 \(S(n,j)\) 表示最小质因子大于 \(p_j\) 处的 \(f_i\) 之和,那么

\[S(n,j)=g(n,\text{maxj})-g(p_{j},j)+\sum_{k>j,e,p_k^e<n}f(p_k^e)\left(S\left(\frac{n}{p_k^e},k\right)+[e\not=1]\right) \]

其中 \(\text{maxj}\)\(g\)\(j\) 的最大取值,也即最大的使 \(p_j\le sqrt(n)\)\(j\)

前面是质数的情况,\(g(n,\text{maxj})\)是所有质数的情况,因为最小质因子 $>\sqrt n $ 的合数并不存在。还要减去\(\le p_j\) 的质数的情况。

后面是合数的情况,枚举最小质因子和次数,然后利用积性函数的性质递归。

边界是\(n\le p_j\),此时\(S(n,j)=0\)

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N = 1e6+1;
#define FOR(i,a,b) for(int i=a;i<=b;++i)
const int mod = 1e9+7,iv2=500000004,iv3=333333336;
ll n,s1[N],s2[N],isp[N],pri[N],top;
ll s;
void init(int len){
	FOR(i,2,len){
		if(!isp[i]){
			pri[++top]=i;
			s1[top]=(s1[top-1]+i)%mod;
			s2[top]=(s2[top-1]+1ll*i*i)%mod;
		}
		for(int j=1;j<=top&&i*pri[j]<=len;j++){
			isp[i*pri[j]]=1;
			if(i%pri[j]==0) break;
		}
	}
}
int pos1[N],pos2[N];
ll g1[N],g2[N],tot,val[N];
ll S(ll x,ll y){
	if(pri[y]>=x) return 0;
	int k=(x<=s)?pos1[x]:pos2[n/x];
	ll res=(2ll*mod-g1[k]+g2[k]+s1[y]-s2[y])%mod;
	for(int i=y+1;pri[i]*pri[i]<=x&&i<=top;i++){
		for(ll j=pri[i],tms=1;j<=x;j*=pri[i],tms++){
			ll wn=j%mod;
			res=(res+wn*(wn-1)%mod*(S(x/j,i)+(tms!=1))%mod)%mod;
		}
	}
	return res;
}
int main(){
	scanf("%lld",&n);
	s=sqrt(n);
	init(s);
	for(ll l=1,r=0;l<=n;l=r+1){
		//printf("%lld %lld\n",l,r);
		r=min(n,n/(n/l));
		ll w=n/l,wn=w%mod;
		val[++tot]=w;
		g1[tot]=(1ll*wn*(wn+1)%mod*iv2%mod+mod-1)%mod;
		g2[tot]=((1ll*wn*(wn+1)%mod*iv2%mod*((2ll*wn+1)%mod)%mod*iv3%mod)+mod-1)%mod;
		if(w<=s) pos1[w]=tot;
		else pos2[n/w]=tot;
	}
	FOR(i,1,top){
		for(int j=1;j<=tot&&pri[i]*pri[i]<=val[j];j++){
			int nexp;if(val[j]/pri[i]<=s) nexp=pos1[val[j]/pri[i]];else nexp=pos2[n/(val[j]/pri[i])];
			g1[j]=(g1[j]-pri[i]*((g1[nexp]-s1[i-1]+mod)%mod)%mod+mod)%mod;
			g2[j]=(g2[j]-pri[i]*pri[i]%mod*((g2[nexp]-s2[i-1]+mod)%mod)%mod+mod)%mod;
		}
	}
	printf("%lld\n",(S(n,0)+1)%mod);
	return 0;
}

用法

\(\mu\)

\(\mu(p^e)=(-1)^e\)\(f'=-1\)

\(\phi\)

\(\phi(p^e)=(p-1)p^{e-1},f'=\text{id}-1\)

ZROI1838

题意:求

\[\sum_{1 \leq a_{1}, a_{2}, \cdots, a k \leq n}\left\lfloor \frac{n}{\operatorname{lcm}\left(a_{1}, a_{2}, \cdots, a_{k}\right)}\right\rfloor \]

推导:

\[=\sum_{1 \leq a_{1}, a_{2}, \cdots, a k \leq n}\sum_{a_1|j,a_2|j,\cdots,a_k|j}^n1 \\=\sum_{j=1}^n\sum_{a_1|j,a_2|j,\cdots,a_k|j}^n1 \\=\sum_{j=1}^n d(j)^k \]

于是使用min25筛,\(d(p^e)^k=(e+1)^k\)\(f'=2^k\)(常数)。

posted @ 2021-04-01 19:27  lcyfrog  阅读(199)  评论(1编辑  收藏  举报