A326844 -id:A326844 - cp's OEIS frontend

A246277 Column index of n in A246278: a(1) = 0, a(2n) = n, a(2n+1) = a(A064989(2n+1)).

0, 1, 1, 2, 1, 3, 1, 4, 2, 5, 1, 6, 1, 7, 3, 8, 1, 9, 1, 10, 5, 11, 1, 12, 2, 13, 4, 14, 1, 15, 1, 16, 7, 17, 3, 18, 1, 19, 11, 20, 1, 21, 1, 22, 6, 23, 1, 24, 2, 25, 13, 26, 1, 27, 5, 28, 17, 29, 1, 30, 1, 31, 10, 32, 7, 33, 1, 34, 19, 35, 1, 36, 1, 37, 9, 38, 3, 39, 1, 40, 8, 41, 1, 42

Offset: 1

Antti Karttunen, Aug 21 2014

nonn

If n >= 2, n occurs in column a(n) of A246278.

By convention, a(1) = 0 because 1 does not occur in A246278.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Index entries for sequences computed from indices in prime factorization

Terms of A348717 halved. A305897 is the restricted growth sequence transform.
Positions of terms 1 .. 8 in this sequence are given by the following sequences: A000040, A001248, A006094, A030078, A090076, A251720, A090090, A030514.
Cf. A078898 (has the same role with array A083221 as this sequence has with A246278).
Cf. A000040, A001222, A001358, A003961, A055396, A064989, A064216, A243055, A246272, A249810, A249820, A249735, A252463.
This sequence is also used in the definition of the following permutations: A246274, A246276, A246675, A246677, A246683, A249815, A249817 (A249818), A249823, A249825, A250244, A250245, A250247, A250249.
Also in the definition of arrays A249821, A251721, A251722.
Sum of prime indices of a(n) is A359358(n) + A001222(n) - 1, cf. A326844.
A112798 lists prime indices, length A001222, sum A056239.
Cf. A005940, A241916, A253565, A356958, A358195.

a246277[n_Integer] := Module[{f, p, a064989, a},
  f[x_] := Transpose@FactorInteger[x];
  p[x_] := Which[
    x == 1, 1,
    x == 2, 1,
    True, NextPrime[x, -1]];
  a064989[x_] := Times @@ Power[p /@ First[f[x]], Last[f[x]]];
  a[1] = 0;
  a[x_] := If[EvenQ[x], x/2, NestWhile[a064989, x, OddQ]/2];
a/@Range[n]]; a246277[84] (* Michael De Vlieger, Dec 19 2014 *)

A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
A246277(n) = { if(1==n, 0, while((n%2), n = A064989(n)); (n/2)); };

A246277(n) = if(1==n, 0, my(f = factor(n), k = primepi(f[1,1])-1); for (i=1, #f~, f[i,1] = prime(primepi(f[i,1])-k)); factorback(f)/2); \\ Antti Karttunen, Apr 30 2022

from sympy import factorint, prevprime
from operator import mul
from functools import reduce
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 a(n): return 0 if n==1 else n//2 if n%2==0 else a(a064989(n))
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 15 2017

;; two different variants, the second one employing memoizing definec-macro)
(define (A246277 n) (if (= 1 n) 0 (let loop ((n n)) (if (even? n) (/ n 2) (loop (A064989 n))))))
(definec (A246277 n) (cond ((= 1 n) 0) ((even? n) (/ n 2)) (else (A246277 (A064989 n)))))

a(1) = 0, a(2n) = n, a(2n+1) = a(A064989(2n+1)) = a(A064216(n+1)). [Cf. the formula for A252463.]
Instead of the equation for a(2n+1) above, we may write a(A003961(n)) = a(n). - Peter Munn, May 21 2022
Other identities. For all n >= 1, the following holds:
For all w >= 0, a(p_{i} * p_{j} * ... * p_{k}) = a(p_{i+w} * p_{j+w} * ... * p_{k+w}).
For all n >= 2, A001222(a(n)) = A001222(n)-1. [a(n) has one less prime factor than n. Thus each semiprime (A001358) is mapped to some prime (A000040), etc.]
For all n >= 2, a(n) = A078898(A249817(n)).
For semiprimes n = p_i * p_j, j >= i, a(n) = A000040(1+A243055(n)) = p_{1+j-i}.
a(n) = floor(A348717(n)/2). - Antti Karttunen, Apr 30 2022
If n has prime factorization Product_{i=1..k} prime(x_i), then a(n) = Product_{i=2..k} prime(x_i-x_1+1). The opposite version is A358195, prime indices A358172, even bisection A241916. - Gus Wiseman, Dec 29 2022

A241916 a(2^k) = 2^k, and for other numbers, if n = 2^e1 * 3^e2 * 5^e3 * ... p_k^e_k, then a(n) = 2^(e_k - 1) * 3^(e_{k-1}) * ... * p_{k-1}^e2 * p_k^(e1+1). Here p_k is the greatest prime factor of n (A006530), and e_k is its exponent (A071178), and the exponents e1, ..., e_{k-1} >= 0.

Original entry on oeis.org

1, 2, 3, 4, 5, 9, 7, 8, 6, 25, 11, 27, 13, 49, 15, 16, 17, 18, 19, 125, 35, 121, 23, 81, 10, 169, 12, 343, 29, 75, 31, 32, 77, 289, 21, 54, 37, 361, 143, 625, 41, 245, 43, 1331, 45, 529, 47, 243, 14, 50, 221, 2197, 53, 36, 55, 2401, 323, 841, 59, 375, 61, 961, 175, 64

Offset: 1

Antti Karttunen, May 03 2014

nonn

For other numbers than the powers of 2 (that are fixed), this permutation reverses the sequence of exponents in the prime factorization of n from the exponent of 2 to that of the largest prime factor, except that the exponents of 2 and the greatest prime factor present are adjusted by one. Note that some of the exponents might be zeros.
Self-inverse permutation of natural numbers, composition of A122111 & A241909 in either order: a(n) = A122111(A241909(n)) = A241909(A122111(n)).
This permutation preserves both bigomega and the (index of) largest prime factor: for all n it holds that A001222(a(n)) = A001222(n) and A006530(a(n)) = A006530(n) [equally: A061395(a(n)) = A061395(n)].
From the above it follows, that this fixes both primes (A000040) and powers of two (A000079), among other numbers.
Even positions from n=4 onward contain only terms of A070003, and the odd positions only the terms of A102750, apart from 1 which is at a(1), and 2 which is at a(2).

Antti Karttunen, Table of n, a(n) for n = 1..8192
A. Karttunen, A few notes on A122111, A241909 & A241916.
Index entries for sequences that are permutations of the natural numbers

A241912 gives the fixed points; A241913 their complement.
Cf. A006530, A137502, A070003, A102750, A278525.
{A000027, A122111, A241909, A241916} form a 4-group.
The sum of prime indices of a(n) is A243503(n).
Even bisection of A358195 = Heinz numbers of rows of A358172.
A112798 lists prime indices, length A001222, sum A056239.
Cf. A005940, A019565, A161511, A246277, A253565, A326844, A356958, A358137.

nn = 65; f[n_] := If[n == 1, {0}, Function[f, ReplacePart[Table[0, {PrimePi[f[[-1, 1]]]}], #] &@ Map[PrimePi@ First@ # -> Last@ # &, f]]@ FactorInteger@ n]; g[w_List] := Times @@ Flatten@ MapIndexed[Prime[#2]^#1 &, w]; Table[If[IntegerQ@ #, n/4, g@ Reverse@(# - Join[{1}, ConstantArray[0, Length@ # - 2], {1}] &@ f@ n)] &@ Log2@ n, {n, 4, 4 nn, 4}] (* Michael De Vlieger, Aug 27 2016 *)

A209229(n) = (n && !bitand(n,n-1));
A241916(n) = if(1==A209229(n), n, my(f = factor(2*n), nbf = #f~, igp = primepi(f[nbf,1]), g = f); for(i=1,nbf,g[i,1] = prime(1+igp-primepi(f[i,1]))); factorback(g)/2); \\ Antti Karttunen, Jul 02 2018

(define (A241916 n) (A122111 (A241909 n)))

a(1)=1, and for n>1, a(n) = A006530(n) * A137502(n)/2.
a(n) = A122111(A241909(n)) = A241909(A122111(n)).
If 2n has prime factorization Product_{i=1..k} prime(x_i), then a(n) = Product_{i=1..k-1} prime(x_k-x_i+1). The opposite version is A000027, even bisection of A246277. - Gus Wiseman, Dec 28 2022

Description clarified by Antti Karttunen, Jul 02 2018

A261079 Sum of index differences between prime factors of n, summed over all unordered pairs of primes present (with multiplicity) in the prime factorization of n.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Sep 23 2015

nonn

			For n = 1 the prime factorization is empty, thus there is nothing to sum, so a(1) = 0.
For n = 6 = 2*3 = prime(1) * prime(2), a(6) = 1 because the (absolute value of) difference between prime indices of 2 and 3 is 1.
For n = 10 = 2*5 = prime(1) * prime(3), a(10) = 2 because the difference between prime indices of 2 and 5 is 2.
For n = 12 = 2*2*3 = prime(1) * prime(1) * prime(2), a(12) = 2 because the difference between prime indices of 2 and 3 is 1, and the pair (2,3) occurs twice as one can pick either one of the two 2's present in the prime factorization to be a pair of a single 3. Note that the index difference between 2 and 2 is 0, thus the pair (2,2) of prime divisors does not contribute to the sum.
For n = 36 = 2*2*3*3, a(36) = 4 because the index difference between 2 and 3 is 1, and the prime factor pair (2,3) occurs 2^2 = four times in total. As the index difference is zero between 2 and 2 as well as between 3 and 3, the pairs (2,2) and (3,3) do not contribute to the sum.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Index entries for sequences computed from indices in prime factorization

Cf. A000720.
Cf. A000961 (positions of zeros), A006094 (positions of ones).
Cf. A243055, A242411.
Cf. also A260737.
A055396 gives minimum prime index, maximum A061395.
A112798 list prime indices, length A001222, sum A056239.
A304818 adds up partial sums of reversed prime indices, row sums of A359361.
A318283 adds up partial sums of prime indices, row sums of A358136.
Cf. A029931, A316413, A325362, A326836, A326844, A355536, A358137, A359362.

Table[Function[p, Total@ Map[Function[b, Times @@ {First@ Differences@ PrimePi@ b, Count[Subsets[p, {2}], c_ /; SameQ[c, b]]}], Subsets[Union@ p, {2}]]][Flatten@ Replace[FactorInteger@ n, {p_, e_} :> ConstantArray[p, e], 2]], {n, 120}] (* Michael De Vlieger, Mar 08 2017 *)

a(n) = A304818(n) - A318283(n). - Gus Wiseman, Jan 09 2023

a(n) = 2*A304818(n) - A359362(n). - Gus Wiseman, Jan 09 2023

A243503 Sums of parts of partitions (i.e., their sizes) as ordered in the table A241918: a(n) = Sum_{i=A203623(n-1)+2..A203623(n)+1} A241918(i).

Original entry on oeis.org

0, 1, 2, 2, 3, 4, 4, 3, 3, 6, 5, 6, 6, 8, 5, 4, 7, 5, 8, 9, 7, 10, 9, 8, 4, 12, 4, 12, 10, 8, 11, 5, 9, 14, 6, 7, 12, 16, 11, 12, 13, 11, 14, 15, 7, 18, 15, 10, 5, 7, 13, 18, 16, 6, 8, 16, 15, 20, 17, 11, 18, 22, 10, 6, 10, 14, 19, 21, 17, 10, 20, 9, 21, 24, 6, 24

Offset: 1

Antti Karttunen, Jun 05 2014

nonn

Each n occurs A000041(n) times in total.

Where are the first and the last occurrence of each n located?

Antti Karttunen, Table of n, a(n) for n = 1..8192

Cf. A243504 (the products of parts), A241918, A000041, A227183, A075158, A056239, A241909.
Sum of prime indices of A241916, the even bisection of A358195.
Sums of even-indexed rows of A358172.
A112798 lists prime indices, length A001222, sum A056239, max A061395.
Cf. A005940, A019565, A246277, A326844, A356958, A359362.

Table[If[n==1,0,With[{y=Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]},Last[y]*Length[y]+Last[y]-Total[y]+Length[y]-1]],{n,100}] (* Gus Wiseman, Jan 09 2023 *)

a(n) = Sum_{i=A203623(n-1)+2..A203623(n)+1} A241918(i).
a(n) = A056239(A241909(n)).
a(n) = A227183(A075158(n-1)).
a(A000040(n)) = a(A000079(n)) = n for all n >= 1.
a(A122111(n)) = a(n) for all n.
a(A243051(n)) = a(n) for all n, and likewise for A243052, A243053 and other rows of A243060.
a(n) = A061395(n) * A001222(n) + A061395(n) - A056239(n) + A001222(n) - 1. - Gus Wiseman, Jan 09 2023
a(n) = A326844(2n) + A001222(n). - Gus Wiseman, Jan 09 2023

A326845 Sum times maximum of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 4, 2, 9, 6, 16, 3, 8, 12, 25, 8, 36, 20, 15, 4, 49, 10, 64, 15, 24, 30, 81, 10, 18, 42, 12, 24, 100, 18, 121, 5, 35, 56, 28, 12, 144, 72, 48, 18, 169, 28, 196, 35, 21, 90, 225, 12, 32, 21, 63, 48, 256, 14, 40, 28, 80, 110, 289, 21, 324, 132, 32, 6, 54, 40

Offset: 1

Gus Wiseman, Jul 26 2019

nonn

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

Cf. A056239, A061395, A106529, A112798, A316413, A326836, A326844, A326846.

Table[If[n==1,0,With[{y=Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]},Total[y]*Max[y]]],{n,100}]

a(n) = A056239(n) * A061395(n).

A361853 Number of integer partitions of n such that (length) * (maximum) = 2n.

Original entry on oeis.org

0, 0, 0, 0, 0, 2, 0, 1, 2, 4, 0, 10, 0, 8, 16, 10, 0, 31, 0, 44, 44, 20, 0, 92, 50, 28, 98, 154, 0, 266, 0, 154, 194, 48, 434, 712, 0, 60, 348, 910, 0, 1198, 0, 1120, 2138, 88, 0, 2428, 1300, 1680, 912, 2506, 0, 4808, 4800, 5968, 1372, 140, 0, 14820, 0, 160

Offset: 1

Gus Wiseman, Mar 29 2023

nonn

Also partitions satisfying (maximum) = 2*(mean).

These are partitions whose diagram has the same size as its complement (see example).

			The a(6) = 2 through a(12) = 10 partitions:
  (411)   .  (4211)  (621)     (5221)   .  (822)
  (3111)             (321111)  (5311)      (831)
                               (42211)     (6222)
                               (43111)     (6321)
                                           (6411)
                                           (422211)
                                           (432111)
                                           (441111)
                                           (32211111)
                                           (33111111)
The partition y = (6,4,1,1) has diagram:
  o o o o o o
  o o o o . .
  o . . . . .
  o . . . . .
Since the partition and its complement (shown in dots) have the same size, y is counted under a(12).

For minimum instead of mean we have A118096.
For length instead of mean we have A237753.
For median instead of mean we have A361849, ranks A361856.
This is the equal case of A361851, unequal case A361852.
The strict case is A361854.
These partitions have ranks A361855.
This is the equal case of A361906, unequal case A361907.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, A058398 by mean.
A051293 counts subsets with integer mean.
A067538 counts partitions with integer mean.
A268192 counts partitions by complement size, ranks A326844.
Cf. A111907, A116608, A188814, A237755, A237824, A237984, A240219, A326849, A327482, A349156, A359894.

Table[Length[Select[IntegerPartitions[n],Length[#]*Max@@#==2n&]],{n,30}]

A361855 Numbers > 1 whose prime indices satisfy (maximum) * (length) = 2*(sum).

Original entry on oeis.org

28, 40, 78, 84, 171, 190, 198, 220, 240, 252, 280, 351, 364, 390, 406, 435, 714, 748, 756, 765, 777, 784, 814, 840, 850, 925, 988, 1118, 1197, 1254, 1330, 1352, 1419, 1425, 1440, 1505, 1564, 1600, 1638, 1716, 1755, 1794, 1802, 1820, 1950, 2067, 2204, 2254

Offset: 1

Gus Wiseman, Mar 29 2023

nonn

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.
Also positive integers whose prime indices satisfy (maximum) = 2*(mean).
Also Heinz numbers of partitions of the same size as their complement (see example).

			The terms together with their prime indices begin:
   28: {1,1,4}
   40: {1,1,1,3}
   78: {1,2,6}
   84: {1,1,2,4}
  171: {2,2,8}
  190: {1,3,8}
  198: {1,2,2,5}
  220: {1,1,3,5}
  240: {1,1,1,1,2,3}
  252: {1,1,2,2,4}
  280: {1,1,1,3,4}
The prime indices of 84 are {1,1,2,4}, with maximum 4, length 4, and sum 8, and 4*4 = 2*8, so 84 is in the sequence.
The prime indices of 120 are {1,1,1,2,3}, with maximum 3, length 5, and sum 8, and 3*5 != 2*8, so 120 is not in the sequence.
The prime indices of 252 are {1,1,2,2,4}, with maximum 4, length 5, and sum 10, and 4*5 = 2*10, so 252 is in the sequence.
The partition (5,2,2,1) with Heinz number 198 has diagram:
  o o o o o
  o o . . .
  o o . . .
  o . . . .
Since the partition and its complement (shown in dots) both have size 10, 198 is in the sequence.

These partitions are counted by A361853, strict A361854.
For median instead of mean we have A361856, counted by A361849.
For minimum instead of mean we have A361908, counted by A118096.
For length instead of mean we have A361909, counted by A237753.
A001222 (bigomega) counts prime factors, distinct A001221 (omega).
A061395 gives greatest prime index.
A112798 lists prime indices, sum A056239.
A326567/A326568 gives mean of prime indices.
Cf. A067801, A237824, A316413, A324521, A326844, A361205, A361851, A361906.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[2,100],Max@@prix[#]*PrimeOmega[#]==2*Total[prix[#]]&]

A061395(a(n)) * A001222(a(n)) = 2*A056239(a(n)).

A268192 Triangle read by rows: T(n,k) is the number of partitions of weight k among the complements of the partitions of n.

Original entry on oeis.org

1, 2, 2, 1, 3, 0, 2, 2, 2, 0, 2, 1, 4, 0, 2, 1, 2, 0, 2, 2, 2, 2, 2, 0, 4, 0, 0, 2, 1, 4, 1, 2, 0, 6, 0, 2, 2, 1, 0, 2, 0, 2, 3, 2, 0, 6, 0, 2, 4, 4, 0, 2, 0, 2, 2, 0, 0, 2, 1, 4, 0, 6, 0, 2, 4, 5, 0, 6, 0, 4, 2, 0, 0, 4, 1, 0, 0, 2, 0, 2, 2, 4, 0, 2, 6, 5, 0, 6, 0, 8

Offset: 1

Emeric Deutsch, Feb 12 2016

The complement of a partition p[1] >= p[2] >=...>= p[k] is p[1]-p[2], p[1]-p[3], ..., p[1]-p[k]. Its Ferrers board emerges naturally from the Ferrers board of the given partition. The weight of a partition of n is n.
Sum of entries in row n is A000041(n) (the partition numbers).
Apparently, number of entries in row n is A033638(n-1) = 1 + floor((n-1)^2/4).
T(n,0) = A000005(n) = number of divisors of n.
T(n,1) = A070824(n+1).
Sum(k*T(n,k),k>0) = A188814(n).

			Row 4 is 3,0,2; indeed, the complements of [4], [3,1], [2,2], [2,1,1], [1,1,1,1] are: empty, [2], empty, [1,1], empty; their weights are 0, 2, 0, 2, 0, respectively.
From _Gus Wiseman_, Sep 24 2019: (Start)
Triangle begins:
  1
  2
  2 1
  3 0 2
  2 2 0 2 1
  4 0 2 1 2 0 2
  2 2 2 2 0 4 0 0 2 1
  4 1 2 0 6 0 2 2 1 0 2 0 2
  3 2 0 6 0 2 4 4 0 2 0 2 2 0 0 2 1
  4 0 6 0 2 4 5 0 6 0 4 2 0 0 4 1 0 0 2 0 2
  2 4 0 2 6 5 0 6 0 8 4 0 0 6 2 0 2 2 0 2 0 2 0 0 2 1
Row  n = 8 counts the following partitions:
  8          332   53      62       71        521     4211   611      5111
  44               22211   422      2111111   32111          311111   41111
  2222                     431
  11111111                 3221
                           3311
                           221111
(End)

Alois P. Heinz, Rows n = 1..70, flattened

Cf. A000005, A000041, A002620, A033638, A070824, A188814, A326844, A326849.

q := 10: with(combinat): a := proc (i, j) options operator, arrow: partition(i)[j] end proc: P[q] := 0: for j to numbpart(q) do P[q] := sort(P[q]+t^(nops(a(q, j))*max(a(q, j))-q)) end do: P[q] := P[q];
# second Maple program:
b:= proc(n, i, l) option remember; expand(`if`(n=0 or i=1,
      x^(`if`(l=0, 0, n*(l-i))), b(n, i-1, l)+`if`(i>n, 0,
      x^(`if`(l=0, 0, l-i))*b(n-i, i, `if`(l=0, i, l)))))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n$2, 0)):
seq(T(n), n=1..15);  # Alois P. Heinz, Feb 12 2016

b[n_, i_, l_] := b[n, i, l] = Expand[If[n == 0 || i == 1, x^(If[l == 0, 0, n*(l - i)]), b[n, i - 1, l] + If[i > n, 0, x^(If[l == 0, 0, l - i])*b[n - i, i, If[l == 0, i, l]]]]]; T[n_] := Function[p, Table[Coefficient[p, x, i], {i, 0, Exponent[p, x]}]][b[n, n, 0]]; Table[T[n], {n, 1, 15}] // Flatten (* Jean-François Alcover, Dec 22 2016, after Alois P. Heinz *)
Table[Length[Select[IntegerPartitions[n],Max[#]*Length[#]-n==k&]],{n,1,11},{k,0,Floor[(n-1)/2]*Ceiling[(n-1)/2]}] (* Gus Wiseman, Sep 24 2019 *)

The weight of the complement of a partition p is (number of parts of p)*(largest part of p) - weight of p.

For a given q, the Maple program yields the generating polynomial of row q.

A326846 Length times maximum of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 2, 2, 3, 4, 4, 3, 4, 6, 5, 6, 6, 8, 6, 4, 7, 6, 8, 9, 8, 10, 9, 8, 6, 12, 6, 12, 10, 9, 11, 5, 10, 14, 8, 8, 12, 16, 12, 12, 13, 12, 14, 15, 9, 18, 15, 10, 8, 9, 14, 18, 16, 8, 10, 16, 16, 20, 17, 12, 18, 22, 12, 6, 12, 15, 19, 21, 18, 12, 20, 10, 21, 24, 9, 24, 10, 18, 22, 15, 8, 26, 23, 16, 14, 28, 20, 20, 24

Offset: 1

Gus Wiseman, Jul 26 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so a(n) is the size of the minimal rectangle containing the Young digram of the integer partition with Heinz number n.

Cf. A047993, A056239, A061395, A106529, A112798.
Cf. A316413, A326836, A326837, A326844, A326845, A326848, A326849.
Cf. also A331297.

Table[PrimeOmega[n]*PrimePi[FactorInteger[n][[-1,1]]],{n,100}]

A326846(n) = if(1==n, 0, bigomega(n)*primepi(vecmax(factor(n)[, 1]))); \\ Antti Karttunen, Jan 18 2020

a(n) = A001222(n) * A061395(n).

More terms from Antti Karttunen, Jan 18 2020

A361906 Number of integer partitions of n such that (length) * (maximum) >= 2*n.

Original entry on oeis.org

0, 0, 0, 0, 0, 2, 3, 5, 9, 15, 19, 36, 43, 68, 96, 125, 171, 232, 297, 418, 529, 676, 853, 1156, 1393, 1786, 2316, 2827, 3477, 4484, 5423, 6677, 8156, 10065, 12538, 15121, 17978, 22091, 26666, 32363, 38176, 46640, 55137, 66895, 79589, 92621, 111485, 133485

Offset: 1

Gus Wiseman, Mar 29 2023

nonn

Also partitions such that (maximum) >= 2*(mean).

These are partitions whose complement (see example) has size >= n.

			The a(6) = 2 through a(10) = 15 partitions:
  (411)   (511)    (611)     (621)      (721)
  (3111)  (4111)   (4211)    (711)      (811)
          (31111)  (5111)    (5211)     (5221)
                   (41111)   (6111)     (5311)
                   (311111)  (42111)    (6211)
                             (51111)    (7111)
                             (321111)   (42211)
                             (411111)   (43111)
                             (3111111)  (52111)
                                        (61111)
                                        (421111)
                                        (511111)
                                        (3211111)
                                        (4111111)
                                        (31111111)
The partition y = (4,2,1,1) has length 4 and maximum 4, and 4*4 >= 2*8, so y is counted under a(8).
The partition y = (3,2,1,1) has length 4 and maximum 3, and 4*3 is not >= 2*7, so y is not counted under a(7).
The partition y = (3,2,1,1) has diagram:
  o o o
  o o .
  o . .
  o . .
with complement (shown in dots) of size 5, and 5 is not >= 7, so y is not counted under a(7).

For length instead of mean we have A237752, reverse A237755.
For minimum instead of mean we have A237821, reverse A237824.
For median instead of mean we have A361859, reverse A361848.
The unequal case is A361907.
The complement is counted by A361852.
The equal case is A361853, ranks A361855.
Reversing the inequality gives A361851.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, A058398 by mean.
A051293 counts subsets with integer mean.
A067538 counts partitions with integer mean, strict A102627, ranks A316413.
A268192 counts partitions by complement size, ranks A326844.
Cf. A027193, A111907, A116608, A237984, A324521, A327482, A349156, A360068.

Table[Length[Select[IntegerPartitions[n],Length[#]*Max@@#>=2n&]],{n,30}]

cp's OEIS Frontend

A246277 Column index of n in A246278: a(1) = 0, a(2n) = n, a(2n+1) = a(A064989(2n+1)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Scheme

Formula

A241916 a(2^k) = 2^k, and for other numbers, if n = 2^e1 * 3^e2 * 5^e3 * ... p_k^e_k, then a(n) = 2^(e_k - 1) * 3^(e_{k-1}) * ... * p_{k-1}^e2 * p_k^(e1+1). Here p_k is the greatest prime factor of n (A006530), and e_k is its exponent (A071178), and the exponents e1, ..., e_{k-1} >= 0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Scheme

Formula

Extensions

A261079 Sum of index differences between prime factors of n, summed over all unordered pairs of primes present (with multiplicity) in the prime factorization of n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A243503 Sums of parts of partitions (i.e., their sizes) as ordered in the table A241918: a(n) = Sum_{i=A203623(n-1)+2..A203623(n)+1} A241918(i).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A326845 Sum times maximum of the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Formula

A361853 Number of integer partitions of n such that (length) * (maximum) = 2n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A361855 Numbers > 1 whose prime indices satisfy (maximum) * (length) = 2*(sum).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A268192 Triangle read by rows: T(n,k) is the number of partitions of weight k among the complements of the partitions of n.

Views

Author

Keywords