mathics.m

(* Some prime stuff in Mathics, an open source mathematica clone with pythonic influence *)
In[1]:= 2^255-19
Out[1]= 57896044618658097711785492504343953926634992332820282019728792003956564819949

In[2]:= 2^255+95
Out[2]= 57896044618658097711785492504343953926634992332820282019728792003956564820063

In[3]:= PrimeQ[2^255-19]
Out[3]= True

In[4]:= PrimeQ[2^255+Range[-19,20]]
Out[4]= {True, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False, False}

In[5]:= Select[Transpose[{PrimeQ[2^255+Range[-19,10000]], Range[-19,10000]}], First]
Out[5]= {{True, -19}, {True, 95}, {True, 141}, {True, 275}, {True, 333}, {True, 443}, {True, 491}, {True, 539}, {True, 611}, {True, 821}, {True, 1073}, {True, 1109}, {True, 1131}, {True, 1155}, {True, 1241}, {True, 1269}, {True, 1271}, {True, 1535}, {True, 1625}, {True, 1661}, {True, 1719}, {True, 2705}, {True, 2819}, {True, 2919}, {True, 3225}, {True, 3885}, {True, 3921}, {True, 4253}, {True, 4355}, {True, 4629}, {True, 4793}, {True, 4809}, {True, 5163}, {True, 5345}, {True, 5409}, {True, 5429}, {True, 5481}, {True, 5663}, {True, 5729}, {True, 5955}, {True, 6135}, {True, 6191}, {True, 6269}, {True, 6501}, {True, 6525}, {True, 6581}, {True, 6701}, {True, 6935}, {True, 7031}, {True, 7395}, {True, 7523}, {True, 7565}, {True, 7821}, {True, 7959}, {True, 8055}, {True, 8363}, {True, 8495}, {True, 8573}, {True, 8739}, {True, 8769}, {True, 8811}, {True, 9189}, {True, 9305}, {True, 9389}, {True, 9405}, {True, 9611}}