hdu4542 小明系列故事——未知剩余系（反素数打表待改！！不会！！）

#include <stdio.h>
#include <algorithm>
#include <string.h>

using namespace std;

typedef long long ll;
const int N = 50005;//互质数的个数不超过k最大值
const ll INF = (((ll)1)<<62)+1;

int p[16] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53};

ll ans, n;//ans表示满足约数为n的反素数，通过多次贪心求的最小值
int num
;//f
代表非约数个数为n的最小x

void table()
{
/*  num[1] = num[2] = 0;
for(int i = 3; i <= N; i++)
{
num[i] = i;
for(int j = 1; j <= (i/2+1); j++)
{
if(i%j == 0) num[i]--;
}
num[i]--;
}
memset(f, 0, sizeof(f));
for(int i = 1; i <= N; i++)
{
if(f[num[i]] == 0) f[num[i]] = i;
}*/
for(int i=1;i<N;i++) num[i] = i;
for(int i=1;i<N;i++)
{
for(int j=i;j<N;j+=i) num[j]--;
if(!num[num[i]]) num[num[i]] = i;
num[i] = 0;
}
}

void dfs(ll step, ll cur, ll cnt, ll Maxnum)//cur表示当前因子个数所对应的反素数，cnt表示当前因子个数
{
if(cnt > n) return;
if(cnt==n && cur<ans) ans = cur;//更新最小反素数
ll tmp = cur;//在原先反素数的基础上乘第step个素数
for(int i = 1; i <= Maxnum; i++)//枚举第step个素数的个数
{
if(ans / p[step] < tmp || cnt*(i+1) > n) break;
tmp *= p[step];
/* tmp*=p[step];
if(tmp>ans || (cnt*(i+1))>n) break;//若乘上这个质因子超过了已知最大反素数，不满足求最小反素数的条件，返回*/
if(n%(cnt*(i+1))==0)
dfs(step+1, tmp, cnt*(i+1), i);//下一个质因子肯定个数要比这个少，因子个数累乘
}
}

int main()
{
//  freopen("in.txt", "r", stdin);
int t, type, Case = 1;
table();
/*  for(int i = 1; i <= 20; i++)
{
printf("%d %lld\n", i, num[i]);
}*/
scanf("%d", &t);
while(t--)
{
//  ans = 0;
scanf("%d%lld", &type, &n);
if(type == 1)
{
ans = num
;
}
else if(type == 0)
{
ans = INF;
dfs(0, 1, 1, 62);
}
printf("Case %d: ", Case++);
if(ans == 0) printf("Illegal\n");
else if(ans >= INF) printf("INF\n");
else printf("%lld\n", ans);
}
return 0;
}