Description
给定两个大整数a和b,请你求a有多少个大于等于b的因子。
Input
第一行输入一个整数t(t<=4000),接下来的t行每行两个整数a,b(1<=b<=a<= 1 0 12 10^{12} 1012) 。
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);
}
}
*/