C++ - 階乗計算(with GMP)！

C++ で任意精度算術ライブラリ GMP(The GNU Multi Precision Arithmetic Library) を使って、階乗の計算をしてみました。

0. 前提条件

• Debian GNU/Linux 10.3 (64bit) での作業を想定。
• GCC 9.2.0 (G++ 9.2.0) (C++17) でのコンパイルを想定。
• GMP 6.2.0 でのコンパイルを想定。
  （GMP がインストール済みであること。インストール方法： GMP - ソースビルドでインストール (on Linux Mint)！）

1. ソースコードの作成

File: factorial.cpp

/***************************************************************
  Factorial of N
  by GMP(The GNU Multi Presicion Arithmetic Library).

    DATE          AUTHOR          VERSION
    2020.08.17    mk-mode.com     1.00 新規作成

  Copyright(C) 2020 mk-mode.com All Rights Reserved.
***************************************************************/
#include <iostream>  // for cout
#include <gmp.h>
#include <gmpxx.h>

class Factorial {

public:
    mpz_class fact(mpz_class);
};

/*
 * Caculation of fatorial
 */
mpz_class Factorial::fact(mpz_class n) {
  try {
    if (n == 0)
      return 1;
    return n * fact(n-1);
  } catch (...) {
    throw;
  }
}

int main(int argc, char* argv[]) {
  Factorial f;
  mpz_class n;
  std::string buf;

  try {
    while (true) {
      std::cout << "n? ";
      getline(std::cin, buf);
      if (buf.empty())
        break;
      n.set_str(buf, 10);
      std::cout << n << "! = "
                << f.fact(n).get_str() << std::endl;
      std::cout << "---" << std::endl;
    }
  } catch (...) {
    std::cerr << "EXCEPTION!" << std::endl;
    return EXIT_FAILURE;
  }

  return EXIT_SUCCESS;
}

• Gist - C++ source code to calculate factorials with GMP.

2. ソースコードのコンパイル

$ g++ -std=c++17 -Wall -O2 --pedantic-errors -lgmp -lgmpxx -o factorial factorial.cpp

3. 動作確認

$ ./factorial
n? 0
0! = 1
---
n? 1
1! = 1
---
n? 2
2! = 2
---
n? 3
3! = 6
---
n? 4
4! = 24
---
n? 5
5! = 120
---
n? 100
100! = 9332621544394415268169923885626670049071596826438162146859296389521759999
32299156089414639761565182862536979208272237582511852109168640000000000000000000
00000
---
n? 1000
1000! = 402387260077093773543702433923003985719374864210714632543799910429938512
39862902059204420848696940480047998861019719605863166687299480855890132382966994
45909974245040870737599188236277271887325197795059509952761208749754624970436014
18278094646496291056393887437886487337119181045825783647849977012476632889835955
73543251318532395846307555740911426241747434934755342864657661166779739666882029
12073791438537195882498081268678383745597317461360853795345242215865932019280908
78297308431392844403281231558611036976801357304216168747609675871348312025478589
32076716913244842623613141250878020800026168315102734182797770478463586817016436
50241536913982812648102130927612448963599287051149649754199093422215668325720808
21333186116811553615836546984046708975602900950537616475847728421889679646244945
16076535340819890138544248798495995331910172335555660213945039973628075013783761
53071277619268490343526252000158885351473316117021039681759215109077880193931781
14194545257223865541461062892187960223838971476088506276862967146674697562911234
08243920816015378088989396451826324367161676217916890977991190375403127462228998
80051954444142820121873617459926429565817466283029555702990243241531816172104658
32036786906117260158783520751516284225540265170483304226143974286933061690897968
48259012545832716822645806652676995865268227280707578139185817888965220816434834
48259932660433676601769996128318607883861502794659551311565520360939881806121385
58600301435694527224206344631797460594682573103790084024432438465657245014402821
88525247093519062092902313649327349756551395872055965422874977401141334696271542
28458623773875382304838656889764619273838149001407673104466402598994902222217659
04339901886018566526485061799702356193897017860040811889729918311021171229845901
64192106888438712185564612496079872290851929681937238864261483965738229112312502
41866493531439701374285319266498753372189406942814341185201580141233448280150513
99694290153483077644569099073152433278288269864602789864321139083506217095002597
38986355427719674282224875758676575234422020757363056949882508796892816275384886
33969099598262809561214509948717012445164612603790293091208890869420285106401821
54399457156805941872748998094254742173582401063677404595741785160829230135358081
84009699637252423056085590370062427124341690900415369010593398383577793941097002
77534720000000000000000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000000000000000000000000000000000
0000000000000000
---
n?

（実際には改行されていない）

---

以上。