%我#21 2023年9月19日12:39:32

%S 7,11,17,41,47,711117117114471471174117471777411141774441，

%电话：44477177741174177477771777411111711171111771141111447，

%电话：1147111717117771411771441114447147171474114747147711711717417

%N仅使用直数位1、4和7的质数。

%C（n）的小数位数永远不能被3整除_罗伯特·伊斯雷尔，2014年5月22日

%C只使用所有三个直数位的最小素数是a（9）=1447（见prime Courios！链接）_伯纳德·肖特，2023年9月8日

%H Robert Israel，n的表，a（n）表示n=1..10000</a>

%H Chris K.Caldwell和G.L.Honaker，Jr.，<a href=“https://t5k.org/curios/page.php？short=1447“>1447，顶级古玩！[古普塔]

%e17是一个术语，因为它是一个素数，仅由直数位1和7组成。

%p f：=过程（x）局部n，d，t，i，a；

%p n：=楼层（对数[3]（（2*x+3））；

%p如果n mod 3=0，则返回0 fi；

%p d：=x-（3^n-3）/2；

%p t：=0；

%i从0到n-1的p do

%pa:=d模块3；

%p t：=t+（3*a+1）*10^i；

%p d:=（d-a）/3；

%日期：

%价格

%p端程序：

%p选择（isprime，map（f，[$1..1000]）；#_罗伯特·伊斯雷尔，2014年5月22日

%t选择[Prime[Range[2000]]，Union[Join[IntegerDigits[#]，{1，4，7}]=={1，4,7}&]

%o（PARI）直线（n）=我的（t）；while（n，t=n%10；如果（t！=1&t！=4&t！=7，返回（0））；n=10）；！！t吨

%o select（直，素数（1000））\\_Charles R Greathouse IV_，2012年9月25日

%Y参考A028373。

%K基，nonn

%O 1,1号机组

%A _Shyam Sunder Gupta，2003年1月23日

%E由_Robert G.Wilson v_修订和扩展，2003年1月24日