#include #include #include #include #include #include using namespace std; typedef __int64 LL; const LL NUM=20; const int maxn=1e6+10; int prime[maxn],c[maxn],tot; LL t,f[1000]; LL mul_mod(LL a,LL b,LL n) { a=a%n; b=b%n; LL s=0; while(b) { if(b&1) s=(s+a)%n; a=(a<<1)%n; b=b>>1; } return s; } LL pow_mod(LL a,LL b,LL n) { a=a%n; LL s=1; while(b) { if(b&1) s=mul_mod(s,a,n); a=mul_mod(a,a,n); b=b>>1; } return s; } bool check(LL a,LL n,LL r,LL s) { LL ans,p,i; ans=pow_mod(a,r,n); p=ans; for(i=1;i<=s;i++) { ans=mul_mod(ans,ans,n); if(ans==1&&p!=1&&p!=n-1)return true; p=ans; } if(ans!=1)return true; return false; } bool Miller_Rabin(LL n) { if(n<2)return false; if(n==2)return true; if(!(n&1))return false; LL i,r,s,a; r=n-1;s=0; while(!(r&1)){r=r>>1;s++;} for(i=0;ix)p=y-x; else p=x-y; d=gcd(p,n); if(d!=1&&d!=n)return d; if(i==k) { y=x; k+=k; } } } void find(LL n) { if(Miller_Rabin(n)) { f[t++]=n; return; } LL p=n; while(p>=n)p=Pollard_rho(p,rand()%(n-1)+1); find(p); find(n/p); } void init() { memset(c,false,sizeof(c)); int i,j; tot=0; for(i=2;i<=1e6;i++) { if(!c[i])prime[tot++]=i; for(j=0;j1e6)break; c[i*prime[j]]=true; if(i%prime[j]==0)break; } } } bool judge(LL n) { int i,j,k; for(i=0;i>n) { if(n==-1)break; if(n==1){cout<<0<=1)ans--; if(k!=1) { if(Miller_Rabin(k))ans++; else { i=(LL)sqrt(k+0.5); if(i*i==k&&Miller_Rabin(i)||judge(k))ans++; else ans+=2; } } if(ans<2)cout<