OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = n + A085072(n). - Amiram Eldar, Dec 26 2020
MATHEMATICA
sig[n_] := Sort @ FactorInteger[n][[;; , 2]]; a[1] = 1; a[n_] := Module[{sign = sig[n], k = n + 1}, While[sig[k] != sign, k++]; k]; Array[a, 70] (* Amiram Eldar, Dec 26 2020 *)
PROG
(PARI) a(n) = {if (n==1, 1, my(k=1, s = vecsort(factor(n)[, 2]~)); while (vecsort(factor(n+k)[, 2]~) != s, k++); n+k; )} \\ Michel Marcus, Nov 02 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Apr 09 2003
STATUS
approved