[搜索 剪枝] 算法马拉松25 E BigPrime

减掉所有只有小质数的数的个数就好了
大概就那么几千万个
首先包含大于根号b的质数只能有一个你可以直接算出来
然后再搜 怎么写都能三四秒出
然后你还要些剪枝
我是把大于三次根号b的质数瞎讨论一波 因为最多出现两个

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;

int n,a,b;

const int N=1000005;


int prime[N],num,vst[N];

inline void Pre(int n){
  for (int i=2;i<=n;i++){
    if (!vst[i]) prime[++num]=i;
    for (int j=1;j<=num && (ll)i*prime[j]<=n;j++){
      vst[i*prime[j]]=1;
      if (i%prime[j]==0) break;
    }
  }
}

int ans=0;
int tag[100005];

int Num,_num;

int sum[N];

inline void dfs(int t,int cur){
  if (a<=cur) ans++;
  int l=(a-1)/cur+1,r=b/cur;
  r=min(r,prime[num]);
  if (l<=r) ans+=sum[r]-sum[l-1];
  int tmp=b/cur;
  for (int i=t;i<=_num && prime[i]<=tmp;i++)
    dfs(i,cur*prime[i]);
}

int main(){
  scanf("%d%d%d",&n,&a,&b);
  Pre(n); 
  int B=sqrt(b);
  Num=num;
  for (int i=1;i<=Num;i++)
    if (prime[i]>B)
      ans+=(b/prime[i])-((a-1)/prime[i]);
    else
      num=i;
  for (int i=1;i<=num;i++)
    if ((ll)prime[i]*prime[i]*prime[i]<=b)
      _num=i;
  for (int i=_num+1;i<=num;i++)
    for (int j=i;j<=num;j++)
      ans+=(b/prime[i]/prime[j])-((a-1)/prime[i]/prime[j]);
  for (int i=_num+1;i<=num;i++)
    sum[prime[i]]++;
  for (int i=1;i<=n;i++) sum[i]+=sum[i-1];
  dfs(1,1);
  printf("%d\n",b-a+1-ans);
  return 0;
}