uestc 1725 吴神数

// 筛选法
// 先求出 sqrt(1<<31)内的素数
// 然后筛选出符合要求的数
// 详情见代码注释
// #pragma comment(linker, "/STACK:1024000000,1024000000")
#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
#include <math.h>
#include <stdio.h>
#include <string.h>
using namespace std;
#define MOD 1000000007
#define maxn 1000010
#define maxm 48010
#define LL  long long
int pr[maxm];
int p;
void getP(){
  int i,j;
  for(i=;i<maxm;i+=)
    pr[i]=;
  for(i=;i*i<maxm;i+=)
    if(!pr[i])
     for(j=i*i;j<maxm;j+=i)
        pr[j]=;
  pr[p++]=; // printf("%d ",p);
  for(i=;i<maxm;i+=)
    if(!pr[i]) pr[p++]=i;//,printf("%d ",i);
}
int fac[maxn],f[maxn],lt[maxn];
int main(){
    getP();
    int T;
    int A,B;
    LL tp;
    scanf("%d",&T);
    int i;
    LL j,k;
    while(T--){
        scanf("%d %d",&A,&B);
        LL len=B-A;
        for(i=;i<=len;i++) fac[i]=,f[i]=false,lt[i]=A+i;
        tp=;
        for(i=;i<p;i++){ // 素数从 3 开始 因为偶数时不可能符合要求的
             tp=tp*pr[i]*pr[i];
             if(tp>B) break;
             for(j=(A-+pr[i])/pr[i];j*pr[i]<=B;j++)
             {
                 k=j*pr[i];
                 if(j%pr[i]==) f[k-A]=true;// 去掉 平方因子
                 if((j*pr[i]-)%(pr[i]-)!=) f[k-A]=true; // 条件3
                 fac[k-A]++; // 条件 2
                 lt[k-A]=lt[k-A]/pr[i];可能出现比sqrt(1<<31)大的素因子，所以需要知道
             }
             tp=;
        }
        int ans=;
        for(i=;i<=len;i++)// 这里就写的比较繁杂了 就是各种情况讨论
            if(!f[i]){
                if(fac[i]==&&lt[i]>){
                    tp=i+A;
                    if((tp-)%(lt[i]-)==)
                        ans++;
                }else if(fac[i]>){
                    tp=i+A;
                    if(lt[i]==) ans++;
                    else
                     if((tp-)%(lt[i]-)==)
                        ans++;
                }
            }
        printf("%d\n",ans);
    }
}