cp's OEIS Frontend

A350509 a(n) = n/A055874(n).

%I A350509 #14 Jan 03 2022 11:08:22
%S A350509 1,1,3,2,5,2,7,4,9,5,11,3,13,7,15,8,17,6,19,10,21,11,23,6,25,13,27,14,
%T A350509 29,10,31,16,33,17,35,9,37,19,39,20,41,14,43,22,45,23,47,12,49,25,51,
%U A350509 26,53,18,55,28,57,29,59,10,61,31,63,32,65,22,67,34,69,35,71,18,73,37,75
%N A350509 a(n) = n/A055874(n).
%H A350509 Michel Marcus, <a href="/A350509/b350509.txt">Table of n, a(n) for n = 1..10000</a>
%t A350509 a[n_] := Module[{k = 1}, While[Divisible[n, k], k++]; k--; n/k]; Array[a, 100] (* _Amiram Eldar_, Jan 02 2022 *)
%o A350509 (PARI) a(n) = my(m = 1); while ((n % m) == 0, m++); n/(m-1);
%o A350509 (Python)
%o A350509 def a(n):
%o A350509     m = 2
%o A350509     while n%m == 0: m += 1
%o A350509     return n//(m-1)
%o A350509 print([a(n) for n in range(1, 80)]) # _Michael S. Branicky_, Jan 02 2022
%Y A350509 Cf. A055874.
%K A350509 nonn
%O A350509 1,3
%A A350509 _Michel Marcus_, Jan 02 2022, after a suggestion from _Charles Kusniec_