A364288 a(n) = n - A243071(n).
1, 1, 0, 2, -2, 0, -8, 4, 4, -4, -20, 0, -50, -16, 2, 8, -110, 8, -236, -8, -8, -40, -488, 0, 14, -100, 18, -32, -994, 4, -2016, 16, -28, -220, 8, 16, -4058, -472, -86, -16, -8150, -16, -16340, -80, 20, -976, -32720, 0, 26, 28, -202, -200, -65482, 36, -4, -64, -452, -1988, -131012, 8, -262082, -4032, 6, 32, -58, -56
Offset: 1
Keywords
Links
- Antti Karttunen, Table of n, a(n) for n = 1..12288
Crossrefs
Programs
-
Mathematica
nn = 60; f[x_] := Times @@ Power[Which[# == 1, 1, # == 2, 1, True, NextPrime[#, -1]] & /@ First[#], Last[#] ] &@ Transpose@ FactorInteger@ x; Do[a[n] = Which[n <= 2, n - 1, OddQ[n], 1 + 2 a[f[n]], True, 2 a[n/2] ], {n, nn}]; Array[# - a[#] &, nn] (* Michael De Vlieger, Jul 25 2023 *)
-
PARI
A064989(n) = { my(f=factor(n>>valuation(n,2))); for(i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f); }; A243071(n) = if(n<=2, n-1, if(!(n%2), 2*A243071(n/2), 1+(2*A243071(A064989(n))))); A364288(n) = (n-A243071(n));