|
|
|
|
1, 2, 6, 17, 67, 166, 676, 1373, 4475, 10446, 30036, 51032, 196386, 315302, 737515, 1654229, 4227565, 6301902, 17975187, 26010425, 70085244, 133337963
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
|
|
MATHEMATICA
|
lps = Cases[Import["https://oeis.org/A025487/b025487.txt", "Table"], {_, _}][[;; , 2]]; s = {}; Do[p = Position[lps, Product[Prime[k]^Floor[n/k], {k, 1, n}]]; If[p == {}, Break[]]; AppendTo[s, p[[1, 1]]], {n, 0, 20}]; s
|
|
PROG
|
(PARI) f(m) = my(c=1, p, q=2, v=vector(logint(m, 2), i, 2^i), w); while(#v, c+=#v; p=q; q=nextprime(q+1); w=List([]); for(i=1, #v, for(j=1, min(valuation(v[i], p), logint(m\v[i], q)), listput(w, v[i]*q^j))); v=w); c;
a(n) = f(prod(k=1, n, prime(k)^(n\k))); \\ Jinyuan Wang, Jul 08 2021
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|