A284002 -id:A284002 - cp's OEIS frontend

A283477 If 2n = 2^e1 + 2^e2 + ... + 2^ek [e1 .. ek distinct], then a(n) = A002110(e1) * A002110(e2) * ... * A002110(ek).

1, 2, 6, 12, 30, 60, 180, 360, 210, 420, 1260, 2520, 6300, 12600, 37800, 75600, 2310, 4620, 13860, 27720, 69300, 138600, 415800, 831600, 485100, 970200, 2910600, 5821200, 14553000, 29106000, 87318000, 174636000, 30030, 60060, 180180, 360360, 900900, 1801800, 5405400, 10810800, 6306300, 12612600, 37837800, 75675600

Offset: 0

Antti Karttunen, Mar 16 2017

nonn

a(n) = Product of distinct primorials larger than one, obtained as Product_{i} A002110(1+i), where i ranges over the zero-based positions of the 1-bits present in the binary representation of n.
This sequence can be represented as a binary tree. Each child to the left is obtained as A283980(k), and each child to the right is obtained as 2*A283980(k), when their parent contains k:
                                      1
                                      |
                   ...................2....................
                  6                                       12
       30......../ \........60                 180......../ \......360
       / \                 / \                 / \                 / \
      /   \               /   \               /   \               /   \
     /     \             /     \             /     \             /     \
  210       420      1260       2520     6300       12600   37800       75600
etc.

Cf. A129912 (same terms, but sorted into ascending order).
Cf. A000120, A001221, A001222, A002110, A005187, A005361, A006068, A007814, A007913, A019565, A029931, A030308, A046660, A048675, A051903, A056169, A065120, A070939, A072411, A090880, A097248, A108951, A124757, A248663, A276075, A276085, A283475, A283483, A283980, A283984, A283985, A284001, A284002, A284003, A284005, A324287, A324289, A324341, A324342, A324343.
Cf. A005940, A052330, A322827, A323505 for other similar trees.
Cf. also A260443.

Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e == 1 :> {Times @@ Prime@ Range@ PrimePi@ p, e}] &[Times @@ Prime@ Flatten@ Position[#, 1] &@ Reverse@ IntegerDigits[n, 2]], {n, 0, 43}] (* Michael De Vlieger, Mar 18 2017 *)

A283477(n) = prod(i=0,exponent(n),if(bittest(n,i),vecprod(primes(1+i)),1)) \\ Edited by M. F. Hasler, Nov 11 2019

from sympy import prime, primerange, factorint
from operator import mul
from functools import reduce
def P(n): return reduce(mul, [i for i in primerange(2, n + 1)])
def a108951(n):
    f = factorint(n)
    return 1 if n==1 else reduce(mul, [P(i)**f[i] for i in f])
def a019565(n): return reduce(mul, (prime(i+1) for i, v in enumerate(bin(n)[:1:-1]) if v == '1')) if n > 0 else 1 # after Chai Wah Wu
def a(n): return a108951(a019565(n))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 22 2017

from sympy import primorial
from math import prod
def A283477(n): return prod(primorial(i) for i, b in enumerate(bin(n)[:1:-1],1) if b =='1') # Chai Wah Wu, Dec 08 2022

(define (A283477 n) (A108951 (A019565 n)))
;; Recursive "binary tree" implementation, using memoization-macro definec:
(definec (A283477 n) (cond ((zero? n) 1) ((even? n) (A283980 (A283477 (/ n 2)))) (else (* 2 (A283980 (A283477 (/ (- n 1) 2)))))))

a(0) = 1; a(2n) = A283980(a(n)), a(2n+1) = 2*A283980(a(n)).
Other identities. For all n >= 0 (or for n >= 1):
a(2n+1) = 2*a(2n).
a(n) = A108951(A019565(n)).
A097248(a(n)) = A283475(n).
A007814(a(n)) = A051903(a(n)) = A000120(n).
A001221(a(n)) = A070939(n).
A001222(a(n)) = A029931(n).
A048675(a(n)) = A005187(n).
A248663(a(n)) = A006068(n).
A090880(a(n)) = A283483(n).
A276075(a(n)) = A283984(n).
A276085(a(n)) = A283985(n).
A046660(a(n)) = A124757(n).
A056169(a(n)) = A065120(n). [seems to be]
A005361(a(n)) = A284001(n).
A072411(a(n)) = A284002(n).
A007913(a(n)) = A284003(n).
A000005(a(n)) = A284005(n).
A324286(a(n)) = A324287(n).
A276086(a(n)) = A324289(n).
A267263(a(n)) = A324341(n).
A276150(a(n)) = A324342(n). [subsequences in the latter are converging towards this sequence]
G.f.: Product_{k>=0} (1 + prime(k + 1)# * x^(2^k)), where prime()# = A002110. - Ilya Gutkovskiy, Aug 19 2019

More formulas and the binary tree illustration added by Antti Karttunen, Mar 19 2017

Four more linking formulas added by Antti Karttunen, Feb 25 2019

A072411 LCM of exponents in prime factorization of n, a(1) = 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 4, 1, 2, 1, 2, 1, 1, 1, 3, 2, 1, 3, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 4, 2, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 2, 2, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 5, 1, 2, 2, 2, 1, 1, 1, 3, 1, 1, 1, 6, 1, 1, 1, 4, 1, 1, 1, 2, 2, 1, 1, 3

Offset: 1

Labos Elemer, Jun 17 2002

The sums of the first 10^k terms, for k = 1, 2, ..., are 14, 168, 1779, 17959, 180665, 1808044, 18084622, 180856637, 1808585068, 18085891506, ... . Apparently, the asymptotic mean of this sequence is limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1.8085... . - Amiram Eldar, Sep 10 2022

			n = 288 = 2*2*2*2*2*3*3; lcm(5,2) = 10; Product(5,2) = 10, max(5,2) = 5;
n = 180 = 2*2*3*3*5; lcm(2,2,1) = 2; Product(2,2,1) = 4; max(2,2,1) = 2; it deviates both from maximum of exponents (A051903, for the first time at n=72), and product of exponents (A005361, for the first time at n=36).
For n = 36 = 2*2*3*3 = 2^2 * 3^2 we have a(36) = lcm(2,2) = 2.
For n = 72 = 2*2*2*3*3 = 2^3 * 3^2 we have a(72) = lcm(2,3) = 6.
For n = 144 = 2^4 * 3^2 we have a(144) = lcm(2,4) = 4.
For n = 360 = 2^3 * 3^2 * 5^1 we have a(360) = lcm(1,2,3) = 6.

Cf. A028234, A067029, A072412-A072414, A273058, A284569, A290103.
Similar sequences: A001222 (sum of exponents), A005361 (product), A051903 (maximal exponent), A051904 (minimal exponent), A052409 (gcd of exponents), A267115 (bitwise-and), A267116 (bitwise-or), A268387 (bitwise-xor).
Cf. also A055092, A060131.
Differs from A290107 for the first time at n=144.
After the initial term, differs from A157754 for the first time at n=360.

Table[LCM @@ Last /@ FactorInteger[n], {n, 2, 100}] (* Ray Chandler, Jan 24 2006 *)

a(n) = lcm(factor(n)[,2]); \\ Michel Marcus, Mar 25 2017

from sympy import lcm, factorint
def a(n):
    l=[]
    f=factorint(n)
    for i in f: l+=[f[i],]
    return lcm(l)
print([a(n) for n in range(1, 151)]) # Indranil Ghosh, Mar 25 2017

a(1) = 1; for n > 1, a(n) = lcm(A067029(n), a(A028234(n))). - Antti Karttunen, Aug 09 2016
From Antti Karttunen, Aug 22 2017: (Start)
a(n) = A284569(A156552(n)).
a(n) = A290103(A181819(n)).
a(A289625(n)) = A002322(n).
a(A290095(n)) = A055092(n).
a(A275725(n)) = A060131(n).
a(A260443(n)) = A277326(n).
a(A283477(n)) = A284002(n). (End)

a(1) = 1 prepended and the data section filled up to 120 terms by Antti Karttunen, Aug 09 2016

A284001 a(n) = A005361(A283477(n)).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 4, 6, 1, 2, 4, 6, 8, 12, 18, 24, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 32, 48, 72, 96, 108, 144, 192, 240, 162, 216, 288, 360, 384, 480, 600, 720, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 32, 48, 72, 96, 108, 144, 192, 240, 162, 216, 288, 360, 384, 480

Offset: 0

Antti Karttunen, Mar 18 2017

a(n) is the product of elements of the multiset that covers an initial interval of positive integers with multiplicities equal to the parts of the n-th composition in standard order (graded reverse-lexicographic, A066099). This composition is obtained by taking the set of positions of 1's in the reversed binary expansion of n, prepending 0, taking first differences, and reversing again. For example, the 13th composition is (1,2,1) giving the multiset {1,2,2,3} with product 12, so a(13) = 12. - Gus Wiseman, Apr 26 2020

Antti Karttunen, Table of n, a(n) for n = 0..8191
Index entries for sequences related to binary expansion of n

Cf. A005361, A283477, A284002, A284005.
Row products of A095684.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Weighted sum is A029931.
- Necklaces are A065609.
- Sum is A070939.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Lyndon words are A275692.
- Distinct parts are counted by A334028.
Cf. A003963, A034691, A053645, A057335, A305936, A318284, A333764, A333940, A333942, A334030.

Table[Times @@ FactorInteger[#][[All, -1]] &[Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e == 1 :> {Times @@ Prime@ Range@ PrimePi@ p, e}] &[Times @@ Prime@ Flatten@ Position[#, 1] &@ Reverse@ IntegerDigits[n, 2]]], {n, 0, 93}] (* Michael De Vlieger, Mar 18 2017 *)

A005361(n) = factorback(factor(n)[, 2]); \\ From A005361
A034386(n) = prod(i=1, primepi(n), prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
A283477(n) = A108951(A019565(n));
A284001(n) = A005361(A283477(n));

(define (A284001 n) (A005361 (A283477 n)))

a(n) = A005361(A283477(n)).
a(n) = A003963(A057335(n)). - Gus Wiseman, Apr 26 2020
a(n) = A284005(A053645(n)) for n > 0 with a(0) = 1. - Mikhail Kurkov, Jun 05 2021 [verification needed]

A329378 Least common multiple of exponents of prime factors of A108951(n), where A108951 is fully multiplicative with a(prime(i)) = prime(i)# = Product_{i=1..i} A000040(i).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 3, 1, 2, 2, 4, 1, 6, 1, 3, 2, 2, 1, 4, 2, 2, 3, 3, 1, 6, 1, 5, 2, 2, 2, 4, 1, 2, 2, 4, 1, 6, 1, 3, 3, 2, 1, 5, 2, 6, 2, 3, 1, 12, 2, 4, 2, 2, 1, 4, 1, 2, 3, 6, 2, 6, 1, 3, 2, 6, 1, 10, 1, 2, 6, 3, 2, 6, 1, 5, 4, 2, 1, 4, 2, 2, 2, 4, 1, 12, 2, 3, 2, 2, 2, 6, 1, 6, 3, 4, 1, 6, 1, 4, 6

Offset: 1

Antti Karttunen, Nov 17 2019

nonn

Cf. A019565, A034386, A072411, A108951, A284002, A329382, A329600, A329605.

Differs from related A329617 for the first time at n=36.

A034386(n) = prod(i=1, primepi(n), prime(i));
A072411(n) = lcm(factor(n)[, 2]); \\ From A072411
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
A329378(n) = A072411(A108951(n));

a(n) = A072411(A108951(n)) = A072411(A329600(n)).
a(n) <= A329617(n) <= A329382(n) <= A329605(n).
a(A019565(n)) = A284002(n).

cp's OEIS Frontend

A283477 If 2n = 2^e1 + 2^e2 + ... + 2^ek [e1 .. ek distinct], then a(n) = A002110(e1) * A002110(e2) * ... * A002110(ek).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Python

Python

Scheme

Formula

Extensions

A072411 LCM of exponents in prime factorization of n, a(1) = 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

Extensions

A284001 a(n) = A005361(A283477(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Formula

A329378 Least common multiple of exponents of prime factors of A108951(n), where A108951 is fully multiplicative with a(prime(i)) = prime(i)# = Product_{i=1..i} A000040(i).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula