实验三

include <stdio.h>

include <stdlib.h>

include <gmp.h>

define MAX_PRIME 10000

int main()
{
mpz_t prime, n;
int i;

// 初始化大数库
mpz_init(prime);
mpz_init(n);
mpz_set_ui(prime, 2);

// 计算10000以内素数
for (i=2; i<=MAX_PRIME; i++) {
    int is_prime = 1;
    int j;
    for (j=2; j*j<=i; j++) {
        if (i % j == 0) {
            is_prime = 0;
            break;
        }
    }
    if (is_prime) {
        mpz_set_ui(n, i);
        mpz_mul(prime, prime, n);
    }
}

// 输出结果
gmp_printf("10000以内素数的乘积为: %Zd\n", prime);

// 释放内存
mpz_clear(n);
mpz_clear(prime);

return 0;

}