cp's OEIS Frontend

A284253 a(n) = n / A284252(n).

%I A284253 #20 Mar 28 2021 07:01:26
%S A284253 1,2,3,4,5,6,7,8,9,2,11,12,13,2,15,16,17,18,19,4,21,2,23,24,25,2,27,4,
%T A284253 29,6,31,32,3,2,35,36,37,2,3,8,41,6,43,4,45,2,47,48,49,10,3,4,53,54,
%U A284253 55,8,3,2,59,12,61,2,63,64,65,6,67,4,3,14,71,72,73,2,75,4,77,6,79,16,81,2,83,12,85,2,3,8,89,18,91,4,3,2,95,96,97,14,9,20,101,6
%N A284253 a(n) = n / A284252(n).
%H A284253 Antti Karttunen, <a href="/A284253/b284253.txt">Table of n, a(n) for n = 1..10001</a>
%F A284253 a(n) = n / A284252(n).
%t A284253 a[n_] := Block[{p = First /@ FactorInteger[n]}, n / SelectFirst[p, # > p[[1]]^2 &, 1]]; Array[a, 120] (* _Giovanni Resta_, Mar 24 2017 *)
%o A284253 (Scheme) (define (A284253 n) (/ n (A284252 n)))
%o A284253 (PARI) a(n) = if(n<2, return(1), my(f=factor(n)[, 1]); for(i=2, #f, if(f[i]>f[1]^2, return(f[i]))); return(1));
%o A284253 for(n=1, 151, print1(n/a(n),", ")) \\ _Indranil Ghosh_, after _David A. Corneth_, Mar 24 2017
%o A284253 (Python)
%o A284253 from sympy import primefactors
%o A284253 def a(n):
%o A284253     pf = primefactors(n)
%o A284253     if pf: min_pf2 = min(pf)**2
%o A284253     for i in pf:
%o A284253         if i > min_pf2: return i
%o A284253     return 1
%o A284253 print([n//a(n) for n in range(1, 151)]) # _Indranil Ghosh_, Mar 24 2017
%Y A284253 Cf. A020639, A284252, A284254, A284255, A284256, A284257, A284258, A284259.
%K A284253 nonn
%O A284253 1,2
%A A284253 _Antti Karttunen_, Mar 24 2017