OFFSET
1,1
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..16384
Michael De Vlieger, 4096 X 4096 pixel raster with origin (0, 0) in the upper left corner and black pixels at (x, y), indicate a number 4096*(y-1) + (x-1) in this sequence. Thus this image contains 7852685 terms of this sequence.
EXAMPLE
MATHEMATICA
f[n_] := f[n] = Which[Abs@ n < 2, 0, PrimeQ[n], 1, True, n Total[#2/#1 & @@@ FactorInteger[Abs@ n]]]; g[n_] := And[n > 0, AnyTrue[FactorInteger[n], #2 >= #1 & @@ # &]]; w = {}; nn = 2^16; k = 1; While[Set[m, #^#] <= nn &[Prime[k]], AppendTo[w, m]; k++]; Reap[Do[If[! g[n], If[g@ NestWhile[f, n, And[! Divisible[#, 4], FreeQ[w, #]] &], Sow[n] ] ], {n, 2, nn}] ][[-1, -1]]
(* or, generate up to 7852685 terms of this sequence from the bitmap by setting y to a number not exceeding 4096: *)
With[{img = https://oeis.org/A359547/a359547.png, y = 2}, Map[4096 (#1 - 1) + #2 - 1 & @@ # &, Position[ImageData[img][[1 ;; y, All]], 0.]] ] (* Michael De Vlieger, Jan 23 2023 *)
PROG
(PARI) isA359547(n) = A359546(n);
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 05 2023
STATUS
approved