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

题干：

给一对数字 a，b 。
a是一个长方形的面积，问有多少种整数的边的组合可以组成面积为a的长方形，要求最短的边不得小于b。
因为给定是长方形的面积，所以就是求所有能整除a的情况。
说白了，就是求a的所有公约数，其最小值不能小于b。

思路：

根据唯一分解定理（一个大于1的合数一定可以写成若干个素数的乘积）
N = p1^a1 * p2^a2 * p3^a3 … * pn^an
根据性质N的因子的个数=(1+a1) * (1+a2) * (1+a3)…(1+an)
因为这样得出的因子一定为素数，所有我们对满足积为n的因子遍历，减去比b小的情况。

#include<cstring> 
#include<cstdio>
#include<iostream>
#include<queue> 
#include <cmath>
#include <map>
#include <algorithm>
using namespace std;
 
const int mx=1e6+10;
int f[mx],p[mx],num;
long long a,b,sum,tmp;

void f1()
{
	for(int i=2;i<mx;i++){
		if(!f[i])  p[num++]=i;
		for(int j=0;j<num&&i*p[j]<mx;j++){
			f[i*p[j]]=1;
			if(!(i%p[j]))  break;
		}
	}
}
void cal()
{
	for(int i=0;i<num&&p[i]<=sqrt(tmp);i++){
		int tt=0;
		while(tmp%p[i]==0)
		{
			tmp/=p[i];
			tt++;
		}
		sum*=(1+tt);
		//printf("%d\n",sum);
	}
	if(tmp>1)
	  sum*=2;
}
int main() 
{
	int t;
	f1();
	f[1]=1;
	scanf("%d",&t);
	for(int i=1;i<=t;i++){
		scanf("%lld%lld",&a,&b);
		if(b>sqrt(a))
		  printf("Case %d: 0\n",i);
		else
		{
			tmp=a;
			sum=1;
			cal();
			//printf("%d\n",sum);
			sum/=2;
			for(int j=1;j<b;j++){                                                                                                   
				if(a%j==0)
				  sum--;
			}
			printf("Case %d: %lld\n",i,sum);
		}
	}
    return 0;
}