2021ccpc威海 I.Distance（min25筛）

传送门：https://codeforces.com/gym/103428/problem/I
题目大意：给定一个数$n$，问你$1$到$n$的整除关系的哈斯图上，相互两点的最短距离和是多少，定义两点$i,j(i>j)$之间的边权是$\frac{i}{j}$。
题解：如果不懂整除关系的哈斯图是什么样子可以看看下图。

不难发现，在图上的每一步行走，相当于乘上一个质数$p$，或者除以一个质数$p$。那么对于$i=p_{i1}^{s_{i1}}{p}_{i2}^{s_{i2}}\cdots p_{ir}^{s_{ir}},j=p_{j1}^{s_{j1}}{p}_{j2}^{s_{j2}}\cdots p_{jr}^{s_{jr}}$，其两点的最短距离一定是
$$\sum _{p\in prime}p\times \left|s_i-s_j\right|$$
稍加变换一下，可以得到
$$\begin{array}{rl}& \sum _{p\in prime}\sum _{i=1}^{p^i\le n}\sum _{j=0}^{i-1}2p(i-j)\times {\mathrm{CNT}}_n(p^i)\times {\mathrm{CNT}}_n(p^j)\end{array}$$
其中定义${\mathrm{CNT}}_n(p^k)$为$1\sim n$中恰好被$p^k$整除且无法被$p^{k+1}$整除的数的个数。
考虑求${\mathrm{CNT}}_n(p^k)$，我们知道$\lfloor \frac{n}{p^k}\rfloor$的意义是$1\sim n$中至少可以被$p^k$整除的数的个数，类比$\lfloor \frac{n}{p^{k+1}}\rfloor$的意义就是$1\sim n$中至少可以被$p^{k+1}$整除的数的个数，用$\lfloor \frac{n}{p^k}\rfloor$减去$\lfloor \frac{n}{p^{k+1}}\rfloor$的话只会剩下只能被$p^k$整除的数的个数，即有${\mathrm{CNT}}_n(p^k)=\lfloor \frac{n}{p^k}\rfloor -\lfloor \frac{n}{p^{k+1}}\rfloor$。
但是$n$的范围在$1\sim 10^{11}$，需要考虑更深一步的优化。不难发现对于每个$\sqrt{n}<p\le n,p\in prime$，其最高幂次只能为$1$，对于它们的贡献，只有$2p\times {\mathrm{CNT}}_n(p)\times {\mathrm{CNT}}_n(1)=2p\lfloor \frac{n}{p}\rfloor (n-\lfloor \frac{n}{p}\rfloor )$。假设我们已经求得每一个整除块内的全部质数和$S_{\lfloor \frac{n}{i}\rfloor }$，就可以整除分块令答案加上$S_{\lfloor \frac{n}{i}\rfloor }\times \lfloor \frac{n}{i}\rfloor (n-\lfloor \frac{n}{i}\rfloor )$。刚好min25筛中的$g(\lfloor \frac{n}{i}\rfloor )$可以求得所有在$1\sim \lfloor \frac{n}{i}\rfloor$中的质数和，只要两两相减便可以得到每一个大于$\sqrt{n}$的整除块的全部质数和了。剩下的$\frac{\sqrt{n}}{\log n}$个质数直接暴力按上式求即可，时间复杂度$O(\frac{n^{0.75}}{\log n}+\frac{\sqrt{n}}{\log n}\times {\log }^{2}n)=O(\frac{n^{0.75}}{\log n})$。
码量不大，就是卡常卡了很久（

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=6e5+5,M=1e9+7;
const ll inv2=M/2+1;
inline ll add(ll x,ll y){return x+y>=M?x+y-M:x+y;}
inline ll sub(ll x,ll y){return x-y<0?x-y+M:x-y;}
inline ll mul(ll x,ll y){return x*y%M;}
int cnt=0,prime[N],sum[N],g[N],id1[N],id2[N],tot=0;
bitset<N> vis;
ll n,sqn,v[N];
int get(ll x){
    if(x>sqn) return id2[n/x];
    return id1[x];
}
void pre(){
    sqn=sqrt(n);
    for(int i=2;i<=sqn;++i){
        if(!vis[i]){
            sum[i]=add(sum[i-1],i);
            prime[++cnt]=i;
        }
        else sum[i]=sum[i-1];
        for(int j=1;i*prime[j]<=sqn;++j){
            vis[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
    for(ll l=1,r;l<=n;l=r+1,++tot){
        ll t=n/l;
        r=n/t;
        g[tot]=sub(mul(mul((t+1)%M,t%M),inv2),1);
        if(t>sqn) id2[r]=tot;
        else id1[t]=tot;
        v[tot]=t;
    }
    for(int i=1;i<=cnt;++i){
        ll s=1ll*prime[i]*prime[i];
        for(int j=0;s<=v[j];++j){
            int id=get(v[j]/prime[i]);
            g[j]=sub(g[j],mul(prime[i],sub(g[id],sum[prime[i-1]])));
        }
    }
}
signed main(){
    cin>>n;
    pre();
    ll res=0,prev=g[get(sqn)];
    for(int i=sqn;i;--i){
        int id=get(n/i);
        res+=mul(sub(g[id],prev),mul(i,n-i));
        prev=g[id];
    }
    for(int i=1;i<=cnt;++i){
        ll t=0;
        for(ll j=1,a=prime[i];a<=n;++j,a*=prime[i]){
            ll t1=0;
            for(ll k=0,b=1;k<=j-1;++k,b*=prime[i])
                t1+=mul(j-k,(n/b-n/b/prime[i])%M);
            t+=mul((n/a-n/a/prime[i])%M,t1%M);
        }
        res+=mul(t%M,prime[i]);
    }
    cout<<mul(res,2)<<endl;
}