在1000如此小的素数判断,在不考虑效率的情况下可以利用素数的定义来判断
//========================================================
由于卡米歇尔数的存在,导致 费马小定理 无法判断一个数是否是素数。
费马小定理: 设p是素数, a是任意整数且 a!三0( mod p ), 则
a^(p-1) 三 1(mod p)
//========================================================
卡米歇尔数:它是合数, 当 1<=a<=n, 都有 a^n 三 a(mod n)
//========================================================
卡米歇尔数的考塞特判别法: 设n是合数,则n是卡米歇尔数当且仅当它是奇数,且整除n的每个素数p满足下述两个条件:
1)p^2 不整除 n
2)p-1 整除 n-1
//========================================================
如果要判断相当大的素数最好使用 合数的拉宾-米勒测试定理
合数的拉宾-米勒测试定理: 设n是奇素数, 记 n-1 = 2^k * q , q 是奇数, 对不被n整除的某个a, 如果下述两个条件都成立,则n是合数.
a) a^q !三 1(mod n);
b) 对所有 i = 0, 1, 2, ...., k-1, a^((2^i)*q) !三 -1(mod n);
//========================================================
这里给出了合数的拉宾-米勒测试定理 a 的取值:
if n < 1,373,653, it is enough to test a = 2 and 3.
if n < 9,080,191, it is enough to test a = 31 and 73.
if n < 4,759,123,141, it is enough to test a = 2, 7, and 61.
if n < 2,152,302,898,747, it is enough to test a = 2, 3, 5, 7, and 11.
//========================================================