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-2 of 2 results.

A108729 Numbers k that divide A076127(k); or k such that A076127(k) is a Niven number.

Original entry on oeis.org

1, 3, 7, 9, 21, 27, 57, 81, 111, 189, 243, 333, 441, 729, 777, 999, 1161, 1701, 2109, 2187, 2529, 2709, 2763, 2793, 2997
Offset: 1

Views

Author

Jason Earls, Jun 21 2005

Keywords

Examples

			a(2)=3 because A076127(3)=101001 and 101001/3 = 33667.
		

Crossrefs

Extensions

More terms from David Wasserman, May 22 2008

A076131 a(n) = 2^n*a(n-1) + 1, a(0) = 0.

Original entry on oeis.org

0, 1, 5, 41, 657, 21025, 1345601, 172236929, 44092653825, 22575438758401, 23117249288602625, 47344126543058176001, 193921542320366288900097, 1588605274688440638669594625, 26027708820495411423962638336001, 852875962629993641540407732994080769
Offset: 0

Views

Author

Kyle Hunter (hunterk(AT)raytheon.com), Oct 31 2002

Keywords

Comments

Base-2 expansion is same as base 10 expansion of A076127.

Crossrefs

Programs

  • Mathematica
    a[0] = 0; a[n_] := 2^n a[n - 1] + 1; Table[ a[n], {n, 0, 13}]
  • PARI
    a(n)=if(n<0,0,subst(Polrev(Vec(sum(k=1,n,x^(k*(k+1)/2)))),x,2))
    
  • PARI
    a(n)=if(n<1,0,1+a(n-1)*2^n)

Formula

a(n) = floor(c*2^((n+1)*(n+2)/2)) where c = sum(k>=1, 1/2^A000217(k))=0.6416325... - Benoit Cloitre, Nov 01 2002

Extensions

Edited and extended by Robert G. Wilson v, Oct 31 2002
Formula corrected by Vaclav Kotesovec, Aug 11 2012
Showing 1-2 of 2 results.