GMP的一个Demo小程序

最新推荐文章于 2024-09-14 14:51:04 发布

原创最新推荐文章于 2024-09-14 14:51:04 发布 · 2.7k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#ui #c #input #文档 #file #测试

本文展示了一个使用GMP数学库的小程序实例，演示了如何进行大数运算，包括计算阶乘、最大公约数、斐波那契数、二项式系数等。此外还实现了浮点数的高精度运算及用户输入的分数加法等功能。

今天上午一口气看完了GMP的文档，边看边整出了一个小程序，贴出来，日后哪天重新编译GMP了，还可以以此作为测试代码

该代码展示了大数的阶乘，大数的二项式数，两大数的最大公约数，较大的fibonacci数，浮点数的高精度乘方和加法运算，以及处理用户输入的任意两分数的加法运算，用户输入的任意大小的整数的下一个素数的计算。

#include<stdio.h> #include<gmp.h> int main (void) { FILE *f=fopen("a.txt","w"); int num=50; mpz_t n,m,c,r; mpz_init_set_str (n, "13703703579", 10); mpz_init_set_str (m, "29012345415", 10); mpz_inits(c,r,'/0'); mpz_gcd(c,n,m); mpz_fac_ui(r,num); mpz_out_str(f,36,r); close(f); gmp_printf ("gcd(%Zd,%Zd) = %Zd/n/n%d! =%Zd:/n in base 36:/t",n,m,c,num,r); mpz_out_str(stdout,36,r); mpz_fib_ui(r,num); gmp_printf("/n/nthe %dth fibonacci number is:/n %Zd/n",num,r); mpz_bin_ui(c,n,num); gmp_printf("/nBinomial number C(%Zd,%d) is:/n%Zd/n",n,num,c); mpz_pow_ui(c,n,num); gmp_printf("/n%Zd^%d =:/n%Zd/n",n,num,c); mpf_set_default_prec(500); mpf_t pi,e,xx,yy; mpf_init_set_str(pi,"3.14159265358979323846264338327950288419716939937510582097494459230781",10); mpf_init_set_str(e ,"2.71828182845904523536028747135266249775724709369995957496696762772407",10); mpf_inits(xx,yy,'/0'); mpf_pow_ui(xx,pi,4);mpf_pow_ui(yy,pi,5); mpf_add(xx,xx,yy); mpf_pow_ui(yy,e,6); gmp_printf("/npi^4+pi^5=%.50Ff/ne^6/t =%.50Ff/n/ninput a number to compute its next prime,",xx,yy); gmp_printf("arbitraryly large:/n"); gmp_scanf("%Zd",c); mpz_nextprime(r,c); gmp_printf("after %Zd,closely followed a prime:/n%Zd/n/n",c,r); mpq_t x,y,z; mpq_inits(x,y,z,'/0'); gmp_printf("you input two fractions with a plus sign(2/4+1/3),I compute their sum:/n"); gmp_scanf ("%Qd + %Qd",x,y); mpq_canonicalize(x); mpq_canonicalize(y); mpq_add(z,x,y); gmp_printf("/nactually,after reduction,your equation is:/n%Qd+%Qd = %Qd/nGoodbye!/n",x,y,z); mpz_clears(m,n,c,r,'/0'); mpq_clears(x,y,z,'/0'); mpf_clears(xx,yy,'/0'); return 0; }