cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A014263 Numbers that contain even digits only.

Original entry on oeis.org

0, 2, 4, 6, 8, 20, 22, 24, 26, 28, 40, 42, 44, 46, 48, 60, 62, 64, 66, 68, 80, 82, 84, 86, 88, 200, 202, 204, 206, 208, 220, 222, 224, 226, 228, 240, 242, 244, 246, 248, 260, 262, 264, 266, 268, 280, 282, 284, 286, 288, 400, 402, 404, 406, 408, 420, 422, 424
Offset: 1

Views

Author

N. J. A. Sloane

Keywords

Comments

The set of real numbers between 0 and 1 that contain no odd digits in their decimal expansion has Hausdorff dimension log 5 / log 10.
Integers written in base 5 and then doubled (in base 10). - Franklin T. Adams-Watters, Mar 15 2006
The carryless mod 10 "even" numbers (cf. A004529) sorted and duplicates removed. - N. J. A. Sloane, Aug 03 2010.
Complement of A007957; A196564(a(n)) = 0; A103181(a(n)) = 0. - Reinhard Zumkeller, Oct 04 2011
If n-1 is represented as a base-5 number (see A007091) according to n-1 = d(m)d(m-1)…d(3)d(2)d(1)d(0) then a(n)= Sum_{j=0..m} c(d(j))*10^j, where c(k)=0,2,4,6,8 for k=0..4. - Hieronymus Fischer, Jun 03 2012

Examples

			a(1000) = 24888.
a(10^4) = 60888.
a(10^5) = 22288888.
a(10^6) = 446888888.

References

  • K. J. Falconer, The Geometry of Fractal Sets, Cambridge, 1985; p. 19.

Links

Crossrefs

Subsequence of A059708.
Cf. A061810, A061811, A007091, A014261, A046034, A052382, A084544, A089581, A084984, A017042, A001743, A202267, A202268, A194182, A196563.

Programs

  • Haskell
    a014263 n = a014263_list !! (n-1)
a014263_list = filter (all (`elem` "02468") . show) [0,2..]
-- Reinhard Zumkeller, Jul 05 2011
  • Magma
    [n: n in [0..424] | Set(Intseq(n)) subset [0..8 by 2]];  // Bruno Berselli, Jul 19 2011
  • Maple
    a:= proc(m) local L,i;
  L:= convert(m-1,base,5);
  2*add(L[i]*10^(i-1),i=1..nops(L))
end proc:
seq(a(i),i=1..100); # Robert Israel, Apr 07 2016
  • Mathematica
    Select[Range[450], And@@EvenQ[IntegerDigits[#]]&] (* Harvey P. Dale, Jan 30 2011 *)
FromDigits/@Tuples[{0,2,4,6,8},3] (* Harvey P. Dale, Jul 07 2025 *)
  • PARI
    a(n) = 2*fromdigits(digits(n-1, 5), 10); \\ Michel Marcus, Nov 04 2022
  • PARI
    is(n)=#setminus(Set(digits(n)), [0,2,4,6,8])==0 \\ Charles R Greathouse IV, Mar 03 2025
  • Python
    from sympy.ntheory.digits import digits
def a(n): return int(''.join(str(2*d) for d in digits(n, 5)[1:]))
print([a(n) for n in range(58)]) # Michael S. Branicky, Jan 13 2022
  • Python
    from itertools import count, islice, product
def agen(): # generator of terms
    yield 0
    for d in count(1):
        for first in "2468":
            for rest in product("02468", repeat=d-1):
                yield int(first + "".join(rest))
print(list(islice(agen(), 58))) # Michael S. Branicky, Jan 13 2022

Formula

A045888(a(n)) = 0. - Reinhard Zumkeller, Aug 25 2009
a(n) = A179082(n) for n <= 25. - Reinhard Zumkeller, Jun 28 2010
From Hieronymus Fischer, Jun 06 2012: (Start)
a(n) = ((2*b_m(n)) mod 8 + 2)*10^m + Sum_{j=0..m-1} ((2*b_j(n)) mod 10)*10^j, where n>1, b_j(n) = floor((n-1-5^m)/5^j), m = floor(log_5(n-1)).
a(1*5^n+1) = 2*10^n.
a(2*5^n+1) = 4*10^n.
a(3*5^n+1) = 6*10^n.
a(4*5^n+1) = 8*10^n.
a(n) = 2*10^log_5(n-1) for n=5^k+1,
a(n) < 2*10^log_5(n-1), else.
a(n) > (8/9)*10^log_5(n-1) n>1.
a(n) = 2*A007091(n-1), iff the digits of A007091(n-1) are 0 or 1.
G.f.: g(x) = (x/(1-x))*Sum_{j>=0} 10^j*x^5^j *(1-x^5^j)* (2+4x^5^j+ 6(x^2)^5^j+ 8(x^3)^5^j)/(1-x^5^(j+1)).
Also: g(x) = 2*(x/(1-x))*Sum_{j>=0} 10^j*x^5^j * (1-4x^(3*5^j)+3x^(4*5^j))/((1-x^5^j)(1-x^5^(j+1))).
Also: g(x) = 2*(x/(1-x))*(h_(5,1)(x) + h_(5,2)(x) + h_(5,3)(x) + h_(5,4)(x) - 4*h_(5,5)(x)), where h_(5,k)(x) = Sum_{j>=0} 10^j*(x^5^j)^k/(1-(x^5^j)^5). (End)
a(5*n+i-4) = 10*a(n) + 2*i for n >= 1, i=0..4. - Robert Israel, Apr 07 2016
Sum_{n>=2} 1/a(n) = A194182. - Bernard Schott, Jan 13 2022

Extensions

Examples and crossrefs added by Hieronymus Fischer, Jun 06 2012

A196564 Number of odd digits in decimal representation of n.

Original entry on oeis.org

0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 0, 1, 0, 1, 0, 1
Offset: 0

Views

Author

Reinhard Zumkeller, Oct 04 2011

Keywords

Links

Crossrefs

Cf. A014261, A014263, A027868, A046034, A055640, A055641, A055642, A061217, A102669-A102685, A122640, A196563.

Programs

  • Haskell
    a196564 n = length [d | d <- show n, d `elem` "13579"]
-- Reinhard Zumkeller, Feb 22 2012, Oct 04 2011
  • Maple
    A196564 := proc(n)
        if n =0 then
                0;
        else
                convert(n,base,10) ;
                add(d mod 2,d=%) ;
        end if:
end proc: # R. J. Mathar, Jul 13 2012
  • Mathematica
    Table[Total[Mod[IntegerDigits[n],2]],{n,0,100}] (* Zak Seidov, Oct 13 2015 *)
  • PARI
    a(n) = #select(x->x%2, digits(n)); \\ Michel Marcus, Oct 14 2015
  • Python
    def a(n): return sum(1 for d in str(n) if d in "13579")
print([a(n) for n in range(100)]) # Michael S. Branicky, May 15 2022

Formula

a(n) = A055642(n) - A196563(n);
a(A014263(n)) = 0; a(A007957(n)) > 0.
From Hieronymus Fischer, May 30 2012: (Start)
a(n) = Sum_{j=0..m} (floor(n/(2*10^j) + (1/2)) - floor(n/(2*10^j))), where m=floor(log_10(n)).
a(10*n+k) = a(n) + a(k), 0<=k<10, n>=0.
a(n) = a(floor(n/10)) + a(n mod 10), n>=0.
a(n) = Sum_{j=0..m} a(floor(n/10^j) mod 10), n>=0.
a(A014261(n)) = floor(log_5(4*n+1)), n>0.
G.f.: g(x) = (1/(1-x))*Sum_{j>=0} x^10^j/(1+x^10^j). (End)

A194377 Numbers m such that Sum_{k=1..m} (<1/2 + k*r> - ) > 0, where r=sqrt(6) and < > denotes fractional part.

Original entry on oeis.org

1, 3, 5, 7, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 21, 23, 25, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 43, 45, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 61, 63, 65, 67, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 81, 83, 85, 87, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100
Offset: 1

Views

Author

Clark Kimberling, Aug 23 2011

Keywords

Comments

See A194368. Although a(n)=A007957(n) for n = 1..70, the number 208, for example, is here but not A007957.

Crossrefs

Cf. A007957, A194368, A194376.

Programs

  • Mathematica
    r = Sqrt[6]; c = 1/2;
x[n_] := Sum[FractionalPart[k*r], {k, 1, n}]
y[n_] := Sum[FractionalPart[c + k*r], {k, 1, n}]
t1 = Table[If[y[n] < x[n], 1, 0], {n, 1, 500}];
Flatten[Position[t1, 1]]   (* empty *)
t2 = Table[If[y[n] == x[n], 1, 0], {n, 1, 400}];
Flatten[Position[t2, 1]]   (* A194376 *)
t3 = Table[If[y[n] > x[n], 1, 0], {n, 1, 100}];
Flatten[Position[t3, 1]]   (* A194377 *)
  • PARI
    is(n)=my(r=sqrt(6),f=x->x-x\1);sum(k=1,n,f(1/2+k*r)-f(k*r))>0 \\ Charles R Greathouse IV, Jul 25 2012
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.