A062457 -id:A062457 - cp's OEIS frontend

A324201 a(n) = A062457(A000043(n)) = prime(A000043(n))^A000043(n), where A000043 gives the exponent of the n-th Mersenne prime.

9, 125, 161051, 410338673, 925103102315013629321, 1271991467017507741703714391419, 49593099428404263766544428188098203, 165163983801975082169196428118414326197216835208154294976154161023

Offset: 1

Antti Karttunen, Feb 18 2019

nonn

If there are no odd perfect numbers, then the terms give all solutions n > 1 to A323244(n) = 0.

Conversely, if these are all numbers k > 1 that satisfy A323244(k) = 0 (which can be proved if one can show, for example, that no number in A007916 can satisfy the equation), then no odd perfect numbers exist. See also A336700. - Antti Karttunen, Jan 12 2024

Amiram Eldar, Table of n, a(n) for n = 1..12

Cf. A000043, A000396, A005940, A007916, A062457, A323244, A324200.
Subsequence of A001597.
Cf. also A336700, A368989.

Prime[#]^#&/@MersennePrimeExponent[Range[8]] (* Harvey P. Dale, Mar 15 2024 *)

a(n) = A062457(A000043(n)).
A323244(a(n)) = 0.
a(n) = A005940(1+A000396(n)). [Provided no odd perfect numbers exist]

A123213 a(n) = A062457(n) concatenated with A014284(n).

Original entry on oeis.org

21, 93, 1256, 240111, 16105118, 482680929, 41033867342, 1698356304159, 180115266146378, 420707233300201101, 25408476896404831130, 6582952005840035281161, 925103102315013629321198, 73885357344138503765449239, 12063348350820368238715343282

Offset: 1

Herman Jamke (hermanjamke(AT)fastmail.fm), Oct 05 2006

			a(1) = 21 = 2 concatenated with 1.
a(2) = 93 = 9 concatenated with 3.
a(3) = 1256 = 125 concatenated with 6.
a(4) = 240111 = 2401 concatenated with 11.

Sergio Silva, Teste Numerico.

Cf. A014284, A062457.

rng=13;A062457=Table[Prime[n]^n, {n, rng}] ;A014284=Accumulate[Join[{1}, Prime[Range[rng-1]]]];Table[FromDigits[Join[IntegerDigits[A062457[[n]]],IntegerDigits[A014284[[n]]]]],{n,rng}] (* James C. McMahon, Nov 17 2024 *)

s=1;print1(21,",");for(n=2,20,s=s+prime(n-1);print1(prime(n)^n,s,","))

A001156 Number of partitions of n into squares.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 2, 3, 4, 4, 4, 5, 6, 6, 6, 8, 9, 10, 10, 12, 13, 14, 14, 16, 19, 20, 21, 23, 26, 27, 28, 31, 34, 37, 38, 43, 46, 49, 50, 55, 60, 63, 66, 71, 78, 81, 84, 90, 98, 104, 107, 116, 124, 132, 135, 144, 154, 163, 169, 178, 192, 201, 209, 220, 235, 247, 256

Offset: 0

N. J. A. Sloane

Number of partitions of n such that number of parts equal to k is multiple of k for all k. - Vladeta Jovovic, Aug 01 2004
Of course p_{4*square}(n)>0. In fact p_{4*square}(32n+28)=3 times p_{4*square}(8n+7) and p_{4*square}(72n+69) is even. These seem to be the only arithmetic properties the function p_{4*square(n)} possesses. Similar results hold for partitions into positive squares, distinct squares and distinct positive squares. - Michael David Hirschhorn, May 05 2005
The Heinz numbers of these partitions are given by A324588. - Gus Wiseman, Mar 09 2019

			p_{4*square}(23)=1 because 23 = 3^2 + 3^2 + 2^2 + 1^2 and there is no other partition of 23 into squares.
G.f.: A(x) = 1 + x + x^2 + x^3 + 2*x^4 + 2*x^5 + 2*x^6 + 2*x^7 +...
such that the g.f. A(x) satisfies the identity [_Paul D. Hanna_]:
A(x) = 1/((1-x)*(1-x^4)*(1-x^9)*(1-x^16)*(1-x^25)*...)
A(x) = 1 + x/(1-x) + x^4/((1-x)*(1-x^4)) + x^9/((1-x)*(1-x^4)*(1-x^9)) + x^16/((1-x)*(1-x^4)*(1-x^9)*(1-x^16)) + ...
From _Gus Wiseman_, Mar 09 2019: (Start)
The a(14) = 6 integer partitions into squares are:
  (941)
  (911111)
  (44411)
  (44111111)
  (41111111111)
  (11111111111111)
while the a(14) = 6 integer partitions in which the multiplicity of k is a multiple of k for all k are:
  (333221)
  (33311111)
  (22222211)
  (2222111111)
  (221111111111)
  (11111111111111)
(End)

N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Alois P. Heinz, Table of n, a(n) for n = 0..10000 (first 1001 terms from T. D. Noe)
J. Bohman et al., Partitions in squares, Nordisk Tidskr. Informationsbehandling (BIT) 19 (1979), 297-301.
J. Bohman et al., Partitions in squares, Nordisk Tidskr. Informationsbehandling (BIT) 19 (1979), 297-301. (Annotated scanned copy)
H. L. Fisher, Letter to N. J. A. Sloane, Mar 16 1989
G. H. Hardy and S. Ramanujan, Asymptotic formulae in combinatory analysis, Proceedings of the London Mathematical Society, 2, XVI, 1917, p. 373.
M. D. Hirschhorn and J. A. Sellers, On a problem of Lehmer on partitions into squares, The Ramanujan Journal 8 (2004), 279-287.
F. Iacobescu, Smarandache Partition Type and Other Sequences, Bull. Pure Appl. Sci. 16E, 237-240, 1997.
Martin Klazar, What is an answer? — remarks, results and problems on PIO formulas in combinatorial enumeration, part I, arXiv:1808.08449 [math.CO], 2018.
Igor Pak, Complexity problems in enumerative combinatorics, arXiv:1803.06636 [math.CO], 2018.
James A. Sellers, Partitions Excluding Specific Polygonal Numbers As Parts, Journal of Integer Sequences, Vol. 7 (2004), Article 04.2.4.
F. Smarandache, Sequences of Numbers Involved in Unsolved Problems, Hexis, Phoenix, 2006.
Florentin Smarandache, Sequences of Numbers Involved in Unsolved Problems, arXiv:math/0604019 [math.GM], 2006.
Eric Weisstein's World of Mathematics, Partition
Eric Weisstein's World of Mathematics, Smarandache Sequences
Eric Weisstein's World of Mathematics, Square Number

Cf. A000041, A000161 (partitions into 2 squares), A000290, A033461, A131799, A218494, A285218, A304046.
Cf. A078134 (first differences).
Cf. A003108, A037444, A046042, A259792, A259793, A294529.
Row sums of A243148.
Euler trans. of A010052 (see also A308297).
Cf. A001462, A003114, A006141, A011757, A039900, A047993, A052335, A062457, A064174, A078135, A109298, A117144.
Cf. A324524, A324572, A324587, A324588.

a001156 = p (tail a000290_list) where
   p _          0 = 1
   p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
-- Reinhard Zumkeller, Oct 31 2012, Aug 14 2011

m:=70; R:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&*[1/(1-x^(k^2)): k in [1..(m+2)]]) )); // G. C. Greubel, Nov 11 2018

b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
      b(n, i-1)+ `if`(i^2>n, 0, b(n-i^2, i))))
    end:
a:= n-> b(n, isqrt(n)):
seq(a(n), n=0..120);  # Alois P. Heinz, May 30 2014

CoefficientList[ Series[Product[1/(1 - x^(m^2)), {m, 70}], {x, 0, 68}], x] (* Or *)
Join[{1}, Table[Length@PowersRepresentations[n, n, 2], {n, 68}]] (* Robert G. Wilson v, Apr 12 2005, revised Sep 27 2011 *)
f[n_] := Length@ IntegerPartitions[n, All, Range@ Sqrt@ n^2]; Array[f, 67] (* Robert G. Wilson v, Apr 14 2013 *)
b[n_, i_] := b[n, i] = If[n==0, 1, If[i<1, 0, b[n, i-1] + If[i^2>n, 0, b[n-i^2, i]]]]; a[n_] := b[n, Sqrt[n]//Floor]; Table[a[n], {n, 0, 120}] (* Jean-François Alcover, Nov 02 2015, after Alois P. Heinz *)

{a(n)=polcoeff(1/prod(k=1, sqrtint(n+1), 1-x^(k^2)+x*O(x^n)), n)} \\ Paul D. Hanna, Mar 09 2012

{a(n)=polcoeff(1+sum(m=1, sqrtint(n+1), x^(m^2)/prod(k=1, m, 1-x^(k^2)+x*O(x^n))), n)} \\ Paul D. Hanna, Mar 09 2012

G.f.: Product_{m>=1} 1/(1-x^(m^2)).
G.f.: Sum_{n>=0} x^(n^2) / Product_{k=1..n} (1 - x^(k^2)). - Paul D. Hanna, Mar 09 2012
a(n) = (1/n)*Sum_{k=1..n} A035316(k)*a(n-k). - Vladeta Jovovic, Nov 20 2002
a(n) = f(n,1,3) with f(x,y,z) = if xReinhard Zumkeller, Nov 08 2009
Conjecture (Jan Bohman, Carl-Erik Fröberg, Hans Riesel, 1979): a(n) ~ c * n^(-alfa) * exp(beta*n^(1/3)), where c = 1/18.79656, beta = 3.30716, alfa = 1.16022. - Vaclav Kotesovec, Aug 19 2015
From Vaclav Kotesovec, Dec 29 2016: (Start)
Correct values of these constants are:
1/c = sqrt(3) * (4*Pi)^(7/6) / Zeta(3/2)^(2/3) = 17.49638865935104978665...
alfa = 7/6 = 1.16666666666666666...
beta = 3/2 * (Pi/2)^(1/3) * Zeta(3/2)^(2/3) = 3.307411783596651987...
a(n) ~ 3^(-1/2) * (4*Pi*n)^(-7/6) * Zeta(3/2)^(2/3) * exp(2^(-4/3) * 3 * Pi^(1/3) * Zeta(3/2)^(2/3) * n^(1/3)). [Hardy & Ramanujan, 1917]
(End)

More terms from Eric W. Weisstein

More terms from Gh. Niculescu (ghniculescu(AT)yahoo.com), Oct 08 2006

A257990 The side-length of the Durfee square of the partition having Heinz number n.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 2, 2, 1, 1, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 2, 2, 1, 2, 1, 1

Offset: 1

Emeric Deutsch, May 18 2015

nonn

The Durfee square of a partition is the largest square that fits inside the Ferrers board of the partition.
We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product(p_j-th prime, j=1...r) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] we get 2*2*3*7*29 = 2436.
In the Maple program the subprogram B yields the partition with Heinz number n.
First appearance of k is a(prime(k)^k) = k. - Gus Wiseman, Apr 12 2019

			a(9)=2; indeed, 9 = 3*3 is the Heinz number of the partition [2,2] and, clearly its Durfee square has side-length =2.

G. E. Andrews, The Theory of Partitions, Addison-Wesley, Reading, Mass. 1976.
G. E. Andrews, K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004, Cambridge.
M. Bona, A Walk Through Combinatorics, World Scientific Publishing Co., 2002.

Alois P. Heinz, Table of n, a(n) for n = 1..20000
Findstat, St000183: The side length of the Durfee square of an integer partition

Positions of 1's are A093641. Positions of 2's are A325164. Positions of 3's are A307386.

Cf. A006918, A056239, A062457, A065770, A112798, A115720, A117485, A215366, A252464, A325163, A325169.

with(numtheory): a := proc (p) local B, S, i: B := proc (n) local nn, j, m: nn := op(2, ifactors(n)): for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: S := {}: for i to nops(B(p)) do if i <= B(p)[nops(B(p))+1-i] then S := `union`(S, {i}) else  end if end do: max(S) end proc: seq(a(n), n = 2 .. 146);
# second Maple program:
a:= proc(n) local l, t;
      l:= sort(map(i-> numtheory[pi](i[1])$i[2], ifactors(n)[2]), `>`);
      for t from nops(l) to 1 by -1 do if l[t]>=t then break fi od; t
    end:
seq(a(n), n=1..120);  # Alois P. Heinz, May 10 2016

a[n_] := a[n] = Module[{l, t}, l = Reverse[Sort[Flatten[Table[PrimePi[ f[[1]] ], {f, FactorInteger[n]}, {f[[2]]}]]]]; For[t = Length[l], t >= 1, t--, If[l[[t]] >= t, Break[]]]; t]; Table[a[n], {n, 1, 120}] (* Jean-François Alcover, Feb 17 2017, after Alois P. Heinz *)

For a partition (p_1 >= p_2 >= ... > = p_r) the side-length of its Durfee square is the largest i such that p_i >=i.

a(1)=0 prepended by Alois P. Heinz, May 10 2016

A109298 Primal codes of finite idempotent functions on positive integers.

Original entry on oeis.org

1, 2, 9, 18, 125, 250, 1125, 2250, 2401, 4802, 21609, 43218, 161051, 300125, 322102, 600250, 1449459, 2701125, 2898918, 4826809, 5402250, 9653618, 20131375, 40262750, 43441281, 86882562, 181182375, 362364750, 386683451, 410338673, 603351125, 773366902, 820677346

Offset: 1

Jon Awbrey, Jul 06 2005

nonn

Finite idempotent functions are identity maps on finite subsets, counting the empty function as the idempotent on the empty set.
From Gus Wiseman, Mar 09 2019: (Start)
Also numbers whose ordered prime signature is equal to the distinct prime indices in increasing order. A prime index of n is a number m such that prime(m) divides n. The ordered prime signature (A124010) is the sequence of multiplicities (or exponents) in a number's prime factorization, taken in order of the prime base. The case where the prime indices are taken in decreasing order is A324571.
Also numbers divisible by prime(k) exactly k times for each prime index k. These are a kind of self-describing numbers (cf. A001462, A304679).
Also Heinz numbers of integer partitions where the multiplicity of m is m for all m in the support (counted by A033461). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
Also products of distinct elements of A062457. For example, 43218 = prime(1)^1 * prime(2)^2 * prime(4)^4.
(End)

			Writing (prime(i))^j as i:j, we have the following table of examples:
Primal Codes of Finite Idempotent Functions on Positive Integers
` ` ` 1 = { }
` ` ` 2 = 1:1
` ` ` 9 = ` ` 2:2
` ` `18 = 1:1 2:2
` ` 125 = ` ` ` ` 3:3
` ` 250 = 1:1 ` ` 3:3
` `1125 = ` ` 2:2 3:3
` `2250 = 1:1 2:2 3:3
` `2401 = ` ` ` ` ` ` 4:4
` `4802 = 1:1 ` ` ` ` 4:4
` 21609 = ` ` 2:2 ` ` 4:4
` 43218 = 1:1 2:2 ` ` 4:4
`161051 = ` ` ` ` ` ` ` ` 5:5
`300125 = ` ` ` ` 3:3 4:4
`322102 = 1:1 ` ` ` ` ` ` 5:5
`600250 = 1:1 ` ` 3:3 4:4
From _Gus Wiseman_, Mar 09 2019: (Start)
The sequence of terms together with their prime indices begins as follows. For example, we have 18: {1,2,2} because 18 = prime(1) * prime(2) * prime(2) has prime signature {1,2} and the distinct prime indices are also {1,2}.
       1: {}
       2: {1}
       9: {2,2}
      18: {1,2,2}
     125: {3,3,3}
     250: {1,3,3,3}
    1125: {2,2,3,3,3}
    2250: {1,2,2,3,3,3}
    2401: {4,4,4,4}
    4802: {1,4,4,4,4}
   21609: {2,2,4,4,4,4}
   43218: {1,2,2,4,4,4,4}
  161051: {5,5,5,5,5}
  300125: {3,3,3,4,4,4,4}
  322102: {1,5,5,5,5,5}
  600250: {1,3,3,3,4,4,4,4}
(End)

David A. Corneth, Table of n, a(n) for n = 1..10000
J. Awbrey, Riffs and Rotes

Cf. A076954, A106177, A108352, A108371, A109297, A109301.
Cf. A001156, A033461, A056239, A062457, A112798, A118914, A124010 (ordered prime signature), A181819, A276078, A304679.
Cf. A324524, A324525, A324570, A324571, A324572, A324587, A324588.
Sequences related to self-description: A000002, A001462, A079000, A079254, A276625, A304360.

Select[Range[10000],And@@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]==k]&]

is(n) = my(f = factor(n)); for(i = 1, #f~, if(prime(f[i, 2]) != f[i, 1], return(0))); 1 \\ David A. Corneth, Mar 09 2019

Sum_{n>=1} 1/a(n) = Product_{n>=1} (1 + 1/prime(n)^n) = 1.6807104966... - Amiram Eldar, Jan 03 2021

Offset set to 1, missing terms inserted and more terms added by Alois P. Heinz, Mar 08 2019

A276078 Numbers n in whose prime factorization no exponent of any prime(k) exceeds k.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 23, 25, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 50, 51, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 85, 86, 87, 89, 90, 91, 93, 94, 95, 97, 98, 99, 101, 102, 103, 105, 106, 107, 109, 110, 111, 113, 114, 115, 117, 118, 119, 121

Offset: 1

Antti Karttunen, Aug 18 2016

nonn

Numbers not divisible by p^(1+A000720(p)) for any prime p, where A000720(p) gives the index of prime p: 1 for 2, 2 for 3, 3 for 5, and so on.
Also Heinz numbers of integer partitions where the multiplicity of i does not exceed i for any i (A052335). Differs from A048103 in lacking {625, 1250, 1875, 3750, 4375, 5625, 6875, 8125, 8750, ...}. - Gus Wiseman, Mar 09 2019
Asymptotic density is Product_{i>=1} 1-prime(i)^(-1-i) = 0.72102334... - Amiram Eldar, Oct 20 2020

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

Positions of zeros in A276077.
Complement: A276079.
Sequence A276076 sorted into ascending order.
Cf. A000040, A000720.
Subsequence of A048103 from which it differs for the first time at n=451, where a(451) = 626, while A048103(451) = 625, a value missing from here.
Cf. A052335, A056239, A062457, A087153, A109298, A112798, A115584, A117144, A118914, A124010, A197987, A276078.
Cf. A324524, A324525, A324571, A324587, A324588.

Select[Range@ 121, Or[# == 1, AllTrue[FactorInteger[#], PrimePi[#1] >= #2 & @@ # &]] &] (* Michael De Vlieger, Jun 24 2017 *)

isok(n) = my(f=factor(n)); for (k=1, #f~, if (f[k, 2] > primepi(f[k, 1]), return(0))); return (1); \\ Michel Marcus, Jun 24 2017

is(n) = {my(t=1);forprime(p = 2, , t++; pp = p^t; if(n%pp==0, return(0)); if(pp > n, return(1)))} \\ David A. Corneth, Jun 24 2017

upto(n) = {my(v = vector(n,i,1), t=1, res=List()); forprime(p=2, , t++; pp = p^t; if(pp>n, break); for(i=1, n\pp, v[pp*i] = 0)); for(i=1, n, if(v[i]==1, listput(res, i))); res} \\ David A. Corneth, Jun 24 2017

from sympy import factorint, primepi
def ok(n):
    f = factorint(n)
    return all(f[i] <= primepi(i) for i in f)
print([n for n in range(1, 151) if ok(n)]) # Indranil Ghosh, Jun 24 2017

A076954 a(n) = Product_{i=1..n} prime(i)^i.

Original entry on oeis.org

1, 2, 18, 2250, 5402250, 870037764750, 4199506113235182750, 1723219765760312626547490750, 29266411525287522788837599332989370750, 52713275010243038997421106186697438702252144407250, 22176856087751973465466098269669474342964368337745368642450857250

Offset: 0

Amarnath Murthy, Oct 20 2002

nonn

For n >= 3, a(n) ends in 250 or 750. - Jianing Song, Jan 27 2019

Gheorghe Coserea, Table of n, a(n) for n = 0..50

Cf. A076955, A006939.

Partial products of A062457. - Michel Marcus, Jan 28 2019

Maple
```
seq(mul(ithprime(i)^i,i=1..n),n=0..13);
```

Table[Times@@Table[Prime[i]^i, {i, n}], {n, 10}] (* Alonso del Arte, Sep 30 2011 *)
Join[{1},FoldList[Times,Table[Prime[i]^i,{i,10}]]] (* Harvey P. Dale, Dec 16 2022 *)

a(n) = { if (n <= 0, return(1)); prod(i = 1, n, prime(i)^i); }
vector(11, i, a(i-1))  \\ Gheorghe Coserea, Aug 24 2015

a(0) = 1, a(n+1) = a(n)*{prime(n+1)^(n+1)}.

More terms from Sascha Kurz, Jan 22 2003

A323244 a(1) = 0; and for n > 1, a(n) = A033879(A156552(n)).

Original entry on oeis.org

0, 1, 1, 2, 1, 4, 1, 6, 0, 5, 1, 10, 1, 16, 2, 6, 1, 12, 1, 18, -3, 18, 1, 22, -4, 46, 4, 22, 1, 10, 1, 30, 14, 82, -2, 14, 1, 256, -12, 22, 1, 36, 1, 66, 8, 226, 1, 46, -12, 19, 8, 130, 1, 28, -19, 70, -12, 748, 1, 42, 1, 1362, 16, 22, 10, 42, 1, 214, 254, 40, 1, 38, 1, 3838, 10, 406, -10, 106, 1, 78, -12, 5458, 1, 26, -72, 12250, -348, 30, 1, 12

Offset: 1

Antti Karttunen, Jan 10 2019

sign

After a(1) = 0, the other zeros occur for k >= 1, at A005940(1+A000396(k)), which, provided no odd perfect numbers exist, is equal to A324201(k) = A062457(A000043(k)): 9, 125, 161051, 410338673, ..., etc.

There are 2321 negative terms among the first 10000 terms.

Antti Karttunen, Table of n, a(n) for n = 1..10000 (based on Hans Havermann's factorization of A156552)
Index entries for sequences related to binary expansion of n
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to sigma(n)

Cf. A000043, A000396, A033879, A064989, A156552, A297112, A323240, A323243, A323245, A323248, A324115, A324051, A324103, A324396, A324398, A324543, A324713.
Cf. A324201 (positions of zeros, conjectured), A324551 (of negative terms), A324720 (of nonnegative terms), A324721 (of positive terms), A324731, A324732.
Cf. A329644 (Möbius transform).
Cf. A323174, A324055, A324185, A324546 for other permutations of deficiency, and also A324574, A324575, A324654.

Array[2 # - If[# == 0, 0, DivisorSigma[1, #]] &@ Floor@ Total@ Flatten@ MapIndexed[#1 2^(#2 - 1) &, Flatten[Table[2^(PrimePi@ #1 - 1), {#2}] & @@@ FactorInteger@ #]] &, 90] (* Michael De Vlieger, Apr 21 2019 *)

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)};
A156552(n) = if(1==n, 0, if(!(n%2), 1+(2*A156552(n/2)), 2*A156552(A064989(n))));
A323244(n) = if(1==n, 0, my(k=A156552(n)); (2*k)-sigma(k));

from sympy import divisor_sigma, primepi, factorint
def A323244(n): return (lambda n: (n<<1)-divisor_sigma(n))(sum((1< 1 else 0 # Chai Wah Wu, Mar 10 2023

a(n) = 2*A156552(n) - A323243(n).
a(1) = 0; and for n > 1, a(n) = A033879(A156552(n)).
a(n) = A323248(n) + A001222(n) = (A323247(n) - A323243(n)) + A001222(n).
From Antti Karttunen, Mar 12 2019 & Nov 23 2019: (Start)
a(n) = Sum_{d|n} (2*A297112(d) - A324543(d)) = Sum_{d|n} A329644(d).
A002487(a(n)) = A324115(n).
a(n) = A329638(n) - A329639(n).
a(n) = A329645(n) - A329646(n).
(End)

A352827 Heinz numbers of integer partitions y with a fixed point y(i) = i. Such a fixed point is unique if it exists.

Original entry on oeis.org

2, 4, 8, 9, 15, 16, 18, 21, 27, 30, 32, 33, 36, 39, 42, 45, 51, 54, 57, 60, 63, 64, 66, 69, 72, 78, 81, 84, 87, 90, 93, 99, 102, 108, 111, 114, 117, 120, 123, 125, 126, 128, 129, 132, 135, 138, 141, 144, 153, 156, 159, 162, 168, 171, 174, 175, 177, 180, 183

Offset: 1

Gus Wiseman, Apr 06 2022

nonn

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

			The terms together with their prime indices begin:
    2: (1)
    4: (1,1)
    8: (1,1,1)
    9: (2,2)
   15: (3,2)
   16: (1,1,1,1)
   18: (2,2,1)
   21: (4,2)
   27: (2,2,2)
   30: (3,2,1)
   32: (1,1,1,1,1)
   33: (5,2)
   36: (2,2,1,1)
   39: (6,2)
   42: (4,2,1)
   45: (3,2,2)
   51: (7,2)
   54: (2,2,2,1)
For example, the partition (3,2,2) with Heinz number 45 has a fixed point at position 2, so 45 is in the sequence.

* = unproved
*These partitions are counted by A001522, strict A352829.
*The complement is A352826, counted by A064428.
The complement reverse version is A352830, counted by A238394.
The reverse version is A352872, counted by A238395
A000700 counts self-conjugate partitions, ranked by A088902.
A001222 counts prime indices, distinct A001221.
A008290 counts permutations by fixed points, unfixed A098825.
A056239 adds up prime indices, row sums of A112798 and A296150.
A115720 and A115994 count partitions by their Durfee square.
A122111 represents partition conjugation using Heinz numbers.
A124010 gives prime signature, sorted A118914, conjugate rank A238745.
A238349 counts compositions by fixed points, complement A352523.
A238352 counts reversed partitions by fixed points, rank statistic A352822.
A352828 counts strict partitions without a fixed point.
A352833 counts partitions by fixed points.
Cf. A062457, A064410, A065770, A093641, A257990, A342192, A352486, A352823, A352824 (characteristic function), A352825, A352831.

pq[y_]:=Length[Select[Range[Length[y]],#==y[[#]]&]];
Select[Range[100],pq[Reverse[Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]]==1&]

A352826 Heinz numbers of integer partitions y without a fixed point y(i) = i. Such a fixed point is unique if it exists.

Original entry on oeis.org

1, 3, 5, 6, 7, 10, 11, 12, 13, 14, 17, 19, 20, 22, 23, 24, 25, 26, 28, 29, 31, 34, 35, 37, 38, 40, 41, 43, 44, 46, 47, 48, 49, 50, 52, 53, 55, 56, 58, 59, 61, 62, 65, 67, 68, 70, 71, 73, 74, 75, 76, 77, 79, 80, 82, 83, 85, 86, 88, 89, 91, 92, 94, 95, 96, 97

Offset: 1

Gus Wiseman, Apr 06 2022

nonn

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.

			The terms together with their prime indices begin:
      1: ()          24: (2,1,1,1)     47: (15)
      3: (2)         25: (3,3)         48: (2,1,1,1,1)
      5: (3)         26: (6,1)         49: (4,4)
      6: (2,1)       28: (4,1,1)       50: (3,3,1)
      7: (4)         29: (10)          52: (6,1,1)
     10: (3,1)       31: (11)          53: (16)
     11: (5)         34: (7,1)         55: (5,3)
     12: (2,1,1)     35: (4,3)         56: (4,1,1,1)
     13: (6)         37: (12)          58: (10,1)
     14: (4,1)       38: (8,1)         59: (17)
     17: (7)         40: (3,1,1,1)     61: (18)
     19: (8)         41: (13)          62: (11,1)
     20: (3,1,1)     43: (14)          65: (6,3)
     22: (5,1)       44: (5,1,1)       67: (19)
     23: (9)         46: (9,1)         68: (7,1,1)

* = unproved
*These partitions are counted by A064428, strict A352828.
The complement is A352827.
The reverse version is A352830, counted by A238394.
A000700 counts self-conjugate partitions, ranked by A088902.
A001222 counts prime indices, distinct A001221.
*A001522 counts partitions with a fixed point.
A008290 counts permutations by fixed points, nonfixed A098825.
A056239 adds up prime indices, row sums of A112798 and A296150.
A115720 and A115994 count partitions by their Durfee square.
A122111 represents partition conjugation using Heinz numbers.
A124010 gives prime signature, sorted A118914.
A238349 counts compositions by fixed points, complement A352523.
A238352 counts reversed partitions by fixed points, rank statistic A352822.
A238395 counts reversed partitions with a fixed point, ranked by A352872.
A352833 counts partitions by fixed points.
Cf. A062457, A064410, A065770, A093641, A257990, A342192, A352486, A352823, A352824, A352825, A352831, A352832.

pq[y_]:=Length[Select[Range[Length[y]],#==y[[#]]&]];
Select[Range[100],pq[Reverse[Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]]==0&]

cp's OEIS Frontend

A324201 a(n) = A062457(A000043(n)) = prime(A000043(n))^A000043(n), where A000043 gives the exponent of the n-th Mersenne prime.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A123213 a(n) = A062457(n) concatenated with A014284(n).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A001156 Number of partitions of n into squares.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Formula

Extensions

A257990 The side-length of the Durfee square of the partition having Heinz number n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A109298 Primal codes of finite idempotent functions on positive integers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A276078 Numbers n in whose prime factorization no exponent of any prime(k) exceeds k.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

PARI

Python

A076954 a(n) = Product_{i=1..n} prime(i)^i.