HDU 6683 Rikka with Geometric Sequence（杜教筛，非常规数论分块）

题目
题意：求$1,2,3,4...n$中有多少个子序列是等比数列,$n<=5e17$

题解：
首先长度为1,2的等比数列我们可以直接算出答案。
对于长度>3的等比数列，
设公比$q=\frac{a}{b}$，那么最后一项被$a^{length-1}$整除。
$a<=\sqrt[3]{n}$可以直接暴力枚举，对于每个最后一项，不同的$<a$且$gcd(a,b)=1$的$b$都可以贡献出一种等比数列，有$\phi (a)$种方案，那么答案就是${\sum }_{a=2}^n{\sum }_{k=4}\phi (a)\lfloor \frac{n}{a^{k-1}}\rfloor$

对于长度=3的等比数列。
数量为${\sum }_{a=2}^n\phi (a)\lfloor \frac{n}{a^2}\rfloor$
对于$\lfloor \frac{n}{a^2}\rfloor$可以数论分块，这个时候可以证明 的取值最多只有$\sqrt[3]{n}$ 种：当$a<=\sqrt[3]{n}$的时候，只有$\sqrt[3]{n}$种值，当$a>\sqrt[3]{n}$的时候，这个商小于等于$\sqrt[3]{n}$,因此最多也就只有$\sqrt[3]{n}$种答案。

对于$\phi (a)$需要求前缀和，那么杜教筛就好了，注意预处理的范围要尽量大，因为这个复杂度很玄，需要常数技巧的支撑。

$\mathrm{A}\mathrm{C}\mathrm{C}\mathrm{o}\mathrm{d}\mathrm{e}:$

#include<bits/stdc++.h>
#define LL long long 
#define maxn 50000005
#define mod 998244353
using namespace std;

int pr[maxn],cnt_pr,phi[maxn],sp[maxn];
bool vis[maxn];

map<LL,int>mp;

int cnt = 0;
int Phi(LL n){
	if(n < maxn) return sp[n];
	if(mp.count(n)) return mp[n];
	int &ret = mp[n] = (n % mod) * (n % mod + 1) / 2 % mod;
	for(LL i=2,nxt;i<=n;i=nxt+1){
		nxt = n/(n/i);
		ret = (ret - (nxt - i + 1ll) % mod * Phi(n/i)) % mod;
	}
	return ret;
}

int main(){
	phi[1] = sp[1] = 1;
	for(int i=2;i<maxn;i++){
		if(!vis[i]) pr[cnt_pr++] = i , phi[i] = i-1;
		for(int j=0;pr[j]*i<maxn;j++){
			vis[i * pr[j]] = 1;
			if(i % pr[j] == 0){
				phi[i * pr[j]] = phi[i] * pr[j];
				break;
			}
			phi[i * pr[j]] = phi[i] * (pr[j] - 1);
		}
		sp[i] = (sp[i-1] + phi[i]) % mod;
	}
	
	int T;
	for(scanf("%d",&T);T--;){
		LL n;scanf("%lld",&n);
		int ans = (n % mod + (n % mod) * (n % mod - 1) / 2) % mod;
		
		for(LL i=2,nxt;i*i<=n;i=nxt+1){
			nxt = floor(sqrt(n/(n/(i*i))));
			ans = (ans + 1ll * (Phi(nxt) - Phi(i-1)) * (n/(i*i)%mod)) % mod;
		}
		
		int lim = floor(pow(n,1.0/3));
		for(int i=2;i<=lim;i++)
			for(LL x=1ll*i*i*i;;x*=i){
				ans = (ans + phi[i] * 1ll * (n / x % mod)) % mod;
				if(x > n / i) break;
			}
		
		printf("%d\n",(ans+mod)%mod);
	}
}