找出所有的生日素数！

无心人

儒略历和格里历转换太复杂，就从1600年开始一直到9999年

1. #include <gmp.h>
   
2. #include <iostream>
   
3. 
   
4. using namespace std;
   
5. 
   
6. int isprime(unsigned long n)
   
7. {
   
8.   if (n % 2 == 0)
   
9.     return 0;
   
10.   if (n % 3 == 0)
    
11.     return 0;
    
12.   if (n % 5 == 0)
    
13.     return 0;
    
14.   if (n % 7 == 0)
    
15.     return 0;
    
16.   if (n % 11 == 0)
    
17.     return 0;
    
18.   if (n % 13 == 0)
    
19.     return 0;
    
20.   if (n % 17 == 0)
    
21.     return 0;
    
22.   if (n % 19 == 0)
    
23.     return 0;
    
24.   mpz_t z;
    
25.   mpz_init_set_ui(z, n);
    
26.   return mpz_probab_prime_p(z, 4);
    
27. }
    
28. 
    
29. int main( void )
    
30. {
    
31.   unsigned long y, m, d, n;
    
32.   for (y = 1600; y <= 9999; y ++)
    
33.     for (m = 1; m <= 12; m ++)
    
34.         for (d = 1; d <= 31; d ++)
    
35.         {
    
36.             if ((m == 4) || (m == 6) || (m == 9) || (m == 11))
    
37.                 if (d > 30)
    
38.                   continue;
    
39.             if (y % 4 != 0)
    
40.             {
    
41.                 if (d > 28)
    
42.                   continue;
    
43.             }       
    
44.             else
    
45.                 if ((y % 100 != 0) || (y % 400 == 0))
    
46.                 {
    
47.                    if (d > 29)
    
48.                      continue;
    
49.                 }
    
50.                 else
    
51.                   if (d > 28)
    
52.                     continue;
    
53. 
    
54.             n = 10000 * y + 100 * m + d;
    
55.             int r = isprime(n);
    
56.             if (r != 0)
    
57.             {
    
58.               if (r == 2)  
    
59.                 cout << "生日素数" << n << endl;
    
60.               else
    
61.                 cout << "可能生日素数" << n << endl;
    
62.             }
    
63.         }
    
64. }
    

复制代码