%I#21 2020年10月15日19:26:45

%S 1,3,2,3,10,3,13,4,5,10,22,3,26,13,10,10,4,34,5,37,10,13,22,46,4,15,26,6，

%电话：13,58,10,61,5,22,34,13,5,73,37,26,10,82,13,86,22,10,46,94,4,14,15,34，

%U 26106,6,22,13,37,58118,10122,61,13,6,26,22134,34,46,13142,5146,73,15,37,22,26157号

%N a（N）是N除以psi（k！）（k>0）的最小k。

%C自2017年9月14日起，以色列：（开始）

%C如果m和n是互质，则a（m*n）=max（a（m），a（n））。

%C a（n）<=2n。

%假设p是一个素数>=5。然后

%如果p在A005382中，则C a（p）=2p-1，否则为2p。

%如果p在A005382中，则C a（p^2）=2p，否则为3p。

%如果p在A005382中，则C a（p^3）=3p；如果p在P062737中，则为4p-1，否则为4p。

%C（结束）

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

%e a（4）=3，因为4除以psi（3！）=12，3是具有此属性的最小值。

%p A:=proc（n）选项记忆；

%p局部F，p，e，t，k；

%p F:=系数（n）[2]；

%p如果nops（F）=1，则

%p p：=F[1]；

%p e：=F[1][2]；

%p如果p＝3，那么

%p t：=1；如果e=1，则返回2 fi

%其他t:=0:

%p fi；

%从2*p到p的k的p

%p如果是素数（k-1），则

%p t：=t+padic：-ordp（k，p）；

%p如果t>=e，则返回（k-1）fi；

%p fi；

%p t：=t+padic：-ordp（k，p）；

%p如果t>=e，则返回k fi；

%日期

%p其他

%p max（seq（进程名（t[1]^t[2]），t=F）

%功率因数

%p端程序：

%p A（1）：=1：

%p映射（A，[1..100]美元）；#_罗伯特·伊斯雷尔，2017年9月14日

%t psi[n_]：=模[{p，e}，乘积[{p、e}=pe；p^e+p^（e-1），{pe，因子整数[n]}]]；

%t a[n_]：=模[{k=1}，而[！可除[psi[k！]，n]，k++]；k] ；a[2]=3；

%t阵列[a，100]（*_Jean-François Alcover_，2020年10月15日，PARI之后*）

%o（PARI）a001615（n）=我的（f=系数（n））；prod（i=1，#f~，f[i，1]^f[i、2]+f[i（1]^（f[i）-1））；

%o a（n）={my（k=1）；while（a001615（k！）%n，k++）；k；}\\在a001615的Charles R Greathouse IV_之后

%Y参见A001615、A005382、A062737、A275985。

%K nonn，看

%O 1,2号机组

%A Altug Alkan，2017年9月7日