בתשובה ליוסי, 13/03/14 12:01
אפרופו אלון עמית 628898
האלגוריתם ''שלי'' (אף פעם לא חשבתי על הבעיה אלא נתקלתי ישר בפתרון שלה, כך שלא מדובר כאן על שום דבר מקורי) אכן טריויאלי.
אפרופו אלון עמית 628899
אם ככה, הטרוויאלי שלי הוא משהו הזה בשורת הפקודה:

perl -e 'sub R(){$a = rand(); $a*$a} sub IN{ R+R < 1} $c = $n = shift; $sum += IN while $c-- ; print 4*$sum/$n, "\n"' 100000

אפרופו אלון עמית 628901
נראה לי ששכ"ג התכוון ל-Buffon's needle.

בקורס בתכנות סטטיסטי שאני מלמד, אני מדגים בדיוק את השיטה שלך (אם כי בשפת R, ולא ב-perl).
אפרופו אלון עמית 628905
לא, עד כמה שידיעתי בפרל מגעת התכוונתי דווקא לפתרונו של יוסי (במלים: היחס בין מס' הנקודות בתוך מעגל היחידה למספר הנקודות בתוך הרבוע החוסם). אבל המחט של בופון יותר מגניבה (והוזכרה כבר באתר זה, נחש בידי מי... תגובה 150004)
אפרופו אלון עמית 629010
בדיוק גדי העלה פוסט על המחט של בופון: http://www.gadial.net/2014/03/14/buffon_needle
אפרופו אלון עמית 629037
ראיתי. בכל פעם שאני נתקל בפתרון הפנטסטי הזה‏1 אני שמח מחדש.
______________
1- אני מתכוון, כמובן, לאינטגרל הכפול :-)

חזרה לעמוד הראשי המאמר המלא

מערכת האייל הקורא אינה אחראית לתוכן תגובות שנכתבו בידי קוראים