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-10 of 17 results. Next

A122111 Self-inverse permutation of the positive integers induced by partition enumeration in A112798 and partition conjugation.

Original entry on oeis.org

1, 2, 4, 3, 8, 6, 16, 5, 9, 12, 32, 10, 64, 24, 18, 7, 128, 15, 256, 20, 36, 48, 512, 14, 27, 96, 25, 40, 1024, 30, 2048, 11, 72, 192, 54, 21, 4096, 384, 144, 28, 8192, 60, 16384, 80, 50, 768, 32768, 22, 81, 45, 288, 160, 65536, 35, 108, 56, 576, 1536, 131072, 42
Offset: 1

Views

Author

Keywords

Comments

Factor n; replace each prime(i) with i, take the conjugate partition, replace parts i with prime(i) and multiply out.
From Antti Karttunen, May 12-19 2014: (Start)
For all n >= 1, A001222(a(n)) = A061395(n), and vice versa, A061395(a(n)) = A001222(n).
Because the partition conjugation doesn't change the partition's total sum, this permutation preserves A056239, i.e., A056239(a(n)) = A056239(n) for all n.
(Similarly, for all n, A001221(a(n)) = A001221(n), because the number of steps in the Ferrers/Young-diagram stays invariant under the conjugation. - Note added Apr 29 2022).
Because this permutation commutes with A241909, in other words, as a(A241909(n)) = A241909(a(n)) for all n, from which follows, because both permutations are self-inverse, that a(n) = A241909(a(A241909(n))), it means that this is also induced when partitions are conjugated in the partition enumeration system A241918. (Not only in A112798.)
(End)
From Antti Karttunen, Jul 31 2014: (Start)
Rows in arrays A243060 and A243070 converge towards this sequence, and also, assuming no surprises at the rate of that convergence, this sequence occurs also as the central diagonal of both.
Each even number is mapped to a unique term of A102750 and vice versa.
Conversely, each odd number (larger than 1) is mapped to a unique term of A070003, and vice versa. The permutation pair A243287-A243288 has the same property. This is also used to induce the permutations A244981-A244984.
Taking the odd bisection and dividing out the largest prime factor results in the permutation A243505.
Shares with A245613 the property that each term of A028260 is mapped to a unique term of A244990 and each term of A026424 is mapped to a unique term of A244991.
Conversely, with A245614 (the inverse of above), shares the property that each term of A244990 is mapped to a unique term of A028260 and each term of A244991 is mapped to a unique term of A026424.
(End)
The Maple program follows the steps described in the first comment. The subprogram C yields the conjugate partition of a given partition. - Emeric Deutsch, May 09 2015
The Heinz number of the partition that is conjugate to the partition with Heinz number n. The Heinz number of a partition p = [p_1, p_2, ..., p_r] is defined as Product(p_j-th prime, j=1...r). Example: a(3) = 4. Indeed, the partition with Heinz number 3 is [2]; its conjugate is [1,1] having Heinz number 4. - Emeric Deutsch, May 19 2015

Crossrefs

Cf. A088902 (fixed points).
Cf. A112798, A241918 (conjugates the partitions listed in these two tables).
Cf. A243060 and A243070. (Limit of rows in these arrays, and also their central diagonal).
Cf. A319988 (parity of this sequence for n > 1), A336124 (a(n) mod 4).
{A000027, A122111, A241909, A241916} form a 4-group.
{A000027, A122111, A153212, A242419} form also a 4-group.
Cf. also array A350066 [A(i, j) = a(a(i)*a(j))].

Programs

  • Maple
    with(numtheory): c := proc (n) local B, C: B := proc (n) local pf: pf := op(2, ifactors(n)): [seq(seq(pi(op(1, op(i, pf))), j = 1 .. op(2, op(i, pf))), i = 1 .. nops(pf))] end proc: C := proc (P) local a: a := proc (j) local c, i: c := 0; for i to nops(P) do if j <= P[i] then c := c+1 else  end if end do: c end proc: [seq(a(k), k = 1 .. max(P))] end proc: mul(ithprime(C(B(n))[q]), q = 1 .. nops(C(B(n)))) end proc: seq(c(n), n = 1 .. 59); # Emeric Deutsch, May 09 2015
    # second Maple program:
    a:= n-> (l-> mul(ithprime(add(`if`(jAlois P. Heinz, Sep 30 2017
  • Mathematica
    A122111[1] = 1; A122111[n_] := Module[{l = #, m = 0}, Times @@ Power @@@ Table[l -= m; l = DeleteCases[l, 0]; {Prime@Length@l, m = Min@l}, Length@Union@l]] &@Catenate[ConstantArray[PrimePi[#1], #2] & @@@ FactorInteger@n]; Array[A122111, 60] (* JungHwan Min, Aug 22 2016 *)
    a[n_] := Function[l, Product[Prime[Sum[If[jJean-François Alcover, Sep 23 2020, after Alois P. Heinz *)
  • PARI
    A122111(n) = if(1==n,n,my(f=factor(n), es=Vecrev(f[,2]),is=concat(apply(primepi,Vecrev(f[,1])),[0]),pri=0,m=1); for(i=1, #es, pri += es[i]; m *= prime(pri)^(is[i]-is[1+i])); (m)); \\ Antti Karttunen, Jul 20 2020
    
  • Python
    from sympy import factorint, prevprime, prime, primefactors
    from operator import mul
    def a001222(n): return 0 if n==1 else a001222(n/primefactors(n)[0]) + 1
    def a064989(n):
        f=factorint(n)
        return 1 if n==1 else reduce(mul, [1 if i==2 else prevprime(i)**f[i] for i in f])
    def a105560(n): return 1 if n==1 else prime(a001222(n))
    def a(n): return 1 if n==1 else a105560(n)*a(a064989(n))
    [a(n) for n in range(1, 101)] # Indranil Ghosh, Jun 15 2017
  • Scheme
    ;; Uses Antti Karttunen's IntSeq-library.
    (definec (A122111 n) (if (<= n 1) n (* (A000040 (A001222 n)) (A122111 (A064989 n)))))
    ;; Antti Karttunen, May 12 2014
    
  • Scheme
    ;; Uses Antti Karttunen's IntSeq-library.
    (definec (A122111 n) (if (<= n 1) n (* (A000079 (A241917 n)) (A003961 (A122111 (A052126 n))))))
    ;; Antti Karttunen, May 12 2014
    
  • Scheme
    ;; Uses Antti Karttunen's IntSeq-library.
    (definec (A122111 n) (if (<= n 1) n (* (expt (A000040 (A071178 n)) (A241919 n)) (A242378bi (A071178 n) (A122111 (A051119 n))))))
    ;; Antti Karttunen, May 12 2014
    

Formula

From Antti Karttunen, May 12-19 2014: (Start)
a(1) = 1, a(p_i) = 2^i, and for other cases, if n = p_i1 * p_i2 * p_i3 * ... * p_{k-1} * p_k, where p's are primes, not necessarily distinct, sorted into nondescending order so that i1 <= i2 <= i3 <= ... <= i_{k-1} <= ik, then a(n) = 2^(ik-i_{k-1}) * 3^(i_{k-1}-i_{k-2}) * ... * p_{i_{k-1}}^(i2-i1) * p_ik^(i1).
This can be implemented as a recurrence, with base case a(1) = 1,
and then using any of the following three alternative formulas:
a(n) = A105560(n) * a(A064989(n)) = A000040(A001222(n)) * a(A064989(n)). [Cf. the formula for A242424.]
a(n) = A000079(A241917(n)) * A003961(a(A052126(n))).
a(n) = (A000040(A071178(n))^A241919(n)) * A242378(A071178(n), a(A051119(n))). [Here ^ stands for the ordinary exponentiation, and the bivariate function A242378(k,n) changes each prime p(i) in the prime factorization of n to p(i+k), i.e., it's the result of A003961 iterated k times starting from n.]
a(n) = 1 + A075157(A129594(A075158(n-1))). [Follows from the commutativity with A241909, please see the comments section.]
(End)
From Antti Karttunen, Jul 31 2014: (Start)
As a composition of related permutations:
a(n) = A153212(A242419(n)) = A242419(A153212(n)).
a(n) = A241909(A241916(n)) = A241916(A241909(n)).
a(n) = A243505(A048673(n)).
a(n) = A064216(A243506(n)).
Other identities. For all n >= 1, the following holds:
A006530(a(n)) = A105560(n). [The latter sequence gives greatest prime factor of the n-th term].
a(2n)/a(n) = A105560(2n)/A105560(n), which is equal to A003961(A105560(n))/A105560(n) when n > 1.
A243505(n) = A052126(a(2n-1)) = A052126(a(4n-2)).
A066829(n) = A244992(a(n)) and vice versa, A244992(n) = A066829(a(n)).
A243503(a(n)) = A243503(n). [Because partition conjugation does not change the partition size.]
A238690(a(n)) = A238690(n). - per Matthew Vandermast's note in that sequence.
A238745(n) = a(A181819(n)) and a(A238745(n)) = A181819(n). - per Matthew Vandermast's note in A238745.
A181815(n) = a(A181820(n)) and a(A181815(n)) = A181820(n). - per Matthew Vandermast's note in A181815.
(End)
a(n) = A181819(A108951(n)). [Prime shadow of the primorial inflation of n] - Antti Karttunen, Apr 29 2022

A317081 Number of integer partitions of n whose multiplicities cover an initial interval of positive integers.

Original entry on oeis.org

1, 1, 1, 2, 3, 5, 5, 9, 11, 16, 20, 30, 34, 50, 58, 79, 96, 129, 152, 203, 243, 307, 375, 474, 563, 707, 850, 1042, 1246, 1532, 1815, 2215, 2632, 3173, 3765, 4525, 5323, 6375, 7519, 8916, 10478, 12414, 14523, 17133, 20034, 23488, 27422, 32090, 37285, 43511, 50559
Offset: 0

Views

Author

Gus Wiseman, Jul 21 2018

Keywords

Comments

Also the number of integer partitions of n with distinct section-sums, where the k-th part of the section-sum partition is the sum of all (distinct) parts that appear at least k times. - Gus Wiseman, Apr 21 2025

Examples

			The a(1) = 1 through a(9) = 16 partitions:
 (1) (2) (3)  (4)   (5)   (6)   (7)    (8)    (9)
         (21) (31)  (32)  (42)  (43)   (53)   (54)
              (211) (41)  (51)  (52)   (62)   (63)
                    (221) (321) (61)   (71)   (72)
                    (311) (411) (322)  (332)  (81)
                                (331)  (422)  (432)
                                (421)  (431)  (441)
                                (511)  (521)  (522)
                                (3211) (611)  (531)
                                       (3221) (621)
                                       (4211) (711)
                                              (3321)
                                              (4221)
                                              (4311)
                                              (5211)
                                              (32211)
		

Crossrefs

The case with parts also covering an initial interval is A317088.
These partitions are ranked by A317090.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A047966 counts partitions with constant section-sums.
A048767 interchanges prime indices and prime multiplicities (Look-and-Say), see A048768.
A055932 lists numbers whose prime indices cover an initial interval.
A116540 counts normal set multipartitions.
A304442 counts partitions with equal run-sums, ranks A353833.
A381436 lists the section-sum partition of prime indices.
A381440 lists the Look-and-Say partition of prime indices.

Programs

  • Mathematica
    normalQ[m_]:=Union[m]==Range[Max[m]];
    Table[Length[Select[IntegerPartitions[n],normalQ[Length/@Split[#]]&]],{n,30}]
  • Python
    from sympy.utilities.iterables import partitions
    def A317081(n):
        if n == 0:
            return 1
        c = 0
        for d in partitions(n):
            s = set(d.values())
            if len(s) == max(s):
                c += 1
        return c # Chai Wah Wu, Jun 22 2020

A317090 Positive integers whose prime multiplicities span an initial interval of positive integers.

Original entry on oeis.org

2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 26, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 82, 83, 84, 85
Offset: 1

Views

Author

Gus Wiseman, Jul 21 2018

Keywords

Comments

The first term in this sequence but absent from A179983 is 180.
The numbers of terms that do not exceed 10^k, for k = 1, 2, ..., are 6, 78, 820, 8379, 84440, 846646, 8473868, 84763404, 847714834, 8477408261, ... . Apparently, the asymptotic density of this sequence exists and equals 0.8477... . - Amiram Eldar, Aug 04 2024

Crossrefs

Subsequences: A129912\{1}, A179983\{1}.
Subsequence of A337533.

Programs

  • Mathematica
    normalQ[m_]:=Union[m]==Range[Max[m]];
    Select[Range[2,100],normalQ[FactorInteger[#][[All,2]]]&]
  • PARI
    is(k) = {my(e = Set(factor(k)[,2])); k > 1 && vecmax(e) == #e;} \\ Amiram Eldar, Aug 04 2024

A242414 Numbers whose prime factorization viewed as a tuple of nonzero powers is palindromic.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 64, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 81, 82, 83, 85, 86, 87, 89, 90, 91, 93, 94, 95, 97, 100
Offset: 1

Views

Author

Antti Karttunen, May 30 2014

Keywords

Comments

The fixed points of permutation A069799.
Differs from its subsequence, A072774, Powers of squarefree numbers, for the first time at n=68, as here a(68) = 90 is included, as 90 = p_1^1 * p_2^2 * p_3^1 has a palindromic tuple of exponents, even although not all of them are identical.
Differs from its another subsequence, A236510, in that, although numbers like 42 = 2^1 * 3^1 * 5^0 * 7^1, with a non-palindromic exponent-tuple (1,1,0,1) are excluded from A236510, it is included in this sequence, because here only the nonzero exponents are considered, and (1,1,1) is a palindrome.
Differs from A085924 in that as that sequence is subtly base-dependent, it excludes 1024 (= 2^10), as then the only exponent present, 10, and thus also its concatenation, "10", is not a palindrome when viewed in decimal base. On the contrary, here a(691) = 1024.

Examples

			As 1 has an empty factorization, (), which also is a palindrome, 1 is present.
As 42 = 2 * 3 * 7 = p_1^1 * p_2^1 * p_4^1, and (1,1,1) is palindrome, 42 is present.
As 90 = 2 * 9 * 5 = p_1^1 * p_2^2 * p_3^1, and (1,2,1) is palindrome, 90 is present.
Any prime power (A000961) is present, as such numbers have a factorization p^e (e >= 1), and any singleton sequence (e) by itself forms a palindrome.
		

Crossrefs

Fixed points of A069799.
Complement: A242416.
A000961, A072774 and A236510 are subsequences.

Programs

  • Mathematica
    Select[Range[100], PalindromeQ[FactorInteger[#][[All, 2]]]&] (* Jean-François Alcover, Feb 09 2025 *)

A242419 Reverse both the exponents and the deltas of the indices of distinct primes present in the prime factorization of n.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 15, 11, 18, 13, 35, 10, 16, 17, 12, 19, 75, 21, 77, 23, 54, 25, 143, 27, 245, 29, 30, 31, 32, 55, 221, 14, 36, 37, 323, 91, 375, 41, 105, 43, 847, 50, 437, 47, 162, 49, 45, 187, 1859, 53, 24, 33, 1715, 247, 667, 59, 150, 61, 899, 147, 64, 65
Offset: 1

Views

Author

Antti Karttunen, May 17 2014

Keywords

Comments

This self-inverse permutation (involution) of natural numbers preserves both the total number of prime divisors and the (index of) largest prime factor of n, i.e. for all n it holds that A001222(a(n)) = A001222(n) and A006530(a(n)) = A006530(n) [equally: A061395(a(n)) = A061395(n)].
Considered as an operation on partitions encoded by the indices of primes in the prime factorization of n (as in table A112798), this implements a bijection which reverses the order of "steps" in Young (or Ferrers) diagram of a partition (but keeps the horizontal line segment of each step horizontal and the vertical line segment vertical). Please see the last example in the example section.
To understand the given recursive formula, it helps to see that in the above context (Young diagrams drawn with French notation), the sequences employed effect the following operations:
A001222: gives the height of whole diagram,
A051119: removes the bottommost step from the diagram,
A241919: gives the length of the horizontal line segment of the bottom step, i.e. its width,
A071178: gives the length of the vertical line segment of the bottom step, i.e. its height,
A242378(k,n): increases the width of whole Young diagram encoded by n by adding a rectangular area A001222(n) squares high and k squares wide to its left,
and finally, multiplying by A000040(a)^b adds a new topmost step whose width is a and height is b. Particularly, multiplying by (A000040(A241919(n))^A071178(n)) transfers the bottommost step to the top.

Examples

			For n = 10 = 2*5 = p_1^1 * p_3^1, we get p_(3-1)^1 * p_3^1 = 3 * 5 = 15, thus a(10) = 15.
For n = 20 = 2*2*5 = p_1^2 * p_3^1, we get p_(3-1)^1 * p_3^2 = 3^1 * 5^2 = 3*25 = 75, thus a(20) = 75.
For n = 84 = 2*2*3*7 = p_1^2 * p_2 * p_4, when we reverse the deltas of indices, and reverse also the order of exponents, we get p_(4-2) * p_(4-1) * p_4^2 = 3 * 5 * 7^2 = 735, thus a(84) = 735.
For n = 2200, we see that it encodes the partition (1,1,1,3,3,5) in A112798 as 2200 = p_1 * p_1 * p_1 * p_3 * p_3 * p_5 = 2^3 * 5^2 * 11. This in turn corresponds to the following Young diagram in French notation:
   _
  | |
  | |
  | |_ _
  |     |
  |     |_ _
  |_ _ _ _ _|
Reversing the order of "steps", so that each horizontal and vertical line segment centered around a "convex corner" moves as a whole, means that the first stair from the top (one unit wide and three units high) is moved to the last position, the second one (two units wide and two units high) stays in the middle, and the original bottom step (two units wide and one unit high) will be the new topmost step, thus we get the following Young diagram:
   _ _
  |   |_ _
  |       |
  |       |_
  |         |
  |         |
  |_ _ _ _ _|
which represents the partition (2,4,4,5,5,5), encoded in A112798 by p_2 * p_4^2 * p_5^3 = 3 * 7^2 * 11^3 = 195657, thus a(2200) = 195657.
		

Crossrefs

Fixed points: A242417.
{A000027, A122111, A153212, A242419} form a 4-group.
{A000027, A069799, A242415, A242419} form also a 4-group.

Formula

If n = p_a^e_a * p_b^e_b * ... * p_h^e_h * p_i^e_i * p_j^e_j * p_k^e_k, where p_a < ... < p_k are distinct primes (sorted into ascending order) in the prime factorization of n, and e_a .. e_k are their nonzero exponents, then a(n) = p_{k-j}^e_k * p_{k-i}^e_j * p_{k-h}^e_i * ... * p_{k-a}^e_b * p_k^e_a.
As a recurrence:
a(1) = 1, and for n>1, a(n) = (A000040(A241919(n))^A071178(n)) * A242378(A241919(n), a(A051119(n))).
By composing related permutations:
a(n) = A122111(A153212(n)) = A153212(A122111(n)).
a(n) = A069799(A242415(n)) = A242415(A069799(n)).
a(n) = A105119(A242420(n)).

A242415 Reverse the deltas of indices of distinct primes in the prime factorization of n.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 15, 11, 12, 13, 35, 10, 16, 17, 18, 19, 45, 21, 77, 23, 24, 25, 143, 27, 175, 29, 30, 31, 32, 55, 221, 14, 36, 37, 323, 91, 135, 41, 105, 43, 539, 20, 437, 47, 48, 49, 75, 187, 1573, 53, 54, 33, 875, 247, 667, 59, 60, 61, 899, 63, 64, 65
Offset: 1

Views

Author

Antti Karttunen, May 24 2014

Keywords

Comments

This self-inverse permutation (involution) of natural numbers preserves both the total number of prime divisors and the (index of) largest prime factor of n, i.e., for all n it holds that A001222(a(n)) = A001222(n) and A006530(a(n)) = A006530(n) [equally: A061395(a(n)) = A061395(n)]. It also preserves the exponent of the largest prime: A053585(a(n)) = A053585(n).
From the above it follows, that this fixes prime powers (A000961), among other numbers.
Considered as a function on partitions encoded by the indices of primes in the prime factorization of n (as in table A112798), this implements an operation which reverses the order of horizontal line segments of the "steps" in Young (or Ferrers) diagram of a partition, but keeps the order of vertical line segments intact. Please see the last example in the example section and compare also to the comments given in A242419.

Examples

			For n = 10 = 2*5 = p_1 * p_3, we get p_(3-1) * p_3 = 3 * 5 = 15, thus a(10) = 15.
For n = 20 = 2*2*5 = p_1^2 * p_3^1, we get p_(3-1)^2 * p_3^1 = 3^2 * 5 = 45, thus a(20) = 45.
For n = 84 = 2*2*3*7 = p_1^2 * p_2 * p_4, when we reverse the deltas of indices, but keep the exponents same, we get p_(4-2)^2 * p_(4-1) * p_4 = p_2^2 * p_3 * p_4 = 3^2 * 5 * 7 = 315, thus a(84) = 315.
For n = 2200, we see that it encodes the partition (1,1,1,3,3,5) in A112798 as 2200 = p_1 * p_1 * p_1 * p_3 * p_3 * p_5 = 2^3 * 5^2 * 11. This in turn corresponds to the following Young diagram in French notation:
   _
  | |
  | |
  | |_ _
  |     |
  |     |_ _
  |_ _ _ _ _|
Reversing the order of horizontal line segment lengths (1,2,2) to (2,2,1), but keeping the order of vertical line segment lengths as (3,2,1), we get a new Young diagram
   _ _
  |   |
  |   |
  |   |_ _
  |       |
  |       |_
  |_ _ _ _ _|
which represents the partition (2,2,2,4,4,5), encoded in A112798 by p_2^3 * p_4^2 * p_5^1 = 3^3 * 7^2 * 11 = 14553, thus a(2200) = 14553.
		

Crossrefs

Formula

If n = p_a^e_a * p_b^e_b * ... * p_h^e_h * p_i^e_i * p_j^e_j * p_k^e_k, where p_a < ... < p_k are distinct primes (sorted into ascending order) in the prime factorization of n, and e_a .. e_k are their nonzero exponents, then a(n) = p_{k-j}^e_a * p_{k-i}^e_b * p_{k-h}^e_c * ... * p_{k-a}^e_j * p_k^e_k.
As a recurrence: a(1) = 1, and for n>1, a(n) = (A000040(A241919(n))^A067029(n)) * A242378(A241919(n), a(A051119(A225891(n)))).
By composing/conjugating related permutations:
a(n) = A069799(A242419(n)) = A242419(A069799(n)).

A085079 Largest number with the prime signature of n using prime divisors of n.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 18, 13, 14, 15, 16, 17, 18, 19, 50, 21, 22, 23, 54, 25, 26, 27, 98, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 250, 41, 42, 43, 242, 75, 46, 47, 162, 49, 50, 51, 338, 53, 54, 55, 686, 57, 58, 59, 150, 61, 62, 147, 64, 65, 66, 67, 578, 69, 70
Offset: 1

Views

Author

Amarnath Murthy, Jul 01 2003

Keywords

Comments

Though a large number of initial terms match, this sequence is different from A069799. Example: a(1500) = a(2^2*3*5^3) = 5^3*3^2*2 = 2250, whereas A069799(1500) = 5^2*3*2^3 = 600.
The first term that is different from A069799 is a(18). - Ivan Neretin, Jul 29 2015

Examples

			20 = 2^2*5, hence a(20) = 5^2*2 = 50.
		

Crossrefs

Programs

  • Haskell
    import Data.List (sort)
    a085079 n = product $ zipWith (^) (a027748_row n) (sort $ a124010_row n)
    -- Reinhard Zumkeller, Apr 27 2013
  • Mathematica
    Table[Times @@ ((ar = Transpose[FactorInteger[n]])[[1]]^Sort[ar[[2]]]), {n, 70}] (* Ivan Neretin, Jul 29 2015*)
    fise[n_]:=Module[{fi=FactorInteger[n]},Times@@(fi[[All,1]]^Sort[ fi[ [All,2]]])]; Array[fise,100] (* Harvey P. Dale, Sep 18 2017 *)

Formula

a(n) >= n. - Michel Marcus, Jul 30 2015

Extensions

Corrected and extended by Ray Chandler, Aug 17 2003

A105119 Numbers obtained by rotating right the indices in the prime signature of n.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 18, 13, 14, 15, 16, 17, 12, 19, 50, 21, 22, 23, 54, 25, 26, 27, 98, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 250, 41, 42, 43, 242, 75, 46, 47, 162, 49, 20, 51, 338, 53, 24, 55, 686, 57, 58, 59, 90, 61, 62, 147, 64, 65, 66, 67, 578, 69, 70
Offset: 1

Views

Author

Yasutoshi Kohmoto, Apr 08 2005

Keywords

Comments

If n = p^a*q^b*r^c, then a(n) = p^c*q^a*r^b.
If n = p^a*q^b*r^c*s^d, then a(n) = p^d*q^a*r^b*s^c.
The sequence is a permutation of the positive integers. The first term which is different from A069799 is a(60).
Inverse permutation to A225891. The fixed points are A072774 (squarefree numbers and their powers). - Ivan Neretin, Jul 26 2015

Examples

			a(60) = a(2^2*3*5) = 2*3^2*5 = 90.
		

Crossrefs

Permutation of A000027(n). Cf. A069799.

Programs

  • Maple
    f:= proc(n) local F,j,m;
          F:= ifactors(n)[2];
          m:= nops(F);
          mul(F[i,1]^F[i-1,2],i=2..m)*F[1,1]^F[m,2] ;
    end proc:
    f(1):= 1:
    map(f, [$1..100]); # Robert Israel, Jul 26 2015
  • Mathematica
    Table[Times @@ ((ar = Transpose[FactorInteger[n]])[[1]]^RotateRight[ar[[2]]]), {n, 70}] (* Ivan Neretin, Jul 26 2015 *)
  • PARI
    a(n)=local(m,s);m=factor(n);s=matsize(m)[1];prod(i=2,s,m[i,1]^m[i-1,2])*m[1,1]^m[s,2] /* Ralf Stephan, Apr 05 2009 */

Extensions

Edited by Ralf Stephan, Apr 05 2009
a(1)=1 prepended by Ivan Neretin, Jul 26 2015

A153212 A permutation of the natural numbers: in the prime factorization of n, swap each prime's index difference (from the previous distinct prime that divides n) and the prime's power.

Original entry on oeis.org

1, 2, 4, 3, 8, 6, 16, 5, 9, 18, 32, 15, 64, 54, 12, 7, 128, 10, 256, 75, 36, 162, 512, 35, 27, 486, 25, 375, 1024, 30, 2048, 11, 108, 1458, 24, 21, 4096, 4374, 324, 245, 8192, 150, 16384, 1875, 45, 13122, 32768, 77, 81, 50, 972, 9375, 65536, 14, 72, 1715, 2916, 39366, 131072, 105, 262144, 118098, 225, 13
Offset: 1

Views

Author

Luchezar Belev (l_belev(AT)yahoo.com), Dec 20 2008

Keywords

Comments

In order for the "index difference" to make sense, we consider the factorization to be sorted with respect to the primes but not the powers to which they are raised; that is, first comes the smallest prime and each subsequent prime is larger than the previous disregarding their powers.
For every n it is true that a(a(n)) = n.
From Antti Karttunen, May 29 2014: (Start)
In other words, this is a self-inverse permutation (involution) of natural numbers.
This permutation maps primes (A000040) to the powers of two larger than one (A000079(n>=1)) and vice versa.
The term a(1) = 1 was added on the grounds that as 1 has an empty prime factorization, there is nothing to swap, thus it stays same. It is also needed as a base case for the given recurrence.
Considered as a function on partitions encoded by the indices of primes in the prime factorization of n (as in table A112798), this implements an operation which exchanges the horizontal and vertical line segment of each "step" in Young (or Ferrers) diagram of a partition. Please see the last example in the example section.
(End)

Examples

			For n = 10 we have 10 = 2^1 * 5^1 = p(1)^1 * p(3)^1 then a(10) = p(1)^1 * p(2)^2 = 2^1 * 3^2 = 18.
For n = 18 we have 18 = 2^1 * 3^2 = p(1)^1 * p(2)^2 then a(18) = p(1)^1 * p(3)^1 = 2^1 * 5^1 = 10.
For n = 19 we have 19 = 19^1 = p(8)^1 then a(19) = p(1)^8 = 2^8 = 256.
For n = 2200, we see that it encodes the partition (1,1,1,3,3,5) in A112798 as 2200 = p_1 * p_1 * p_1 * p_3 * p_3 * p_5 = 2^3 * 5^2 * 11. This in turn corresponds to the following Young diagram in French notation:
   _
  | |
  | |
  | |_ _
  |     |
  |     |_ _
  |_ _ _ _ _|
Exchanging the order of the horizontal and vertical line segment of each "step", results the following Young diagram:
   _ _ _
  |     |_ _
  |         |
  |         |_
  |           |
  |_ _ _ _ _ _|
which represents the partition (3,5,5,6,6), encoded in A112798 by p_3 * p_5^2 * p_6^2 = 5 * 11^2 * 13^2 = 102245, thus a(2200) = 102245.
		

Crossrefs

Fixed points: A242421.
{A000027, A122111, A153212, A242419} form a 4-group.

Programs

  • PARI
    a(n) = {my(f = factor(n)); my(g = f); for (i=1, #f~, if (i==1, g[i,1] = prime(f[i,2]), g[i,1] = prime(f[i,2]+ primepi(g[i-1,1]))); if (i==1, g[i,2] = primepi(f[i,1]), g[i,2] = primepi(f[i,1]) - primepi(f[i-1,1]));); factorback(g);} \\ Michel Marcus, Dec 16 2014

Formula

Denote the i-th prime with p(i): p(1)=2, p(2)=3, p(3)=5, p(4)=7, etc. Let n = p(a1)^b1 * p(a2)^b2 * ... * p(ak)^bk is the factorization of n where p(i)^j is the i-th prime raised to power j. As mentioned above, we assume that the primes are sorted, i.e., a1 < a2 < a3 ... < ak. Then a(n) = p(c1)^d1 * p(c2)^d2 * ... * p(ck)^dk where c1 = b1 and c(i) = b(i) + c(i-1) for i > 1 d1 = a1 and d(i) = a(i) - a(i-1) for i > 1.
From Antti Karttunen, May 16 2014: (Start)
a(1) = 1 and for n>1, let r = a(A051119(n)). Then a(n) = r * (A000040(A061395(r)+A071178(n)) ^ A241919(n)).
a(n) = A122111(A242419(n)) = A242419(A122111(n)).
(End)

Extensions

Term a(1)=1 prepended, and also more terms computed by Antti Karttunen, May 16 2014

A242416 Numbers whose prime factorization viewed as a tuple of nonzero powers is not palindromic.

Original entry on oeis.org

12, 18, 20, 24, 28, 40, 44, 45, 48, 50, 52, 54, 56, 60, 63, 68, 72, 75, 76, 80, 84, 88, 92, 96, 98, 99, 104, 108, 112, 116, 117, 120, 124, 132, 135, 136, 140, 144, 147, 148, 150, 152, 153, 156, 160, 162, 164, 168, 171, 172, 175, 176, 180, 184, 188, 189, 192, 200
Offset: 1

Views

Author

Antti Karttunen, May 29 2014

Keywords

Comments

These are terms that appear in 2-cycles of permutation A069799.
Complement of A242414.

Examples

			12 = p_1^2 * p_2^1 is present, as (2,1) is not a palindrome.
		

Crossrefs

Complement: A242414.
A subsequence of A059404, from which this differs for the first at n=23, as 90 = A059404(23) is not member of this sequence, as the exponents in the prime factorization of 90 = 2^1 * 3^2 * 5^1 form a palindrome, even though 90 is not a power of a squarefree number.
Cf. A069799.

Programs

  • Maple
    q:= n-> (l-> is(n<>mul(l[i, 1]^l[-i, 2], i=1..nops(l))))(sort(ifactors(n)[2])):
    select(q, [$1..200])[];  # Alois P. Heinz, Feb 04 2022
  • Mathematica
    Select[Range[200], !PalindromeQ[FactorInteger[#][[All, 2]]]&] (* Jean-François Alcover, Feb 09 2025 *)
Showing 1-10 of 17 results. Next