login
A365465
a(n) = A356867(n) / gcd(n, A356867(n)).
5
1, 1, 1, 5, 4, 1, 10, 1, 1, 7, 14, 5, 25, 10, 4, 25, 16, 1, 35, 7, 10, 125, 40, 1, 4, 16, 1, 11, 22, 7, 55, 11, 14, 35, 8, 5, 49, 49, 25, 35, 140, 10, 250, 20, 4, 245, 196, 25, 625, 4, 16, 125, 64, 1, 7, 55, 35, 275, 88, 7, 350, 56, 10, 343, 98, 125, 875, 70, 40, 125, 160, 1, 1225, 196, 4, 3125, 400, 16, 1000, 8, 1, 13, 26
OFFSET
1,4
COMMENTS
Denominator of n / A356867(n).
LINKS
FORMULA
a(n) = A356867(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];
A365465(n) = (A356867(n)/gcd(n, A356867(n)));
CROSSREFS
Cf. A356867, A365462, A365463, A365464 (numerators).
Cf. also A364502.
Sequence in context: A213055 A005752 A364355 * A098494 A008955 A182824
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 15 2023
STATUS
approved