%I #32 May 25 2024 10:44:36
%S 1,2,3,2,5,3,7,4,3,5,11,6,13,7,5,8,17,9,19,10,7,11,23,12,5,13,9,14,29,
%T 15,31,16,11,17,7,18,37,19,13,20,41,21,43,22,15,23,47,24,7,25,17,26,
%U 53,27,11,28,19,29,59,30,61,31,21,32,13,33,67,34,23,35,71,36,73,37,25,38
%N a(n) = n if n is 1 or a prime, otherwise a(n) = n divided by the least prime factor of n (A032742(n)).
%C A026741 generalized to give either a prime or the largest proper divisor of a nonprime.
%C Sometimes called "Conway's subprime function", although it surely predates John Conway. - _N. J. A. Sloane_, Sep 29 2017
%H Reinhard Zumkeller, <a href="/A117818/b117818.txt">Table of n, a(n) for n = 1..10000</a>
%p A117818 := proc(n)
%p local a,d;
%p if isprime(n) or n =1 then
%p return n;
%p end if;
%p a := -1 ;
%p for d in numtheory[divisors](n) do
%p if d < n and d> a then
%p a := d ;
%p end if;
%p end do:
%p a ;
%p end proc:
%p seq(A117818(n),n=1..100) ; # _R. J. Mathar_, Apr 30 2024
%t Table[If[PrimeQ[n], n, If[n == 1, 1, n/FactorInteger[n][[1, 1]]]], {n, 1, 76}]
%t Table[Which[n==1,1,PrimeQ[n],1,True,Divisors[n][[-2]]],{n,80}] (* _Harvey P. Dale_, Feb 02 2022 *)
%o (Haskell)
%o a117818 n = if a010051 n == 1 then n else a032742 n
%o -- _Reinhard Zumkeller_, Jun 24 2013
%o (Python)
%o import sympy
%o def A117818(n):
%o if n == 1:
%o return 1
%o else:
%o _=sympy.ntheory.factor_.primefactors(n)
%o return _[-1]
%o print([A117818(n) for n in range(1,100)])
%o # _R. J. Mathar_, May 24 2024
%Y Cf. A026741, A032742, A292772.
%K nonn,easy
%O 1,2
%A _Roger L. Bagula_, Apr 30 2006
%E Edited by _Stefan Steinerberger_, Jul 22 2007
%E Extended by _Charles R Greathouse IV_, Jul 28 2010