OFFSET
1,2
COMMENTS
Sequence is prime-free.
LINKS
K. Strassburger, Plot of shotgun numbers
EXAMPLE
B(1) = 1,2,3,4,5,6,7,8, 9,10,11,12,13,14,...
B(2) = 1,4,3,2,5,8,7,6, 9,12,11,10,13,16,...
B(3) = 1,4,8,2,5,3,7,6,10,12,11, 9,13,16,...
B(4) = 1,4,8,6,5,3,7,2,10,12,11,14,13,16,...
MATHEMATICA
max = 66; b[1, j_] := j; b[k_, j_] := b[k, j] = b[k-1, j]; Do[b[k, 2j*k-k] = b[k-1, 2j*k]; b[k, 2j*k] = b[k-1, 2j*k-k], {k, 2, max}, {j, 1, max}]; a[n_] := b[max, n]; Table[a[n], {n, 1, max}] (* Jean-François Alcover, Oct 11 2012 *)
PROG
(SageMath)
def divsign(s, k):
if not k.divides(s): return 0
return (-1)^(s//k)*k
def A064494(n):
s = n
for k in srange(n, 1, -1):
s -= divsign(s, k)
return s
print([A064494(n) for n in (1..66)]) # Peter Luschny, Sep 16 2019
CROSSREFS
KEYWORD
nice,nonn
AUTHOR
Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Oct 16 2001
STATUS
approved