Min_25筛

前言：

因为我不会州阁筛，所以就只有学习一个他的一个相对简单的但是功能大致相同的筛法：Min_25 筛。
首先，其时间复杂度是$\frac{n^{\frac{3}{4}}}{\ln n}$的。其功能是求出某个积性函数$f(x)$的前缀和，
并且需要满足如下几个条件：
$$\{\begin{array}{l}1.f(p)是一个多项式(\text{Let p be a prime})\\ 2.f(p^c)容易被算出\end{array}$$

做法

先规定一些符号以及一些定义：
$p_i$表示第$i$个质数，比如$p_1=2,p_2=3,p_3=5...$

首先考虑如下函数：
$${\pi }_k(n)=\sum _{i=1}^n{i}^k[\text{i is a prime}]$$
那么考虑如何求${\pi }_k(n)$。定义一个新的函数
$${\pi }_{k,j}(n)=\sum _{i=1}^n{i}^k[\text{i is a prime or i’s minimal prime factor}\ge p_j]$$首先，我们不妨一开始把所有合数都看成质数，那么此时${\pi }_{k,0}(n)=({\sum }_{i=1}^n{i}^k)-1$，之所以减一是因为一既不是合数也不是质数。
那么现在应该要去掉重复的部分，考虑如何去掉，我们枚举某个合数的最小质因数$p_i$，然后此时有$p_i^k\times {\pi }_{k,j-1}(\lfloor \frac{n}{p_i}\rfloor )$的贡献，但是注意到这会算重，因为剩下的中有可能有最小质因数比$p_i$还小的数，是哪一部分呢？明显就是前$i-1$小的质数，所以此时应当减去的是：$p_i^k\times ({\pi }_{k,j-1}(\lfloor \frac{n}{p_i}\rfloor )-{\pi }_{k,j-1}(p_i-1))$。
那么既然都可以枚举$p_i$了，那么${\pi }_{k,0}(p_i-1)$就可以直接预处理了。
然后我们注意到，这个$\lfloor \frac{n}{p_i}\rfloor$只有$\sqrt{n}$种，所以可以离散化掉。然后考虑更新顺序，我们可以看到，大的质数要在小的质数之后被筛掉，所以直接$\text{for}$一遍把质数从小到大就行了。代码如下（这里以${\pi }_0$为例）：

int val , pi0 , vcnt;
//离散化
	for (ll i = 1;i <= n; i++) {
		val[++vcnt] = n / i;
		pi0[vcnt] = val[vcnt] - 1;
		if (val[vcnt] < N) id1[val[vcnt]] = vcnt;
		else id2[n / val[vcnt]] = vcnt;//离散化的数组
		i = n / (n / i);
	}
//从这开始筛	
	for (int j = 1;j <= pcnt; j++)
	for (int i = 1;i <= vcnt && p[j] * p[j] <= val[i]; i++) {
		int k = (val[i] / p[j]) < N ? id1[val[i] / p[j]] : id2[n / (val[i] / p[j] ) ];
		pi0[i] -= pi0[k] - ( j - 1 );
	}

可以发现，实现很简单。复杂度？就是$\frac{n^{\frac{3}{4}}}{\ln n}$，然而我不会证明。

下面继续，我们设$F(x,j)$表示$$\sum _{i=1}^{x}f(i)[\text{i is a prime or i’s minimal prime factor >}{p}_j]$$，那么首先$F(x,j)$应当有所有质数的$f$值的和，这个之前已经处理过了。

然后还是看合数，同理枚举最小的质因数$p_i(i\ge j)$，枚举其次幂$e$,满足$p_i^e\le x$，那么答案便可以加上$F(\frac{x}{p_i^e},i+1)\times f(p_i^e)$，然而值得注意的是，无论$j$是几，$F$中都不包含$1$，所以我们还要另外加上$f(p_i^e)$。
这部分很好理解，相当于就是暴力枚举约数，然后递归下去。
代码如下：

ll F(ll x ,int j){
	if (x <= 1 || pri[j] > x) return 0;
	int k = x < N ? id1[x] : id2[n / x];
	ll ans = //sum of primes;
	
	for (int i = j;i <= pcnt && pri[i] * pri[i] <= x; i++ ) {
		ll t = pri[i];
		for (int e = 1; t <= x ; e++ , t = t * pri[i] ){
			ans += F(x / t1 , i + 1) * f(t) + f(t);
		}
			ans -= f(pri[i]);
	}
	
	return ans;
}

最后，我们注意到$F(n,1)+f(1)$就是答案，因为$1$没有质因子，所以$F(n,1)$算不到他。
然后就完了。

当然，我们也就注意到了这个筛可以算许多东西，比如算$\pi (x)$这种。
例如：

Loj 6028(6027)

求$[1,n]$中模$m$为$p$的质数个数。
#Solution
注意这个初始，就只把$tm+p$的数算起来，然后讨论一下就完了。就是说，如果你用质数$p_k$来更新$pi(i)$，然后$p_k==b(\mathrm{m}\mathrm{o}\mathrm{d}m)$，然后你就得用$p_k\times inv(b)$来更新一下。