2011-08-18から1日間の記事一覧

ユークリッドの互除法

最大公約数を求める。 #! usr/bin/env python # -*- coding: utf-8 -*- #euclidian1.py m = int(raw_input("Set m : ")) n = int(raw_input("Set n : ")) while m != n: if m >= n: m = m-n else: n = n-m print "GCC = ",m n>nの時有効。 #! usr/bin/env py…

モンテカルロ法

モンテカルロ法で円周率πを求める。 このページとかが参考になる。 http://www.f.waseda.jp/takezawa/math/number/Pi/monte.html #!usr/bin/env python # -*- coding: utf-8 -*- #montecarlo.py import sys import random N=raw_input("input trial number :…

ランダムな順列

乱数を使ってランダムな順列を得る。 その前にpythonで擬似乱数を得る方法について。 >>> random.randint(1, 10) # Integer from 1 to 10, endpoints included 7 はじめに効率の悪い方法。 アルゴリズムは、 (1)1~Nの乱数を1つ得る。これを順列の1番目のデ…