（哈斯克尔）
导入数据。列表（排列，（\\））
a258706 n=a258706列表！！（n-1）
a258706_list=f a000040_list，其中
f ps'@（p:ps）|any（==0）（映射a010051'dps）=f ps
|否则=p:f（ps'\\dps）
其中dps=映射读取$permutations$show p
（PARI）
{A=[2，5]；对于（n=1317，my（D=[1，3，7，9]，r=10^n\9）；对于（A=1，4，对于（b=A^（n<3），4，针对（j=0，if（b！=A，n-1），ispseudoprime（D[A]*r+（D[b]-D[A]）*10^j）||next（2））；A=集合并<A，n-1）]）；A}
是（n）={（n=数字（n））[#n]>=n[1]&#select（d->d，n[^1]-n[^-1]）<2&&！对于（i=1，（#n）^（n[#n]>1），isprime（fromdigits（n=concat（n[^1]，n[1]）））||return）}\\根据约翰逊定理和这里所要求的极小性，数字必须是ab.b或a…ab（=>数字的第一个差至多有1个非零分量）然后就足以考虑数字的旋转。
