（*OEIS序列A088628的Mathematica程序，Ray Chandler，2014年4月8日*）nmax=40；a={11227}；ds={1,2}；总和=3；做[id=整数位数[n]；ds=连接[ds，id]；sumd+=总计[id]；wds=静止[ds]；（*删除首字母1*）如果[Mod[sumd，3]==0，AppendTo[wds，1]]；wds=前缀[Sort[wds]，1]；（*恢复初始值1*）il=长度[wds]；找到=0；k=计数[wds，Last[wds]]+1；While期间[pre=取[wds，il-k]；post=取[wds，-k]；ps=排列[后]；执行[c=起始数字[Join[pre，p]]；如果[PrimeQ[c]，则找到=1；中断[]]；，{p，ps}]；如果[found==1，AppendTo[a，c]]；找到==0，k++；];，{n，3，nmax}]；一