基础数学1003 LightOJ 1341

题意:
有一个矩形(不能使正方形)的魔毯,
给你面积和能接受的最小的边长(最短的边不能小于这个的大小
问你有多少种魔毯可以被接受
思路:
因为长*宽=面积
所以长和宽必然是面积的质因子分解里的质因子组成的
所以我们可以暴力枚举最短的边的质因子
因为最短的边是<=sqrt(1e12)<=1e6的
因为是质因子组成的,即使是MinSide很接近sqrt(Area)
我们的时间也是比O(1e6)小很多的

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
#define lowbit(x) (x&(-x))
typedef long long LL;
const int maxn = 100005;
const int inf=(1<<28)-1;
#define maxp 2000005
bool notprime[maxp];
int primes[maxp];
void get_prime()
{
    notprime[1]=true;
    for(int i=2;i<maxp;++i)
    if(!notprime[i])
    {
        primes[++primes[0]]=i;
        for(LL j=(LL)i*i;j<maxp;j+=i)
        notprime[j]=true;
    }
}
int Fac[maxn],Num[maxn],tot;
void GetFac(LL Area)
{
    tot=0;
    for(int i=1;i<=primes[0]&&primes[i]*primes[i]<=Area;++i)
    if(Area%primes[i]==0)
    {
        int Cnt=1;
        Fac[++tot]=primes[i];
        Area/=primes[i];
        while(Area%primes[i]==0)
            Area/=primes[i],Cnt++;
        Num[tot]=Cnt;
    }
    if(Area>1)
    {
        Fac[++tot]=Area;
        Num[tot]=1;
    }
}
LL MaxSide,MinSide;
LL Ans=0;
void dfs(int pos,LL tmp)
{
    if(pos==tot+1)
    {
        if(tmp>=MinSide) Ans++;
        return ;
    }
    dfs(pos+1,tmp);
    for(int i=1;i<=Num[pos];++i)
    {
        tmp*=Fac[pos];
        if(tmp>=MaxSide) return ;
        dfs(pos+1,tmp);
    }
}
int main()
{
    get_prime();
    int Case=0,T;
    scanf("%d",&T);
    while(T--)
    {
        LL Area;
        scanf("%lld%lld",&Area,&MinSide);
        MaxSide=sqrt(Area);
        if(MaxSide*MaxSide!=Area) MaxSide++;
        if(MinSide>=MaxSide)
        {
            printf("Case %d: 0\n",++Case);
            continue;
        }
        GetFac(Area);
        Ans=0;
        dfs(1,1);
        printf("Case %d: %lld\n",++Case,Ans);
    }
    return 0;
}