python
最急降下法とは、目的関数の最小化元(或いは最大化元)を求める問題に対する反復法と呼ばれる解法の一つです。 目的関数fと勾配、初期点を問題に合わせて入力する必要があります。
要素が一つのタプルを作りたいときに a = (1) としても, >>> a 1 >>> type(a) <class 'int'> でaは整数型になっていることが分かります。 じゃあ、どうすればいいかというと a = (1, ) もしくは a = 1, としてやれば、いいです。 なんか気持ち悪いですね。</class>
pythonのpulpモジュールを用いて、数独を線形計画法で解きます。 pulpはpip3 install pulpでインストールできます。
pythonの標準モジュールmultiprocessingを使用することで、簡単にプロセス並列処理を行うことができます。 Pool map 基本的な使い方は from multiprocessing import Pool with Pool(processes=None) as pool: pool.map(func, iter) Poolの引数であるprocesse…
普通にリストの内包表記の感覚でタプルの内包表記を書こうとすると a = (i for i in range(10)) type(a) > <class 'generator'> と、ジェネレータが作成されてしまいます。 タプルが欲しいのに….と言うときは、 b = tuple(a) とtupleでジェネレーターをラップさせればokです。 </class>…
今までpickleファイルを作成するときは、 with open(file, 'wb') as f: pickle.dump(obj=obj, file=f) としていましたが、プロトコルバージョンを高いものに指定してやると、より速く読み込めるpickleファイルが作れるらしい。。 最新のプロトコルバージョン…
functoolsのreduce関数を使って、リストの中身を全て足したり掛けたり、といった操作を行えます。 from functools import reduce from operator import add, mul a = [1, 2, 3, 4] # 適当なリスト reduce(add, a) # 総和 > 10 reduce(mul, a) # 総乗 > 24 こ…
pythonでopen street mapのデータを読み込めるosmnxという便利なモジュールがあるのですが、 pip3 install osmnxとインストールしようとすると、 fiona/_transform.cpp:473:10: fatal error: 'cpl_conv.h' file not found #include "cpl_conv.h" ^ 1 error g…
pythonのformat形式で例えば{ゴリラ} という文字列を作りたいときは animal = 'ゴリラ' print('{{}}'.format(animal)) print('{{animal}}'.format(animal=animal)) と{を2つ重ねます。
Cython その2 - サブロウ丸の続きです。 今回はCythonでの配列の型付けの話をします。 配列の要素が全て同一の型の場合、(当然ですが)そのことを宣言した方がより高速になります。 Cythonで上記を実現するにはメモリービューと呼ばれるものを使用します。 nu…
Cython その1 - サブロウ丸 前回からの続きです。 cython tsp.pyx --annotate もしくは cython tsp.pyx -a を実行すると tsp.html というファイルが作られます。 tsp.html · GitHub これをブラウザ等で開くと という、画面が出てきます。 簡単に言うと、濃い…
Cythonを使ってPythonのコードを高速化します。 Cythonのインストール pip3 install cython コンパイル cythonのインストールが終わったら、スクリプトをc言語にコンパイルするために、以下のsetup.pyを用意します。 (ext_modules = [Extension('tsp', ['tsp…
巡回セールスマン問題をpythonで解きます。 01整数計画法を使っても解けますが、今回は動的計画法を使用します。 解説 解説用のpdfを文書を作成したので是非ご覧ください。 drive.google.com pythonコード DPは配列ではなく、辞書を使用しています。 次回か…
全く同じ画像だけではなく、より幅広く'似てそう'な画像を探します。 OpenCVのインストール OpenCVのインストールに関するページがネット上に多数あることから、OpenCVのインストールの難しさ(というより、やっかいさ?)が伺えますが、例に違わず僕もOpenCVの…
下記のような、集合と自然数の対応付けを考えます。 {0, 1, 2, 3} <—> 15 基本的な考えは2進数表記です。 {0, 1, 2, 3} <—> 20 + 21 + 22 + 23 = 15 {0, 1, 3} <—> 20 + 21 + 23 = 11 これをそのまま関数にすると、 def set2int(_iter): res = 0 for i in _i…
subprocess シェルの中で他のプログラムを起動し、そのプログラムが生成した出力を知りたいだけなら、getoutput()関数を使う。 import subprocess as sp ls = sp.getoutput('ls') ls > 'Applications\nDesktop\nDocuments\n...(略)' オプションとパイプ、リ…
下記の問題を解きます。 kamipeipaa君は新しいものが大好きです。 kamipeipaa君はある日N個の整数A1,A2,A3,....,ANを見つけました。 整数Aiに対して,Ai=Ajとなるjがi以外に存在しなければ,Aiは新規性があるのではないかとkamipeipaa君は考えました。 上記…
プログラムのループ文の中で計算の途中結果を出力したいけど、 画面につらつらと出力されるのは鬱陶しいなー...という時に便利な、出力を上書きしていく方法です。 サンプルコード from time import sleep for i in range(10): print('\r', i, end='') sleep…
PythonパッケージのNumbaのインストールに手こずったので、記録。 とりあえず、やったこと numbaのインストールにはllvmとllvmliteが必要とのことなので e-1. llvmのインストール brwe install llvm e-2. export LLVM_CONFIG=/usr/local/Cellar/llvm/3.9.1/b…
文字列の中に様々な形式で値を入れる方法です。 %を使う方法 a = 'I am %d years old.' % 20 > a 'I am 20 years old. ' 様々な言語でお馴染みですね。 指定 データ型 %s 文字列 %d 10進整数 %x 16進整数 %o 8進整数 %f 10進float %e 指数形式float %g 10進f…
pythonプログラムのパフォーマンス、ここでは 1. 行単位の実行時間と実行頻度 2. 行単位のメモリ使用量 を調べる方法を紹介します。 参考ページ(Pythonスクリプトのパフォーマンス計測ガイド | Yakst) 行単位の実行時間と実行頻度 line_profilerモジュールを…
lru_cache lru_cacheは関数の引数と返り値を保存する。 from functools import lru_cache @lru_cache(maxsize=None) def fib(n): if n < 2: return n return fib(n-1) + fib(n-2) 公式の例。再帰的にフィボナッチ数を計算していますが、そのまま実行すると呆…
エラトステネスの篩を使った素数判定法をpythonで実装しました
pythonスクリプト内に定数やその他データを入力する方法で、僕が知っているものをまとめました。 大きく分けて以下の3つです。 1. sys(引数入力) 2. input(標準入力) 3. fileからの入力 sysを使う(引数入力) python3 file名.py [引数1] [引数2] ... という形…
前3回ほど、pythonのlist操作についての記事を書きましたが、 Python コーディング(基) list編 - サブロウ丸 Pythonコーディング(注) list編 - サブロウ丸 Python コーディング(豆) list編 - サブロウ丸 その応用として、 みなさんお馴染みのsortingをコーデ…
pythonのlist操作に関する豆達です。 enumerate: index付きforループ a = ['I', 'like', 'cats'] for i, word in enumerate(a): print(i, word) > 0 I > 1 like > 2 cats index: listにおける指定した要素のindexを取得 a = [1, 2, 3] a.index(2) > 1 #要素…
copy pythonではすべてのオブジェクトは参照渡しになっています。だからlist等をcopyするときには注意が必要です。 a = [0, 1, 2] b = a b[0] = 'Hello' b > ['Hello', 1, 2] a > ['Hello', 1, 2] # bを変更するとaも変更されてしまう これを避けるには copy…
list操作(基)sortの話です. sort sortとsortedの2種類. 使い方その他は下記. a = [2, 3, 4, 1] sorted(a) > [1, 2, 3, 4] #aのsort後のlistが出力される a > [2, 3, 4, 1] #aに変化はなし a.sort() a > [1, 2, 3, 4] #a自体がsortされている. sorted(a, reve…
あえてターミナルから辞書を開く必要はないと思いますが、今日はその方法を紹介します。 単純な方法 ホームディレクトリの.bashrcを開いて dict(){ open dict://$1 } を追記して保存する。 すると、例えばdict appleと入力すると、辞書appのappleのページが…
前回からの続きです pipの導入 - サブロウ丸 pip 基本ではないかもコマンド pip freeze インストールしたパッケージリストを出力します。 例えばpip freeze > hoge.txtでパッケージリストをリダイレクトで保存しておき、pip install -r hoge.txtとすることで…