ÊÖ¿®¤¢¤ê¤¬¤È¤¦¤´¤¶¤¤¤Þ¤¹¡£
³Î¤«¤ËÁ²²½¼°¤Ë¤è¤ë£´¼¡¼ý«¤Î¼°¤ò»È¤¦¤Î¤¬¿ô³Ø³¦¤Ç¤Ï°ìÈÌŪ¤Ê¤è¤¦¤Ç¤¹¡£
Borwein¤Î¸ø¼°¤Ç¤ÎÊýË¡¤Ê¤éÅìµþÂç³Ø¤Î¶â¤ÀÀèÀ¸¤Îºî¤é¤ì¤¿¥×¥í¥°¥é¥à¤¬ÌµÎÁ¤ÇÇÛÉÛ¤µ¤ì¤Æ¤¤¤Þ¤¹¤¬¡¢¤¢¤Þ¤ê¤Ë¤âÆñ²ò¤Ê¤Î¤Ç¤³¤³¤Î¥Û¡¼¥à¥Ú¡¼¥¸¤Î°Õ¿Þ¤È¤·¤Æ¡¢¤½¤°¤ï¤Ê¤¤¤È»×¤¤¤Þ¤¹¡£
Ãæ³ØÀ¸¤Ç¤â¡¢¶½Ì£¤ò»ý¤Æ¤ëÈϰϤǡ¢¤Ç¤¤ì¤ÐÍý²ò¤Ç¤¤ëÆâÍƤò¡¢¤È¹Í¤¨¤Þ¤·¤¿¡£
»×¤¤¤Ä¤¤¤¿ÊýË¡¤Ê¤Î¤Ç¤¹¤¬¡¢ÀµÊý·Á¤ÎÆâÉô¤ËÆâÀܤ¹¤ë±ß¤òÉÁ¤¡¢¥é¥ó¥À¥à¤Î¤½¤ÎÀµÊý·Á¤ËÅÀ¤òÃÖ¤¤Þ¤¹¡£¤½¤ÎÅÀ¤¬±ß¤Ë¾è¤ë¤«¾è¤é¤Ê¤¤¤«¤ÎÈ椫¤é±ß¼þΨ¤ò·×»»¤Ç¤¤Þ¤¹¡£
¤³¤ì¤ÏÍð¿ô¤ÎÀºÌ©¤µ¤Ç²ò¤Ë±Æ¶Á¤¬½Ð¤Þ¤¹¤Î¤Ç¡¢¤³¤ÎÊÕ¤ê¤Î²ÝÂê¤ò¥É¥é¤¨¤Þ¤ó¤µ¤ó¤¬¸¡Æ¤¤·¤Æ¤ß¤Æ¤¯¤À¤µ¤¤¡£
»ä¤¬»î¤·¤¿¤Î¤Ï°Ê²¼¤Ç¤¹
#include <stdlib.h>
int main(){
int n=0,m=0;
long double j=0,k=0;
while(1){
n++;
j=double(rand())/32767;
k=double(rand())/32767;
if(j*j+k*k<=1)m++;
printf("%.10f %d\n",4*double(m)/double(n),n);
}
return 0;
}
193,700,000¸Ä¤Î¥é¥ó¥À¥à¤ÊÅÀ¤òÃÖ¤¤¤¿¤È¤¡¢3.14153207753¡¡¤È¤Ê¤ê¤Þ¤·¤¿¡£
¤Þ¤¿ÌÌÇò¤½¤¦¤Ê¤â¤Î¤ò»×¤¤¤Ä¤¤¤¿¤é½ñ¤¹þ¤Þ¤»¤Æ¤¤¤¿¤À¤¤Þ¤¹¡£¿ô³Ø¡¢¾ðÊó¤Ë´Ø¤¹¤ëÌÌÇò¤¤¤ªÏ䬤¢¤Ã¤¿¤é¥á¡¼¥ë¤¯¤À¤µ¤¤¤Í¡£