hdu4542 求反素数

题意Problem Description
　　“今有物不知其数，三三数之有二，五五数之有三，七七数之有二，问物几何？”

　　这个简单的谜题就是中国剩余定理的来历。

　　在艰难地弄懂了这个定理之后，小明开始设计一些复杂的同余方程组X mod ai = bi 来调戏别人，结果是必然的，都失败了。

　　可是在这个过程中，小明发现有时并不一定要把ai和bi告诉你。他只需要告诉你，ai在区间 [1, X] 范围内每个值取一次时，有K个ai使bi等于0，或有K个ai使bi不等于0，最小的X就可以求出来了。

　　你来试试看吧！

Input
输入第一行为T，表示有T组测试数据。
每组数据包含两个整数Type和K，表示小明给出的条件。Type为0表示“有K个ai使bi等于0”，为1表示“有K个ai使bi不等于0”。

[Technical Specification]

1. 1 <= T <= 477
2. 1 <= K <= 47777, Type = 0 | 1

Output
对每组数据，先输出为第几组数据，如果没有这样的数，输出“Illegal”，否则输出满足条件的最小的X，如果答案大于2^62, 则输出“INF”。

Sample Input
3
0 3
1 3
0 10

Sample Output
Case 1: 4
Case 2: 5
Case 3: 48

思路：当type为1时，可以打表，ip数组表示k对应的最小数ip[k]；当type为0时，即求反素数，用一个dfs求出；
代码：

#include<cstdio>  
#include<cstring>  
#include<cmath>  
#include<algorithm>  
using namespace std;  
typedef __int64 LL;
#define ns 50100  
int top=15;  
const LL INF=((LL)1<<62)+1;  
__int64  ans,cnt;  
int m,k,type;  
__int64 ip[ns+10];
int p[100]={1,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};  
void init()
{
    int i,j;
    for(i=0;i<ns;i++)
        ip[i]=i;  //表示i最多的质数为i;
    for(i=1;i<ns;i++)
    {
        for(j=i;j<ns;j=j+i)
            ip[j]--;  //因子有i就相应减一;
        if(!ip[ip[i]])//ip[i]表示与i对应的质数数目tt,ip[tt]表示质数数目为tt对应最小的数;
            ip[ip[i]]=i;
        ip[i]=0;//与i对应的质数数目不能为i；
    }
}
void dfs(__int64 t,int sum,int pos,int limit)  //可以去看求反素数的模板
{
    if(sum==k&&ans>t)
        ans=t;
    if(sum>=k||pos>top)
        return;
    for(int i=1;i<=limit;i++)
    {
        if(t>ans/p[pos]||sum*(i+1)>k)
            break;
        t=t*p[pos];
        if(k%(sum*(i+1))==0)
        dfs(t,sum*(i+1),pos+1,i);
    }
}
int main()
{
    int n,s=1;
    init();
    scanf("%d",&n);
    while(n--)
    {
        scanf("%d%d",&type,&k);
        if(type)
        {
            ans=ip[k];
        }
        else
        {
            ans=INF;
            dfs(1,1,1,62);
        }
        printf("Case %d: ",s++);
        if(ans>=INF)
            printf("INF\n");
        else if(ans==0)
            printf("Illegal\n");
        else
            printf("%I64d\n",ans);
    }
}