GCPC2015 F Divisions (约数个数定理)

 

题意很简单，输入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;
}