[go: up one dir, main page]

login
A342918
a(n) = (1+n) / A143771(n).
3
1, 1, 1, 5, 1, 7, 1, 3, 5, 11, 1, 13, 1, 5, 2, 17, 1, 19, 1, 7, 11, 23, 1, 25, 13, 9, 7, 29, 1, 31, 1, 11, 17, 35, 3, 37, 1, 13, 5, 41, 1, 43, 1, 15, 23, 47, 1, 49, 25, 17, 13, 53, 1, 55, 7, 19, 29, 59, 1, 61, 1, 21, 8, 65, 11, 67, 1, 23, 35, 71, 1, 73, 1, 25, 19, 77, 13, 79, 1, 27, 41, 83, 1, 85, 43, 29, 11, 89, 1, 91, 23
OFFSET
1,4
LINKS
FORMULA
a(n) = (1+n) / A143771(n).
PROG
(PARI)
A143771(n) = my(d = divisors(n)); gcd(vector(#d, k, d[k]+n/d[k])); \\ From A143771.
A342918(n) = ((1+n)/A143771(n));
CROSSREFS
After n=1 differs from A342916 for the first time at n=44, where a(44) = 15, while A342916(44) = 5.
Sequence in context: A364491 A190643 A068328 * A109375 A051712 A086892
KEYWORD
nonn
AUTHOR
Antti Karttunen, Mar 29 2021
STATUS
approved