A365464 a(n) = n / gcd(n, A356867(n)).

1, 1, 1, 4, 5, 1, 7, 1, 1, 10, 11, 4, 13, 7, 5, 8, 17, 1, 19, 5, 7, 22, 23, 1, 1, 13, 1, 28, 29, 10, 31, 8, 11, 17, 5, 4, 37, 19, 13, 8, 41, 7, 43, 11, 5, 46, 47, 8, 49, 1, 17, 13, 53, 1, 5, 28, 19, 58, 59, 5, 61, 31, 7, 64, 13, 22, 67, 17, 23, 7, 71, 1, 73, 37, 1, 76, 77, 13, 79, 5, 1, 82, 83, 28, 17, 43, 29, 44, 89, 10

Offset

1,4

Comments

Numerator of n / A356867(n).

Links

Antti Karttunen, Table of n, a(n) for n = 1..19683

Formula

a(n) = n / A365463(n).

Prog

(PARI)
up_to = 19683;
A356867list(up_to) = { my(v=vector(up_to), met=Map(), h=0, ak); for(i=1, #v, if(1==vecsum(digits(i, 3)), v[i] = i; h = i, ak = v[i-h]; forprime(p=2, , if(3!=p && !mapisdefined(met, p*ak), v[i] = p*ak; break))); mapput(met, v[i], i)); (v); };
v356867 = A356867list(up_to);
A356867(n) = v356867[n];
A365464(n) = n/gcd(n, A356867(n));

Crossrefs

Cf. A356867, A365462, A365463, A365465 (denominators).

Cf. also A364501.

Sequence in context: A094642 A069284 A272638 * A299630 A068447 A375822

Adjacent sequences: A365461 A365462 A365463 * A365465 A365466 A365467

Keyword

nonn,frac

Author

Antti Karttunen, Sep 15 2023

Status

approved