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.
2^5 = 32 in base 3 = "1012" (= A007089(32)) as 1*27 + 1*3 + 2*1 = 32. 2^6 = 64 in base 3 = "2101" as 2*27 + 1*9 + 1*1 = 64. "1012" reversed is "2101" and vice versa, thus a(5) = 64 and a(6) = 32.
base(n) = {my(a=[n%3]); while(0Altug Alkan, Dec 29 2015
from sympy import factorint
from sympy.ntheory.factor_ import digits
from operator import mul
def a030102(n): return 0 if n==0 else int(''.join(map(str, digits(n, 3)[1:][::-1])), 3)
def a038502(n):
f=factorint(n)
return 1 if n==1 else reduce(mul, [1 if i==3 else i**f[i] for i in f])
def a038500(n): return n/a038502(n)
def a263273(n): return 0 if n==0 else a030102(a038502(n))*a038500(n)
def a263272(n): return a263273(2*n)/2
def a(n): return 2*a263272(a263273(n)) # Indranil Ghosh, May 25 2017
(define (A265342 n) (* 2 (A265351 n)))
from sympy import factorint
from sympy.ntheory.factor_ import digits
from operator import mul
def a030102(n): return 0 if n==0 else int(''.join(map(str, digits(n, 3)[1:][::-1])), 3)
def a038502(n):
f=factorint(n)
return 1 if n==1 else reduce(mul, [1 if i==3 else i**f[i] for i in f])
def a038500(n): return n/a038502(n)
def a263273(n): return 0 if n==0 else a030102(a038502(n))*a038500(n)
def a(n): return a263273(4*n)/2 # Indranil Ghosh, May 23 2017
(define (A264986 n) (A263272 (+ n n)))
from sympy import factorint
from sympy.ntheory.factor_ import digits
from operator import mul
def a030102(n): return 0 if n==0 else int(''.join(map(str, digits(n, 3)[1:][::-1])), 3)
def a038502(n):
f=factorint(n)
return 1 if n==1 else reduce(mul, [1 if i==3 else i**f[i] for i in f])
def a038500(n): return n/a038502(n)
def a263273(n): return 0 if n==0 else a030102(a038502(n))*a038500(n)
def a(n): return a263273(2*(2*n + 1))/2 # Indranil Ghosh, May 23 2017
(define (A264987 n) (A263272 (+ 1 n n)))
Comments