cp's OEIS Frontend

A264987 Odd bisection of A263272.

%I A264987 #8 May 23 2017 20:59:38
%S A264987 1,3,5,11,9,7,13,15,23,29,33,17,35,27,19,37,21,25,31,39,41,95,45,59,
%T A264987 113,69,77,83,87,47,101,99,65,119,51,71,89,105,53,107,81,55,109,57,73,
%U A264987 91,111,43,97,63,61,115,75,79,85,93,49,103,117,67,121,123,203,257,285,149,311,135,167,329,177,221,275
%N A264987 Odd bisection of A263272.
%F A264987 a(n) = A263272((2*n)+1).
%o A264987 (Scheme) (define (A264987 n) (A263272 (+ 1 n n)))
%o A264987 (Python)
%o A264987 from sympy import factorint
%o A264987 from sympy.ntheory.factor_ import digits
%o A264987 from operator import mul
%o A264987 def a030102(n): return 0 if n==0 else int(''.join(map(str, digits(n, 3)[1:][::-1])), 3)
%o A264987 def a038502(n):
%o A264987     f=factorint(n)
%o A264987     return 1 if n==1 else reduce(mul, [1 if i==3 else i**f[i] for i in f])
%o A264987 def a038500(n): return n/a038502(n)
%o A264987 def a263273(n): return 0 if n==0 else a030102(a038502(n))*a038500(n)
%o A264987 def a(n): return a263273(2*(2*n + 1))/2 # _Indranil Ghosh_, May 23 2017
%Y A264987 Cf. A263272, A264986, A264989.
%K A264987 nonn
%O A264987 0,2
%A A264987 _Antti Karttunen_, Dec 05 2015