Pi günü

Üçüncü ayın ondördü, Monte Carlo'dan selamlar...

python2 -c "
from random import random
l = lambda (a, b), c: (a + c, b + 1)
g = (1 if random()**2 + random()**2 <= 1 else 0 for i in range(10**6))
r, t = reduce(l, g, (0.0, 0))
print r / t * 4.0
"

Görüşler

FZ
7 yıl önce

A-aaa, Buffon! :)

auselen
7 yıl önce

tras olunca taniyamadim.

auselen
7 yıl önce

Genel sıkıntım, bu tarz basit programları bir satırda yazacak bir programlama dili bulamamış olmam. Herhangi bir öneriye açığım. Aynı şekilde birisi python ile daha düzgününü yazabiliyorsa, görmeyi çok isterim.

FZ
7 yıl önce

Valla her birkac senede bir ogrenmeye heves ettigim APL turevi dillerden biri care olabilir zira uzaktan bakinca asiri 'kompakt' gorunuyorlar. Misal:

auselen
7 yıl önce (güncellenmiş)

https://news.ycombinator.com/item?id=13874165

tongucyumruk
7 yıl önce (güncellenmiş)

Hmm, tek satir mi? Challenge accepted:

print (lambda (r, t): r / t * 4.0)((lambda random: reduce(lambda (a, b), c: (a + c, b + 1), (1 if random()**2 + random()**2 <= 1 else 0 for i in range(10**6)), (0.0, 0)))(__import__("random").random))

Daha duzgunu mu? Tabi ki hayir, ama saniyorum APL ile yazilmis olsa da benzer seviyede okunamaz olurdu.

tongucyumruk
7 yıl önce

Denemek isteyenler icin: https://repl.it/GW3m

auselen
7 yıl önce (güncellenmiş)

from random import random; print 4 * sum(1.0 if random()**2 + random()**2 <= 1 else 0 for i in range(10**6)) / 10**6

10**6 yi tekrar etmeden, import kismi kaldirilsa tatmin olurdum, herhalde. "Running average" tarzi bir yöntem bekli de?

FZ
7 yıl önce

J kodu:

pi=: (0.25&* %~ +/@(1 >: [: +/&.:*: _1 2 p. 0 ?@$~ 2&,))"0

Komut satirindan denemek icin:

$ ijconsole
   pi=: (0.25&* %~ +/@(1 >: [: +/&.:*: _1 2 p. 0 ?@$~ 2&,))"0
   pi 1e6
3.14104

FZ
7 yıl önce

Octave kodu:

>> sum(rand(1,100000).^2+rand(1,100000).^2<1)/100000*4
ans =  3.1435

FZ
7 yıl önce

Biraz daha sadelestirirsek:

>> t=10000000;sum(rand(1,t).^2+rand(1,t).^2<1)/t*4
ans =  3.1415

auselen
7 yıl önce (güncellenmiş)

araya noktali virgul koyunca, ya da komut satirindan bir seferde calistiramayinca pek tadi yok.

Moving average versiyonu. import olmasa bu yeterdi bana.

from random import random
print 4 * reduce(
    lambda a, (c, d): (a * (d - 1) + c) / d,
    ((1.0 if random()**2 + random()**2 <= 1 else 0, i) for i in range(1, 10**6)),
   0)

FZ
7 yıl önce (güncellenmiş)

Komut satirindan calistirmakta bir sıkıntı yok Octave icin de, noktali virgulu de kaldiralim gitsin:

$  octave --eval "sum(rand(1,10^6).^2+rand(1,10^6).^2<1)/10^6*4"
ans =  3.1488
$

auselen
7 yıl önce (güncellenmiş)

Şimdi de 10**6 yi azaltalım.

auselen
7 yıl önce

Vektör çarpımı ile olabilir belki.

FZ
7 yıl önce

Requirement: "Şimdi de 10**6 yi azaltalım."

Deliverable:

octave --eval 'feval(@(x) sum(eval(x) + eval(x) < 1), "rand(1,10^6).^2") / 10^6*4'
ans =  3.1421

Bu kadar "golf" yeter bugun icin. Herkese bol vektorel islemli gunler :)

onur
7 yıl önce

Python tek satırlıklarıyla pek güçlü bir dil değil, bu işin ehli elbette Perl. Farklı bir algoritma ve perl ile yazılmış bir Pi hesaplayıcı¹:

perl -e '()while(printf"%.16f\n",$,+=((-1/++$_/(($?||=1)*=($_-$_++)*-1)))*4)'

butch
7 yıl önce (güncellenmiş)

Klavyenizden Pi sayısına piyano ile eşlik etmek için buraya.

Görüş belirtmek için giriş yapın...