每日一题（01.16）

Description

给定两个大整数a和b，请你求a有多少个大于等于b的因子。

Input

第一行输入一个整数t(t<=4000)，接下来的t行每行两个整数a,b(1<=b<=a<=$10^{12}$) 。

Output

a大于等于b的因子数。格式详见Sample。

Sample Input

2
10 2
12 2

Sample Output

Case 1: 1
Case 2: 2

Reference Code

#include <cstdio>
#include <cstring>
#define ll long long
using namespace std;
int const maxn=1e6+3;
bool np[maxn];
int p[maxn];
int res,primes=0;
int main(){
    memset(np,false,sizeof(np));
    for(int i=2;i*i<=maxn;i++) if(!np[i])
        for(int j=i*i;j<=maxn;j+=i)
            np[j]=true;
    for(int i=2;i<=maxn;i++)
        if(!np[i]) p[primes++]=i;
    int t;
    scanf("%d",&t);
    ll a,b,buf;
    for(int cnt=1;cnt<=t;cnt++){
        scanf("%lld%lld",&a,&b);
        if(a<b*b) printf("Case %d: 0\n",cnt);
        else{
            res=1,buf=a;
            for(int i=0;i<primes&&p[i]*p[i]<=buf;i++){
                int sum=1;
                while(buf%p[i]==0){
                    sum++;
                    buf/=p[i];
                }
                res*=sum;
            }
            if(buf>1) res<<=1;
            res>>=1;
            for(int i=1;i<b;i++) if(!(a%i)) res--;
            printf("Case %d: %d\n",cnt,res);
        }
    }
}

/*暴力法，会超时
#include <cstdio>
#include <cmath>
#define ll long long
int main(){
    int t,res;
    ll a,b,high;
    scanf("%d",&t);
    for(int cnt=1;cnt<=t;cnt++){
        res=0;
        scanf("%lld%lld",&a,&b);
        high=sqrt(a);
        for (ll i=b;i<=high;i++){
            if (!(a%i)) res+=1;
        }
        printf("Case %d: %d\n",cnt,res);
    }
}
*/