%I#38 2019年7月27日00:40:23

%S 1,1,2,1,3,1,2,3,5,1,2,7,5,2,1,3,2,7,11,1,2,5,13,3,2,1,1,2,11，

%第17、7、2、1、19、13、2、1、3、1、2、3、23、1、2、7、5、17、2、1、3、11、19、29、1、2、1、3、13，

%U 2,13,3,1,2,23,5,1,2,17,5,2,1,1,1,1,2,1,1,2,3,41,12,17,43,29,2,1,3,13,13,2,31,47

%N的最大真除数的最小素因子。

%C对于n>1:a（n）=1，当n是素数时；a（A001358（n））=A084127（n）；a（A025475（n））=A020639（A02547（n）。【由_Peter Munn于2017年2月19日更正】

%当n是复合的时，这是第二个因子，当n被写成素数的乘积时，按非递减顺序。例如，12=2*2*3，那么a（12）=2.-_Peter Munn_，2017年2月19日

%C对于所有足够大的n，a（1）、a（2）、……的中值。。。a（n）为A281889（2）=7_Peter Munn，2019年7月12日

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

%F a（n）=A020639（A032742（n））。

%F A117357（n）=A020639（A054576（n））；A117358（n）=A032742（A054576（n））=A054575（n）/A117357（n）。-_Reinhard Zumkeller，2006年3月10日

%F如果A001222（n）>=2，a（n）=A027746（n，2），否则a（n_Peter Munn，2019年7月13日

%t PrimeFactors[n_]：=扁平[表[#[[1]]，{1}]&/@FactorInteger[n]]；f[n_]：=块[{gpd=Divisors[n][[-2]]}，如果[gpd==1，1，素数因子[gpd][[1]]]；表[如果[n==1，1，f[n]]，{n，1，95}]

%t（*第二个程序：*）

%t表[If[Or[PrimeQ@n，n==1]，1，FactorInteger[n/SelectFirst[Prime@Range@PrimePi[Sqrt@n]，Divisible[n，#]&]][[1，1]]，{n，94}]（*_Michael De Vlieger_，2017年8月14日*）

%o（PARI）lpf（n）=如果（n>1，系数（n）[1,1]，1）

%o a（n）=lpf（n/lpf（n））\\查尔斯·格里特豪斯IV，2013年5月9日

%o（PARI）a（n）=如果（n<4||isprime（n），return（1））；my（f=系数（n））；如果（f[1,2]>1，f[1,1]，f[2,1]）\\-Charles R Greathouse IV_，2013年5月9日

%o（方案）（定义（A014673n）（A020639（/n（A020639n）））；；该条目下给出的A020639代码-Antti Karttune_，2017年8月12日

%Y参见A001222、A001358、A020639、A025475、A027746、A032742、A054576、A084127、A085392、A08539、A115561、A117357、A117358、A281889。

%K nonn公司

%O 1,4型

%A _Reinhard Zumkeller，2003年6月24日