GYM 100796 K.Profact（dfs）

Description

判断一个数字是否可以表示成几个数字阶乘的乘积

Input

第一行一整数$T$表示用例组数，每组用例输入一整数$n(1\le T\le 10^5,1\le n\le 10^{18})$

Output

如果$n$可以表示为几个数字阶乘的乘积则输出$YES$，否则输出$NO$

Sample Input

7
1
2
3
8
12
24
25

Sample Output

YES
YES
NO
YES
YES
YES
NO

Solution

$n$最多以$19!$为因子，直接$dfs$，从大到小搜，一个重要剪枝是如果当前是一个素数的阶乘且 n <script type="math/tex" id="MathJax-Element-9">n</script>可以被其整除则一定除，如果除完无解说明无解

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=100001;
ll f[maxn],n;
int T,res,flag;
int p[]={2,3,5,7,11,13,17,19};
void init()
{
    f[0]=1;
    res=1;
    while(1)
    {
        f[res]=f[res-1]*res;
        if(f[res]>=1e18)break;
        res++;
    }
}
void dfs(ll n,int k)
{
    if(n==1)flag=1;
    if(flag)return ;
    int mark=0;
    for(int i=k;i>1;i--)
        if(n%f[i]==0)
        {
            dfs(n/f[i],i);
            if(!flag)
            {
                for(int j=0;j<8;j++)
                {
                    if(p[j]==i)return ;
                    if(p[j]>i)break;
                }
            }
            if(flag)return ;
        }
}
int main()
{
    init();
    scanf("%d",&T);
    while(T--)
    {
        scanf("%I64d",&n);
        flag=0;
        dfs(n,19);
        printf("%s\n",flag?"YES":"NO");
    }
    return 0;
}