login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A351436
a(n) = n - A351168(n).
0
0, 1, 1, 3, 1, 2, 1, 7, 5, 2, 1, 4, 1, 2, 3, 15, 1, 10, 1, 4, 3, 2, 1, 8, 9, 2, 19, 4, 1, 6, 1, 31, 3, 2, 5, 20, 1, 2, 3, 8, 1, 6, 1, 4, 9, 2, 1, 16, 13, 18, 3, 4, 1, 38, 5, 8, 3, 2, 1, 12, 1, 2, 9, 63, 5, 6, 1, 4, 3, 10, 1, 40, 1, 2, 27, 4, 7, 6, 1, 16, 65, 2
OFFSET
1,4
FORMULA
a(n) = n * (1 - ((A006530(n) - 1)/A006530(n))^A071178).
MATHEMATICA
a[n_] := n - Module[{f = FactorInteger[n]}, n*(1 - 1/f[[-1, 1]])^f[[-1, 2]]]; a[1] = 0; Table[a[n], {n, 2, 83}] (* Robert P. P. McKone, Feb 11 2022, from Amiram Eldar in A351168 *)
CROSSREFS
Cf. A351168, A006530 (largest prime factor), A071178 (exponent).
Sequence in context: A106790 A078897 A322034 * A226629 A349620 A349380
KEYWORD
nonn
AUTHOR
Ben Polson, Feb 11 2022
STATUS
approved