学习笔记（数论）

CF1139D

最后一个整除分块 + 前缀和 复杂度O( sqrt(m) )
预处理O(m)，所以没必要整除分块了~

#include <bits/stdc++.h>
#define ll long long
#define endl '\n'
using namespace std;
const int maxn=1e5+5;
const int mod=1e9+7;
int prime[maxn],inv[maxn],mo[maxn],cnt=0;
bool vis[maxn];
void init(){
	inv[1]=1;mo[1]=1;
    for(int i=2;i<maxn;i++)
        inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
	for(int i=2;i<maxn;i++){
		if(!vis[i]){
			prime[cnt++]=i;
			mo[i]=-1;
		}
		for(int j=0;j<cnt&&i*prime[j]<maxn;j++){
			vis[i*prime[j]]=1;
			if(i%prime[j]==0){
				mo[i*prime[j]]=0;
				break;
			}
			mo[i*prime[j]]=-mo[i];
		}
	}
	for(int i=1;i<maxn;i++) mo[i]=mo[i-1]+mo[i];
}
int main(){
	ios::sync_with_stdio(false);cin.tie(0);
	init();
	int n,ans=1;cin>>n;
	for(int i=2,j;i<=n;i=j+1){
		j=n/(n/i);
		ans=(ans-1ll*(mo[j]-mo[i-1])*(n/i)%mod*inv[n-n/i]%mod)%mod;
	}
	cout<<(ans+mod)%mod<<endl;
}