数学 - Python - 素数について(SymPy モジュール、isprime 関数、itertools、count メソッド) - ( @Servantprime )MathPoppy⌬ さんのツイートより

Pythonからはじめる数学入門
原著

学習環境

• Surface 3 (4G LTE)、Surface 3 タイプ カバー、Surface ペン(端末)
• Windows 10 Pro (OS)
• Nebo(Windows アプリ)
• iPad Pro + Apple Pencil
• MyScript Nebo(iPad アプリ)
• 参考書籍
  • Pythonからはじめる数学入門
  • JavaScriptリファレンス 第6版
  • インタラクティブ・データビジュアライゼーション

【今日のトリビア】
1,2,3,4…と数字を順に並べていくと、217で素数になる。
ほえー！そすーそすー！(546桁の素数) pic.twitter.com/osz8VSXi00

— MathPoppy⌬ (@Servantprime) 2017年12月23日

SymPy(Python)で1、2、3、4…と数字を順に並べていくと、217で素数になるか確認してみた。(ついでに桁数も指定できるようにしてみた。)

コード(Emacs)

Python 3

#!/usr/bin/env python3
from sympy import pprint, isprime
import sys
import itertools

if len(sys.argv) > 1:
    m = int(sys.argv[1])
else:
    m = 547
s = ''
for n in itertools.count(1):
    s = f'{s}{n}'
    if isprime(int(s)):
        print(f'{s}: {len(s)}桁')
    if len(s) >= m:
        print(f'{len(s)}桁で終了\n{s}')
        break

入出力結果(Terminal, Jupyter(IPython))

$ ./sample.py 550
552桁で終了
123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220
$

ということで、MathPoppy⌬さんに言うような素数にはならなかった。上記のコードに何か間違いがあるのか、SymPyのisprime 関数の問題なのか、MathPoppy⌬さんの計算ミスなのか…