A066323 -id:A066323 - cp's OEIS frontend

A342729 Self numbers in base i-1: numbers not of the form k + A066323(k).

1, 3, 5, 7, 9, 22, 24, 26, 39, 41, 43, 56, 58, 60, 73, 75, 77, 90, 92, 94, 107, 109, 111, 136, 138, 140, 153, 155, 157, 170, 172, 174, 199, 201, 203, 216, 218, 220, 233, 235, 237, 262, 264, 266, 279, 281, 283, 296, 298, 300, 313, 315, 317, 330, 332, 334, 347, 349

Offset: 1

Amiram Eldar, Mar 19 2021

Equivalently, self numbers in base -4, since A066323(k) is also the sum of the digits of k in base -4.
Analogous to self numbers (A003052) using base i-1 representation (A271472) instead of decimal expansion.
The number of terms not exceeding 10^k, for k=1,2,..., is 5, 20, 155, 1507, 15008, 150007, 1500014, 15000011. Is the asymptotic density of this sequence exactly 3/20?

József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 4, p. 384-386.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Walter Penney, A "binary" system for complex numbers, Journal of the ACM, Vol. 12, No. 2 (1965), pp. 247-248.
Eric Weisstein's World of Mathematics, Self Number.
Wikipedia, Self number.

Cf. A007608, A066323, A066321, A271472, A342725, A342726, A342727, A342728.

Similar sequences: A003052 (decimal), A010061 (binary), A010064 (base 4), A010067 (base 6), A010070 (base 8), A339211 (Zeckendorf), A339212 (dual Zeckendorf), A339213 (base phi), A339214 (factorial base), A339215 (primorial base).

s[n_] := Module[{v = {{0, 0, 0, 0}, {0, 0, 0, 1}, {1, 1, 0, 0}, {1, 1, 0, 1}}}, Plus @@ Flatten @ v[[1 + Reverse @ Most[Mod[NestWhileList[(# - Mod[#, 4])/-4 &, n, # != 0 &], 4]]]]]; f[n_] := n + s[n]; m = 1000; Complement[Range[m], Select[Union@Array[f, m], # <= m &]]

A342728 a(n) is the least number k such that A066323(k) = n.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 23, 39, 55, 71, 87, 103, 359, 615, 871, 1127, 1383, 1639, 5735, 9831, 13927, 18023, 22119, 26215, 91751, 157287, 222823, 288359, 353895, 419431, 1468007, 2516583, 3565159, 4613735, 5662311, 6710887, 23488103, 40265319, 57042535, 73819751

Offset: 0

Amiram Eldar, Mar 19 2021

a(n) is the least number k whose sum of digits in base i-1 (or in base -4) is n.

Amiram Eldar, Table of n, a(n) for n = 0..1000
Walter Penney, A "binary" system for complex numbers, Journal of the ACM, Vol. 12, No. 2 (1965), pp. 247-248.
Index entries for linear recurrences with constant coefficients, signature (1,0,0,0,0,16,-16).

Cf. A007608, A066321, A066323, A271472, A342725, A342726, A342727, A342729.

Join[{0}, LinearRecurrence[{1,0,0,0,0,16,-16}, Range[7], 50]]

a(n) = n for n <= 7, and a(n) = a(n-1) + 16*a(n-6) - 16*a(n-7) for n > 7.
G.f.: x*(1 + x + x^2 + x^3 + x^4 + x^5 - 15*x^6)/(1 - x - 16*x^6 + 16*x^7). - Stefano Spezia, Mar 20 2021
From Greg Dresden, Jun 21 2021: (Start)
a(3*n+1) = (24 + (4^n)*(25 -  9*(-1)^n))/40.
a(3*n+2) = (24 + (4^n)*(50 +  6*(-1)^n))/40.
a(3*n+3) = (24 + (4^n)*(75 + 21*(-1)^n))/40. (End)

A007608 Nonnegative integers in base -4.

Original entry on oeis.org

0, 1, 2, 3, 130, 131, 132, 133, 120, 121, 122, 123, 110, 111, 112, 113, 100, 101, 102, 103, 230, 231, 232, 233, 220, 221, 222, 223, 210, 211, 212, 213, 200, 201, 202, 203, 330, 331, 332, 333, 320, 321, 322, 323, 310, 311, 312, 313, 300, 301, 302, 303, 13030

Offset: 0

N. J. A. Sloane, Robert G. Wilson v, Mira Bernstein

The base 2i representation (quater-imaginary representation) of nonnegative integers is obtained by interleaving with zeros, cf. A212494.

More precisely, a(n) is the number n written in base -4; numbers [which represent some nonnegative integer] in base -4 are 0, 1, 2, 3, 100, 101, 102, 103, 110, 111, 112, 113, 120, 121, 122, 123, 130, 131, 132, 133, ... (A212556) - M. F. Hasler, May 20 2012

D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, 1969, Vol. 2, p. 189.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Joerg Arndt, Table of n, a(n) for n = 0..1000
Matthew Szudzik, A Mathematica programming contest
Eric Weisstein's World of Mathematics, Negabinary
Wikipedia, Negative base

Cf. A212556 (sorted), A066323 (sum of digits), A212526 (negative integers in base -4).

Other bases: A007090, A039724, A073785, A073786, A073787, A073788, A073789, A073790 & A039723.

a007608 0 = 0
a007608 n = a007608 n' * 10 + m where
   (n', m) = if r < 0 then (q + 1, r + 4) else (q, r)
             where (q, r) = quotRem n (negate 4)
-- Reinhard Zumkeller, Jul 15 2012

ToNegaBases[i_Integer, b_Integer] := FromDigits[ Rest[ Reverse[ Mod[ NestWhileList[(#1 - Mod[ #1, b])/-b &, i, #1 != 0 &], b]]]]; Table[ ToNegaBases[n, 4], {n, 0, 55}]

A007608(n,s="")={until(!n\=-4,s=Str(n%-4,s));eval(s)}  \\ M. F. Hasler, May 20 2012

def A007608(n):
    s, q = '', n
    while q >= 4 or q < 0:
        q, r = divmod(q, -4)
        if r < 0:
            q += 1
            r += 4
        s += str(r)
    return int(str(q)+s[::-1]) # Chai Wah Wu, Apr 09 2016

A066321 Binary representation of base-(i-1) expansion of n: replace i-1 with 2 in base-(i-1) expansion of n.

Original entry on oeis.org

0, 1, 12, 13, 464, 465, 476, 477, 448, 449, 460, 461, 272, 273, 284, 285, 256, 257, 268, 269, 3280, 3281, 3292, 3293, 3264, 3265, 3276, 3277, 3088, 3089, 3100, 3101, 3072, 3073, 3084, 3085, 3536, 3537, 3548, 3549, 3520, 3521, 3532, 3533, 3344, 3345, 3356

Offset: 0

Marc LeBrun, Dec 14 2001

Here i = sqrt(-1).
First differences follow a strange period-16 pattern: 1 11 1 XXX 1 11 1 -29 1 11 1 -189 1 11 1 -29 where XXX is given by A066322. Number of one-bits is A066323.
From Andrey Zabolotskiy, Feb 06 2017: (Start)
(Observations.)
Actually, the sequence of the first differences can be split into blocks of size of any power of 2, and there will be only one position in the block that does not repeat. In this sense, one may say that the first differences follow (almost-)period-2^s pattern for any s > 0.
Specifically, the first differences are given by the formula: a(n+1)-a(n) = A282137(A007814((n xor ...110011001100) + 1)). Here binary representation of n is bitwise-xored with the period-4 bit sequence (A021913 written right-to-left) which is infinite or simply long enough; A007814(m) does not depend on the bits of m other than the least significant 1.
A282137 gives all first differences in the order of decreasing occurrence frequency.
(End)
Penney shows that since (i-1)^4 = -4, the representation a(n) of a real integer n is found by writing n in base -4 using digits 0 to 3 (A007608), changing those digits to bit strings 0000, 0001, 1100, 1101 respectively, and interpreting as binary. - Kevin Ryde, Sep 07 2019

			a(4) = 464 = 2^8 + 2^7 + 2^6 + 2^4 since (i-1)^8 + (i-1)^7 + (i-1)^6 + (i-1)^4 = 4.

D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, 1969, Vol. 2, p. 172. (See also exercise 16, p. 177; answer, p. 494.)

Paul Tek, Table of n, a(n) for n = 0..10000
Joerg Arndt, The fxt demos: bit wizardry, radix(-1+i), C++ radix-m1pi.h function bin_real_to_radm1pi().
Solomon I. Khmelnik, Specialized Digital Computer for Operations with Complex Numbers, Questions of Radio Electronics, 12 (1964), 60-82 [in Russian].
W. J. Penney, A "binary" system for complex numbers, JACM 12 (1965), 247-248.
N. J. A. Sloane, Table of n, (I-1)^n for n = 0..100
Paul Tek, Perl program for this sequence
Andrey Zabolotskiy, negation & addition of Gaussian integers written in base i-1 [Python script].

Cf. A066322, A066323, A282137.
See A271472 for the conversion of these decimal numbers to binary.
See A009116 and A009545 for real and imaginary parts of (i-1)^n (except for signs).
See A256441 for expansions of -n.

f:= proc(n) option remember; local t,m;
   t:= n mod 4;
   procname(t) + 16*procname((t-n)/4)
end proc:
f(0):= 0: f(1):= 1: f(2):= 12: f(3):= 13:
seq(f(i),i=0..100); # Robert Israel, Oct 21 2016

a(n) = my(ret=0,p=0); while(n, ret+=[0,1,12,13][n%4+1]<Kevin Ryde, Sep 07 2019

Perl
```
See Links section.
```

from gmpy2 import c_divmod
u = ('0000','1000','0011','1011')
def A066321(n):
    if n == 0:
        return 0
    else:
        s, q = '', n
        while q:
            q, r = c_divmod(q, -4)
            s += u[r]
        return int(s[::-1],2) # Chai Wah Wu, Apr 09 2016

In "rebase notation" a(n) = (i-1)[n]2.

G.f. g(z) satisfies g(z) = z*(1+12*z+13*z^2)/(1-z^4) + 16*z^4*(13+12*z^4+z^8)/((1-z)*(1+z^4)*(1+z^8)) + 256*(1-z^16)*g(z^16)/(z^12-z^13). - Robert Israel, Oct 21 2016

A342726 Niven numbers in base i-1: numbers that are divisible by the sum of their digits in base i-1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 10, 12, 15, 16, 18, 20, 24, 25, 30, 32, 33, 35, 36, 40, 42, 44, 45, 48, 50, 54, 60, 64, 65, 66, 70, 77, 80, 88, 90, 96, 99, 100, 110, 112, 120, 124, 125, 126, 130, 140, 144, 145, 147, 150, 156, 160, 168, 170, 180, 182, 184, 185, 186, 190, 192

Offset: 1

Amiram Eldar, Mar 19 2021

Numbers k that are divisible by A066323(k).

Equivalently, Niven numbers in base -4, since A066323(k) is also the sum of the digits of k in base -4.

			2 is a term since its representation in base i-1 is 1100 and 1+1+0+0 = 2 is a divisor of 2.
10 is a term since its representation in base i-1 is 111001100 and 1+1+1+0+0+1+1+0+0 = 5 is a divisor of 10.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Walter Penney, A "binary" system for complex numbers, Journal of the ACM, Vol. 12, No. 2 (1965), pp. 247-248.

Cf. A007608, A066321, A066323, A271472, A342725, A342727, A342728, A342729.

Similar sequences: A005349 (decimal), A049445 (binary), A064150 (ternary), A064438 (quaternary), A064481 (base 5), A118363 (factorial), A328208 (Zeckendorf), A328212 (lazy Fibonacci), A331085 (negaFibonacci), A333426 (primorial), A334308 (base phi), A331728 (negabinary), A342426 (base 3/2).

v = {{0, 0, 0, 0}, {0, 0, 0, 1}, {1, 1, 0, 0}, {1, 1, 0, 1}}; q[n_] := Divisible[n, Total[Flatten @ v[[1 + Reverse @ Most[Mod[NestWhileList[(# - Mod[#, 4])/-4 &, n, # != 0 &], 4]]]]]]; Select[Range[200], q]

A342727 Digitally balanced numbers in base i-1: numbers that in base i-1 have the same number of 0's as 1's.

Original entry on oeis.org

2, 21, 26, 31, 36, 41, 46, 51, 310, 315, 325, 330, 335, 340, 345, 350, 355, 360, 365, 370, 375, 390, 395, 405, 410, 415, 420, 425, 430, 435, 455, 470, 475, 485, 490, 495, 535, 550, 555, 565, 570, 575, 580, 585, 590, 595, 600, 605, 610, 620, 625, 630, 635, 645

Offset: 1

Amiram Eldar, Mar 19 2021

			2 is a term since its representation in base i-1, 1100, has 2 0's and 2 1's.
21 is a term since its representation in base i-1, 110011010001, has 6 0's and 6 1's.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Walter Penney, A "binary" system for complex numbers, Journal of the ACM, Vol. 12, No. 2 (1965), pp. 247-248.

Cf. A066321, A066323, A271472, A342725, A342726, A342728, A342729.

Similar sequences: A031443 (binary), A210619 (Zeckendorf).

v = {{0, 0, 0, 0}, {0, 0, 0, 1}, {1, 1, 0, 0}, {1, 1, 0, 1}}; balQ[n_] := Plus @@ (d = IntegerDigits[n]) == Length[d]/2; q[n_] := balQ @ FromDigits[Flatten@v[[1 + Reverse @ Most[Mod[NestWhileList[(# - Mod[#, 4])/-4 &, n, # != 0 &], 4]]]]]; Select[Range[1000], q]

A342725 Numbers that are palindromic in base i-1.

Original entry on oeis.org

0, 1, 13, 17, 189, 205, 257, 273, 3005, 3069, 3277, 3341, 4033, 4097, 4305, 4369, 48061, 48317, 49149, 49405, 52173, 52429, 53261, 53517, 64449, 64705, 65537, 65793, 68561, 68817, 69649, 69905, 768957, 769981, 773309, 774333, 785405, 786429, 789757, 790781, 834509

Offset: 1

Amiram Eldar, Mar 19 2021

Amiram Eldar, Table of n, a(n) for n = 1..512
Walter Penney, A "binary" system for complex numbers, Journal of the ACM, Vol. 12, No. 2 (1965), pp. 247-248.

Cf. A066321, A066323, A271472, A342726, A342727, A342728, A342729.

Similar sequences: A002113 (decimal), A006995 (binary), A014190 (base 3), A014192 (base 4), A029952 (base 5), A029953 (base 6), A029954 (base 7), A029803 (base 8), A029955 (base 9), A046807 (factorial base), A094202 (Zeckendorf), A331191 (dual Zeckendorf), A331891 (negabinary), A333423 (primorial base).

v = {{0, 0, 0, 0}, {0, 0, 0, 1}, {1, 1, 0, 0}, {1, 1, 0, 1}}; q[n_] := PalindromeQ @ FromDigits[Flatten @ v[[1 + Reverse @ Most[Mod[NestWhileList[(# - Mod[#, 4])/-4 &, n, # != 0 &], 4]]]]]; Select[Range[0, 10^4], q]

13 is a term since its base-(i-1) presentation is 100010001 which is palindromic.

cp's OEIS Frontend

A342729 Self numbers in base i-1: numbers not of the form k + A066323(k).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

A342728 a(n) is the least number k such that A066323(k) = n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A007608 Nonnegative integers in base -4.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

A066321 Binary representation of base-(i-1) expansion of n: replace i-1 with 2 in base-(i-1) expansion of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

PARI

Perl

Python

Formula

A342726 Niven numbers in base i-1: numbers that are divisible by the sum of their digits in base i-1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A342727 Digitally balanced numbers in base i-1: numbers that in base i-1 have the same number of 0's as 1's.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

A342725 Numbers that are palindromic in base i-1.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula