|
%I #41 Mar 12 2021 07:53:17
%S 2,6,18,30,90,270,450,210,630,6750,20250,9450,15750,47250,22050,2310,
%T 6930,330750,3543750,1653750,4961250,53156250,24806250,727650,1212750,
%U 57881250,173643750,18191250,8489250,25467750,2668050,30030,90090,40020750,1910081250,891371250,9550406250,455814843750,212713593750
%N Odd bisection of A260443 (the even terms): a(n) = A260443((2*n)+1).
%C From _David A. Corneth_, Oct 22 2016: (Start)
%C The exponents of the prime factorization of a(n) are first nondecreasing, then nonincreasing.
%C The exponent of 2 in the prime factorization of a(n) is 1. (End)
%H Antti Karttunen, <a href="/A277324/b277324.txt">Table of n, a(n) for n = 0..1024</a>
%F a(n) = A260443((2*n)+1).
%F a(0) = 2; for n >= 1, a(n) = A260443(n) * A260443(n+1).
%F Other identities. For all n >= 0:
%F A007949(a(n)) = A005811(n). [See comments in A125184.]
%F A156552(a(n)) = A277189(n), a(n) = A005940(1+A277189(n)).
%F A048675(a(n)) = 2n + 1. - _David A. Corneth_, Oct 22 2016
%F A001222(a(n)) = A007306(1+n).
%F A056169(a(n)) = A284267(n).
%F A275812(a(n)) = A284268(n).
%F A248663(a(n)) = A283975(n).
%F A000188(a(n)) = A283484(n).
%F A247503(a(n)) = A284563(n).
%F A248101(a(n)) = A284564(n).
%F A046523(a(n)) = A284573(n).
%F a(n) = A277198(n) * A284008(n).
%F a(n) = A284576(n) * A284578(n) = A284577(n) * A000290(A284578(n)).
%e A method to find terms of this sequence, explained by an example to find a(7). To find k = a(7), we find k such that A048675(k) = 2*7+1 = 15. 7 has the binary partitions: {[7, 0, 0], [5, 1, 0], [3, 2, 0], [1, 3, 0], [3, 0, 1], [1, 1, 1]}. To each of those, we prepend a 1. This gives the binary partitions of 15 starting with a 1. For example, for the first we get [1, 7, 0, 0]. We see that only [1, 5, 1, 0], [1, 3, 2, 0] and [1, 1, 1, 1] start nondecreasing, then nonincreasing, so we only check those. These numbers will be the exponents in a prime factorization. [1, 5, 1, 0] corresponds to prime(1)^1 * prime(2)^5 * prime(3)^1 * prime(4)^0 = 2430. We find that [1, 1, 1, 1] gives k = 210 for which A048675(k) = 15 so a(7) = 210. - _David A. Corneth_, Oct 22 2016
%t a[n_] := a[n] = Which[n < 2, n + 1, EvenQ@ n, Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[# == 1] &@ a[n/2], True, a[#] a[# + 1] &[(n - 1)/2]]; Table[a[2 n + 1], {n, 0, 38}] (* _Michael De Vlieger_, Apr 05 2017 *)
%o (Scheme, two versions)
%o (define (A277324 n) (A260443 (+ 1 n n)))
%o (define (A277324 n) (if (zero? n) 2 (* (A260443 n) (A260443 (+ 1 n)))))
%o (Python)
%o from sympy import factorint, prime, primepi
%o from operator import mul
%o def a003961(n):
%o F=factorint(n)
%o return 1 if n==1 else reduce(mul, [prime(primepi(i) + 1)**F[i] for i in F])
%o def a260443(n): return n + 1 if n<2 else a003961(a260443(n//2)) if n%2==0 else a260443((n - 1)//2)*a260443((n + 1)//2)
%o def a(n): return a260443(2*n + 1)
%o print([a(n) for n in range(101)]) # _Indranil Ghosh_, Jun 21 2017
%Y Cf. A005811, A005940, A007306, A007949, A156552, A260443, A277189, A277323, A283484, A283975, A284267, A284268, A284563, A284564, A284573.
%Y Cf. A277200 (same sequence sorted into ascending order).
%K nonn
%O 0,1
%A _Antti Karttunen_, Oct 10 2016
%E More linking formulas added by _Antti Karttunen_, Apr 16 2017
|
