BZOJ-5244 最大真因数（min25筛）

题意：一个数的真因数指不包括其本身的所有因数，给定L，R，求这个区间的所有数的最大真因数之和。

思路：min25筛可以求出所有最小因子为p的数的个数，有可以求出最小因子为p的所有数之和。

那么此题就是对于所有素数因子，求它对应的和。

#include<bits/stdc++.h>
using namespace std;
#define ll unsigned long long
const int maxn=;
ll Sqr,vis[maxn],pri[maxn],sp[maxn],tot,m,id1[maxn],id2[maxn],g[maxn],h[maxn];
ll w[maxn];
void Sieve(int n)
{
    tot=; vis[]=;
    for(int i=;i<=n;i++){
        if(!vis[i]) pri[++tot]=i,sp[tot]=sp[tot-]+i;
        for(int j=;i*pri[j]<=n;j++){
               vis[i*pri[j]]=;
            if(i%pri[j]==) break;
        }
    }
}
ll solve(ll n)
{
    Sqr=sqrt(n); Sieve(Sqr); ll res=; m=;
    for(ll i=,j;i<=n;i=j+){
        j=n/(n/i); w[++m]=n/i;
        if(w[m]<=Sqr) id1[w[m]]=m;
        else id2[n/w[m]]=m;
       if(w[m]&) g[m]=(w[m]+)/*w[m]-;
       else g[m]=w[m]/*(w[m]+)-;
    }
    for(int j=;j<=tot;j++)
      for(int i=;i<=m&&pri[j]*pri[j]<=w[i];i++){
        int k=(w[i]/pri[j]<=Sqr)?id1[w[i]/pri[j]]:id2[n/(w[i]/pri[j])];
        g[i]=g[i]-pri[j]*(g[k]-sp[j-]);
        if(i==) res+=g[k]-sp[j-];
    }
    return res;
}
int main()
{
    ll l,r;
    scanf("%llu%llu",&l,&r);
    printf("%llu\n",solve(r)-solve(l-));
    return ;
}