问题来源
Timus Online Judge 网站上有这么一道题目:1356.
Something Easier。这道题目的输入是一组 2 到 109 之间整数,对于每个输入的整数,要求用最少个数的素数的和来表示。这道题目的时间限制是 1 秒。
问题解答
我们知道著名的哥德巴赫猜想是:
任何一个充分大的偶数都可以表示为两个素数之和
于是我们有以下的 C 语言程序(1356.c):
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 |
|
上述程序分析如下:
- 根据哥德巴赫猜想,充分大的偶数 n = p + q,这里 p <= q 是素数。我们猜测当 n <= 109 时,p < 104。第 8 行就是定义 p 的最大值。
- 根据素数定理,我们知道 104 以内的素数有 1229 个。第 9 行就是定义程序中要用到的素数的个数。
- 第 17 到 26 行的 getSieve 函数用埃拉托斯特尼筛法筛选出素数。
- 第 28 到 36 行的 getPrimes 函数从筛中取出这些素数。
- 第 38 到 79 行的一系列函数最终是为了 probablyPrime 函数,用于检测素数。请参见我在2010年7月写的随笔:【算法】米勒-拉宾素性检验。
- 第 81 到 84 行的 isPrime 函数调用 probablyPrime 函数来检测素数。
- 第 86 到 93 行的 outEven 函数对大于 2 的偶数验证哥德巴赫猜想,即输出一对素数 p 和 q。
- 第 95 到 110 行是 main 函数。其中:
- 第 103 行处理 n 是素数的情况,直接输出该素数(包括素数 2,所以 outEven 函数处理的偶数肯定大于 2)。
- 第 104 行对大于 2 的偶数输出一对素数(通过调用 outEven 函数,强哥德巴赫猜想)。
- 第 105 行处理大于 5 的奇数能够分解为 2 和另外一个素数的和的情况(注意不要遗漏这个情形!)。
- 第 106 行处理大于 5 的奇数的其他情况,首先输出一个 3,然后调用 outEven 函数处理偶数 n - 3 (弱哥德巴赫猜想)。
上述程序在 Timus Online Judge 网站的运行时间是 0.015 秒。
版权声明:本文为博主http://www.zuiniusn.com 原创文章,未经博主允许不得转载。