费马二平方定理(pascal编程)

除了2这个特殊的素数外，所有的素数都可以分两类：第一类是被4除余1的素数，如5，13，17，29，37，41；第二类则是被4除余3的素数如3，7，11，19，23，31.很奇怪的是：第一类素数都能表示为两个整数的平方和，第二类都不能。如：

　　5=1*1+2*2;

　　13=2*2+3*3;

　　17=1*1+4*4;

29=2*2+5*5;

这就是著名的费马二平方定理。有趣的是：上述有的等式右侧的数又恰恰是两个素数，如上面的13和29两个数的等号右侧就都是素数，我们把这样的素数取名为“费马二平方素数”。即：如果一个素数F能够表示成两个素数的平方和形式F=X*X+Y*Y，其中X、Y都是素数，那么它就是费马二平方素数。

编程打印出N以内的所有费马二平方素数（N<=2*10⁹）。

输出格式：no.：F=X*X+Y*Y

注意，由于F为质数(F>2)，所以必定为奇数，所以X*X和Y*Y中必定一奇一偶，由于X和Y也是质数,所以X和Y必定有一个是2，假定X=2，所以问题转化为F=4+Y*Y，F和Y为质数。程序如下：