题目链接: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;
}