f[0]=0;f[n_Integer]:=10^(楼层[Log[10,n]]+1)*f[n-1]+n;Do[If[PrimeQ[FromDigits[Take[IntegerDigits[f[n],n]]],打印[n]],{n,1,3000}]
案例[FromDigits/@Rest[FoldList[Append,{},RealDigits[N[ChampernowneNumber[],1000]][[1]]],p_?PrimeQ:>整数长度[p]](*埃里克·韦斯特因2015年11月4日*)