%I#19 2015年8月6日04:17:24

%S 10000103100001211000022310000330310000511100060910000721，

%电话：100007231000081310008191000083110001009100011071119，

%电话：1000120310001207100012091000121310001221100012271000123110010101100113100102131001227100110311100103171001032310010401

%N素数天数：可以解释为ISO日期的素数。

%C ISO日期的格式为YYYY-MM-DD。

%C此序列的最后一项是99991207。

%H Christian Perfect，n表，n的a（n）=1..24974</a>

%H<a href=“http://www.mei.org.uk/month_item网站“>MEI每月数学项目。

%H<a href=“http://www.iso.org/iso/home/standards/iso8601.htm“>日期和时间格式-ISO 8601。

%e 20150821属于此序列，因为它是首要的，代表2015年8月21日。

%t d=排序@Join[Flatten[Plus[{100，300，500，700，800，1000，1200}，#]&/@选择[Range[1，31，2]！Mod[#，5]==0&]]，Flatten[Plus[{200，400，600，900，1100}，#]&/@选择[Range[1，29，2]！型号[#，5]==0&]]]；选择[Flatten[Plus[If[Or[Mod[#，10000000]==0，Mod[#,40000]==0&&！Mod[#，1000000]==0]，d，Delete[d，25]]，#]&/@（10000 Range[10001001]）]，PrimeQ]（*_Michael De Vlieger_，2015年8月4日*）

%o（Python）

%o导入日期时间

%o、。

%o定义prime_number_days（）：

%o。。。j=0

%o。。。date=datetime.date（1000,1,1）

%o。。。当日期年份<2100时：

%o。。。。。。x=int（date.strftime（“%Y%m%d”））

%o。。。。。。而素数[j]<x:

%o。。。。。。。。。j+=1

%哦。。。。。。如果素数[j]==x：

%o。。。。。。。。。产量x

%哦。。。。。。date+=datetime.timedelta（days=1）

%Y参考A000040（质数）。

%K nonn、easy、fini、base

%O 1,1号机组

%《完美基督》，2015年8月4日