莉雁ｹｴ縺ｮ縲�13譌･縺ｮ驥第屆譌･縲阪�菴募屓��
==============================
.. container::
窶� 縺ゅｋ縺��縲∝ヵ縺後�13譌･縺ｮ驥第屆譌･縲阪↓縺､縺�※遏･縺｣縺ｦ縺�ｋ2縲�3縺ｮ縺薙→ 窶�
縺ｾ縺溘�13譌･縺ｮ驥第屆譌･縲阪′繧�▲縺ｦ縺上ｋ縲�
----------------------------------
2015蟷ｴ縺ｮ縺ゅｋ譌･縲∵ｬ｡縺ｮ謇薙■蜷医ｏ縺帙�遒ｺ隱阪〒繧ｫ繝ｬ繝ｳ繝繝ｼ繧堤惻繧√ｋ縺ｨ縲√◎縺ｮ譌･縺後�13譌･縺ｮ驥第屆譌･縲阪↓縺ｪ繧九％縺ｨ縺ｫ豌励′縺､縺�◆縲�
縺ｪ繧薙□縺九▽縺�％縺ｮ髢薙�謇薙■蜷医ｏ縺帙ｂ縲�13譌･縺ｮ驥第屆譌･縲阪□縺｣縺溘ｈ縺�′豌励′縺吶ｋ縲�
縺ｪ縺懊□繧阪≧縺ｨ豌励↓縺ｪ縺｣縺溘�縺ｧ縲√■繧�▲縺ｨ隱ｿ縺ｹ縺ｦ縺ｿ繧医≧縺ｨPython縺ｮNotebook繧帝幕縺�※縺ｿ縺溘�
python繝励Ο繧ｰ繝ｩ繝�繧偵▽縺九▲縺ｦ縲�13譌･縲阪�驥第屆譌･縺ｯ莉悶�譖懈律縺ｫ豈斐∋縺ｦ螟壹＞縺ｮ縺句ｰ代↑縺��縺具ｼ溘ｒ隱ｿ縺ｹ縺ｦ縺ｿ繧九�
----------------------------------------------------------------------------------------------
謇句ｧ九ａ縺ｫ縲≫昴�13譌･縺ｮ驥第屆譌･縲阪′莉悶�譖懈律縺ｫ豈斐∋螟壹＞縺ｨ縺九∝ｰ代↑縺�→縺九＞縺�％縺ｨ縺後≠繧九�縺具ｼ溪�
繧堤｢ｺ縺九ａ縺ｦ縺ｿ繧九％縺ｨ縺九ｉ蛻昴ａ縺ｦ縺ｿ繧医≧縲�
縺ｾ縺壹�縲√＞縺､繧ゅ�縺翫∪縺倥↑縺�→縲∵律莉倥ョ繝ｼ繧ｿ繧貞叙繧頑桶縺�◆繧√�\ ``datetime``\ 繝｢繧ｸ繝･繝ｼ繝ｫ繧偵う繝ｳ繝昴�繝医＠縺ｦ縺翫％縺��
.. code:: ipython3
#!python
#-*- coding:utf-8 -*-
# python2/python3縺ｮ莠呈鋤諤ｧ縺ｮ縺溘ａ縺ｫ
from __future__ import print_function
# 譌･莉倥ョ繝ｼ繧ｿ繧貞叙繧頑桶縺�◆繧√�繝｢繧ｸ繝･繝ｼ繝ｫ繧偵う繝ｳ繝昴�繝医＠縺ｦ縺翫￥縲�
import datetime
貅門ｙ縺悟�譚･縺溘→縺薙ｍ縺ｧm2021蟷ｴ縺九ｉ2121蟷ｴ縺ｾ縺ｧ縺ｮ100蟷ｴ髢薙↓縺ゅｋ豈取怦縺ｮ13縺ｮ譖懈律縺ｮ鬆ｻ蠎ｦ繧呈焚縺井ｸ翫￡縺ｦ縺ｿ繧後�縺�＞縺�繧阪≧縲�
.. code:: ipython3
#縲譛域屆縺九ｉ譌･譖懊�譖懈律縺ｮ蜷榊燕繧貞ｮ夂ｾｩ縺励※縺翫″縺ｾ縺吶�
wdname_e=("mon","tue","wed","thu","fri","sat","sun")
# 譌･譛ｬ隱槭〒縺ｮ蜷榊燕
wdname_j=(u"譛�",u"轣ｫ",u"豌ｴ",u"譛ｨ", u"驥�", u"蝨�", u"譌･")
wdname=wdname_j
wdcount={} #譖懈律豈弱�蝗樊焚繧定ｦ壹∴縺ｦ縺翫￥縺溘ａ縺ｮ霎樊嶌蝙九ョ繝ｼ繧ｿ繧堤畑諢上☆繧九�
ny=100 #100蟷ｴ蛻�ｒ隱ｿ縺ｹ縺ｦ縺ｿ縺ｾ縺吶�
for y in range(2021,2021+ny): #range(2021,2022)->(2021,2022,...,2120)
for m in range(1,13,1): #range(1,13,1) -> (1,2,...,12)
d=datetime.datetime(y,m,13) # y蟷ｴm譛�13譌･縺ｮ譌･莉倥ョ繝ｼ繧ｿ繧剃ｽ懈�縲�
wd=d.weekday() # 縺昴�譌･縺ｮ譖懈律繧貞叙繧雁�縺吶�
wdcount[wd] = wdcount.get(wd,0) + 1 #wdcount縺ｮkey縺詣d縺ｮ繝��繧ｿ縺ｫ1繧貞刈縺医ｋ縲�
print("譖懈律","蝗樊焚" )
for k in range(7): # 譖懈律豈弱�逋ｻ蝣ｴ蝗樊焚繧貞魂蛻ｷ縺吶ｋ縲�
print (wdname[k], wdcount[k])
.. parsed-literal::
譖懈律 蝗樊焚
譛� 171
轣ｫ 172
豌ｴ 172
譛ｨ 170
驥� 171
蝨� 172
譌･ 172
縺薙ｌ繧峨�蜊倡ｴ斐↑蟷ｳ蝮�､ $ 171.4=
100:raw-latex:`\times`:raw-latex:`\frac{12}{7}`$ 縺ｫ荳閾ｴ縺励※縺�ｋ縲�
蠖薙◆繧雁燕縺�縺代←縲∫音縺ｫ13譌･縺ｯ驥第屆譌･縺悟､壹＞縺ｨ縺九√◎縺ｮ騾�→縺�≧繧上￠縺ｧ縺ｯ縺ｪ縺�ゅ％
縺ｮ蛻�ｸ�ｒ譖懈律豈弱�蝗樊焚繧堤､ｺ縺吶ヰ繝ｼ繧ｰ繝ｩ繝輔↓縺吶ｌ縺ｰ縲√ｏ縺九ｊ繧�☆縺�□繧阪≧縺九�
繧ｰ繝ｩ繝輔ｒ譖ｸ縺�※縺ｿ繧九�
--------------------
繧ｰ繝ｩ繝輔ｒ謠上￥縺溘ａ縺ｮ繝｢繧ｸ繝･繝ｼ繝ｫ縺ｯ縺�￥縺､縺九≠繧九′縲√％繧後ｂ縺�▽繧ゅ�繧医≧縺ｫ\ ``matplotlib``\ 繧剃ｽｿ縺�％縺ｨ縺ｫ縺吶ｋ縲�
``matplotlib``\ 縺ｮ繧ｵ繝悶Δ繧ｸ繝･繝ｼ繝ｫ\ ``pyplot``\ 縺ｯ邁｡蜊倥↑繧ｰ繝ｩ繝輔�菴懈�縺ｫ縺ｯ繝斐ャ繧ｿ繝ｪ縺�縲�
邏ｰ驛ｨ縺ｫ縺ｾ縺ｧ蜃昴ｋ縺ｪ繧噂 ``matplotlib``\ 縺ｮ繧ｪ繝悶ず繧ｧ繧ｯ繝医ｒ菴ｿ縺�％縺ｪ縺吶％縺ｨ縺ｫ縺ｪ繧九�
.. code:: ipython3
import matplotlib.pyplot as pyplot #matplotli.pyplot 繧� pyplot縺ｨ縺�≧蜷榊燕縺ｧ菴ｿ縺�％縺ｨ縺ｫ縺吶ｋ縲�
pyplot.bar(wdcount.keys(), wdcount.values(),align="center") #繝舌�繧ｰ繝ｩ繝輔ｒ菴ｿ縺｣縺ｦ繧ｰ繝ｩ繝募喧縺吶ｋ縲�
# pyplot.bar(*zip(*wdcount.items()),align="center") #縺薙≧縺�≧譖ｸ縺肴婿繧ゅ≠繧九�
pyplot.xticks(range(7), wdname_j, font="IPAexGothic") # x霆ｸ縺ｮ繝ｩ繝吶Ν繧呈屆譌･縺ｫ縺励※縺翫％縺��
pyplot.draw()
.. image:: Friday13th_files/Friday13th_7_0.png
豈主ｹｴ縺ｮ縲�13譌･縺ｮ驥第屆譌･縲阪�蝗樊焚繧呈焚縺医※縺ｿ繧医≧縲�
--------------------------------------------
縺薙�繧医≧縺ｫ縲∵ｯ主ｹｴ縺ｮ縲�13譌･縺ｮ驥第屆譌･縲阪�豈主ｹｴ 12/7=1.7蝗�
縺ゅｋ繧上￠縺�縺後�1蝗槭�蟷ｴ繧ゅ≠繧後�縲�2蝗槭�蟷ｴ繧ゅ≠繧九□繧阪≧縲� 繧ゅ≧荳蠎ｦPython
繧ｹ繧ｯ繝ｪ繝励ヨ
繧偵▽縺九▲縺ｦ縲∵ｯ主ｹｴ縺ｮ縲�13譌･縺ｮ驥第屆譌･縲阪�蝗樊焚繧呈焚縺医※遒ｺ縺九ａ縺ｪ縺代ｌ縺ｰ縺ｪ繧峨↑縺��
.. code:: ipython3
F13inY={} #蟷ｴ豈弱�縲�13譌･縺ｮ驥第屆譌･縲阪�蝗樊焚繧定ｨ倬鹸縺吶ｋ霎樊嶌蝙九ョ繝ｼ繧ｿ`F13inY`繧堤畑諢上☆繧九�
def count_F13inY(year):
count=0 # 繧ｫ繧ｦ繝ｳ繧ｿ count繧�0繧ｯ繝ｪ繧｢縺吶ｋ縲�
for m in range(1,13,1): # 豈取怦縺ｮ
d=datetime.datetime(year,m,13) #13譌･縺�
wd=d.weekday() # 譖懈律縺�
if wd == 4:# 驥第屆譌･(4)縺�縺｣縺溘ｉ
count +=1 # 繧ｫ繧ｦ繝ｳ繧ｿ繧�1蠅励ｄ縺吶�
return count
for y in range(2015,2115):
count=count_F13inY(y)
F13inY[count] = F13inY.get(count,0)+1 # 荳蟷ｴ縺ゅ◆繧翫�縲�13譌･縺ｮ驥第屆譌･縲阪�蝗樊焚縺慶ount縺ｨ縺ｪ繧句ｹｴ謨ｰ繧�1蠅励ｄ縺�
for k, v in F13inY.items(): # 荳蟷ｴ縺ゅ◆繧翫�縲�13譌･縺ｮ驥第屆譌･縲阪�蝗樊焚縺渓縺ｮ縺ｪ繧句ｹｴ謨ｰ繧貞魂蛻ｷ縺吶ｋ縲�
print (k,"occurense ", v,"times")
.. parsed-literal::
3 occurense 16 times
1 occurense 43 times
2 occurense 41 times
繝励Ο繧ｰ繝ｭ繝ｩ繝�縺ｮ邨先棡縺ｧ縺ｯ縲�2015-2114縺ｮ100蟷ｴ髢薙〒縲√�13譌･縺ｮ驥第屆譌･縲阪′1譌･縺励°辟｡縺�ｹｴ縺�43蝗槭�2譌･縺ゅｋ蟷ｴ縺�42蝗槭�
縺昴＠縺ｦ縺ｪ繧薙→縲悟ｹｴ縺ｫ3譌･縺ｮ縲�13譌･縺ｮ驥第屆譌･縲阪ｒ謖√▽縲榊ｽ薙◆繧雁ｹｴ縺�16蝗槭≠繧九�
縺ｾ縺ゅ√％縺ｮ蟷ｳ蝮��遒ｺ縺九↓縲�
.. code:: ipython3
(43+2*42+3*15)/100.
.. parsed-literal::
1.72
縺ｨ縺ｪ繧翫� 蟷ｳ蝮�☆繧後�蝮�ｭ峨↑遒ｺ邇�↓繧医ｋ邨先棡騾壹ｊ縺ｫ縺ｪ縺｣縺ｦ縺ｯ縺�ｋ繧薙□縺代ｌ縺ｩ縲�
蟷ｴ縺ｫ3蝗槭�縲�13譌･縺ｮ驥第屆譌･縲阪′縺ゅｋ蟷ｴ繧呈爾縺励※縲�
-------------------------------------------
縺昴ｌ縺ｧ縺ｯ縺ｨ縺薙�100蟷ｴ縺ｮ縺�■縲�13譌･縺ｮ驥第屆譌･縲阪′3蝗樊怏繧九�縺ｯ縺ｩ縺ｮ蟷ｴ縺九ｒ隱ｿ縺ｹ縺ｦ縺ｿ縺溘￥縺ｪ繧九�
谺｡縺ｮ繝励Ο繧ｰ繝ｩ繝�繧剃ｽｿ縺�→縲�2015-1214蟷ｴ縺ｮ100蟷ｴ髢薙↓縺ｯ縲�
2015, 2026, 2037, 2040, 2043, 2054, 2065, 2068, 2071, 2082, 2093, 2096,
2099, 2105, 2108, 縺昴＠縺ｦ 2111
縺ｮ16蝗槭〒縺ゅｋ縺薙→縺後ｏ縺九ｋ縲�
遒ｺ縺九↓莉雁ｹｴ(2015)縺ｯ3蝗槭�縲�13譌･縺ｮ驥第屆譌･縲阪′縺ゅｋ蟷ｴ縺ｧ縲√後∪縺溘°縲阪→諤昴▲縺溘�繧ゅ◎縺ｮ縺溘ａ縺�繧阪≧縲�
縺薙�繝励Ο繧ｰ繝ｩ繝�繧偵▽縺九∴縺ｰ縲∵ｬ｡縺ｫ縲悟ｹｴ縺ｫ3蝗槭�縲�13譌･縺ｮ驥第屆譌･縲上阪′縺ゅｋ谺｡縺ｮ蟷ｴ縺後ｏ縺九▲縺ｦ縺励∪縺��
縺昴ｌ縺ｯ縺ｪ繧薙→2026蟷ｴ縺�縺ｨ繝励Ο繧ｰ繝ｩ繝�縺ｯ謨吶∴縺ｦ縺上ｌ繧九�
蠖灘�縺ｯ縲後∴縺｣縲√∪縺溘�13譌･縺ｮ驥第屆譌･縲上阪→諢溘§繧区ｧ倥↑蟷ｴ縺ｯ蟾｡縺｣縺ｦ縺薙↑縺�ｏ縺代□縲�
.. code:: ipython3
years=[]
for y in range(2015,2115):
count=count_F13inY(y)
if count >=3:
print ("year",y, "has ",count, "day(s) of \"Friday 13th\"")
years.append(y)
print(len(years), years)
.. parsed-literal::
year 2015 has 3 day(s) of "Friday 13th"
year 2026 has 3 day(s) of "Friday 13th"
year 2037 has 3 day(s) of "Friday 13th"
year 2040 has 3 day(s) of "Friday 13th"
year 2043 has 3 day(s) of "Friday 13th"
year 2054 has 3 day(s) of "Friday 13th"
year 2065 has 3 day(s) of "Friday 13th"
year 2068 has 3 day(s) of "Friday 13th"
year 2071 has 3 day(s) of "Friday 13th"
year 2082 has 3 day(s) of "Friday 13th"
year 2093 has 3 day(s) of "Friday 13th"
year 2096 has 3 day(s) of "Friday 13th"
year 2099 has 3 day(s) of "Friday 13th"
year 2105 has 3 day(s) of "Friday 13th"
year 2108 has 3 day(s) of "Friday 13th"
year 2111 has 3 day(s) of "Friday 13th"
16 [2015, 2026, 2037, 2040, 2043, 2054, 2065, 2068, 2071, 2082, 2093, 2096, 2099, 2105, 2108, 2111]
繝励Ο繧ｰ繝ｩ繝�繧堤洒縺上☆繧九◆繧√↓縺薙ｓ縺ｪ譖ｸ縺肴婿繧ゅ〒縺阪ｋ繧峨＠縺��
.. code:: ipython3
years=[y for y in range(2015,2115) if count_F13inY(y)>=3 ] # 3蝗�(莉･荳�)縺ｮ 縲�13譌･縺ｮ驥第屆譌･縲阪′縺ゅｋ蟷ｴ縺ｮ繝ｪ繧ｹ繝�
for y in years:
print ("year",y, "has ","more than 3 days of \"Friday 13th\"")
print()
print(len(years), years)
.. parsed-literal::
year 2015 has more than 3 days of "Friday 13th"
year 2026 has more than 3 days of "Friday 13th"
year 2037 has more than 3 days of "Friday 13th"
year 2040 has more than 3 days of "Friday 13th"
year 2043 has more than 3 days of "Friday 13th"
year 2054 has more than 3 days of "Friday 13th"
year 2065 has more than 3 days of "Friday 13th"
year 2068 has more than 3 days of "Friday 13th"
year 2071 has more than 3 days of "Friday 13th"
year 2082 has more than 3 days of "Friday 13th"
year 2093 has more than 3 days of "Friday 13th"
year 2096 has more than 3 days of "Friday 13th"
year 2099 has more than 3 days of "Friday 13th"
year 2105 has more than 3 days of "Friday 13th"
year 2108 has more than 3 days of "Friday 13th"
year 2111 has more than 3 days of "Friday 13th"
16 [2015, 2026, 2037, 2040, 2043, 2054, 2065, 2068, 2071, 2082, 2093, 2096, 2099, 2105, 2108, 2111]
縺薙�邨先棡繧偵ｈ縺上∩繧九→縲�2015-2065縺ｾ縺ｧ縺ｯ縲√�13譌･縺ｮ驥第屆譌･縺悟ｹｴ縺ｫ3蝗槭≠繧句ｹｴ縲阪�螟ｧ菴�11蟷ｴ縺ゅｋ縺��13蟷ｴ縺ｮ髢薙ｒ縺翫＞縺ｦ迴ｾ縺ｦ縺�ｋ縲�
縺励°縺�(2065,2068,2071) 縺翫ｈ縺ｳ(2093,2096, 2099)
縺ｯ荳牙ｹｴ縺翫″縺ｮ髢馴囈縺ｧ窶懊�13譌･縺ｮ驥第屆譌･縺悟ｹｴ縺ｫ3蝗槭≠繧句ｹｴ窶昴′郢ｰ繧願ｿ斐＠縺ｦ縺�ｋ縲�
縺薙ｌ縺ｯ驕主悉縺ｮ郢ｰ繧願ｿ斐＠繧りｪｿ縺ｹ縺ｦ縺ｿ繧句ｿ�ｦ√′縺ゅｋ縺�繧阪≧縲�
.. code:: ipython3
years=[]
for y in range(1915,2016):
count=count_F13inY(y)
if count >=3:
print ("year",y, "has ",count, "day(s) of \"Friday 13th\"")
years.append(y)
print(len(years), years)
.. parsed-literal::
year 1925 has 3 day(s) of "Friday 13th"
year 1928 has 3 day(s) of "Friday 13th"
year 1931 has 3 day(s) of "Friday 13th"
year 1942 has 3 day(s) of "Friday 13th"
year 1953 has 3 day(s) of "Friday 13th"
year 1956 has 3 day(s) of "Friday 13th"
year 1959 has 3 day(s) of "Friday 13th"
year 1970 has 3 day(s) of "Friday 13th"
year 1981 has 3 day(s) of "Friday 13th"
year 1984 has 3 day(s) of "Friday 13th"
year 1987 has 3 day(s) of "Friday 13th"
year 1998 has 3 day(s) of "Friday 13th"
year 2009 has 3 day(s) of "Friday 13th"
year 2012 has 3 day(s) of "Friday 13th"
year 2015 has 3 day(s) of "Friday 13th"
15 [1925, 1928, 1931, 1942, 1953, 1956, 1959, 1970, 1981, 1984, 1987, 1998, 2009, 2012, 2015]
1915-2015縺ｮ髢薙〒縺ｯ縲�(1925, 1928, 1931), (1953, 1956,
1959)縺翫ｈ縺ｳ(2009,2012,2015)繧ゅ御ｸ牙ｹｴ縺翫″縺ｮ髢馴囈縺ｧ縲�13譌･縺ｮ驥第屆譌･縺悟ｹｴ縺ｫ3蝗槭≠繧句ｹｴ縺檎ｹｰ繧願ｿ斐☆縲肴凾譛溘↓蠖薙◆縺｣縺ｦ縺�ｋ縺薙→縺悟愛譏弱☆繧九�
縺昴�諢丞袖縺ｧ繧�2015蟷ｴ縺ｯ縺｡繧�▲縺ｨ縺�縺醍音蛻･縺ｪ蟷ｴ縺ｧ縺ゅ▲縺溘ｏ縺代□縲�
谺｡縺ｮ縲御ｸ牙ｹｴ縺翫″縺ｮ髢馴囈縺ｧ縲�13譌･縺ｮ驥第屆譌･縺悟ｹｴ縺ｫ3蝗槭≠繧句ｹｴ縺檎ｹｰ繧願ｿ斐☆縲肴悄髢薙�2037蟷ｴ繝ｼ2043蟷ｴ縺ｨ縺�≧縺薙→縺ｫ縺ｪ繧九�
縺昴�鬆�↓縺ｯ縺ｾ縺溘≫晄怙霑�13譌･縺ｮ驥第屆譌･縺悟､壹＞縺ｪ窶昴→諤昴▲縺溘ｊ縺吶ｋ縺ｮ縺�繧阪≧縺九�
�育ｵ�)
縺翫∪縺�
------
縲�1蟷ｴ髢薙↓縺ゅｋ13譌･縺ｮ驥第屆譌･縺ｮ譌･謨ｰ縲上�蛻�ｸ��
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
繧ｰ繝ｬ繧ｴ繝ｪ繧ｪ豁ｴ縺ｯ400蟷ｴ縺ｧ縺｡繧�≧縺ｩ縲�1460972譌･縲�20870騾ｱ縺ｨ縺ｪ繧九�縺ｧ縲∝�蟶��繝代ち繝ｼ繝ｳ縺ｯ400蟷ｴ縺ｮ蜻ｨ譛溘ｒ謖√▽縲�
繧ゅ＠縲√％縺ｮ譌･謨ｰ縺後■繧�≧縺ｩ7縺ｧ蜑ｲ繧雁�繧後ｋ縺薙％縺ｧ閠�∴縺ｦ縺�ｋ繧医≧縺ｪ譖懈律縺ｮ郢ｰ繧願ｿ斐＠縺ｯ2800蟷ｴ縺ｮ蜻ｨ譛溘ｒ繧ゅ▲縺ｦ縺�◆縺�繧阪≧縲�
縺薙�400蟷ｴ縺ｮ蜻ｨ譛溘�髢薙↓縲�1蟷ｴ髢薙↓縺ゅｋ13譌･縺ｮ驥第屆譌･縺ｮ譌･謨ｰ縲上�蛻�ｸ�ｒ繧ｰ繝ｩ繝輔↓縺励※縺ｿ繧医≧縲�
.. code:: ipython3
## %%timeit # 螳溯｡梧凾髢捺ｸｬ螳壹�縺溘ａ縺ｮ縺翫∪縺倥↑縺�(cell magic)
ymax=401
# 478 ms ﾂｱ 7.68 ms per loop (mean ﾂｱ std. dev. of 7 runs, 1 loop each) on Intel Mac with ymax=9999
# 220 ms ﾂｱ 8.01 ms per loop (mean ﾂｱ std. dev. of 7 runs, 1 loop each) on M1 Mac with arc -x86_64 and ymax=9999
F13inY={}
for y in range(1,ymax):
count=count_F13inY(y)
F13inY[count] = F13inY.get(count,0)+1
for k in F13inY:
print (k, "occurense ", F13inY[k], "times")
years=range(1,ymax)
pyplot.figure(figsize=(16,9))
pyplot.bar(years,[count_F13inY(y) for y in years])
pyplot.draw()
.. parsed-literal::
2 occurense 170 times
1 occurense 171 times
3 occurense 59 times
.. image:: Friday13th_files/Friday13th_21_1.png
繧ｰ繝ｬ繧ｴ繝ｪ繧ｪ豁ｴ縺ｧ縺ｯ縲�4蟷ｴ縺ｫ荳蝗槭≠繧矩柾蟷ｴ縺後�100蟷ｴ縺ｫ荳蝗槭↑縺上↑縺｣縺ｦ縺励∪縺��
縺薙�蜉ｹ譫懊ｒ隕九ｋ縺溘ａ縺ｫ縲�100蟷ｴ豈弱↓蛹ｺ蛻�▲縺ｦ蛻�ｸ�ｒ繝励Ο繝�ヨ縺励※縺ｿ繧医≧縲�
.. code:: ipython3
ny=100
years=range(1,1+ny)
pyplot.figure(figsize=(20,5))
for n in range(1,401,100):
pyplot.bar(years,[count_F13inY(y) for y in range(n,n+ny)])
pyplot.draw()
.. image:: Friday13th_files/Friday13th_23_0.png
100蟷ｴ
蜷後§譖懈律繧呈戟縺､譛�
~~~~~~~~~~~~~~~~
豈取怦縺ｮ譌･謨ｰ縺ｯ縲�柾蟷ｴ縺ｮ2譛医ｒ髯､縺�※縺ｯ縲√＞縺､繧ょ酔縺倥ょｾ薙▲縺ｦ縲�1譛医�縺ゅｋ譌･莉倥�譖懈律縺ｨ莉悶�譛医�蜷後§譌･莉倥�譖懈律縺ｯ荳螳壹�髢｢菫ゅ′縺ゅｋ縲�
髢丞ｹｴ縺ｧ縺ｯ縺ｪ縺�√≠繧句ｹｴ縺ｮ豈取怦縺ｮ13譌･縺ｮ譖懈律繧定ｪｿ縺ｹ縺ｦ隕九ｋ縲�
.. code:: ipython3
y=2021;print([(i,wdname_j[datetime.date(y,i,13).weekday()]) for i in range(1,13)])
.. parsed-literal::
[(1, '豌ｴ'), (2, '蝨�'), (3, '蝨�'), (4, '轣ｫ'), (5, '譛ｨ'), (6, '譌･'), (7, '轣ｫ'), (8, '驥�'), (9, '譛�'), (10, '豌ｴ'), (11, '蝨�'), (12, '譛�')]
縺薙�繧医≧縺ｫ縲�柾蟷ｴ縺ｧ縺ｪ縺�ｹｴ縺ｯ ﾂ���1譛医�10譛茨ｼ�, ��2譛�, 3譛�, 11譛茨ｼ峨�
(4譛医�7譛茨ｼ峨� (9譛�,12譛茨ｼ� 縺後◎繧後◇繧悟酔縺俶屆譌･縺ｫ縺ｪ繧九�
2譛医′縺｡繧�≧縺ｩ蝗幃ｱ髢難ｼ�28譌･縺ｪ縺ｮ縺ｧ縲�2譛医→3譛医�1-28譌･縺ｮ譖懈律縺檎ｹｰ繧願ｿ斐☆縺ｮ縺ｯ繧医￥繧上°繧九�
.. code:: ipython3
pyplot.hist([datetime.date(y,i,13).weekday() for i in range(1,13)],bins=range(0,8));
.. parsed-literal::
(array([2., 2., 2., 1., 1., 3., 1.]),
array([0, 1, 2, 3, 4, 5, 6, 7]),
<BarContainer object of 7 artists>)
.. image:: Friday13th_files/Friday13th_28_1.png
縺ｧ縺ｯ髢丞ｹｴ縺ｫ縺ｯ縺ｩ縺�↑繧九□繧阪≧縺具ｼ溘2譛医′29譌･縺ゅｋ縺ｮ縺ｧ縲�3譛医�譖懈律縺ｯ2譛医�蜷後§譌･莉倥°繧我ｸ縺､繧ｺ繝ｬ繧九�縺壹□縲�
.. code:: ipython3
y=2020;print([(i,wdname_j[datetime.date(y,i,13).weekday()]) for i in range(1,13)])
.. parsed-literal::
[(1, '譛�'), (2, '譛ｨ'), (3, '驥�'), (4, '譛�'), (5, '豌ｴ'), (6, '蝨�'), (7, '譛�'), (8, '譛ｨ'), (9, '譌･'), (10, '轣ｫ'), (11, '驥�'), (12, '譌･')]
縺ｨ縺�≧繧上￠縺ｧ縲�柾蟷ｴ縺ｫ縺ｯ縲�(1譛医�4譛�,7譛茨ｼ�, (2譛医�8譛茨ｼ�, (3譛医�11譛�),
(9譛医�12譛�) 縺後◎繧後◇繧悟酔縺俶屆譌･縺ｨ縺ｪ繧九％繧上￠縺�縲�
縺､縺�〒縺ｫ縲√←縺ｮ蟷ｴ縺ｧ繧�
3譛医→11譛医�4譛医→7譛医�9譛医→12譛医�縺昴ｌ縺槭ｌ,31譌･繧帝勁縺代�縲∝酔縺俶律莉倥�譌･縺ｯ蜷後§譖懈律縺ｨ縺ｪ繧九�
.. code:: ipython3
pyplot.hist([datetime.date(y,i,13).weekday() for i in range(1,13)],bins=range(0,8));
.. parsed-literal::
(array([3., 1., 1., 2., 2., 1., 2.]),
array([0, 1, 2, 3, 4, 5, 6, 7]),
<BarContainer object of 7 artists>)
.. image:: Friday13th_files/Friday13th_32_1.png
譛蠕後↓譛医�31譌･縺ｮ譖懈律縺ｮ髢｢菫ゅｒ縺ｿ縺ｦ縺ｿ繧九→縲�
髢丞ｹｴ縺ｮ蟷ｴ縺ｫ縺ｯ縲�
.. code:: ipython3
y=2020;print([(i,wdname_j[datetime.date(y,i,31).weekday()]) for i in (1,3,5,7,8,10,12)])
.. parsed-literal::
[(1, '驥�'), (3, '轣ｫ'), (5, '譌･'), (7, '驥�'), (8, '譛�'), (10, '蝨�'), (12, '譛ｨ')]
髢丞ｹｴ縺ｧ縺ｪ縺�ｹｴ縺ｫ縺ｯ縲�
.. code:: ipython3
y=2021;print([(i,wdname_j[datetime.date(y,i,31).weekday()]) for i in (1,3,5,7,8,10,12)])
.. parsed-literal::
[(1, '譌･'), (3, '豌ｴ'), (5, '譛�'), (7, '蝨�'), (8, '轣ｫ'), (10, '譌･'), (12, '驥�')]
縺ｨ縺ｪ繧九�縺ｧ縲�31譌･縺ｯ縲√朱柾蟷ｴ縺ｫ縺ｯ縲∽ｸ譛医→荳�怦縲√◎繧御ｻ･螟悶�蟷ｴ縺ｫ縺ｯ荳譛医→蜊∵怦縺悟酔縺俶屆譌･縺ｨ縺ｪ繧九ゅ上％縺ｨ縺後ｏ縺九ｋ縲�
蜈ｷ菴鍋噪縺ｫ縺ｯ縲�1螟ｪ髯ｽ蟷ｴ��365.24219040 - 6.15ﾃ�10-6 T
[譌･]縲ゅ◆縺�縺励ゝ縺ｯ繝ｦ繝ｪ繧ｦ繧ｹ荳也ｴ謨ｰ.繧ｰ繝ｬ繧ｴ繝ｪ繧ｪ證ｦ縺九ｉ縺ｮ縺壹ｌ縺ｯ10000蟷ｴ縺ｫ3譌･蛻�洒縺��
::
>>> 365+1/4-1/100+1/400-1/3200+1/80000
365.2422
縺�縺九ｉ縲�3200=8*400
蟷ｴ縺ｫ荳蝗槭＆繧峨↓髢丞ｹｴ繧呈ｸ帙ｉ縺�,80000=25*3200=200*400蟷ｴ縺ｮ蛟肴焚縺ｯ髢丞ｹｴ縺ｨ縺吶ｌ縺ｰ縲√■繧�≧縺ｩ縺�＞縲�
繧ゅ▲縺ｨ繧ゅ◎繧後∪縺ｧ蟷ｳ蝮�､ｪ髯ｽ譌･縺檎樟蝨ｨ縺ｮ縺ｾ縺ｾ縺ｨ縺�≧縺薙→縺ｯ縺ｪ縺��縺ｧ縺薙�繝ｫ繝ｼ繝ｫ繧ゅ≠縺ｾ繧頑э蜻ｳ縺ｯ縺ｪ縺�□繧阪≧縲�