题意很简单,输入n,输出n的约数的个数。
求约数个数两种方式,暴力sqrt(n)或者约数个数定理。
很明显1e18的范围,暴力1e9次循环会T。
那么,只能用约数个数定理,分解质因数。
但是,一般是1e12的范围,预处理1e6的素数筛,或者1e14,处理1e7的素数筛。
这道题怎么处理呢?可以用miller_rabbin大素数判定。
先跑一个1e7的素数筛,那么跑完后我们可以确定的是n的因子没有小于1e7的素数。
跑到最后的n,如果小于等于1e14,那么一定是素数,ans*2即可。
如果n大于1e14,那么用大素数判定,如果是素数,那么ans*2.
如果不是素数,那么n只能由两个大于1e7的素数构成。
如果这两个素数相等,ans*3,如果不等ans*2*2,即ans*4。
#include <bits/stdc++.h>
using namespace std;
typedef unsigned long long ll;
ll prime[10000100];
bool IsPrime[10000100];
int num;
typedef unsigned long long LL;
LL modular_multi(LL x,LL y,LL mo)
{
LL t;
x%=mo;
for(t=0;y;x=(x<<1)%mo,y>>=1)
if (y&1)
t=(t+x)%mo;
return t;
}
LL modular_exp(LL num,LL t,LL mo)
{
LL ret=1,temp=num%mo;
for(;t;t>>=1,temp=modular_multi(temp,temp,mo))
if (t&1)
ret=modular_multi(ret,temp,mo);
return ret;
}
bool miller_rabbin(LL n)
{
if (n==2)return true;
if (n<2||!(n&1))return false;
int t=0;
LL a,x,y,u=n-1;
while((u&1)==0) t++,u>>=1;
for(int i=0;i<100;i++)
{
a=rand()%(n-1)+1;
x=modular_exp(a,u,n);
for(int j=0;j<t;j++)
{
y=modular_multi(x,x,n);
if (y==1&&x!=1&&x!=n-1)
return false;
x=y;
}
if (x!=1)
return false;
}
return true;
}
void init_prime()
{
num=0;
memset(IsPrime,false,sizeof(IsPrime));
for(ll i=2; i<=10000000; i++)
{
if(IsPrime[i]) continue;
prime[num++]=i;
for(ll j=i+i; j<=10000000; j+=i)
IsPrime[j]=true;
}
}
ll solve(ll n)
{
ll ans=1,sum;
for(int i=0; i<num&&prime[i]*prime[i]<=n; i++)
{
ll tmp=prime[i];
sum=0;
if(n%tmp==0)
{
while(n%tmp==0)
{
n/=tmp;
sum++;
}
ans*=(sum+1ll);
}
}
if(n>=1e14)
{
if(miller_rabbin(n)) ans*=2;
else
{
ll w=sqrt(n);
if(w*w==n) ans*=3;
else ans*=4;
}
n=1;
}
if(n>1) ans*=2;
return ans;
}
int main()
{
int cas=1,T;
ll n;
init_prime();
while(~scanf("%llu",&n))
{
printf("%llu\n",solve(n));
}
return 0;
}