HDU4542 未知剩余系 反素数

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4542

 

题目大意：两种操作：type为0是让求出约数为k的最小的整数n；type为1是让求有k个数不是n的约数的最小的n。

 

 

分析：数据范围为0到47777，对于type为0的情况，就是找出约数为k的最小的反素数；对于type为1，打表就行。

 

实现代码如下：

#include<stdio.h>
typedef long long ll;
const ll INF=((ll)1<<62)+1;
const int maxn=48000;
const int prime[16]= {1,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};
ll bestnum;
int d[maxn],n;
void init()
{
    for(int i=1;i<maxn;i++) d[i] = i;
    for(int i=1;i<maxn;i++)
    {
        for(int j=i;j<maxn;j+=i) d[j]--;
        if(!d[d[i]]) d[d[i]] = i;
        d[i] = 0;
    }
}
void getantiprime(ll num, int k,int sum,int limit)
{//num：当前枚举到的数，k：枚举到的第k大的质因子；sum：该数的约数个数；limit :质因数上限；
    if(sum>n) return ;
    if(sum==n && bestnum > num)
        bestnum = num; //如果约数个数相同，将最优解更新为较小的数；
    for(int i=1; i<=limit; i++) //开始枚举每个质因子的个数；
    {
        if(bestnum/prime[k] < num ||sum*(i+1)>n)  break;
        num *= prime[k]; //累乘到当前数；
        if( n% (sum*(i+1) )==0)
          getantiprime(num, k+1, sum*(i+1),i); //继续下一步搜索；
    }
}
int main()
{
    init();
    int t,T=1,type;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d", &type,&n);
        if(type) bestnum=d[n];
        else
        {
            bestnum=INF;
            getantiprime(1,1,1,62);
        }
        printf("Case %d: ", T++);
        if(bestnum==0) printf("Illegal\n");
        else if(bestnum>=INF) printf("INF\n");
        else printf("%lld\n",bestnum);
    }
    return 0;
}