LightOJ - 1341 Aladdin and the Flying Carpet 唯一分解定理

https://vjudge.net/contest/307652#problem/C
题意
t组数据（t<=4000）,每组两个数a,b，1 ≤ b ≤ a ≤ 1e12。a代表矩形的面积，b代表矩形的可能的最短的一条边。求矩形的长宽有多少种可能的情形，矩形不能是正方形。
题解
首先不考虑b的限制，这题就是求a的因子数除以2。题目要求b是最小可能边长，当b>=sqrt(a)时，没有符合条件的边（=是正方形），暴力枚举[1,n]，求出a的因子减掉就可以了。
但我咋看这题4000*1e6也是超时的啊。。。
代码

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxn=(1e6+10);
bool isprime[maxn];
ll su[maxn];
int cnt;
void init()
{
    cnt=0;
    isprime[0]=isprime[1]=1;
    for(int i=2; i<maxn; i++)
    {
        if(!isprime[i])
            su[cnt++]=i;
        for(int j=0; j<cnt&&su[j]*i<maxn; j++)
        {
            isprime[su[j]*i]=1;
            if(i%su[j]==0) break;
        }
    }
}
int fj(ll x)
{
    ll tmp=sqrt(x);
    int fact=1;
    for(int i=0;i<cnt&&su[i] <= sqrt(x);i++)//更迷的是把这里的sqrt(x)换成tmp竟然就超时...这题可真的是..,
    {
        int num=0;
        while(x%su[i]==0)
        {
            num++;
            x/=su[i];
        }
        fact*=(num+1);
    }
    if(x>1) fact*=2;
    return fact;
}
int main()
{
    int t;
    int cmt=0;
    ll a,b;
    int ans;
    init();
    scanf("%d",&t);
    while(t--)
    {
        cmt++;
        scanf("%lld%lld",&a,&b);
        if(b>sqrt(a)) ans=0;
        else
        {
            int fact=fj(a)/2;
            int num=0;
            for(ll i=1; i<b; i++)
            {
                if(a%i==0) num++;
            }
            ans=fact-num;
        }
        printf("Case %d: %d\n",cmt,ans);
    }
    return 0;
}