推荐：HugeCalc 超大整数完全精度快速计算器/算法库

forcal

HugeCalc 是一款高精度算法库（同时支持 MBCS + UNICODE 版），适合于大规模科学计算，尤其适用于数论、密码学等领域研究，其核心算法耗费作者十余年的心血。具有占用资源少、效率高、使用便捷、易二次开发、可移植性强、可扩展性好等特点。关键文件 HugeCalc.dll 虽然很小，却提供了公共函数接口 709 个（标准C++接口 473 个；标准C接口 236 个），且其计算速度完全可与大型专业数学工具软件媲美！
    例如：精确计算出5千万的阶乘，在一台工作站上，HugeCalc 耗时654.3秒，内存耗费峰值1105MB；普通程序一般都会因寻址空间不足等原因而中途崩溃。
7 v* Q5 D% v3 _$ P& g+ @( K　　现已提供了如下功能：
2 K1 Q7 G- H: }  l9 @　　⊙ 高精度快速加法
　　⊙ 高精度快速减法
" i2 r% s5 d8 _+ _7 J　　⊙ 高精度快速乘法
　　⊙ 高精度快速除法
　　⊙ 高精度快速同余
. |4 S9 r# A7 X0 {' E8 m) D5 k　　⊙ 高精度快速位运算
　　⊙ 高精度快速乘方
" G5 {5 s) H8 E$ ^' R: D　　⊙ 高精度快速开方
　　⊙ 超大整数快速取对数
　　⊙ 高精度快速求排列
　　⊙ 高精度快速求组合
　　⊙ 高精度快速阶乘、双阶乘、素数阶乘
　　⊙ 高精度快速计算 Fibonacci、Lucas 数列
　　⊙ 高精度快速乘积取模
: i5 O: p) s( ?8 m1 K3 L　　⊙ 高精度快速数论倒数取模运算
4 h+ u- `9 j+ G( F* u  n1 e, X+ G　　⊙ 高精度快速乘方取模（支持负指数）
# ]. }" N2 i% b3 X  n$ p　　⊙ 高精度快速求最大公约数（支持群组运算）
　　⊙ 高精度快速计算扩展最大公约数
　　⊙ 高精度快速求最小公倍数（支持群组运算）
　　⊙ 高精度快速“等幂和”（支持群组运算）
, h6 o8 B+ Z5 a; q/ ?5 a　　⊙ 高精度快速任意进制转换
# Y( u2 B% H/ p0 a1 R1 s7 @" B% R　　⊙ 超大整数素性快速检测
5 s4 Z2 N1 G! ^6 S% u" C　　⊙ 生成随机超大（素）整数、快速生成最邻近素数
4 ?( W6 N+ Z" W　　⊙ 自由指定有效位运算
1 C3 ^. z6 I- F6 {2 m! ?　　⊙ 强大而灵活的输出

    官方网站：http://www.emath.ac.cn/

forcal

HugeCalc的确很强大，计算：

1：1! + 2! + 3! + ... ... + 1000!
2：1! *  2! *  3! * ... ... *  1000!

' e# Y2 z1 S9 `  F! s8 D9 ~" H+ jForcal+HugeCalc代码：
5 {/ }! b, T6 ]. I2 A3 X% B+ i9 l

2. !using["HugeCalc"];3 i5 ]& b6 E\\" o* u* l3 A# K1 N0 Y
   
4. i: main(:i,a,sum,pro)=
   
5. ; j\\" d2 v; O* _( w: k7 L7 s
6.    oo{a=HI(0),sum=HI(0),pro=HI(1)},: Q* Y) V/ v, R/ U' z
   
8.    i=0,(++i<=1000).while{oo{sum.=sum+a.Fac(i), pro.=pro*a}},* V9 r3 R; Z# n2 j4 T- k
   
10.    printff["\r\nsum=\r\n"], sum.Show(),
    
11. ' N* _# p6 b\\" X/ W; g) r, Z# [  _$ x
12.    printff["\r\npro=\r\n"], pro.Show();

!using["HugeCalc"];<font class=\\"jammer\\">3 i5 ]& b6 E\\" o* u* l3 A# K1 N0 Y</font><br /> i: main(:i,a,sum,pro)=<br /> <span style=\\"display:none\\">; j\\" d2 v; O* _( w: k7 L7 s</span>    oo{a=HI(0),sum=HI(0),pro=HI(1)},<font class=\\"jammer\\">: Q* Y) V/ v, R/ U' z</font><br />    i=0,(++i<=1000).while{oo{sum.=sum+a.Fac(i), pro.=pro*a}},<font class=\\"jammer\\">* V9 r3 R; Z# n2 j4 T- k</font><br />    printff["\r\nsum=\r\n"], sum.Show(),<br /> <span style=\\"display:none\\">' N* _# p6 b\\" X/ W; g) r, Z# [  _$ x</span>    printff["\r\npro=\r\n"], pro.Show();

耗时约11～12秒。大家可用其他软件试一下。

结果：
sum=
$ T$ H) ?/ a6 T2 Q8 k: ^402 790 050 531 223 457 680 670 558 292 183 044 930 441 094 551 528 966 419 334 642 920 075 902 370 082 550 849 878 938 011 122 874 840 322 664 395 727 767 319 128 483 846 892 726 398 968 790 560 901 035 018 977 865 543 415 443 376 965 144 616 549 510 538 851 077 440 410 165 819 782 073 475 583 082 079 396 504 111 380 884 916 810 768 467 471 317 168 913 183 948 275 897 490 803 646 698 793 846 190 852 594 417 739 671 786 914 563 472 162 093 180 956 284 720 843 607 311 204 509 397 095 907 530 114 810 381 703 803 860 648 328 656 344 389 023 077 975 647 748 628 886 817 466 804 724 858 587 747 840 936 784 462 436 173 977 429 991 182 753 207 302 348 016 538 494 642 001 301 432 395 279 726 745 538 178 727 708 503 090 870 745 742 249 581 124 437 084 158 031 806 471 602 854 863 650 302 741 825 130 214 652 788 342 621 283 870 175 833 906 572 656 890 104 432 413 095 171 095 723 550 851 047 360 182 700 817 666 447 833 450 097 928 957 386 876 271 765 554 137 407 111 747 666 512 309 980 678 539 886 786 878 262 902 584 539 118 231 519 316 835 523 885 400 137 767 203 426 696 953 126 303 292 828 224 235 349 134 630 280 084 480 064 650 497 764 755 473 824 102 259 453 127 730 077 777 531 063 115 138 764 952 527 597 098 445 380 046 443 997 540 925 880 086 936 317 933 421 768 424 925 084 707 414 812 932 329 166 745 886 234 137 226 923 024 813 461 900 363 225 999 464 524 935 227 947 647 167 947 434 411 637 139 203 517 404 371 214 232 998 239 493 954 562 323 327 076 903 915 931 443 305 369 962 856 388 265 027 346 155 623 332 543 302 943 446 674 383 336 900 484 885 947 609 695 640 139 797 208 838 438 493 998 636 482 572 950 930 777 919 644 406 154 624 084 228 042 985 270 573 774 543 689 565 298 092 301 062 689 807 689 071 120 317 578 510 142 181 892 798 247 199 071 452 022 029 172 707 400 688 758 980 305 702 923 813 291 544 851 703 914 211 657 597 582 497 807 598 793 035 565 731 652 297 925 769 886 549 168 450 867 091 163 868 402 442 018 371 341 022 143 349 040 941 897 422 377 931 540 131 712 376 088 451 926 646 802 255 265 721 820 332 564 453 961 091 792 448 710 886 674 897 251 899 281 432 897 966 817 352 381 951 694 370 184 383 299 109 654 078 410 356 489 144 960 977 274 480 320 655 205 474 842 571 310 858 767 154 579 079 644 215 514 199 003 832 599 232 631 020 684 528 203 955 813 111 140 473 074 686 810 426 051 322 892 063 650 122 841 004 997 271 652 977 702 547 274 337 834 459 757 112 850 838 135 183 696 311 712 170 990 224 739 928 697 694 772 853 140 476 802 033 735 594 574 707 189 941 820 656 793 611 572 800 671 948 131 128 109 558 012 358 679 446 259 496 244 143 021 997 936 801 273 888 393 102 412 918 463 505 832 701 613 704 464 047 394 891 693 567 824 946 365 095 976 820 092 895 872 874 833 067 345 472 399 549 213 736 568 277 881 574 159 120 670 765 215 962 463 587 184 455 433 922 702 506 813 953 637 742 383 768 964 226 648 365 610 070 168 839 324 749 896 910 992 390 804 658 988 479 107 678 562 550 844 858 514 512 762 455 637 677 863 067 821 690 083 222 834 405 613 388 709 119 188 719 479 181 528 190 836 888 908 152 845 729 021 513 491 658 586 509 889 455 613 570 273 121 541 638 995 145 894 561 079 999 213 500 076 101 248 781 500 642 979 091 198 044 736 325 729 107 144 579 768 273 776 743 716 439 557 581 557 253 878 204 647 391 176 582 851 167 312 182 026 188 795 156 620 056 856 503 340 092 247 479 478 684 738 621 107 994 804 323 593 105 039 052 556 442 336 528 920 420 940 313
pro=
+ U" p6 N/ w7 s3 245 708 184 223 678 ... ...

结果太长，放到了一个sumpro.txt文件中。
sumpro.txt (1.5 MB, 下载次数: 0)

2010-11-6 16:09 上传

点击文件名下载附件
下载积分: 体力 -2 点

无痕无泪

路过，呵呵。。。。。

xxt19890711

飘过。。。。。。