UOJ #188 Sanrd —— min_25筛

题目：http://uoj.ac/problem/188

参考博客：https://www.cnblogs.com/cjoieryl/p/10149748.html

关键是枚举最小质因子...所以构造的 S 与最小质因子有关。

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
int const xn=1e6+;
int pri[xn],cnt,sqr;
ll n,w[xn],m,f[xn];
bool vis[xn];
void init(ll mx)
{
  for(int i=;i<=mx;i++)
    {
      if(!vis[i])pri[++cnt]=i;
      for(int j=;j<=cnt&&(ll)i*pri[j]<=mx;j++)
    {
      vis[i*pri[j]]=;
      if(i%pri[j]==)break;
    }
    }
}
int Id(ll x)
{
  if(x>sqr)return n/x;
  return m-x+;
}
ll S(ll x,int y)
{
  if(pri[y]>x)return ;
  ll ret=;
  for(int i=y;i<=cnt&&(ll)pri[i]*pri[i]<=x;i++)
      for(ll p0=pri[i];p0*pri[i]<=x;p0*=pri[i])
    ret+=S(x/p0,i+)+(ll)pri[i]*(f[Id(x/p0)]-i+);
  return ret;
}
ll solve(ll nw)
{
  m=; n=nw; sqr=sqrt(n);
  for(ll i=,j;i<=n;i=j+)
    {w[++m]=n/i; j=n/(n/i); f[m]=w[m]-;}
  for(int j=;j<=cnt;j++)
    for(int i=;i<=m&&(ll)pri[j]*pri[j]<=w[i];i++)
    f[i]=f[i]-(f[Id(w[i]/pri[j])]-j+);
  return S(n,);
}
int main()
{
  ll L,R; scanf("%lld%lld",&L,&R); init(sqrt(R));
  printf("%lld\n",solve(R)-solve(L-));
  return ;
}