A242414 -id:A242414 - cp's OEIS frontend

A025065 Number of palindromic partitions of n.

1, 1, 2, 2, 4, 4, 7, 7, 12, 12, 19, 19, 30, 30, 45, 45, 67, 67, 97, 97, 139, 139, 195, 195, 272, 272, 373, 373, 508, 508, 684, 684, 915, 915, 1212, 1212, 1597, 1597, 2087, 2087, 2714, 2714, 3506, 3506, 4508, 4508, 5763, 5763, 7338, 7338, 9296, 9296, 11732, 11732, 14742, 14742, 18460, 18460, 23025, 23025, 28629, 28629

Offset: 0

Clark Kimberling

nonn

That is, the number of partitions of n into parts which can be listed in palindromic order.
Alternatively, number of partitions of n into parts from the set {1,2,4,6,8,10,12,...}. - T. D. Noe, Aug 05 2005
Also, partial sums of A035363.
Also number of partitions of n with at most one part occurring an odd number of times. - Reinhard Zumkeller, Dec 18 2013
The first Mathematica program computes terms of A025065; the second computes the k palindromic partitions of user-chosen n. - Clark Kimberling, Jan 20 2014
a(n) is the number of partitions p of n+1 such that 2*max(p) > n+1. - Clark Kimberling, Apr 20 2014.
From Gus Wiseman, Nov 28 2018: (Start)
Also the number of integer partitions of n + 2 that are the vertex-degrees of some hypertree. For example, the a(6) = 7 partitions of 8 that are the vertex-degrees of some hypertree, together with a realizing hypertree are:
     (41111): {{1,2},{1,3},{1,4},{1,5}}
     (32111): {{1,2},{1,3},{1,4},{2,5}}
     (22211): {{1,2},{1,3},{2,4},{3,5}}
    (311111): {{1,2},{1,3},{1,4,5,6}}
    (221111): {{1,2},{1,3},{2,4,5,6}}
   (2111111): {{1,2},{1,3,4,5,6,7}}
  (11111111): {{1,2,3,4,5,6,7,8}}
(End)
Conjecture: a(n) is the length of maximal initial segment of A308355(n-1) that is identical to row n of A128628, for n >= 2. - Clark Kimberling, May 24 2019
From Gus Wiseman, May 21 2021: (Start)
The Heinz numbers of palindromic partitions are given by A265640. The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.
Also the number of integer partitions of n with a part greater than or equal to n/2. This is equivalent to Clark Kimberling's final comment above. The Heinz numbers of these partitions are given by A344414. For example, the a(1) = 1 through a(8) = 12 partitions are:
  (1)  (2)   (3)   (4)    (5)    (6)     (7)     (8)
       (11)  (21)  (22)   (32)   (33)    (43)    (44)
                   (31)   (41)   (42)    (52)    (53)
                   (211)  (311)  (51)    (61)    (62)
                                 (321)   (421)   (71)
                                 (411)   (511)   (422)
                                 (3111)  (4111)  (431)
                                                 (521)
                                                 (611)
                                                 (4211)
                                                 (5111)
                                                 (41111)
Also the number of integer partitions of n with at least n/2 parts. The Heinz numbers of these partitions are given by A344296. For example, the a(1) = 1 through a(8) = 12 partitions are:
  (1)  (2)   (21)   (22)    (221)    (222)     (2221)     (2222)
       (11)  (111)  (31)    (311)    (321)     (3211)     (3221)
                    (211)   (2111)   (411)     (4111)     (3311)
                    (1111)  (11111)  (2211)    (22111)    (4211)
                                     (3111)    (31111)    (5111)
                                     (21111)   (211111)   (22211)
                                     (111111)  (1111111)  (32111)
                                                          (41111)
                                                          (221111)
                                                          (311111)
                                                          (2111111)
                                                          (11111111)
(End)

			The partitions for the first few values of n are as follows:
n: partitions .......................... number
1: 1 ................................... 1
2: 2 11 ................................ 2
3: 3 111 ............................... 2
4: 4 22 121 1111 ....................... 4
5: 5 131 212 11111 ..................... 4
6: 6 141 33 222 1221 11211 111111 ...... 7
7: 7 151 313 11311 232 21112 1111111 ... 7
From _Reinhard Zumkeller_, Jan 23 2010: (Start)
Partitions into 1,2,4,6,... for the first values of n:
1: 1 ....................................... 1
2: 2 11 .................................... 2
3: 21 111 .................................. 2
4: 4 22 211 1111 ........................... 4
5: 41 221 2111 11111 ....................... 4
6: 6 42 4211 222 2211 21111 111111.......... 7
7: 61 421 42111 2221 22111 211111 1111111 .. 7. (End)

David A. Corneth, Table of n, a(n) for n = 0..10000 (first 101 terms from Reinhard Zumkeller that were corrected by _Georg Fischer_, Jan 20 2019)

Cf. A172033, A004277. - Reinhard Zumkeller, Jan 23 2010
Cf. A004526, A030019, A056503, A147878, A320921, A322136.
The bisections are both A000070.
The ordered version (palindromic compositions) is A016116.
The complement is counted by A233771 and A210249.
The case of palindromic prime signature is A242414.
Palindromic partitions are ranked by A265640, with complement A229153.
The case of palindromic plane trees is A319436.
The multiplicative version (palindromic factorizations) is A344417.
A000569 counts graphical partitions.
A027187 counts partitions of even length, ranked by A028260.
A035363 counts partitions into even parts, ranked by A066207.
A058696 counts partitions of even numbers, ranked by A300061.
A110618 counts partitions with length <= half sum, ranked by A344291.
Cf. A000041, A067538, A143773, A209816, A338914, A338915, A340387, A344296, A344414, A344415, A344416.

a025065 = p (1:[2,4..]) where
   p [] _ = 0
   p _  0 = 1
   p ks'@(k:ks) m | m < k     = 0
                  | otherwise = p ks' (m - k) + p ks m
-- Reinhard Zumkeller, Aug 12 2011

import Data.List (group)
a025065 = length . filter (<= 1) .
                   map (sum . map ((`mod` 2) . length) . group) . ps 1
   where ps x 0 = [[]]
         ps x y = [t:ts | t <- [x..y], ts <- ps t (y - t)]
-- Reinhard Zumkeller, Dec 18 2013

Map[Length[Select[IntegerPartitions[#], Count[OddQ[Transpose[Tally[#]][[2]]], True] <= 1 &]] &, Range[40]] (* Peter J. C. Moses, Jan 20 2014 *)
n = 8; Select[IntegerPartitions[n], Count[OddQ[Transpose[Tally[#]][[2]]], True] <= 1 &] (* Peter J. C. Moses, Jan 20 2014 *)
CoefficientList[Series[1/((1 - x) Product[1 - x^(2 n), {n, 1, 50}]), {x, 0, 60}], x] (* Clark Kimberling, Mar 14 2014 *)

N=66; q='q+O('q^N); Vec( 1/((1-q)*eta(q^2)) ) \\ Joerg Arndt, Mar 11 2014

a(n) = A000070(A004526(n)). - Reinhard Zumkeller, Jan 23 2010
G.f.: 1/((1-q)*prod(n>=1, 1-q^(2*n))). [Joerg Arndt, Mar 11 2014]
a(2*k+2) = a(2*k) + A000041(k + 1). - David A. Corneth, May 29 2021
a(n) ~ exp(Pi*sqrt(n/3)) / (2*Pi*sqrt(n)). - Vaclav Kotesovec, Nov 16 2021

Edited by N. J. A. Sloane, Dec 29 2007

Prepended a(0)=1, added more terms, Joerg Arndt, Mar 11 2014

A072774 Powers of squarefree numbers.

Original entry on oeis.org

Offset: 1

Reinhard Zumkeller, Jul 10 2002

nonn

Essentially the same as A062770. - R. J. Mathar, Sep 25 2008
Numbers m such that in canonical prime factorization all prime exponents are identical: A124010(m,k) = A124010(m,1) for k = 2..A000005(m). - Reinhard Zumkeller, Apr 06 2014
Heinz numbers of uniform partitions. An integer partition is uniform if all parts appear with the same multiplicity. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 16 2018

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Complement of A059404.
Cf. A072775, A072776, A072777 (subsequence), A005117, A072778, A124010, A329332 (tabular arrangement), A384667 (characteristic function).
A subsequence of A242414.
Cf. A000005, A000009, A000837, A007916, A013661, A047966, A052409, A052410, A072774, A078374, A289023, A289509, A300486, A302491, A302796, A302979.

import Data.Map (empty, findMin, deleteMin, insert)
import qualified Data.Map.Lazy as Map (null)
a072774 n = a072774_list !! (n-1)
(a072774_list, a072775_list, a072776_list) = unzip3 $
   (1, 1, 1) : f (tail a005117_list) empty where
   f vs'@(v:vs) m
    | Map.null m || xx > v = (v, v, 1) :
                             f vs (insert (v^2) (v, 2) m)
    | otherwise = (xx, bx, ex) :
                  f vs' (insert (bx*xx) (bx, ex+1) $ deleteMin m)
    where (xx, (bx, ex)) = findMin m
-- Reinhard Zumkeller, Apr 06 2014

isA := n -> n=1 or is(1 = nops({seq(p[2], p in ifactors(n)[2])})):
select(isA, [seq(1..97)]);  # Peter Luschny, Jun 10 2025

Select[Range[100], Length[Union[FactorInteger[#][[All, 2]]]] == 1 &] (* Geoffrey Critzer, Mar 30 2015 *)

is(n)=ispower(n,,&n); issquarefree(n) \\ Charles R Greathouse IV, Oct 16 2015

from math import isqrt
from sympy import mobius, integer_nthroot
def A072774(n):
    def g(x): return int(sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1)))-1
    def f(x): return n-2+x-sum(g(integer_nthroot(x,k)[0]) for k in range(1,x.bit_length()))
    kmin, kmax = 1,2
    while f(kmax) >= kmax:
        kmax <<= 1
    while True:
        kmid = kmax+kmin>>1
        if f(kmid) < kmid:
            kmax = kmid
        else:
            kmin = kmid
        if kmax-kmin <= 1:
            break
    return kmax # Chai Wah Wu, Aug 19 2024

a(n) = A072775(n)^A072776(n).
Sum_{n>=1} 1/a(n)^s = 1 + Sum_{k>=1} (zeta(k*s)/zeta(2*k*s)-1) for s > 1. - Amiram Eldar, Mar 20 2025
a(n)/n ~ Pi^2/6 (A013661). - Friedjof Tellkamp, Jun 09 2025

A317089 Numbers whose prime factors span an initial interval of prime numbers and whose prime multiplicities span an initial interval of positive integers.

Original entry on oeis.org

2, 6, 12, 18, 30, 60, 90, 150, 180, 210, 300, 360, 420, 450, 540, 600, 630, 1050, 1260, 1350, 1470, 1500, 2100, 2250, 2310, 2520, 2940, 3150, 3780, 4200, 4410, 4620, 5880, 6300, 6930, 7350, 8820, 9450, 10500, 11550, 12600, 13230, 13860, 14700, 15750, 16170

Offset: 1

Gus Wiseman, Jul 21 2018

nonn

			The sequence of rows of A296150 indexed by the terms of this sequence begins: (1), (21), (211), (221), (321), (3211), (3221), (3321), (32211), (4321), (33211), (322111), (43211).

Andrew Howroyd, Table of n, a(n) for n = 1..300

Cf. A001222, A007916, A055932, A056239, A124010, A133808, A242414, A296150.

Cf. A317087, A317088, A317090, A317091, A317092.

normalQ[m_]:=Union[m]==Range[Max[m]];
Select[Range[10000],And[normalQ[PrimePi/@FactorInteger[#][[All,1]]],normalQ[FactorInteger[#][[All,2]]]]&]

ok(n)={my(f=factor(n), p=f[,1], e=vecsort(f[,2],,8)); n > 1 && #p==primepi(p[#p]) && #e==e[#e]} \\ Andrew Howroyd, Aug 26 2018

A069799 The number obtained by reversing the sequence of nonzero exponents in the prime factorization of n with respect to distinct primes present, as ordered by their indices.

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, 150, 61, 62, 147, 64, 65

Offset: 1

Amarnath Murthy, Apr 13 2002

nonn

Equivalent description nearer to the old name: a(n) is a number obtained by reversing the indices of the primes present in the prime factorization of n, from the smallest to the largest, while keeping the nonzero exponents of those same primes at their old positions.
This self-inverse permutation of natural numbers fixes the numbers in whose prime factorization the sequence of nonzero exponents form a palindrome: A242414.
Integers which are changed are A242416.
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 vertical line segments of the "steps" in Young (or Ferrers) diagram of a partition, but keeps the order of horizontal line segments intact. Please see the last example in the example section.

			a(24) = 54 as 24 = p_1^3 * p_2^1 = 2^3 * 3^1 and 54 = p_1^1 * p_2^3 = 2 * 3^3.
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 vertical line segment lengths (3,2,1)  to (1,2,3), but keeping the order of horizontal line segment lengths as (1,2,2), we get a new Young diagram
   _
  | |_ _
  |     |
  |     |_ _
  |         |
  |         |
  |_ _ _ _ _|
which represents the partition (1,3,3,5,5,5), encoded in A112798 by p_1 * p_3^2 * p_5^3 = 2 * 5^2 * 11^3 = 66550, thus a(2200) = 66550.

Alois P. Heinz (first 1000 terms) & Antti Karttunen, Table of n, a(n) for n = 1..8192
Wikipedia, Young diagram
Index entries for sequences that are permutations of the natural numbers

A242414 gives the fixed points and A242416 is their complement.
Cf. A112798, A059404, A122111, A153212, A137502, A241916, A242418.
{A000027, A069799, A242415, A242419} form a 4-group.
The set of permutations {A069799, A105119, A225891} generates an infinite dihedral group.

a069799 n = product $
            zipWith (^) (a027748_row n) (reverse $ a124010_row n)
-- Reinhard Zumkeller, Apr 27 2013
(MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library)
(require 'factor)
(define (A069799 n) (let ((pf (ifactor n))) (apply * (map expt (uniq pf) (reverse (multiplicities pf))))))
(define (ifactor n) (cond ((< n 2) (list)) (else (sort (factor n) <))))
(define (uniq lista) (let loop ((lista lista) (z (list))) (cond ((null? lista) (reverse! z)) ((and (pair? z) (equal? (car z) (car lista))) (loop (cdr lista) z)) (else (loop (cdr lista) (cons (car lista) z))))))
(define (multiplicities lista) (let loop ((mults (list)) (lista lista) (prev #f)) (cond ((not (pair? lista)) (reverse! mults)) ((equal? (car lista) prev) (set-car! mults (+ 1 (car mults))) (loop mults (cdr lista) prev)) (else (loop (cons 1 mults) (cdr lista) (car lista))))))
;; Antti Karttunen, May 24 2014

A069799 := proc(n) local e,j; e := ifactors(n)[2]:
mul (e[j][1]^e[nops(e)-j+1][2], j=1..nops(e)) end:
seq (A069799(i), i=1..40);
# Peter Luschny, Jan 17 2011

f[n_] := Block[{a = Transpose[ FactorInteger[n]], m = n}, If[ Length[a] == 2, Apply[ Times, a[[1]]^Reverse[a[[2]] ]], m]]; Table[ f[n], {n, 1, 65}]

a(n) = {my(f = factor(n)); my(g = f); my(nbf = #f~); for (i=1, nbf, g[i, 1] = f[nbf-i+1, 1];); factorback(g);} \\ Michel Marcus, Jul 02 2015

If n = p_a^e_a * p_b^e_b * ... * p_j^e_j * p_k^e_k, where p_a < ... < p_k are distinct primes of the prime factorization of n (sorted into ascending order), and e_a, ..., e_k are their nonzero exponents, then a(n) = p_a^e_k * p_b^e_j * ... * p_j^e_b * p_k^e_a.
a(n) = product(A027748(o(n)+1-k)^A124010(k): k=1..o(n)) = product(A027748(k)^A124010(o(n)+1-k): k=1..o(n)), where o(n) = A001221(n). - Reinhard Zumkeller, Apr 27 2013
From Antti Karttunen, Jun 01 2014: (Start)
Can be obtained also by composing/conjugating related permutations:
a(n) = A242415(A242419(n)) = A242419(A242415(n)).
a(n) = A122111(A242415(A122111(n))) = A153212(A242415(A153212(n))).
(End)

Edited, corrected and extended by Robert G. Wilson v and Vladeta Jovovic, Apr 15 2002
Definition corrected by Reinhard Zumkeller, Apr 27 2013
Definition again reworded, Comments section edited and Young diagram examples added by Antti Karttunen, May 30 2014

A329395 Numbers whose binary expansion without the most significant (first) digit has Lyndon and co-Lyndon factorizations of equal lengths.

Original entry on oeis.org

1, 2, 3, 4, 7, 8, 10, 13, 15, 16, 22, 25, 31, 32, 36, 42, 46, 49, 53, 59, 63, 64, 76, 82, 94, 97, 109, 115, 127, 128, 136, 148, 156, 162, 166, 169, 170, 172, 181, 182, 190, 193, 201, 202, 211, 213, 214, 217, 221, 227, 235, 247, 255, 256, 280, 292, 306, 308

Offset: 1

Gus Wiseman, Nov 13 2019

nonn

We define the Lyndon product of two or more finite sequences to be the lexicographically maximal sequence obtainable by shuffling the sequences together. For example, the Lyndon product of (231) with (213) is (232131), the product of (221) with (213) is (222131), and the product of (122) with (2121) is (2122121). A Lyndon word is a finite sequence that is prime with respect to the Lyndon product. Equivalently, a Lyndon word is a finite sequence that is lexicographically strictly less than all of its cyclic rotations. Every finite sequence has a unique (orderless) factorization into Lyndon words, and if these factors are arranged in lexicographically decreasing order, their concatenation is equal to their Lyndon product. For example, (1001) has sorted Lyndon factorization (001)(1).
Similarly, the co-Lyndon product is the lexicographically minimal sequence obtainable by shuffling the sequences together, and a co-Lyndon word is a finite sequence that is prime with respect to the co-Lyndon product, or, equivalently, a finite sequence that is lexicographically strictly greater than all of its cyclic rotations. For example, (1001) has sorted co-Lyndon factorization (1)(100).
Conjecture: also numbers k such that the k-th composition in standard order (A066099) is a palindrome, cf. A025065, A242414, A317085, A317086, A317087, A335373. - Gus Wiseman, Jun 06 2020

			The sequence of terms together with their trimmed binary expansions and their co-Lyndon and Lyndon factorizations begins:
   1:      () =               0 = 0
   2:     (0) =             (0) = (0)
   3:     (1) =             (1) = (1)
   4:    (00) =          (0)(0) = (0)(0)
   7:    (11) =          (1)(1) = (1)(1)
   8:   (000) =       (0)(0)(0) = (0)(0)(0)
  10:   (010) =         (0)(10) = (01)(0)
  13:   (101) =         (10)(1) = (1)(01)
  15:   (111) =       (1)(1)(1) = (1)(1)(1)
  16:  (0000) =    (0)(0)(0)(0) = (0)(0)(0)(0)
  22:  (0110) =        (0)(110) = (011)(0)
  25:  (1001) =        (100)(1) = (1)(001)
  31:  (1111) =    (1)(1)(1)(1) = (1)(1)(1)(1)
  32: (00000) = (0)(0)(0)(0)(0) = (0)(0)(0)(0)(0)
  36: (00100) =     (0)(0)(100) = (001)(0)(0)
  42: (01010) =     (0)(10)(10) = (01)(01)(0)
  46: (01110) =       (0)(1110) = (0111)(0)
  49: (10001) =       (1000)(1) = (1)(0001)
  53: (10101) =     (10)(10)(1) = (1)(01)(01)
  59: (11011) =     (110)(1)(1) = (1)(1)(011)
  63: (11111) = (1)(1)(1)(1)(1) = (1)(1)(1)(1)(1)

Lyndon and co-Lyndon compositions are (both) counted by A059966.
Numbers whose reversed binary expansion is Lyndon are A328596.
Numbers whose binary expansion is co-Lyndon are A275692.
Cf. A001037, A060223, A102659, A211100, A329131, A329312, A329313, A329318, A329326, A329394, A329398.
Cf. A329358, A329399, A329400, A329401.

lynQ[q_]:=Array[Union[{q, RotateRight[q, #]}]=={q, RotateRight[q, #]}&, Length[q]-1, 1, And];
lynfac[q_]:=If[Length[q]==0, {}, Function[i, Prepend[lynfac[Drop[q, i]], Take[q, i]]][Last[Select[Range[Length[q]], lynQ[Take[q, #]]&]]]];
colynQ[q_]:=Array[Union[{RotateRight[q, #], q}]=={RotateRight[q, #], q}&, Length[q]-1, 1, And];
colynfac[q_]:=If[Length[q]==0, {}, Function[i, Prepend[colynfac[Drop[q, i]], Take[q, i]]]@Last[Select[Range[Length[q]], colynQ[Take[q, #]]&]]];
Select[Range[100],Length[lynfac[Rest[IntegerDigits[#,2]]]]==Length[colynfac[Rest[IntegerDigits[#,2]]]]&]

A317087 Numbers whose prime factors span an initial interval of prime numbers and whose sequence of prime multiplicities is a palindrome.

Original entry on oeis.org

1, 2, 4, 6, 8, 16, 30, 32, 36, 64, 90, 128, 210, 216, 256, 270, 300, 512, 810, 900, 1024, 1296, 2048, 2310, 2430, 2700, 2940, 3000, 3150, 4096, 7290, 7776, 8100, 8192, 9000, 11550, 16384, 21870, 24300, 27000, 30000, 30030, 32768, 41160, 44100, 46656, 47250, 48510

Offset: 1

Gus Wiseman, Jul 21 2018

nonn

3^m*10^k for k, m > 0 are terms of this sequence. - Chai Wah Wu, Jun 23 2020

			The sequence of rows of A296150 indexed by the terms of this sequence begins: (1), (11), (21), (111), (1111), (321), (11111), (2211), (111111), (3221), (1111111), (4321), (222111), (11111111), (32221), (33211), (111111111), (322221), (332211).

Giovanni Resta, Table of n, a(n) for n = 1..10000
Wikipedia, Palindrome

Cf. A025065, A055932, A124010, A133808, A242414, A296150, A317085, A317086.

nrmpalQ[n_]:=With[{f=If[n==1,{},FactorInteger[n]]}, And[PrimePi/@ Sort[First/@f] == Range[ Length[f]], Reverse[Last/@f] == Last/@f]]; Select[Range[100],nrmpalQ]
upto = 10^20; pL[n_] := Block[{p = Prime@Range@n, h = Ceiling[n/2]}, Take[p, h] Reverse@ If[n == 2 h, Take[p, -h], Prepend[ Take[p, 1-h], 1]]]; ric[v_, p_] := If[p == {}, AppendTo[L, v], Block[{w = v}, While[w <= upto, ric[w, Rest@ p]; w *= First@ p]]]; np = 1; L = {1}; While[(b = Times @@ Prime[Range@ np]) <= upto, ric[b, pL[np++]]]; Sort[L] (* Giovanni Resta, Jun 23 2020 *)

from sympy import factorint, primepi
A317087_list = [1]
for n in range(1,10**5):
    d = factorint(n)
    k, l = sorted(d.keys()), len(d)
    if l > 0 and l == primepi(max(d)):
        for i in range(l//2):
            if d[k[i]] != d[k[l-i-1]]:
                break
        else:
            A317087_list.append(n) # Chai Wah Wu, Jun 23 2020

A317086 Number of normal integer partitions of n whose sequence of multiplicities is a palindrome.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 3, 1, 2, 2, 4, 1, 4, 1, 4, 5, 4, 1, 7, 1, 8, 6, 6, 1, 10, 5, 7, 8, 11, 1, 20, 1, 9, 12, 9, 13, 25, 1, 10, 17, 21, 1, 37, 1, 21, 36, 12, 1, 44, 16, 23, 30, 33, 1, 53, 17, 55, 38, 15, 1, 103, 1, 16, 95, 51, 28, 69, 1, 73, 57, 82

Offset: 0

Gus Wiseman, Jul 21 2018

A partition is normal if its parts span an initial interval of positive integers.
a(n) = 1 if and only if n = 0, 1, 2, 4 or a prime > 3. - Chai Wah Wu, Jun 22 2020
From David A. Corneth, Jul 08 2020: (Start)
Let [f_1, f_2, ,..., f_i, ..., f_m] be the multiplicities of parts i in a partition of Sum_{i=1..m} (f_i * i). Then, as the sequence of multiplicities is a palindrome, we have f_1 = f_m, ..., f_i = f_(m+1-i). So the sum is f_1 * (1 + m) + f_2 * (2 + m-1) + ... + f_(floor(m/2)) * m/2 (the last term depending on the parity of m.). This way it becomes a list of Diophantine equations for which we look for the number of solutions.
For example, for m = 4 we look for solutions to the Diophantine equation 5 * (c + d) = n where c, d are positive integers >= 1. A similar technique is used in A254524. (End)

			The a(20) = 8 partitions:
  (44432111), (44332211), (43332221),
  (3333221111), (3332222111), (3322222211), (3222222221),
  (11111111111111111111).

David A. Corneth, Table of n, a(n) for n = 0..9999 (first 215 terms from Chai Wah Wu)
Wikipedia, Palindrome

Cf. A000009, A000041, A000837, A025065, A055932, A242414, A254524, A317085, A317087.

Table[Length[Select[IntegerPartitions[n],And[Union[#]==Range[First[#]],Length/@Split[#]==Reverse[Length/@Split[#]]]&]],{n,30}]

from sympy.utilities.iterables import partitions
from sympy import integer_nthroot, isprime
def A317086(n):
    if n > 3 and isprime(n):
        return 1
    else:
        c = 1
        for d in partitions(n,k=integer_nthroot(2*n,2)[0],m=n*2//3):
            l = len(d)
            if l > 0:
                k = max(d)
                if l == k:
                    for i in range(k//2):
                        if d[i+1] != d[k-i]:
                            break
                    else:
                        c += 1
        return c # Chai Wah Wu, Jun 22 2020

A317085 Number of integer partitions of n whose sequence of multiplicities is a palindrome.

Original entry on oeis.org

1, 1, 2, 3, 4, 4, 8, 6, 11, 12, 18, 16, 31, 25, 40, 47, 60, 58, 92, 85, 125, 135, 165, 173, 248, 246, 310, 351, 435, 450, 602, 608, 766, 846, 997, 1098, 1382, 1421, 1713, 1912, 2272, 2413, 2958, 3118, 3732, 4135, 4718, 5127, 6170, 6550, 7638, 8396, 9667, 10433

Offset: 0

Gus Wiseman, Jul 21 2018

nonn

			The a(10) = 18 partitions:
(ten),
(91), (82), (73), (64), (55),
(721), (631), (541), (532),
(5221), (4411), (4321), (3322),
(33211), (32221), (22222),
(1111111111).

Chai Wah Wu, Table of n, a(n) for n = 0..166
Wikipedia, Palindrome

Cf. A000041, A000837, A025065, A124010, A242414, A317086, A317087.

Table[Length[Select[IntegerPartitions[n],Length/@Split[#]==Reverse[Length/@Split[#]]&]],{n,30}]

from sympy.utilities.iterables import partitions
def A317085(n):
    c = 1
    for d in partitions(n,m=n*2//3):
        l = len(d)
        if l > 0:
            k = sorted(d.keys())
            for i in range(l//2):
                if d[k[i]] != d[k[l-i-1]]:
                    break
            else:
                c += 1
    return c # Chai Wah Wu, Jun 22 2020

A360247 Numbers for which the prime indices have the same mean as the distinct prime indices.

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, 101, 102, 103, 105, 106, 107, 109, 110, 111, 113, 114, 115, 118, 119, 121, 122, 123, 125, 127, 128, 129, 130

Offset: 1

Gus Wiseman, Feb 07 2023

nonn

First differs from A072774 in having 90.
First differs from A242414 in lacking 126.
Includes all squarefree numbers and perfect powers.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

			The prime indices of 900 are {3,3,2,2,1,1} with mean 2, and the distinct prime indices are {1,2,3} also with mean 2, so 900 is in the sequence.

Signature instead of parts: A324570, counted by A114638.
Signature instead of distinct parts: A359903, counted by A360068.
These partitions are counted by A360243.
The complement is A360246, counted by A360242.
For median instead of mean the complement is A360248, counted by A360244.
For median instead of mean we have A360249, counted by A360245.
For greater instead of equal mean we have A360252, counted by A360250.
For lesser instead of equal mean we have A360253, counted by A360251.
A008284 counts partitions by number of parts, distinct A116608.
A058398 counts partitions by mean, also A327482.
A088529/A088530 gives mean of prime signature (A124010).
A112798 lists prime indices, length A001222, sum A056239.
A316413 = numbers whose prime indices have integer mean, distinct A326621.
A326567/A326568 gives mean of prime indices.
A326619/A326620 gives mean of distinct prime indices.
Cf. A000975, A051293, A067340, A067538, A360005, A360241.

isA360247 := proc(n)
    local ifs,pidx,pe,meanAll,meanDist ;
    if n = 1 then
        return true ;
    end if ;
    ifs := ifactors(n)[2] ;
    # list of prime indices with multiplicity
    pidx := [] ;
    for pe in ifs do
        [numtheory[pi](op(1,pe)),op(2,pe)] ;
        pidx := [op(pidx),%] ;
    end do:
    meanAll := add(op(1,pe)*op(2,pe),pe=pidx) / add(op(2,pe),pe=pidx) ;
    meanDist := add(op(1,pe),pe=pidx) / nops(pidx) ;
    if meanAll = meanDist then
        true;
    else
        false;
    end if;
end proc:
for n from 1 to 130 do
    if isA360247(n) then
        printf("%d,",n) ;
    end if;
end do: # R. J. Mathar, May 22 2023

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],Mean[prix[#]]==Mean[Union[prix[#]]]&]

A332643 Neither the unsorted prime signature of a(n) nor the negated unsorted prime signature of a(n) is unimodal.

Original entry on oeis.org

2100, 3300, 3900, 4200, 4410, 5100, 5700, 6468, 6600, 6900, 7644, 7800, 8400, 8700, 9300, 9996, 10200, 10500, 10780, 10890, 11100, 11172, 11400, 12300, 12740, 12900, 12936, 13200, 13230, 13524, 13800, 14100, 15210, 15246, 15288, 15600, 15900, 16500, 16660

Offset: 1

Gus Wiseman, Feb 28 2020

nonn

A sequence of positive integers is unimodal if it is the concatenation of a weakly increasing and a weakly decreasing sequence.

A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization.

			The sequence of terms together with their prime indices begins:
   2100: {1,1,2,3,3,4}
   3300: {1,1,2,3,3,5}
   3900: {1,1,2,3,3,6}
   4200: {1,1,1,2,3,3,4}
   4410: {1,2,2,3,4,4}
   5100: {1,1,2,3,3,7}
   5700: {1,1,2,3,3,8}
   6468: {1,1,2,4,4,5}
   6600: {1,1,1,2,3,3,5}
   6900: {1,1,2,3,3,9}
   7644: {1,1,2,4,4,6}
   7800: {1,1,1,2,3,3,6}
   8400: {1,1,1,1,2,3,3,4}
   8700: {1,1,2,3,3,10}
   9300: {1,1,2,3,3,11}
   9996: {1,1,2,4,4,7}
  10200: {1,1,1,2,3,3,7}
  10500: {1,1,2,3,3,3,4}
  10780: {1,1,3,4,4,5}
  10890: {1,2,2,3,5,5}

MathWorld, Unimodal Sequence

Not requiring non-unimodal negation gives A332282.
These are the Heinz numbers of the partitions counted by A332640.
Not requiring non-unimodality gives A332642.
The case of compositions is A332870.
Unimodal compositions are A001523.
Non-unimodal permutations are A059204.
Non-unimodal compositions are A115981.
Unsorted prime signature is A124010.
Non-unimodal normal sequences are A328509.
Partitions whose 0-appended first differences are unimodal are A332283, with Heinz numbers the complement of A332287.
Compositions whose negation is unimodal are A332578.
Compositions whose negation is not unimodal are A332669.
Partitions whose 0-appended first differences are not unimodal are A332744, with Heinz numbers A332832.
Numbers whose signature is neither increasing nor decreasing are A332831.
Cf. A007052, A056239, A072704, A112798, A242031, A242414, A332280, A332281, A332288, A332294, A332639, A332728, A332742.

unimodQ[q_]:=Or[Length[q]<=1,If[q[[1]]<=q[[2]],unimodQ[Rest[q]],OrderedQ[Reverse[q]]]]
Select[Range[10000],!unimodQ[Last/@FactorInteger[#]]&&!unimodQ[-Last/@FactorInteger[#]]&]

Intersection of A332282 and A332642.

cp's OEIS Frontend

A025065 Number of palindromic partitions of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Haskell

Mathematica

PARI

Formula

Extensions

A072774 Powers of squarefree numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

A317089 Numbers whose prime factors span an initial interval of prime numbers and whose prime multiplicities span an initial interval of positive integers.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A069799 The number obtained by reversing the sequence of nonzero exponents in the prime factorization of n with respect to distinct primes present, as ordered by their indices.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A329395 Numbers whose binary expansion without the most significant (first) digit has Lyndon and co-Lyndon factorizations of equal lengths.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A317087 Numbers whose prime factors span an initial interval of prime numbers and whose sequence of prime multiplicities is a palindrome.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

A317086 Number of normal integer partitions of n whose sequence of multiplicities is a palindrome.

Views

Author

Keywords

Comments