一个开源的高精度运算库——gmp

全称是GNU Multiple Precision Arithmetic Library，即GNU高精度算术运算库，官方网站是： http://gmplib.org/
它的功能非常强大，接口很简单，文档详尽，有C风格的接口也有C++的精心封装后的接口，其中不但有普通的整数、实数、浮点数的高精度运算，还有随机数生成，尤其是提供了非常完备的数论中的运算接口，比如Miller-Rabin素数测试算法，大素数生成，欧几里德算法，求域中元素的逆，Jacobi符号，legendre符号等。
它本身提供了很多例子程序，学习过程非常快，很容易将它们集成到自己的代码中去。
下面是最近做一个加密算法的时候使用gmp提供的大素数生成的一段程序，功能是找到给定数量的不小于2^127的素数，其逻辑就像打印一个字符串那样简单：
#include <stdio.h>
#include <stdlib.h>
#include <gmp.h>
 
int main(int argc, char **argv) {
    mpz_t begin, m1, m2;
    int count;
 
    /* set begin to 2^127 */
    mpz_init_set_str(begin, "170141183460469231731687303715884105728", 0);
 
    count = (argc==1)?10:atoi(argv[1]);
 
 
    while(count--) {
        mpz_nextprime(begin, begin);
        gmp_printf("%Zd/n", begin);
    }
 
    mpz_clear(begin);
    return 0;
}
生成的素数如下：
170141183460469231731687303715884105757
170141183460469231731687303715884105773
170141183460469231731687303715884105793
170141183460469231731687303715884105829
170141183460469231731687303715884105851
170141183460469231731687303715884105979
170141183460469231731687303715884106001
170141183460469231731687303715884106031
170141183460469231731687303715884106123
170141183460469231731687303715884106207
170141183460469231731687303715884106213
170141183460469231731687303715884106231
170141183460469231731687303715884106273
170141183460469231731687303715884106303
170141183460469231731687303715884106309