|
|
A015053
|
|
Smallest positive integer for which n divides a(n)^6.
|
|
9
|
|
|
1, 2, 3, 2, 5, 6, 7, 2, 3, 10, 11, 6, 13, 14, 15, 2, 17, 6, 19, 10, 21, 22, 23, 6, 5, 26, 3, 14, 29, 30, 31, 2, 33, 34, 35, 6, 37, 38, 39, 10, 41, 42, 43, 22, 15, 46, 47, 6, 7, 10, 51, 26, 53, 6, 55, 14, 57, 58, 59, 30, 61, 62, 21, 2, 65, 66, 67, 34, 69, 70, 71, 6, 73, 74, 15, 38, 77, 78
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
Sum_{k=1..n} a(k) ~ c * n^2, where c = (zeta(11)/2) * Product_{p prime} (1 - 1/p^2 + 1/p^3 - 1/p^4 + 1/p^5 - 1/p^6 + 1/p^7 - 1/p^8 + 1/p^9 - 1/p^10) = 0.3522558764... . - Amiram Eldar, Oct 27 2022
|
|
MATHEMATICA
|
spi[n_]:=Module[{k=1}, While[PowerMod[k, 6, n]!=0, k++]; k]; Array[spi, 80] (* Harvey P. Dale, Feb 29 2020 *)
f[p_, e_] := p^Ceiling[e/6]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 18 2020 *)
|
|
PROG
|
(PARI) a(n) = my(f=factor(n)); for (i=1, #f~, f[i, 2] = ceil(f[i, 2]/6)); factorback(f); \\ Michel Marcus, Feb 15 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,mult,easy
|
|
AUTHOR
|
R. Muller (Research37(AT)aol.com)
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|