三浦と窮理とブログ

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

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

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

剰余の変域

問題 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…

Qiskit で量子フーリエ変換を実装する.

量子フーリエ変換について勉強してQiskitで実装してみたのでまとめていきたいと思います. 量子フーリエ変換(QFT) については日本語で解説しているサイトがたくさんあるのでとても助かりました. Qiskit による実装 python3.7 で実行しています. IBM Q Expe…

Mathjax version 3 で数式表示が速くなった

このブログでは数式の描画にMathjaxを使っています。 Mathjaxのversion 3 がでてたので今までの version 2.7.5 から移行しました。 設定の仕方が少し変更されてました。 デモが用意されていたのですぐにできました。 今まで書いた記事全てを確認したわけでは…

量子コンピュータで量子もつれ状態を生成する(Qiskitチュートリアル)

Qiskit をインストールして,チュートリアルの 量子もつれ状態の生成を実機で行うまでやってみたのでまとめていきます.公式のドキュメントがとても詳しいので簡単でした. windows10 ,python 3.7.3 で行いました. Qiskit のインストール Qiskit の APIド…

TOSHIBA dynabook R734/K ノートパソコンのバッテリー交換

ノートパソコンのバッテリーの持ちが悪くなってきたので交換しました. 仕様 ノートパソコンの機種は TOSHIBA dynabook R734/K PR734KEF637AD31 です. 標準で装着されているバッテリーは バッテリパック31AQ,PA5161U-1BRS,PABAS276 という型番です. 10.8…

C++ で2次元配列の最大値を求める.

C++ で1次元配列の最大値を求めるならばfor文を使わずに max_element を使うと楽です. 2次元配列の最大値を一度に求めるSTL関数は私は見つけられませんでしたが,各行に対して max_element を繰り返せばいいでしょう. max_element のおかげで for を書く手…

C++ で std::queue を定数時間で空にする方法

C++ でstd::queue の内容を全て削除したくなったのですが,std::queue には clear 関数がありませんでした. なので,私は直観的には queue が空になるまで pop を繰り返せばいいと思ったのですが,これはqueueのサイズに比例した時間がかかります. もっと…

nvidia/cudaのdockerイメージをrunする時にNVidiaドライバーのバージョンが古いため出るエラーメッセージ

動作環境 グラボ GeForce GTX1080Ti NVidiaドライババージョン 384.66 Docker バージョン 18.09.7, build 2d0083d nvidia-dockerバージョン 2.0.3 エラー発生 上述の環境でCUDA10.1が入ったdockerイメージを次のようにrunしようとしました. $ docker run --…

numpy.gradientで配列の微分をするとき,両端の微分の精度を上げるには

numpy.gradientで配列を微分するとき,配列の境界部分(始端終端)の値だけ目的の正しい微分値が出ないことがあった. すぐに解決したのでメモです. 境界部分での微分はedge_orderというオプションで精度を選べる. numpy.gradient — NumPy v1.16 Manual 何…

TOSHIBA dynabook R734/K ノートパソコンのmSATA SSDの換装(128GB→240GB)

ノートパソコンに入っていた128GBのmSATA-SSDを240GBのものに換装しましたのでその記録です.Cドライブが入ったSSDなのでクローン作業も必要でした.mSATAを扱うのは初めてだったので,面白かったです. 目次 仕様 作業工程 新SSDをmSATA-USB変換アダプタに…

【熱力学】気体の定圧比熱と定積比熱の差とエンタルピー

気体の定圧比熱と定積比熱の差はで与えられる.また,温度変化に伴う体積変化が $ dV = (\partial V / \partial T) _ {過程} d T $ であるような準静的過程での比熱 $ C _ {過程} $ はエンタルピーを $ H=U+p V $ としたとき,比熱は熱力学第一法則より,

【熱力学】体膨張率と等温圧縮率が一定な気体の状態方程式

気体の体膨張率を $ \alpha $ ,等温圧縮率を $ \kappa _ T $ とするとき,次の関係が成り立つ. \begin{equation} \frac{1}{V} d V=\alpha d T-\kappa _ {T} d p \end{equation} $ \alpha $ と $ \beta $ が温度や圧力に依らぬ定数であったとし,温度 $ T _…

【熱力学】ファン・デル・ワースル気体のビリアル展開とジュール温度の計算

ファン・デル・ワースルの状態方程式で表される気体について.等温圧縮率熱圧力係数体膨張率である.次に,式\eqref{eq:netu5p}よりビリアル展開することを考える.となるので,第2ビリアル係数は $ b-\frac{a}{RT} $ である.この係数は高温では正の値 $ b $…

【熱力学】気体の状態方程式といろいろな係数

気体の状態方程式を $ f(p,T,V)=0 $ とする.関数 $ f $ を $ p,T,V $ で微分して得られる導関数をそれぞれ $ f _ p,f _ T,f _ V $ とする.これらを用いると\begin{alignat}{2}&等温圧縮率 \quad& \kappa _ T &=-\frac{1}{V}\left ( \frac{\partial V}{\par…

python3 のループ文のelse節についてのチュートリアル

python3.6のチュートリアルでループ制御でのelse節について知ったのでメモです. docs.python.jp ループのelse節の動作についてはチュートリアルに書いてある通りでたぶん理解できたと思います. チュートリアルにあったサンプルプログラムについて考えてみ…

【熱力学】高所ほど気圧が低くなることについて

地球上で空気が高さ方向に静止していて熱平衡であり,温度が $ T $ とする.単位面積の底面をもつ空気柱を考える.高さ $ h $ と $ h+dh $ の間の圧力差 $ dp $ は \begin{equation} d p = - \rho ( h ) g d h \end{equation} となる. $ \rho(h) $ [質量/体…

TikZ Wolke を使ってみる (TikZをsvgで表示してくれるjavascript)

TikZのソースコードを使ってsvg画像を生成してwebページ上で表示してくれるjavascriptが4か月くらい前に公開されたらしいので使ってみる.mathjaxはtikzに対応していなかったので,実際に使うかわからないが少なくとも方法があることを知れてよかった. 以下…

温度計の表示温度と物体の温度の関係

接触式温度計で物体の温度を計る時,温度計の示度は温度計自体の温度を表しているため,温度計と熱平衡になった物体の温度を示している.温度計と接触させる前の物体の温度 $ T $ を求めよう. 物体の熱容量を $ C $ ,接触前の温度計の示度を $ T _ i $ ,…

位相空間論の連結性についての例題

位相空間論の連結性についての例題を解きましたので,解答を載せていきます.連結な5点位相空間の例連結ではない4点位相空間の例互いに交わる2つの連結集合の共通部分が連結集合にならない例互いに交わる2つの連結集合の和集合は連結である.一般線形群は連…

位相空間論のコンパクト性についての例題

位相空間論のコンパクト性についての例題を解きましたので,解答を載せていきます.コンパクトではない集合の例n次元ユークリッド空間 はコンパクトではない密着空間 はすべての部分集合がコンパクトである2つのコンパクト集合の和集合もコンパクト集合であ…

位相空間論の積空間・商空間についての例題

位相空間論の積空間や商空間についての例題を解きましたので解答を載せていきます. 定義などの細かいことは次の本を参考にしています. 集合と位相 (数学シリーズ)作者: 内田伏一出版社/メーカー: 裳華房発売日: 1986/11/01メディア: 単行本購入: 4人 クリ…

位相的性質・ハウスドルフ性 についての例題

位相空間のハウスドルフ性についての例題を解いたので解答を載せていきます.位相空間 $ X $ がハウスドルフ空間である.#ハウスドルフ空間ではない位相空間の例#ハウスドルフ空間の部分位相空間はハウスドルフ空間である.#ハウスドルフ空間の3点を分離する…

百均のスマホ充電ACアダプタが変形したので分解した

百均で買ったスマホの充電ACアダプタの端子の付け根が以下の写真のように変形していた. 左の端子の根本が黄色く焼け焦げたようになっていて,端子が浮き出ている. こういう電源周りのものは発火する危険があるので,いろいろ無頓着な私でも(だからこのア…

TOSHIBA dynabook R734/K ノートパソコンのメモリ増設(12GB→16GB)

以前ノートパソコンのメモリを 4GB から 4 + 8 = 12GB に増設したのですが,今回は 8 + 8 = 16GB に増設します. 目次 仕様 作業内容・動画 メモリの取りはずし メモリの取り付け 動作確認 あと一年はハードに使っていきたいです. 仕様 前回と同じノートパ…