UVa Problem Solution: 10006 - Carmichael Numbers

Compute the remainder with the following equation:
    (x*y)%n = (x%n)*(y%n)%n

Code:

1. /***************************************************************************
2.  *   Copyright (C) 2008 by Liu Kaipeng                                     *
3.  *   LiuKaipeng at gmail dot com                                           *
4.  ***************************************************************************/
6. /* @JUDGE_ID 00000 10006 C++ "Carmichael Numbers" */
8. #include <algorithm>
9. #include <cstdio>
10. #include <cstring>
11. #include <deque>
12. #include <fstream>
13. #include <iostream>
14. #include <list>
15. #include <map>
16. #include <queue>
17. #include <set>
18. #include <stack>
19. #include <string>
20. #include <vector>
22. using namespace std;
23.      
24. int const max_num = 65001;
25. bool is_prime[max_num];
27. void gen_primes()
28. {
29.   for (int i = 0; i < max_num; ++i)
30.     is_prime[i] = true;
31.   for (int i = 2; i < max_num; ++i) {
32.     for (; i < max_num && !is_prime[i]; ++i) {}
33.     for (int j = 2; i*j < max_num; ++j)
34.       is_prime[i*j] = false;
35.   } 
36. }
38. /*
39.  * Calculate a^n % n with a divide and conquer method
40.  */
41. int power_mod(int a, int p, int n)
42. {
43.   if (p == 1) return a;
44.   unsigned int r = power_mod(a, p/2, n);
45.   if (p % 2 == 0) {
46.     return r * r % n;
47.   } else {
48.     r = r * r % n;
49.     return r * a % n;
50.   }
51. }
53. bool is_carmichael(int n)
54. {
55.   if (is_prime[n]) return false;
56.   for (int a = 2; a < n; ++a) {
57.     if (power_mod(a, n, n) != a) return false;
58.   }
59.   return true;
60. }
61.           
62. int main(int argc, char *argv[])
63. {
64. #ifndef ONLINE_JUDGE
65.   freopen((string(argv[0]) + ".in").c_str(), "r", stdin);
66.   freopen((string(argv[0]) + ".out").c_str(), "w", stdout);
67. #endif
68.   
69.   gen_primes();
70.   for (int n; cin >> n && n != 0; )
71.     if (is_carmichael(n)) 
72.       cout << "The number " << n << " is a Carmichael number./n";
73.     else
74.       cout << n << " is normal./n";
76.   return 0;
77. }
78.