|
|
A346691
|
|
Replace 7^k with (-1)^k in base-7 expansion of n.
|
|
5
|
|
|
0, 1, 2, 3, 4, 5, 6, -1, 0, 1, 2, 3, 4, 5, -2, -1, 0, 1, 2, 3, 4, -3, -2, -1, 0, 1, 2, 3, -4, -3, -2, -1, 0, 1, 2, -5, -4, -3, -2, -1, 0, 1, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 0, 1, 2, 3, 4, 5, 6, -1, 0, 1, 2, 3, 4, 5, -2, -1, 0, 1, 2, 3, 4, -3, -2, -1, 0, 1, 2, 3, -4, -3, -2, -1, 0, 1, 2, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
If n has base-7 expansion abc..xyz with least significant digit z, a(n) = z - y + x - w + ...
|
|
LINKS
|
|
|
FORMULA
|
G.f. A(x) satisfies: A(x) = x * (1 + 2*x + 3*x^2 + 4*x^3 + 5*x^4 + 6*x^5) / (1 - x^7) - (1 + x + x^2 + x^3 + x^4 + x^5 + x^6) * A(x^7).
a(n) = n + 8 * Sum_{k>=1} (-1)^k * floor(n/7^k).
|
|
EXAMPLE
|
83 = 146_7, 6 - 4 + 1 = 3, so a(83) = 3.
|
|
MATHEMATICA
|
nmax = 104; A[_] = 0; Do[A[x_] = x (1 + 2 x + 3 x^2 + 4 x^3 + 5 x^4 + 6 x^5)/(1 - x^7) - (1 + x + x^2 + x^3 + x^4 + x^5 + x^6) A[x^7] + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]
Table[n + 8 Sum[(-1)^k Floor[n/7^k], {k, 1, Floor[Log[7, n]]}], {n, 0, 104}]
|
|
PROG
|
(Python)
from sympy.ntheory.digits import digits
def a(n):
return sum(bi*(-1)**k for k, bi in enumerate(digits(n, 7)[1:][::-1]))
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign,base
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|