%I#23 2017年11月16日15:52:25

%序号1,6,26102510321022890153690150729015618090

%N a（N）是使k的素因式分解中的不同位数为N的最小数k（将形式p^1的项计数为p）。

%将k写成素数的乘积；则a（n）是最小的数字k，使得k的乘积表示中的非重复数字总数（所有素数中的非循环数字数和所有指数中大于1的非循环数）等于n。按照惯例，第一项a（1）=1。序列已完成。

%C属性：所有指数都等于1（参见下面的示例）。

%e a（1）=1；

%e a（2）=6=2*3和A254315（6）=2；

%e a（3）=26=2*13和A254315（26）=3；

%e a（4）=102=2*3*17和A254315（102）=4；

%e a（5）=510=2*3*5*17和A254315（510）=5；

%e a（6）=3210=2*3*5*107和A254315（3210）=6；

%e a（7）=22890=2*3*5*7*109和A254315（22890）=7；

%e a（8）=153690=2*3*5*47*109和A254315（153690）=8；

%e a（9）=1507290=2*3*5*47*1069和A254315（1507290）=9；

%e a（10）=15618090=2*3*5*487*1069和A254315（15618090）=10。

%p使用（ListTools）：

%p代表n从2到10 do：

%p ii：=0：

%p代表k从2到10^9，而（ii=0）do：

%p n0:=长度（k）：lst:={}:x0:=系数（k）

%p y：=压扁（x0[2]）：z：=转换（y，集）：

%pz1:=z减去{1}:nn0:=nops（z1）：

%p代表从1到nn0的m do：

%p t1:=转换（z1[m]，基数，10）：z2:=转换（t1，集合）：

%p lst:=lst联合z2：

%日期：

%p nn1：=无（lst）：

%p如果nn1=n，则ii:=1:printf（“%d%d\n”，n，k）：

%p其他

%p菲：

%日期：

%日期：

%t f[n_]：=块[{pf=FactorInteger@n，i}，长度@DeleteDuplicates@Flatten@IntegerDigits@Rest@Flattan@Reap@Do[If[Last[pf[i]]]==1，母猪@First@pf[i]，母猪@FromDigits@Flatten[IntegerDigits/@pf[[i]]，{i，长度@pf}]]；b=-1；压扁@Last@Reap@Do[a=f[n]；如果[a>b，母猪[n]；b=a]，{n，10^6}]（*迈克尔·德弗里格，2015年1月29日*）

%t使用[{s=Array[CountDistinct@Flatten@IntegerDigits[FactorInteger[#]/.{p_，e_}/；e==1:>{p}]&，10^6]}，地图[FirstPosition[s，#][[1]]&，Range@Max@s]]（*_Michael De Vlieger_，2017年11月3日*）

%o（PARI）a（n）=对于（k=1,10^5，s=[]；F=系数（k）；对于（i=1，#F[，1]，s=concat（s，digits（F[i，1]））；如果（F[i，2]>1，s=concat（s，digits（F[i，2]））；如果（#vecsort（s，，8）==n，返回（k））

%o打印1（1，“，”）；（n=2,7，打印1（a（n），“，”）\\_Derek Orr_2015年1月30日

%Y参考A043537，A254315。

%K nonn、base、fini、full

%O 1,2号机组

%2015年1月28日，拉格瑙市

%E a（10）由Giovanni Resta更正，2017年11月3日