|
|
|
|
1, 2, 3, 8, 5, 12, 7, 32, 27, 20, 11, 48, 13, 28, 45, 128, 17, 108, 19, 80, 63, 44, 23, 192, 125, 52, 243, 112, 29, 180, 31, 512, 99, 68, 175, 432, 37, 76, 117, 320, 41, 252, 43, 176, 405, 92, 47, 768, 343, 500, 153, 208, 53, 972, 275, 448, 171, 116, 59, 720
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
This sequence has similarities with A087019.
These are the positions of first appearances of each positive integer in A346701, and also in A346703. - Gus Wiseman, Aug 09 2021
|
|
LINKS
|
|
|
FORMULA
|
a(n) = n iff n = 1 or n is a prime number.
a(p^k) = p^(2*k-1) for any k > 0 and any prime number p.
If g = A006530(n) is the greatest prime factor of n, then a(n) = n^2/g.
(End)
|
|
EXAMPLE
|
For n = 42:
- 42 = 2 * 3 * 7, so:
2 3 7
x 2 3 7
-------
2 3 7
2 3 3
+ 2 2 2
-----------
2 2 3 3 7
- hence a(42) = 2 * 2 * 3 * 3 * 7 = 252.
|
|
MATHEMATICA
|
Table[n^2/FactorInteger[n][[-1, 1]], {n, 100}] (* Gus Wiseman, Aug 09 2021 *)
|
|
PROG
|
(PARI) See Links section.
|
|
CROSSREFS
|
The version for even indices is A129597(n) = 2*a(n) for n > 1.
These are the positions of first appearances in A346701 and in A346703.
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A209281 adds up the odd bisection of standard compositions (even: A346633).
A346697 adds up the odd bisection of prime indices (reverse: A346699).
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|