UVa 10837 (欧拉函数 搜索) A Research Problem

发现自己搜索真的很弱，也许做题太少了吧。代码大部分是参考别人的，=_=||

题意：

给出一个phi(n)，求最小的n

分析：

回顾一下欧拉函数的公式：，注意这里的Pi是互不相同的素数，所以后面搜索的时候要进行标记。

先找出所有的素数p，满足(p - 1)整除题目中所给的phi(n)

然后暴搜。。

素数打表打到1e4就够了，如果最后剩下一个大素数单独进行判断。

 1 #include <cstdio>
 2 #include <cmath>
 3 #include <cstring>
 4 #include <algorithm>
 5 using namespace std;
 6
 7 const int maxn = 10000;
 8 const int INF = 200000000;
 9
10 int phi_n, ans;
11
12 bool vis[maxn + 10];
13 int prime[1500], pcnt = 0;
14 int fac[1500], tot;
15
16 void prime_table()
17 {
18     int m = sqrt(maxn + 0.5);
19     for(int i = 2; i <= m; i++) if(!vis[i])
20         for(int j = i * i; j <= maxn; j += i) vis[j] = true;
21     for(int i = 2; i <= maxn; i++) if(!vis[i]) prime[pcnt++] = i;
22 }
23
24 void factor(int n)
25 {
26     tot = 0;
27     for(int i = 0; i < pcnt && (prime[i]-1) * (prime[i]-1) <= n; i++) if(n % (prime[i] - 1) == 0)
28         fac[tot++] = prime[i];
29 }
30
31 int judge(int n)
32 {
33     if(n == 1) return 1;
34     n++;
35     //判断剩余的phin中加1后是否为素数
36     for(int i = 0; i < pcnt && prime[i] * prime[i] <= n; i++) if(n % prime[i] == 0)
37         return -1;
38     //如果为素数的话，判断是否标记过
39     for(int i = 0; i < tot; i++) if(vis[i] && n == fac[i])
40         return -1;
41     return n;
42 }
43
44 void dfs(int n, int phin, int d)
45 {
46     if(d == tot)
47     {
48         int t = judge(phin);
49         if(t > 0) ans = min(ans, n * t);
50         return;
51     }
52
53     dfs(n, phin, d+1);
54     if(phin % (fac[d] - 1) == 0)
55     {
56         vis[d] = true;
57         n *= fac[d];
58         phin /= (fac[d] - 1);
59         for(;;)
60         {
61             dfs(n, phin, d+1);
62             if(phin % fac[d] != 0)
63                 return;
64             phin /= fac[d]; n *= fac[d];
65         }
66     }
67     vis[d] = false;
68 }
69
70 int main()
71 {
72     freopen("in.txt", "r", stdin);
73
74     int kase = 0;
75     prime_table();
76     while(scanf("%d", &phi_n) == 1 && phi_n)
77     {
78         ans = INF;
79         memset(vis, false, sizeof(vis));
80         factor(phi_n);
81         dfs(1, phi_n, 0);
82         printf("Case %d: %d %d\n", ++kase, phi_n, ans);
83     }
84
85     return 0;
86 }

代码君