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

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, 29, 10, 31, 16, 33, 17, 35, 9, 37, 19, 39, 20, 41, 14, 43, 22, 45, 23, 47, 12, 49, 25, 51, 26, 53, 18, 55, 28, 57, 29, 59, 10, 61, 31, 63, 32, 65, 22, 67, 34, 69, 35, 71, 18, 73, 37, 75

Offset

1,3

Links

Michel Marcus, Table of n, a(n) for n = 1..10000

Mathematica

a[n_] := Module[{k = 1}, While[Divisible[n, k], k++]; k--; n/k]; Array[a, 100] (* Amiram Eldar, Jan 02 2022 *)

Prog

(PARI) a(n) = my(m = 1); while ((n % m) == 0, m++); n/(m-1);
(Python)
def a(n):
    m = 2
    while n%m == 0: m += 1
    return n//(m-1)
print([a(n) for n in range(1, 80)]) # Michael S. Branicky, Jan 02 2022

Crossrefs

Cf. A055874.

Sequence in context: A328579 A046524 A086571 * A133945 A124219 A103833

Adjacent sequences: A350506 A350507 A350508 * A350510 A350511 A350512

Keyword

nonn

Author

Michel Marcus, Jan 02 2022, after a suggestion from Charles Kusniec

Status

approved