Lightoj 1341(质因子分解)

题意：给出一个长方形的面积，然后求长和宽都不能小于m的方案数。

思路：先对这个面积进行质因子分解，得到这个数的所有的因子数，然后因为是不能是正方形，所以他能组成的方案书为因子数/2，要得到大于等于m的话，直接减去小于m的数。

注意刚开始要判断m * m 是否是小于等于n，这样的话m就小于等于1e6了，不会超时。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
#define fi first
#define se second
#define INF 0x3f3f3f3f
#define clr(x,y) memset(x,y,sizeof x)
#define PI acos(-1.0)


const int maxn = 2e6 + 20;
const int Mod = 1e9 + 7;

int prime[maxn],is_[maxn];
ll n,m;int len;
void Init()
{
    len = 0;clr(is_,true);
    for(int i = 2; i < maxn; i ++)
    {
        if(is_[i])
        {
            prime[len ++] = i;for(int j = i * 2; j < maxn; j += i)is_[j] = false;
        }
    }
}
int main()
{
    Init();
    int Tcase;scanf("%d",&Tcase);
    for(int ii = 1; ii <= Tcase; ii ++)
    {
        scanf("%lld%lld",&n,&m);
        if(m * m > n){printf("Case %d: 0\n",ii);continue;}
        ll ans = 1,t = n;
        for(int i = 0; i < len && prime[i] * prime[i] <= n; i ++)
        {
            if(n % prime[i] == 0)
            {
                int cnt = 0;while(n % prime[i] == 0){n /= prime[i];cnt ++;}ans = ans * (cnt + 1);
            }
        }
        if(n > 1)ans <<= 1;
//        cout << ans << endl;
        ans /= 2;
        for(int i = 1; i < m; i ++)if(t % i == 0)ans --;
        printf("Case %d: %lld\n",ii,ans);
    }
    return 0;
}