OFFSET
0,9
COMMENTS
The indices of this sequence are sorted by evaluation at x=2, i.e. 0, 1, x, x+1, x^2, x^2+1, etc. Generalization of arithmetic derivative is as described by Victor Ufnarovski and Bo Ahlander. For a general UFD this requires a choice of canon primes for which p' = 1. However, GF(2)[x] has only one unit, so there is a unique arithmetic derivative over this UFD.
The arithmetic derivative of a square polynomial in GF(2)[x] is 0 so square polynomials act like units: given f, g in GF(2)[x], (f*f*g)' = f*f*(g'). The arithmetic derivative of a positive degree squarefree polynomial in GF(2)[x] is always nonzero.
LINKS
Victor Ufnarovski and Bo Ahlander, How to Differentiate a Number, J. Integer Seqs., Vol. 6, 2003.
EXAMPLE
To find a(9), first convert 9 into its corresponding GF(2) polynomial x^3 + 1. Then find its arithmetic derivative, x^2. Finally, convert to an integer via evaluation at x=2, giving a(9) = 4.
PROG
(SageMath)
P.<x> = GF(2)[]
def AD(p):
F = list(p.factor())
f = [f[0] for f in F for _ in range(f[1])]
return SymmetricFunctions(P).e()([len(f)-1, 0]).expand(len(f))(f)
def more(l):
return [x*p for p in l]+[x*p+1 for p in l]
L = [x, x+1]
L = L + more(L) + more(more(L)) + more(more(more(L))) + more(more(more(more(L))))
L.sort()
', '.join(map(str, ([0, 0]+[AD(p).change_ring(ZZ)(2) for p in L])))
CROSSREFS
KEYWORD
nonn
AUTHOR
Keith J. Bauer, Feb 27 2024
STATUS
approved