#はじめに
今回はVisual Studio 2017でMPIR(GMP)の導入方法について解説していきます。
そろそろMPIR(GMP)に手を出してみようかな
と思っている方や、
Visual StudioでプログラミングしてるけどMPIR(GMP)の入れ方がわからない
という方に向けて書いてます。
これを見れば誰でもMPIR(GMP)が導入できる！はず
⇒こちらからMPIR〇〇source.zip をダウンロードします。
MPIR〇〇source.zip を選択します。
この解説では、 MPIR 3.0.0 source.zip を使います。
#ダウンロード後
ダウンロードしたMPIR〇〇source.zipを すべて展開(解凍) します。
今回、展開場所は Cドライブ直下(Windows (C:)) に置きます。
展開はそこそこ時間がかかります。
展開が終わったら次へ。
#dll_mpir_gcをビルドする
展開したフォルダ、C:\mpir-〇〇を開くと
build.vcXX(XXは数字)というフォルダがあると思います。
１番数字が大きいものを開いてください。
mpir.slnファイルがあるのでVisual Studio 2017で起動します。
この解説ではC:\mpir-3.0.0\build.vc15\mpir.slnにありました。
起動したら最初に左上にある
DebugをReleaseに
Win32をx64に変えてください。
次にソリューションエクスプローラーから
dll_mpir_gcを選択し、右クリックすると出てくる
プロジェクトの再ターゲットを選択します。
OKを選択します。
次もソリューションエクスプローラーから
dll_mpir_gcを選択し、右クリックすると出てくる
ビルドを選択します。
無事に出力できました。
mpir.slnを閉じましょう。
#Visual Studioでの設定
さて、ここからが本番です。
画面上のメニューの プロジェクト(P) から 〇〇のプロパティ を押します。
左の C/C++ ⇒ 全般 ⇒ 追加のインクルードディレクトリ に下記を追加します。
※ダウンロードしたMPIR(GMP)のバージョンによって適宜、名前を変えてください。
C:\mpir-3.0.0
左の リンカー ⇒ 全般 ⇒ 追加のライブラリディレクトリ に下記のどちらかを追加します。
※ダウンロードしたMPIR(GMP)のバージョンによって適宜、名前を変えてください。
C:\mpir-3.0.0\build.vc15\dll_mpir_gc\Win32\Release
C:\mpir-3.0.0\build.vc15\dll_mpir_gc\x64\Release
左の リンカー ⇒ 入力 ⇒ 追加の依存ファイル に下記を追加します。
※ダウンロードしたMPIR(GMP)のバージョンによって適宜、名前を変えてください。
mpir.lib
今回は、コンソールのプロジェクトとして使用するので、
左の リンカー ⇒ システム ⇒ サブシステム を
コンソール (/SUBSYSTEM:CONSOLE)
にしましょう。
最後に作成したプロジェクトの
ソースコードが入っているフォルダに
**mpir.dll**を入れたら準備完了です。
mpir.dll は C:\mpir-〇〇\dll\(Win32 or x64)\Releaseにあります。
プロジェクト名を右クリックした後、
"エクスプローラーでフォルダを開く"
を押すと該当するフォルダが開きます。
この解説では、
C:\mpir-3.0.0\dll\x64\Release\mpir.dll
にありました。
#MPIR(GMP)が使えるかどうかのテスト
今度はメルセンヌ素数を求めてみましょう。
#include <cstdint>
#include <iostream>
#include <mpirxx.h>
#include <chrono>
bool lucasLehmerTest(const int_fast32_t p_) {
mpz_class m = 1;
m <<= p_;
--m;
mpz_class s = 4;
for (int_fast32_t i = 2; i < p_; ++i)
s = (s*s - 2) % m;
if (s == 0) return true;
return false;
}
bool lucasLehmerTestFast(const int_fast32_t p_) {
mpz_class m{ 1 };
m <<= p_;
--m;
mpz_class s{ 4 };
mpz_class sqrt;
for (int_fast32_t i{ 2 }; i < p_; ++i) {
sqrt = s * s;
s = (sqrt & m) + (sqrt >> p_);
if (s >= m) s -= m;
s -= 2;
}
if (s == 0) return true;
return false;
}
int main() {
std::chrono::system_clock::time_point start, end;
start = std::chrono::system_clock::now();
if (lucasLehmerTestFast(9689)) std::cout << "TRUE" << std::endl;
else std::cout << "FALSE" << std::endl;
end = std::chrono::system_clock::now();
std::cout << std::chrono::duration_cast<std::chrono::milliseconds>(end - start).count() << std::endl;
return 0;
}
TRUE
232
実行できました！
##ソースコードのライセンス
These codes are licensed under CC0.
ソースコードは自由に使用してください。