A007089 -id:A007089 - cp's OEIS frontend

A074940 Numbers having at least one 2 in their ternary representation.

2, 5, 6, 7, 8, 11, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 29, 32, 33, 34, 35, 38, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 83, 86, 87, 88, 89, 92

Offset: 1

Benoit Cloitre and Reinhard Zumkeller, Oct 04 2002; revised Dec 03 2003

Also, numbers m such that 3 divides C(2m,m).
Also, numbers m such that the central trinomial coefficient A002426(m) == 0 (mod 3). - Emeric Deutsch and Bruce E. Sagan, Dec 04 2003
Also, numbers m such that A092255(m) == 0 (mod 3). - Benoit Cloitre, Mar 22 2004
Also, numbers m such that the coefficient of x^m equals 0 in Product_{k>=0} (1-x^(3^k)). - N. J. A. Sloane, Jun 01 2010

			12 is not in the sequence since it is 110_3, but 11 is in the sequence since it is 102_3. - _Michael B. Porter_, Jun 30 2016

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Emeric Deutsch and Bruce E. Sagan, Congruences for Catalan and Motzkin numbers and related sequences, arXiv:math/0407326 [math.CO], 2004.
Emeric Deutsch and Bruce E. Sagan, Congruences for Catalan and Motzkin numbers and related sequences, J. Num. Theory 117 (2006), 191-215.

Complement of A005836.
Cf. A006996, A007089, A081603, A081610, A081605, A081606.
A039966(a(n)) = 0.

a074940 n = a074940_list !! (n-1)
a074940_list = filter ((== 0) . a039966) [0..]
-- Reinhard Zumkeller, Jun 06 2012, Sep 29 2011

Select[Range@ 120, MemberQ[IntegerDigits[#, 3], 2] &] (* or *)
Select[Range@ 120, Divisible[Binomial[2 #, #], 3] &] (* Michael De Vlieger, Jun 29 2016 *)
Select[Range[100],DigitCount[#,3,2]>0&] (* Harvey P. Dale, Aug 25 2019 *)

is(n)=while(n,if(n%3==2,return(1));n\=3);0 \\ Charles R Greathouse IV, Aug 21 2011

from gmpy2 import digits
def A074940(n):
    def f(x):
        s = digits(x,3)
        for i in range(l:=len(s)):
            if s[i]>'1':
                break
        else:
            return n+int(s,2)
        return n+int(s[:i]+'1'*(l-i),2)
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Oct 29 2024

a(n) = n + O(n^0.631). - Charles R Greathouse IV, Aug 21 2011

More terms from Emeric Deutsch and Bruce E. Sagan, Dec 04 2003

A081604 Number of digits in ternary representation of n.

Original entry on oeis.org

1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5

Offset: 0

Reinhard Zumkeller, Mar 23 2003

a(n) is the length of row n in table A054635. - Reinhard Zumkeller, Sep 05 2014

			a(8) = 2 because 8 = 22_3, having 2 digits.
a(9) = 3 because 9 = 100_3, having 3 digits.

Reinhard Zumkeller, Table of n, a(n) for n = 0..10000
Eric Weisstein's World of Mathematics, Ternary.

Cf. A003137, A007089, A030341, A049803, A054635, A062153, A070939, A080342.
Cf. A062756, A077267, A081603.
Cf. A134021, A134025, A134026, A246435.

a081604 n = if n < 3 then 1 else a081604 (div n 3) + 1
-- Reinhard Zumkeller, Sep 05 2014, Feb 21 2013

A081604 := proc(n)
    max(1,1+ilog[3](n)) ;
end proc: # R. J. Mathar, Jul 12 2016

Table[Length[IntegerDigits[n, 3]], {n, 0, 99}] (* Alonso del Arte, Dec 30 2012 *)
Join[{1},IntegerLength[Range[120],3]] (* Harvey P. Dale, Apr 07 2019 *)

a(n) = A062153(n) + 1 for n >= 1.
a(n) = A077267(n) + A062756(n) + A081603(n);
From Reinhard Zumkeller, Oct 19 2007: (Start)
0 <= A134021(n) - a(n) <= 1;
a(A134025(n)) = A134021(A134025(n));
a(A134026(n)) = A134021(A134026(n)) - 1. (End)
a(n+1) = -Sum_{k=1..n} mu(3*k)*floor(n/k). - Benoit Cloitre, Oct 21 2009
a(n) = floor(log_3(n)) + 1. - Can Atilgan and Murat Erşen Berberler, Dec 05 2012
a(n) = if n < 3 then 1 else a(floor(n/3)) + 1. - Reinhard Zumkeller, Sep 05 2014
G.f.: 1 + (1/(1 - x))*Sum_{k>=0} x^(3^k). - Ilya Gutkovskiy, Jan 08 2017

A003137 Write n in base 3 and juxtapose.

Original entry on oeis.org

1, 2, 1, 0, 1, 1, 1, 2, 2, 0, 2, 1, 2, 2, 1, 0, 0, 1, 0, 1, 1, 0, 2, 1, 1, 0, 1, 1, 1, 1, 1, 2, 1, 2, 0, 1, 2, 1, 1, 2, 2, 2, 0, 0, 2, 0, 1, 2, 0, 2, 2, 1, 0, 2, 1, 1, 2, 1, 2, 2, 2, 0, 2, 2, 1, 2, 2, 2, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 2, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 2, 1, 0, 2, 0, 1, 0, 2, 1, 1, 0, 2, 2, 1

Offset: 1

N. J. A. Sloane

An irregular table in which the n-th row lists the base-3 digits of n, see A007089. - Jason Kimberley, Dec 07 2012

The base-3 Champernowne constant (A077771): it is normal in base 3. - Jason Kimberley, Dec 07 2012

			1,
2,
1,0,
1,1,
1,2,
2,0,
2,1,
2,2,
1,0,0,
1,0,1,.... _R. J. Mathar_, Aug 16 2021

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Reinhard Zumkeller, Rows n = 1..1000 of triangle, flattened
Eric Weisstein's World of Mathematics, Ternary.
Wikipedia, Ternary numeral system

Tables in which the n-th row lists the base b digits of n: A030190 and A030302 (b=2), this sequence and A054635 (b=3), A030373 (b=4), A031219 (b=5), A030548 (b=6), A030998 (b=7), A031035 and A054634 (b=8), A031076 (b=9), A007376 and A033307 (b=10). - Jason Kimberley, Dec 06 2012

Cf. A081604 (row lengths), A053735 (row sums), A030341 (rows reversed), A077771, A007089.

a003137 n k = a003137_tabf !! (n-1) !! k
a003137_row n = a003137_tabf !! (n-1)
a003137_tabf = map reverse $ tail a030341_tabf
a003137_list = concat a003137_tabf
-- Reinhard Zumkeller, Feb 21 2013

&cat[Reverse(IntegerToSequence(n,3)):n in[1..31]]; // Jason Kimberley, Dec 07 2012

Flatten@ IntegerDigits[ Range@ 40, 3] (* or *)
almostNatural[n_, b_] := Block[{m = 0, d = n, i = 1, l, p}, While[m <= d, l = m; m = (b - 1) i*b^(i - 1) + l; i++]; i--; p = Mod[d - l, i]; q = Floor[(d - l)/i] + b^(i - 1); If[p != 0, IntegerDigits[q, b][[p]], Mod[q - 1, b]]]; Array[ a[#, 3] &, 105] (* Robert G. Wilson v, Jul 01 2014 *)

from itertools import count, islice
from sympy.ntheory.factor_ import digits
def A003137_gen(): return (d for m in count(1) for d in digits(m,3)[1:])
A003137_list = list(islice(A003137_gen(),30)) # Chai Wah Wu, Jan 07 2022

More terms from Larry Reeves (larryr(AT)acm.org), Sep 25 2000

A062153 a(n) = floor(log_3(n)).

Original entry on oeis.org

0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4

Offset: 1

Henry Bottomley, Jun 06 2001

Harry J. Smith, Table of n, a(n) for n = 1..1000

Cf. A000244, A000523, A007089, A102572.

[ Ilog(3,n) : n in [1..100] ]; // N. J. A. Sloane, Dec 23 2006

Floor[Log[3,Range[120]]]  (* Harvey P. Dale, Apr 30 2011 *)

a(n)=logint(n,3) \\ Charles R Greathouse IV, Nov 15 2022

a(n) = (number of digits of n when written in base 3) - 1.
a(n) = if n > 2 then a(floor(n / 3)) + 1 else 0. - Reinhard Zumkeller, Oct 29 2001
G.f.: (1/(1 - x))*Sum_{k>=1} x^(3^k). - Ilya Gutkovskiy, Jan 08 2017

A134451 Ternary digital root of n.

Original entry on oeis.org

0, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2

Offset: 0

Reinhard Zumkeller, Oct 27 2007

Continued fraction expansion of sqrt(3) - 1. - N. J. A. Sloane, Dec 17 2007. Cf. A040001, A048878/A002530.
Minimum number of terms required to express n as a sum of odd numbers.
Shadow transform of even numbers A005843. - Michel Marcus, Jun 06 2013
From Jianing Song, Nov 01 2022: (Start)
For n > 0, a(n) is the minimal gap of distinct numbers coprime to n. Proof: denote the minimal gap by b(n). For odd n we have A058026(n) > 0, hence b(n) = 1. For even n, since 1 and -1 are both coprime to n we have b(n) <= 2, and that b(n) >= 2 is obvious.
The maximal gap is given by A048669. (End)

			n=42: A007089(42) = '1120', A053735(42) = 1+1+2+0 = 4,
A007089(4)='11', A053735(4)=1+1=2: therefore a(42) = 2.
0.732050807568877293527446341... = 0 + 1/(1 + 1/(2 + 1/(1 + 1/(2 + ...)))). - _Harry J. Smith_, May 31 2009

Harry J. Smith, Table of n, a(n) for n = 0..20000
Lorenz Halbeisen and Norbert Hungerbuehler, Number theoretic aspects of a combinatorial function, Notes on Number Theory and Discrete Mathematics 5(4) (1999), 138-150; see Definition 7 for the shadow transform.
N. J. A. Sloane, Transforms.
Eric Weisstein's World of Mathematics, Digital Root.
Eric Weisstein's World of Mathematics, Ternary.

Cf. A000010, A055034, A134452, A160390 (decimal expansion).
Apart from a(0) the same as A040001.
Related base-3 sequences: A053735, A134451, A230641, A230642, A230643, A230853, A230854, A230855, A230856, A230639, A230640, A010063 (trajectory of 1).

a134451 = until (< 3) a053735
-- Reinhard Zumkeller, May 12 2011

A134451:=n->((n+1) mod 2)+2*signum(n)-1; seq(A134451(n), n=0..100); # Wesley Ivan Hurt, Dec 06 2013

Table[Mod[n + 1, 2] + 2 Sign[n] - 1, {n, 0, 100}] (* Wesley Ivan Hurt, Dec 06 2013 *)

a(n) = if(!n, 0, 2 - n%2) \\ Andrew Howroyd, Dec 08 2024

a(n) = n if n <= 2, otherwise a(A053735(n)).
a(A005408(n)) = 1; a(A005843(n)) = 2 for n>0;
a(n) = 0 if n=0, otherwise A000034(n-1).
a(n) = ((n+1) mod 2) + 2*sign(n) - 1. - Wesley Ivan Hurt, Dec 06 2013
Multiplicative with a(2^e) = 2, a(p^e) = 1 for odd prime p. - Andrew Howroyd, Aug 06 2018
a(0) = A055034(1) / A000010(1), a(n) = A000010(n+1) / A055034(n+1), n>1. - Torlach Rush, Oct 29 2019
Dirichlet g.f.: zeta(s)*(1+1/2^s). - Amiram Eldar, Jan 01 2023

A356867 For n >= 1, write n = 3^m + k, where m >= 0 is the greatest power of 3 <= n, and k is in the range 0 <= k < 3^(m+1) - 3^m, then for n such that k=0, a(n)=n, and for n such that k > 0, a(n) is the smallest prime multiple p*a(k), p != 3, that is not already a term.

Original entry on oeis.org

1, 2, 3, 5, 4, 6, 10, 8, 9, 7, 14, 15, 25, 20, 12, 50, 16, 18, 35, 28, 30, 125, 40, 24, 100, 32, 27, 11, 22, 21, 55, 44, 42, 70, 56, 45, 49, 98, 75, 175, 140, 60, 250, 80, 36, 245, 196, 150, 625, 200, 48, 500, 64, 54, 77, 110, 105, 275, 88, 84, 350, 112, 90, 343

Offset: 1

David James Sycamore, Sep 01 2022

Any prime p may be used to generate a sequence D(p) of this kind. The present sequence is D(3), and D(2) is the Doudna sequence, A005940.
Conjectured to be a permutation of the positive integers in which the primes appear in order.
From Antti Karttunen, Sep 16 2023: (Start)
The conjecture is true: Sequence is a permutation of natural numbers. By definition it is injective, and the surjectivity is guaranteed by the fact that there are infinitely many such n > k encountered by the greedy algorithm that a(n) will be a multiple of a(k), and "the smallest prime multiple" condition guarantees that all multiples of a(k) will eventually appear. That the primes and A100484 appear in order follows from the formulas a(3^m + 1) = prime(m+2), and a(3^m + 2) = 2*prime(m+2).
If the base-3 representation of n-1 has the base-3 representation of k-1 as its suffix, then a(n) is a multiple of a(k). For example, A007089(16-1) = 120, and A007089(43-1) = 1120, thus the former is the suffix of the latter, and a(16) = 50 indeed divides a(43) = 250.
(End)

			n=1=3^0+0 so a(1)=1. n=2=3^0+1 so k=1 and a(2)=2. Similarly a(3)=3 and a(9)=9.
n=10=3^2+1, therefore k=1 and a(1)=1 so a(10)=1*7=7 (since 2 and 5 have already occurred).

Michael De Vlieger, Table of n, a(n) for n = 1..19683 (19683 = 3^9)
Michael De Vlieger, Fan style ternary tree showing a(n) for n = 1..3^9, with a heat map color function for level m where 3^m is blue, smaller values are bluer, and larger are yellow-green. The smallest value in level m is shown in purple and largest is shown in red.
Index entries for sequences that are permutations of the natural numbers

Cf. A007089, A007949, A011655, A048473, A100484, A053735, A364958 (fixed points), A365390 (inverse permutation), A365424, A365459, A365462 [= a(n)-n], A365463 [= gcd(a(n),n)], A365464, A365465, A365717 [= A348717(a(1+n))], A365719 [= A046523(a(1+n))], A365721 [= omega(a(1+n))], A365722 [= bigomega(a(1+n))].

Cf. also A005940, A364611, A364628 for variants D(2), D(5) and D(7).

Block[{a, c, i, j, k, m, t, nn}, nn = 64; m = 1; i = 2; p = Prime[i]; c[] = False; Do[Set[{m, k}, {1, n - p^Floor[Log[p, n]]}]; If[k == 0, Set[{a[n], c[n]}, {n, True}], While[Set[t, Prime[m] a[k]]; Or[m == i, c[t]], m++]; Set[{a[n], c[t]}, {t, True}]], {n, nn}]; Array[a, nn] ] (* _Michael De Vlieger, Sep 01 2022 *)

up_to = 19683;
A356867list(up_to) = { my(v=vector(up_to),met=Map(),h=0,ak); for(i=1,#v,if(1==vecsum(digits(i,3)), v[i] = i; h = i, ak = v[i-h]; forprime(p=2,,if(3!=p && !mapisdefined(met,p*ak), v[i] = p*ak; break))); mapput(met,v[i],i)); (v); };
v356867 = A356867list(up_to);
A356867(n) = v356867[n]; \\ Antti Karttunen, Sep 15 2023

from sympy import nextprime
from sympy.ntheory import digits
from itertools import count, islice
def b(n): return n - 3**(len(digits(n,3)) - 2)
def agen():
    aset, alst = set(), [None]
    for n in count(1):
        k = b(n)
        if k == 0: an = n
        else:
            ak, p = alst[k], 2
            while p == 3 or p*ak in aset: p = nextprime(p)
            an = p*ak
        yield an; aset.add(an); alst.append(an)
print(list(islice(agen(), 64))) # Michael S. Branicky, Sep 02 2022

a(3^m + 1) = prime(m+2) for m >= 1.
Conjectures from Jianing Song, Nov 23 2022: (Start)
(1) a(3^m+2) = 2*prime(m+2) for m >= 2. - [The conjecture is true because a(2) = 2 and 3^m + 2 < 3^(1+m) + (3^m) + 1 for all m - Antti Karttunen, Sep 16 2023]
(2) For n > m >= 1, a(3^n+3^m+1) = prime(m+2)^2 for n = m+1; prime(n+2)*prime(m+2)^2 for n >= m+2.
(3) For n > m >= 1, a(3^n+3^m+2) = 4*prime(n+2) for n >= 3, m = 1; 2*prime(m+2)^2 for n = m+1, m >= 2; 2*prime(m+2)*prime(m+3) for n = m+2, m >= 2; 2*prime(n+2)*prime(m+2)^2 for n >= m+3, m >= 2. (End)
From Antti Karttunen, Sep 17 2023: (Start)
If A053735(n) = 1, then a(n) = n, otherwise a(n) = A365424(n) * a(A365459(n)).
For all n >= 1, A007949(a(n)) = A007949(n) and a(3*n) = 3*a(n).
For n >= 1, a(3^n - 1) = 2^(2n - 1), a(A048473(n)) = 2^(2*(n-1)).
These are conjectures so far:
For n >= 1, a(3^n - 2) = 10^(n-1).
For n >= 2, a(3^n - 3) = A002023(n-2) = 6*4^(n-2).
(End)

More terms from Michael De Vlieger, Sep 01 2022

A036953 Primes having only {0, 1, 2} as digits.

Original entry on oeis.org

2, 11, 101, 211, 1021, 1201, 2011, 2111, 2221, 10111, 10211, 12011, 12101, 12211, 20011, 20021, 20101, 20201, 21001, 21011, 21101, 21121, 21211, 21221, 22111, 101021, 101111, 101221, 102001, 102101, 102121, 110221, 111121, 111211, 112111

Offset: 1

Patrick De Geest, Jan 04 1999

Number of n-digit terms d(n) = (1, 1, 2, 5, 16, 34, 76, 194, 543, 1469, 4094, 11017, ...); e.g., there are five 4-digit terms: 1021, 1201, 2011, 2111, 2221, hence d(4) = 5. - Zak Seidov, Jun 30 2013

Also, primes in A007089. - M. F. Hasler, Jul 25 2015

Zak Seidov, Table of n, a(n) for n = 1..10000
James Maynard and Brady Haran, Primes without a 7, Numberphile video (2019)
Index to entries for primes with digits in a given set

Subsequence of A107715.
Cf. A036952-A036964.
Cf. A020450 - A020472, A260044, A260267 - A260271, A199325 - A199329, A061247, A199340 - A199349, A217039, A079651.

Select[FromDigits/@Tuples[{0,1,2},6],PrimeQ] (* Harvey P. Dale, Jul 11 2017 *)

lista(n) = {forprime(p=2, n, if (vecmax(digits(p)) <= 2, print1(p, ", ")))} \\ Michel Marcus, Aug 02 2014

A036953={(n,show=0)->for(d=1,1e9,my(u=vector(d,i,10^(d-i))~);forvec(v=vector(d,i,if(i>1,if(iM. F. Hasler, Jul 25 2015

from gmpy2 import digits
from sympy import isprime
[int(digits(n,3)) for n in range(1000) if isprime(int(digits(n,3)))] # Chai Wah Wu, Jul 31 2014

Edited by M. F. Hasler, Jul 25 2015

A054635 Champernowne sequence: write n in base 3 and juxtapose.

Original entry on oeis.org

0, 1, 2, 1, 0, 1, 1, 1, 2, 2, 0, 2, 1, 2, 2, 1, 0, 0, 1, 0, 1, 1, 0, 2, 1, 1, 0, 1, 1, 1, 1, 1, 2, 1, 2, 0, 1, 2, 1, 1, 2, 2, 2, 0, 0, 2, 0, 1, 2, 0, 2, 2, 1, 0, 2, 1, 1, 2, 1, 2, 2, 2, 0, 2, 2, 1, 2, 2, 2, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 2, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 2, 1, 0, 2, 0, 1, 0, 2, 1

Offset: 0

N. J. A. Sloane, Apr 16 2000

Essentially the same as A003137. - R. J. Mathar, Aug 29 2009
An irregular table in which the n-th row lists the base-3 digits of n. - Jason Kimberley, Dec 07 2012
The base-3 Champernowne constant (A077771): it is normal in base 3. - Jason Kimberley, Dec 07 2012

Reinhard Zumkeller, Rows n = 0..1000 of triangle, flattened
Eric Weisstein's World of Mathematics, Ternary Champernowne Constant
Wikipedia, Ternary numeral system

Cf. A054637 (partial sums).
Cf. A081604 (row lengths), A053735 (row sums), A030341 (rows reversed), A007089, A077771.
Table in which the n-th row lists the base b digits of n: A030190 and A030302 (b=2), A003137 and this sequence (b=3), A030373 (b=4), A031219 (b=5), A030548 (b=6), A030998 (b=7), A031035 and A054634 (b=8), A031076 (b=9), A007376 and A033307 (b=10). - Jason Kimberley, Dec 06 2012

a054635 n k = a054635_tabf !! n !! k
a054635_row n = a054635_tabf !! n
a054635_tabf = map reverse a030341_tabf
a054635_list = concat a054635_tabf
-- Reinhard Zumkeller, Feb 21 2013

[0]cat &cat[Reverse(IntegerToSequence(n,3)):n in[1..31]]; // Jason Kimberley, Dec 07 2012

almostNatural[n_, b_] := Block[{m = 0, d = n, i = 1, l, p}, While[m <= d, l = m; m = (b - 1) i*b^(i - 1) + l; i++]; i--; p = Mod[d - l, i]; q = Floor[(d - l)/i] + b^(i - 1); If[p != 0, IntegerDigits[q, b][[p]], Mod[q - 1, b]]]; Array[ almostNatural[#, 3] &, 105, 0] (* Robert G. Wilson v, Jun 29 2014 *)
First[RealDigits[ChampernowneNumber[3], 3, 100, 0]] (* Paolo Xausa, Jun 19 2024 *)

from sympy.ntheory.digits import digits
def agen(limit):
    for n in range(limit):
        yield from digits(n, 3)[1:]
print([an for an in agen(35)]) # Michael S. Branicky, Sep 01 2021

A140267 Nonnegative integers in balanced ternary representation (with 2 standing for -1 digit).

Original entry on oeis.org

0, 1, 12, 10, 11, 122, 120, 121, 102, 100, 101, 112, 110, 111, 1222, 1220, 1221, 1202, 1200, 1201, 1212, 1210, 1211, 1022, 1020, 1021, 1002, 1000, 1001, 1012, 1010, 1011, 1122, 1120, 1121, 1102, 1100, 1101, 1112, 1110, 1111, 12222, 12220, 12221

Offset: 0

Antti Karttunen, May 19 2008, prompted by Eric Angelini's posting on SeqFan mailing list on Sep 15 2005

Sequence A117967 in ternary. (See there for more references.)
From Daniel Forgues, Mar 22 2010: (Start)
The balanced ternary digits {-1, 0, +1} (balanced trits) of a(n) are being represented by {2, 0, 1} respectively in this sequence.
The sign of a(n) is given by the sign of its leading trit.
The number k, k >= 0, of trailing "0"s of a(n) indicates that a(n) is divisible by 3^k.
a(n) is even/odd if it has an even/odd count of nonzero trits. (End)

			For example a(2) = 12, as 1*3 + -1*1 = 2. Similarly, a(19) = 1201, as 1*27 + -1*9 + 0*3 + 1*1 = 19.

Daniel Forgues, Table of n, a(n) for n = 0..100000
Jeff Connelly, Ternary Computing Testbed 3-Trit Computer Architecture, 2008. - _Daniel Forgues_, Mar 23 2010
Brian Hayes, Third Base, American Scientist, November-December 2001. - _Daniel Forgues_, Mar 23 2010
Ternary.info Forum, Balanced ternary arithmetics. - _Daniel Forgues_, Mar 23 2010

a(n) = A007089(A117967(n)). Cf. A140268.

from sympy.ntheory.factor_ import digits
def a004488(n): return int("".join([str((3 - i)%3) for i in digits(n, 3)[1:]]), 3)
def a117968(n):
    if n==1: return 2
    if n%3==0: return 3*a117968(n/3)
    elif n%3==1: return 3*a117968((n - 1)/3) + 2
    else: return 3*a117968((n + 1)/3) + 1
def a117967(n): return 0 if n==0 else a004488(a117968(n))
def a(n): return int("".join(map(str, digits(a117967(n), 3)[1:]))) # Indranil Ghosh, Jun 06 2017

Definition edited by Daniel Forgues, Mar 24 2010

A029984 Numbers k such that k^2 is palindromic in base 3.

Original entry on oeis.org

0, 1, 2, 4, 10, 11, 20, 22, 28, 34, 56, 82, 89, 113, 154, 164, 244, 262, 488, 524, 730, 755, 802, 862, 1021, 1342, 1358, 1460, 2188, 2242, 2684, 2716, 3046, 4276, 4376, 4484, 6562, 6641, 6778, 8030, 8215, 8350, 8887, 12482, 13124, 14810, 19684

Offset: 1

Patrick De Geest

Giovanni Resta, Table of n, a(n) for n = 1..400 (first 100 terms from Harvey P. Dale)
Patrick De Geest, Palindromic Squares in bases 2 to 17
G. J. Simmons, On palindromic squares of non-palindromic numbers, J. Rec. Math., 5 (No. 1, 1972), 11-19. [Annotated scanned copy]

Cf. A002778, A003166, A007089, A029985, A263607, A263608.

pal3Q[n_]:=Module[{idn3=IntegerDigits[n^2,3]},idn3==Reverse[idn3]]; Select[Range[0,20000],pal3Q]  (* Harvey P. Dale, May 22 2012 *)

cp's OEIS Frontend

A074940 Numbers having at least one 2 in their ternary representation.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Python

Formula

Extensions

A081604 Number of digits in ternary representation of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Formula

A003137 Write n in base 3 and juxtapose.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

Python

Extensions

A062153 a(n) = floor(log_3(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A134451 Ternary digital root of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A356867 For n >= 1, write n = 3^m + k, where m >= 0 is the greatest power of 3 <= n, and k is in the range 0 <= k < 3^(m+1) - 3^m, then for n such that k=0, a(n)=n, and for n such that k > 0, a(n) is the smallest prime multiple p*a(k), p != 3, that is not already a term.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI