%I#27 2023年3月6日02:59:40

%S 2,3,5,7199919991299992911177117171173133333171171，

%电话：717111616161119949911113331131131131331111311313311，

%电话：311313311113113111313113113113111311311311313133377733773773773377773773737733377737337773773333333333

%N个素数，其小数位数之和与素数具有相同的小数位数集。

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

%ea（5）=199是一个项，因为199是质数，1+9+9=19与199具有相同的十进制数字集{1,9}。

%p dmax:=7:#用于dmax数字以内的术语

%p dsets:=proc（s，s）选项记忆；

%p#带有总和s和元素集s的非递减列表[x_1，…，x_n]

%p局部i，x1；

%p如果S={}，那么如果S=0，则返回{[]}，否则返回{}fi-fi；

%px1:=最小值（S）；

%p`union`（seq（map（t->[x1$i，op（t）]，procname（s-i*x1，s减去{x1}）），i=1..`if`（x1=0，dmax，floor（s/x1）））

%p端程序：

%p R：={2,3,5,7}：计数：=4：

%s从2到9的p*dmax-1 do

%p如果s mod 3＝0，则下一个fi；

%p ds：=转换（转换（s，基数，10），集合）；

%p DS：=选择（t->nops（t）>1和nops（t）<=dmax，数据集（s，DS））；

%DS do中r的p

%p代表移除（t->成员（t[1]，[0,2,4,5,6,8]）或t[-1]=0，组合：-置换（r））do中的v

%p p：=加（v[i]*10^（i-1），i=1..nops（v））；

%p如果是素数（p），则R:=R并{p}；计数：=计数+1；

%功率因数

%日期日期：

%p排序（转换（R，列表））；

%o（PARI）isok（p）=如果（i素数（p），my（d=数字（p））；集合（d）==集合（数字（vecsum（d）））；\\_米歇尔·马库斯，2023年2月28日

%Y参考A000040，A158473。

%A249515中的Y引物。

%K nonn，基础

%O 1，1

%以色列罗贝尔特，2023年2月27日