三浦と窮理とブログ

自分の経験したことを検索可能にしていくブログ.誰かの役に立ってくれれば嬉しいです.

2次元イジングモデルのオンサーガー解の数値計算

無限体積の2次元イジングモデルの厳密解であるオンサーガー解を数値的に計算する. 1スピン当たりの自由エネルギー \begin{align} -\beta f &= \log(2\cosh 2K) + \frac{1}{2\pi ^ 2} \int ^ \pi _ 0 d\omega _ 1 \int ^ \pi _ 0 d\omega _ 2 \log(1 - k ^ 2…

Qiskit で平方根ゲート を使う方法

Qiskit で平方根ゲートを作るには power メソッドを使います. 例えば $\sqrt{NOT}$ ゲートは XGate().powre(.5) です. power メソッドは Gate クラスに対して使えます. 参考 Gate — Qiskit 0.13.0 ドキュメント 使用例 実際に $\sqrt{NOT}$ ゲートをつく…

Qiskit で一度追加した量子ゲートを削除する方法

Qiski で量子回路を組んでいて,一度 QuantumCircuit に追加した Gate (およびInstruction )を削除する方法を説明します. ドキュメント QuantumCircuit に追加した Instruction の情報は,data という属性にリストになって保存されています *1. ですので…

量子プログラミング入門【物理学科生向け】

前置きの前置き 量子コンピュータについて少し勉強してみたので,量子プログラミング入門という題でまとめていきたいと思います. 「現在利用可能な量子コンピュータを触ってみよう」というようなモチベーションです. ハードウェア自体は現在も世界各地で絶…

剰余の変域

問題 L ≦ x ≦ R な自然数 x に対して, 自然数 m を法とする剰余の値 y = x%m の変域は if (R - L ≥ m) 0 ≦ y ≦ m - 1 else if ($\lfloor$L/m$\rfloor$ = $\lfloor$R/m$\rfloor$) L%m ≦ y ≦ R%m else if ($\lfloor$L/m$\rfloor$ + 1 = $\lfloor$R/m$\rfloor$…

A - Diverse Word:AGC22

問題 A - Diverse Word 入力された多彩な単語に対し,辞書順で次の多彩な単語を求める. 理解 多彩な単語の辞書がすぐには想像できないので書き下すことにする.「見たものは想像できる」の精神. a から始まる多彩な単語について,次の図のように木の形に単…

sundiskのusbメモリが書き込み禁止になったので,sumsungのメモリを買った.

1年前に買った Sundiskのusbメモリが突然書き込み禁止になりました. ネットで調べた書き込み禁止解除の方法を試しても治らなかったので新しいusbメモリを買いました. 今回壊れたのは SanDisk Ultra Fit 128GB です. 【 サンディスク 正規品 】5年保証 USB…

qiskit.aqua で shor アルゴリズムの実行

qiskit.aqua に実装されている shor アルゴリズムを実行してみます. 環境は windows10, python3.7 です. 実行例 jupyter で実行してみました. シミュレーションだと,とても遅いです.N=33のときには3時間かかった上,素因数を見つけられませんでした. …

冪剰余と剰余の逆元を求めるC++/Boostライブラリ

冪剰余 冪剰余 $a ^ b \bmod m$ を求めます.たぶん二分累乗法だと思います. #include <boost/multiprecision/integer.hpp> boost::multiprecision::powm(a, b, m) 剰余の逆元(拡張ユークリッド互除法) 剰余の逆元 $a^{-1} \bmod m$ が存在すれば求めます. (2019年10月5日追記:これはboos</boost/multiprecision/integer.hpp>…

Qiskit.Aqua を使おうとしたら,numpy+mklを要求されたこと.

環境 win10,python3.7 >>>import qiskit >>>qiskit.__qiskit_version__ {'qiskit-terra': '0.9.0', 'qiskit-ignis': '0.2.0', 'qiskit': '0.12.0', 'qiskit-aer': '0.3.0', 'qiskit-ibmq-provider': '0.3.2', 'qiskit-aqua': '0.6.0'} エラー内容 Qiskit.Aq…