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

A348717 a(n) is the least k such that A003961^i(k) = n for some i >= 0 (where A003961^i denotes the i-th iterate of A003961).

Original entry on oeis.org

1, 2, 2, 4, 2, 6, 2, 8, 4, 10, 2, 12, 2, 14, 6, 16, 2, 18, 2, 20, 10, 22, 2, 24, 4, 26, 8, 28, 2, 30, 2, 32, 14, 34, 6, 36, 2, 38, 22, 40, 2, 42, 2, 44, 12, 46, 2, 48, 4, 50, 26, 52, 2, 54, 10, 56, 34, 58, 2, 60, 2, 62, 20, 64, 14, 66, 2, 68, 38, 70, 2, 72, 2
Offset: 1

Views

Author

Rémy Sigrist, Oct 31 2021

Keywords

Comments

All terms except a(1) = 1 are even.
To compute a(n) for n > 1:
- if n = Product_{j = 1..o} prime(p_j)^e_j (where prime(i) denotes the i-th prime number, p_1 < ... < p_o and e_1 > 0)
- then a(n) = Product_{j = 1..o} prime(p_j + 1 - p_1)^e_j.
This sequence has similarities with A304776: here we shift down prime indexes, there prime exponents.
Smallest number generated by uniformly decrementing the indices of the prime factors of n. Thus, for n > 1, the smallest m > 1 such that the first differences of the indices of the ordered prime factors (including repetitions) are the same for m and n. As a function, a(.) preserves properties such as prime signature. - Peter Munn, May 12 2022

Links

Crossrefs

Positions of particular values (see formula section): A000040, A001248, A006094, A030078, A030514, A046301, A050997, A090076, A090090, A166329, A251720.
Also see formula section for the relationship to: A000265, A003961, A004277, A005940, A020639, A046523, A055396, A071364, A122111, A156552, A243055, A243074, A297845, A322993.
Sequences with comparable definitions: A304776, A316437.
Cf. A246277 (terms halved), A305897 (restricted growth sequence transform), A354185 (Möbius transform), A354186 (Dirichlet inverse), A354187 (sum with it).

Programs

  • Mathematica
    a[1] = 1; a[n_] := Module[{f = FactorInteger[n], d}, d = PrimePi[f[[1, 1]]] - 1; Times @@ ((Prime[PrimePi[#[[1]]] - d]^#[[2]]) & /@ f)]; Array[a, 100] (* Amiram Eldar, Oct 31 2021 *)
  • PARI
    a(n) = { my (f=factor(n)); if (#f~>0, my (pi1=primepi(f[1,1])); for (k=1, #f~, f[k,1] = prime(primepi(f[k,1])-pi1+1))); factorback(f) }

Formula

a(n) = n iff n belongs to A004277.
A003961^(A055396(n)-1)(a(n)) = n for any n > 1.
a(n) = 2 iff n belongs to A000040 (prime numbers).
a(n) = 4 iff n belongs to A001248 (squares of prime numbers).
a(n) = 6 iff n belongs to A006094 (products of two successive prime numbers).
a(n) = 8 iff n belongs to A030078 (cubes of prime numbers).
a(n) = 10 iff n belongs to A090076.
a(n) = 12 iff n belongs to A251720.
a(n) = 14 iff n belongs to A090090.
a(n) = 16 iff n belongs to A030514.
a(n) = 30 iff n belongs to A046301.
a(n) = 32 iff n belongs to A050997.
a(n) = 36 iff n belongs to A166329.
a(1) = 1, for n > 1, a(n) = 2*A246277(n). - Antti Karttunen, Feb 23 2022
a(n) = A122111(A243074(A122111(n))). - Peter Munn, Feb 23 2022
From Peter Munn and Antti Karttunen, May 12 2022: (Start)
a(1) = 1; a(2n) = 2n; a(A003961(n)) = a(n). [complete definition]
a(n) = A005940(1+A322993(n)) = A005940(1+A000265(A156552(n))).
Equivalently, A156552(a(n)) = A000265(A156552(n)).
A297845(a(n), A020639(n)) = n.
A046523(a(n)) = A046523(n).
A071364(a(n)) = A071364(n).
a(n) >= A071364(n).
A243055(a(n)) = A243055(n).
(End)

A316431 Least common multiple divided by greatest common divisor of the integer partition with Heinz number n > 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 4, 6, 1, 1, 2, 1, 3, 2, 5, 1, 2, 1, 6, 1, 4, 1, 6, 1, 1, 10, 7, 12, 2, 1, 8, 3, 3, 1, 4, 1, 5, 6, 9, 1, 2, 1, 3, 14, 6, 1, 2, 15, 4, 4, 10, 1, 6, 1, 11, 2, 1, 2, 10, 1, 7, 18, 12, 1, 2, 1, 12, 6, 8, 20, 6, 1, 3, 1, 13, 1, 4, 21, 14, 5, 5, 1, 6, 6, 9, 22, 15, 24, 2, 1, 4, 10, 3, 1, 14, 1, 6, 12
Offset: 2

Views

Author

Gus Wiseman, Jul 02 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			63 is the Heinz number of (4,2,2), which has LCM 4 and GCD 2, so a(63) = 4/2 = 2.
91 is the Heinz number of (6,4), which has LCM 12 and GCD 2, so a(91) = 12/2 = 6.

Links

Crossrefs

Cf. A005117, A056239, A074761, A289508, A289509, A290103, A296150, A316429, A316430, A316437.

Programs

  • Mathematica
    Table[With[{pms=Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]]},LCM@@pms/GCD@@pms],{n,2,100}]
  • PARI
    A316431(n) = if(1==n,1,my(pis = apply(p -> primepi(p), factor(n)[, 1]~)); lcm(pis)/gcd(pis)); \\ Antti Karttunen, Sep 06 2018

Formula

a(n) = A290103(n)/A289508(n).
a(n) = a(A005117(n)). - David A. Corneth, Sep 06 2018

Extensions

More terms from Antti Karttunen, Sep 06 2018

A316438 Heinz numbers of integer partitions whose product is strictly greater than the LCM of the parts.

Original entry on oeis.org

9, 18, 21, 25, 27, 36, 39, 42, 45, 49, 50, 54, 57, 63, 65, 72, 75, 78, 81, 84, 87, 90, 91, 98, 99, 100, 105, 108, 111, 114, 115, 117, 121, 125, 126, 129, 130, 133, 135, 144, 147, 150, 153, 156, 159, 162, 168, 169, 171, 174, 175, 180, 182, 183, 185, 189, 195
Offset: 1

Views

Author

Gus Wiseman, Jul 03 2018

Keywords

Comments

Also numbers n > 1 such that A290104(n) > 1.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Sequence of partitions whose product is greater than their LCM begins: (22), (221), (42), (33), (222), (2211), (62), (421), (322), (44), (331), (2221), (82), (422), (63), (22111), (332), (621), (2222), (4211).

Crossrefs

Cf. A056239, A074761, A143773, A285572, A289508, A290103, A296150, A316429, A316431, A316433, A316437.

Programs

  • Mathematica
    Select[Range[2,300],With[{pms=Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]},Times@@pms/LCM@@pms>1]&]

A316436 Sum divided by GCD of the integer partition with Heinz number n > 1.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 3, 2, 4, 1, 4, 1, 5, 5, 4, 1, 5, 1, 5, 3, 6, 1, 5, 2, 7, 3, 6, 1, 6, 1, 5, 7, 8, 7, 6, 1, 9, 4, 6, 1, 7, 1, 7, 7, 10, 1, 6, 2, 7, 9, 8, 1, 7, 8, 7, 5, 11, 1, 7, 1, 12, 4, 6, 3, 8, 1, 9, 11, 8, 1, 7, 1, 13, 8, 10, 9, 9, 1, 7, 4, 14, 1, 8, 10, 15, 6, 8, 1, 8, 5, 11, 13, 16, 11, 7, 1, 9, 9, 8, 1, 10, 1, 9, 9
Offset: 2

Views

Author

Gus Wiseman, Jul 03 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Links

Crossrefs

Cf. A056239, A289508, A289509, A290103, A290104, A296150, A316430, A316431, A316432, A316437.

Programs

  • Maple
    a:= n-> (l-> add(i, i=l)/igcd(l[]))(map(i->
      numtheory[pi](i[1])$i[2], ifactors(n)[2])):
seq(a(n), n=2..100);  # Alois P. Heinz, Jul 03 2018
  • Mathematica
    Table[With[{pms=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]},Total[pms]/GCD@@pms],{n,2,100}]
  • PARI
    A316436(n) = { my(f = factor(n), pis = apply(p -> primepi(p), f[, 1]~), es = f[, 2]~, g = gcd(pis)); sum(i=1, #f~, pis[i]*es[i])/g; }; \\ Antti Karttunen, Sep 10 2018

Extensions

More terms from Antti Karttunen, Sep 10 2018
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.