どこから左を切っても素数となる日をpythonで探索する

2019年5月23日になりました。

20190523→素数
0190523→素数
190523→素数
90523→素数
0523→素数
523→素数
23→素数
3→素数

どこから左を切っても素数となる珍しい日です！

— もりしー@カレーパンドラ/素数大富豪 (@9973_prm) May 22, 2019

こちらのツイートでそんな考え方もあるんだなぁと感動しました｡ということで次のそんな日をpythonで探索してみます｡

D:\Documents\program\python> python .\prime_yyyymmdd.py
20190823 is full prime day!
20300317 is full prime day!
20360317 is full prime day!
20400307 is full prime day!
20400823 is full prime day!
------------------------------
From 2019/05/24 to 2046/10/09
We find 5 full prime days which means 0.05 %

次が8月23日で案外近くてびっくりしました｡でもその次は10年以上後ですね｡

コードはgithubにもあげておきます｡

github.com

import datetime
import math
import sys

MAX_DAY = 10000


def is_prime(n):
    '''ref to https://qiita.com/srtk86/items/874639e361917e5016d4'''
    if n == 1:
        return False

    for k in range(2, int(math.sqrt(n)) + 1):
        if n % k == 0:
            return False

    return True


def increment_day(dt):
    return dt + datetime.timedelta(days=1)


def find_full_prime_day():
    d = datetime.datetime.now()
    start_day = d.strftime('%Y/%m/%d')
    full_prime_day_cnt = 0

    for _ in range(MAX_DAY):
        n = d.strftime('%Y%m%d')
        while n != '':
            if not is_prime(int(n)):
                break

            n = n[1:]
        else:
            print(f'{d.strftime("%Y%m%d")} is full prime day!')
            full_prime_day_cnt += 1

        d = increment_day(d)

    end_day = d.strftime('%Y/%m/%d')
    full_prime_day_ratio = full_prime_day_cnt / MAX_DAY * 100

    print('-' * 30)
    print(f'From {start_day} to {end_day}')
    print(f'We find {full_prime_day_cnt} full prime days '
          f'which means {full_prime_day_ratio} %')


if __name__ == "__main__":
    args = sys.argv
    if len(args) == 1:
        find_full_prime_day()
    elif len(args) == 2 and str.isdecimal(args[1]):
        MAX_DAY = int(args[1])
        find_full_prime_day()
    else:
        print('Usage: $ python prime_yyymmdd [MAX_DAY=10000]')