A289507 -id:A289507 - cp's OEIS frontend

A000027 The positive integers. Also called the natural numbers, the whole numbers or the counting numbers, but these terms are ambiguous.

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77

Offset: 1

N. J. A. Sloane

For some authors, the terms "natural numbers" and "counting numbers" include 0, i.e., refer to the nonnegative integers A001477; the term "whole numbers" frequently also designates the whole set of (signed) integers A001057.
a(n) is smallest positive integer which is consistent with sequence being monotonically increasing and satisfying a(a(n)) = n (cf. A007378).
Inverse Euler transform of A000219.
The rectangular array having A000027 as antidiagonals is the dispersion of the complement of the triangular numbers, A000217 (which triangularly form column 1 of this array). The array is also the transpose of A038722. - Clark Kimberling, Apr 05 2003
For nonzero x, define f(n) = floor(nx) - floor(n/x). Then f=A000027 if and only if x=tau or x=-tau. - Clark Kimberling, Jan 09 2005
Numbers of form (2^i)*k for odd k (i.e., n = A006519(n)*A000265(n)); thus n corresponds uniquely to an ordered pair (i,k) where i=A007814, k=A000265 (with A007814(2n)=A001511(n), A007814(2n+1)=0). - Lekraj Beedassy, Apr 22 2006
If the offset were changed to 0, we would have the following pattern: a(n)=binomial(n,0) + binomial(n,1) for the present sequence (number of regions in 1-space defined by n points), A000124 (number of regions in 2-space defined by n straight lines), A000125 (number of regions in 3-space defined by n planes), A000127 (number of regions in 4-space defined by n hyperplanes), A006261, A008859, A008860, A008861, A008862 and A008863, where the last six sequences are interpreted analogously and in each "... by n ..." clause an offset of 0 has been assumed, resulting in a(0)=1 for all of them, which corresponds to the case of not cutting with a hyperplane at all and therefore having one region. - Peter C. Heinig (algorithms(AT)gmx.de), Oct 19 2006
Define a number of points on a straight line to be in general arrangement when no two points coincide. Then these are the numbers of regions defined by n points in general arrangement on a straight line, when an offset of 0 is assumed. For instance, a(0)=1, since using no point at all leaves one region. The sequence satisfies the recursion a(n) = a(n-1) + 1. This has the following geometrical interpretation: Suppose there are already n-1 points in general arrangement, thus defining the maximal number of regions on a straight line obtainable by n-1 points, and now one more point is added in general arrangement. Then it will coincide with no other point and act as a dividing wall thereby creating one new region in addition to the a(n-1)=(n-1)+1=n regions already there, hence a(n)=a(n-1)+1. Cf. the comments on A000124 for an analogous interpretation. - Peter C. Heinig (algorithms(AT)gmx.de), Oct 19 2006
The sequence a(n)=n (for n=1,2,3) and a(n)=n+1 (for n=4,5,...) gives to the rank (minimal cardinality of a generating set) for the semigroup I_n\S_n, where I_n and S_n denote the symmetric inverse semigroup and symmetric group on [n]. - James East, May 03 2007
The sequence a(n)=n (for n=1,2), a(n)=n+1 (for n=3) and a(n)=n+2 (for n=4,5,...) gives the rank (minimal cardinality of a generating set) for the semigroup PT_n\T_n, where PT_n and T_n denote the partial transformation semigroup and transformation semigroup on [n]. - James East, May 03 2007
"God made the integers; all else is the work of man." This famous quotation is a translation of "Die ganzen Zahlen hat der liebe Gott gemacht, alles andere ist Menschenwerk," spoken by Leopold Kronecker in a lecture at the Berliner Naturforscher-Versammlung in 1886. Possibly the first publication of the statement is in Heinrich Weber's "Leopold Kronecker," Jahresberichte D.M.V. 2 (1893) 5-31. - Clark Kimberling, Jul 07 2007
Binomial transform of A019590, inverse binomial transform of A001792. - Philippe Deléham, Oct 24 2007
Writing A000027 as N, perhaps the simplest one-to-one correspondence between N X N and N is this: f(m,n) = ((m+n)^2 - m - 3n + 2)/2. Its inverse is given by I(k)=(g,h), where g = k - J(J-1)/2, h = J + 1 - g, J = floor((1 + sqrt(8k - 7))/2). Thus I(1)=(1,1), I(2)=(1,2), I(3)=(2,1) and so on; the mapping I fills the first-quadrant lattice by successive antidiagonals. - Clark Kimberling, Sep 11 2008
a(n) is also the mean of the first n odd integers. - Ian Kent, Dec 23 2008
Equals INVERTi transform of A001906, the even-indexed Fibonacci numbers starting (1, 3, 8, 21, 55, ...). - Gary W. Adamson, Jun 05 2009
These are also the 2-rough numbers: positive integers that have no prime factors less than 2. - Michael B. Porter, Oct 08 2009
Totally multiplicative sequence with a(p) = p for prime p. Totally multiplicative sequence with a(p) = a(p-1) + 1 for prime p. - Jaroslav Krizek, Oct 18 2009
Triangle T(k,j) of natural numbers, read by rows, with T(k,j) = binomial(k,2) + j = (k^2-k)/2 + j where 1 <= j <= k. In other words, a(n) = n = binomial(k,2) + j where k is the largest integer such that binomial(k,2) < n and j = n - binomial(k,2). For example, T(4,1)=7, T(4,2)=8, T(4,3)=9, and T(4,4)=10. Note that T(n,n)=A000217(n), the n-th triangular number. - Dennis P. Walsh, Nov 19 2009
Hofstadter-Conway-like sequence (see A004001): a(n) = a(a(n-1)) + a(n-a(n-1)) with a(1) = 1, a(2) = 2. - Jaroslav Krizek, Dec 11 2009
a(n) is also the dimension of the irreducible representations of the Lie algebra sl(2). - Leonid Bedratyuk, Jan 04 2010
Floyd's triangle read by rows. - Paul Muljadi, Jan 25 2010
Number of numbers between k and 2k where k is an integer. - Giovanni Teofilatto, Mar 26 2010
Generated from a(2n) = r*a(n), a(2n+1) = a(n) + a(n+1), r = 2; in an infinite set, row 2 of the array shown in A178568. - Gary W. Adamson, May 29 2010
1/n = continued fraction [n]. Let barover[n] = [n,n,n,...] = 1/k. Then k - 1/k = n. Example: [2,2,2,...] = (sqrt(2) - 1) = 1/k, with k = (sqrt(2) + 1). Then 2 = k - 1/k. - Gary W. Adamson, Jul 15 2010
Number of n-digit numbers the binary expansion of which contains one run of 1's. - Vladimir Shevelev, Jul 30 2010
From Clark Kimberling, Jan 29 2011: (Start)
Let T denote the "natural number array A000027":
   1    2    4    7 ...
   3    5    8   12 ...
   6    9   13   18 ...
  10   14   19   25 ...
T(n,k) = n+(n+k-2)*(n+k-1)/2. See A185787 for a list of sequences based on T, such as rows, columns, diagonals, and sub-arrays. (End)
The Stern polynomial B(n,x) evaluated at x=2. See A125184. - T. D. Noe, Feb 28 2011
The denominator in the Maclaurin series of log(2), which is 1 - 1/2 + 1/3 - 1/4 + .... - Mohammad K. Azarian, Oct 13 2011
As a function of Bernoulli numbers B_n (cf. A027641: (1, -1/2, 1/6, 0, -1/30, 0, 1/42, ...)): let V = a variant of B_n changing the (-1/2) to (1/2). Then triangle A074909 (the beheaded Pascal's triangle) * [1, 1/2, 1/6, 0, -1/30, ...] = the vector [1, 2, 3, 4, 5, ...]. - Gary W. Adamson, Mar 05 2012
Number of partitions of 2n+1 into exactly two parts. - Wesley Ivan Hurt, Jul 15 2013
Integers n dividing u(n) = 2u(n-1) - u(n-2); u(0)=0, u(1)=1 (Lucas sequence A001477). - Thomas M. Bridge, Nov 03 2013
For this sequence, the generalized continued fraction a(1)+a(1)/(a(2)+a(2)/(a(3)+a(3)/(a(4)+...))), evaluates to 1/(e-2) = A194807. - Stanislav Sykora, Jan 20 2014
Engel expansion of e-1 (A091131 = 1.71828...). - Jaroslav Krizek, Jan 23 2014
a(n) is the number of permutations of length n simultaneously avoiding 213, 231 and 321 in the classical sense which are breadth-first search reading words of increasing unary-binary trees. For more details, see the entry for permutations avoiding 231 at A245898. - Manda Riehl, Aug 05 2014
a(n) is also the number of permutations simultaneously avoiding 213, 231 and 321 in the classical sense which can be realized as labels on an increasing strict binary tree with 2n-1 nodes. See A245904 for more information on increasing strict binary trees. - Manda Riehl, Aug 07 2014
a(n) = least k such that 2*Pi - Sum_{h=1..k} 1/(h^2 - h + 3/16) < 1/n. - Clark Kimberling, Sep 28 2014
a(n) = least k such that Pi^2/6 - Sum_{h=1..k} 1/h^2 < 1/n. - Clark Kimberling, Oct 02 2014
Determinants of the spiral knots S(2,k,(1)). a(k) = det(S(2,k,(1))). These knots are also the torus knots T(2,k). - Ryan Stees, Dec 15 2014
As a function, the restriction of the identity map on the nonnegative integers {0,1,2,3...}, A001477, to the positive integers {1,2,3,...}. - M. F. Hasler, Jan 18 2015
See also A131685(k) = smallest positive number m such that c(i) = m (i^1 + 1) (i^2 + 2) ... (i^k+ k) / k! takes integral values for all i>=0: For k=1, A131685(k)=1, which implies that this is a well defined integer sequence. - Alexander R. Povolotsky, Apr 24 2015
a(n) is the number of compositions of n+2 into n parts avoiding the part 2. - Milan Janjic, Jan 07 2016
Does not satisfy Benford's law [Berger-Hill, 2017] - N. J. A. Sloane, Feb 07 2017
Parametrization for the finite multisubsets of the positive integers, where, for p_j the j-th prime, n = Product_{j} p_j^(e_j) corresponds to the multiset containing e_j copies of j ('Heinz encoding' -- see A056239, A003963, A289506, A289507, A289508, A289509). - Christopher J. Smyth, Jul 31 2017
The arithmetic function v_1(n,1) as defined in A289197. - Robert Price, Aug 22 2017
For n >= 3, a(n)=n is the least area that can be obtained for an irregular octagon drawn in a square of n units side, whose sides are parallel to the axes, with 4 vertices that coincide with the 4 vertices of the square, and the 4 remaining vertices having integer coordinates. See Affaire de Logique link. - Michel Marcus, Apr 28 2018
a(n+1) is the order of rowmotion on a poset defined by a disjoint union of chains of length n. - Nick Mayers, Jun 08 2018
Number of 1's in n-th generation of 1-D Cellular Automata using Rules 50, 58, 114, 122, 178, 186, 206, 220, 238, 242, 250 or 252 in the Wolfram numbering scheme, started with a single 1. - Frank Hollstein, Mar 25 2019
(1, 2, 3, 4, 5, ...) is the fourth INVERT transform of (1, -2, 3, -4, 5, ...). - Gary W. Adamson, Jul 15 2019

T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 1.
T. M. Apostol, Modular Functions and Dirichlet Series in Number Theory, Springer-Verlag, 1990, page 25.
John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 22.
W. Fulton and J. Harris, Representation theory: a first course, (1991), page 149. [From Leonid Bedratyuk, Jan 04 2010]
I. S. Gradstein and I. M. Ryshik, Tables of series, products, and integrals, Volume 1, Verlag Harri Deutsch, 1981.
R. E. Schwartz, You Can Count on Monsters: The First 100 numbers and Their Characters, A. K. Peters and MAA, 2010.
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).

N. J. A. Sloane, Table of n, a(n) for n = 1..500000 [a large file]
Archimedes Laboratory, What's special about this number?
Affaire de Logique, Pick et Pick et Colegram (in French), No. 1051, 18-04-2018.
Paul Barry, A Catalan Transform and Related Transformations on Integer Sequences, Journal of Integer Sequences, Vol. 8 (2005), Article 05.4.5.
James Barton, The Numbers
A. Berger and T. P. Hill, What is Benford's Law?, Notices, Amer. Math. Soc., 64:2 (2017), 132-134.
A. Breiland, L. Oesper, and L. Taalman, p-Coloring classes of torus knots, Online Missouri J. Math. Sci., 21 (2009), 120-126.
N. Brothers, S. Evans, L. Taalman, L. Van Wyk, D. Witczak, and C. Yarnall, Spiral knots, Missouri J. of Math. Sci., 22 (2010).
C. K. Caldwell, Prime Curios
Case and Abiessu, interesting number
S. Crandall, notes on interesting digital ephemera
O. Curtis, Interesting Numbers
M. DeLong, M. Russell, and J. Schrock, Colorability and determinants of T(m,n,r,s) twisted torus knots for n equiv. +/-1(mod m), Involve, Vol. 8 (2015), No. 3, 361-384.
Walter Felscher, Historia Matematica Mailing List Archive.
P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009; see page 371.
Robert R. Forslund, A logical alternative to the existing positional number system, Southwest Journal of Pure and Applied Mathematics, Vol. 1 1995 pp. 27-29.
E. Friedman, What's Special About This Number?
R. K. Guy, Letter to N. J. A. Sloane
Milan Janjic, Enumerative Formulas for Some Functions on Finite Sets
Kival Ngaokrajang, Illustration about relation to many other sequences, when the sequence is considered as a triangular table read by its antidiagonals. Additional illustrations when the sequence is considered as a centered triangular table read by rows.
Mike Keith, All Numbers Are Interesting: A Constructive Approach
Leonardo of Pisa [Leonardo Pisano], Illustration of initial terms, from Liber Abaci [The Book of Calculation], 1202 (photo by David Singmaster).
Feihu Liu, Guoce Xin, and Chen Zhang, Ehrhart Polynomials of Order Polytopes: Interpreting Combinatorial Sequences on the OEIS, arXiv:2412.18744 [math.CO], 2024. See pp. 15, 24.
Robert Munafo, Notable Properties of Specific Numbers.
G. Pfeiffer, Counting Transitive Relations, Journal of Integer Sequences, Vol. 7 (2004), Article 04.3.2.
R. Phillips, Numbers from one to thirty-one
J. Striker, Dynamical Algebraic Combinatorics: Promotion, Rowmotion, and Resonance, Notices of the AMS, June/July 2017, pp. 543-549.
G. Villemin's Almanac of Numbers, NOMBRES en BREF (in French)
Eric Weisstein's World of Mathematics, Natural Number, Positive Integer, Counting Number Composition, Davenport-Schinzel Sequence, Idempotent Number, N, Smarandache Ceil Function, Whole Number, Engel Expansion, and Trinomial Coefficient
Wikipedia, List of numbers, Interesting number paradox, and Floyd's triangle
Robert G. Wilson v, English names for the numbers from 0 to 11159 without spaces or hyphens
Robert G. Wilson v, American English names for the numbers from 0 to 100999 without spaces or hyphens
Index entries for "core" sequences
Index entries for sequences of the a(a(n)) = 2n family
Index entries for sequences that are permutations of the natural numbers
Index entries for related partition-counting sequences
Index entries for linear recurrences with constant coefficients, signature (2,-1).
Index to divisibility sequences
Index entries for sequences related to Benford's law

A001477 = nonnegative numbers.
Partial sums of A000012.
Cf. A001478, A001906, A007931, A007932, A027641, A074909, A089353 (multisets), A178568, A194807.
Cf. A026081 = integers in reverse alphabetical order in U.S. English, A107322 = English name for number and its reverse have the same number of letters, A119796 = zero through ten in alphabetical order of English reverse spelling, A005589, etc. Cf. A185787 (includes a list of sequences based on the natural number array A000027).
Cf. Boustrophedon transforms: A000737, A231179;
Cf. A038722 (mirrored when seen as triangle), A056011 (boustrophedon).
Cf. A048993, A048994, A000110 (see the Feb 03 2015 formula).
Cf. A289187, A000010, A008683, A000203, A049310.

a000027 = id
a000027_list = [1..]  -- Reinhard Zumkeller, May 07 2012

print([n for n in 1:280]) # Paul Muljadi, Apr 09 2024

Magma
```
[ n : n in [1..100]];
```

A000027 := n->n; seq(A000027(n), n=1..100);

Range@ 77 (* Robert G. Wilson v, Mar 31 2015 *)

makelist(n, n, 1, 30); /* Martin Ettl, Nov 07 2012 */

PARI
```
{a(n) = n};
```

print join(", ", 1..280) # Paul Muljadi, May 29 2024

def A000027(n): return n # Chai Wah Wu, May 09 2022

R
```
1:100
```
Shell
```
seq 1 100
```

a(2k+1) = A005408(k), k >= 0, a(2k) = A005843(k), k >= 1.
Multiplicative with a(p^e) = p^e. - David W. Wilson, Aug 01 2001
Another g.f.: Sum_{n>0} phi(n)*x^n/(1-x^n) (Apostol).
When seen as an array: T(k, n) = n+1 + (k+n)*(k+n+1)/2. Main diagonal is 2n*(n+1)+1 (A001844), antidiagonal sums are n*(n^2+1)/2 (A006003). - Ralf Stephan, Oct 17 2004
Dirichlet generating function: zeta(s-1). - Franklin T. Adams-Watters, Sep 11 2005
G.f.: x/(1-x)^2. E.g.f.: x*exp(x). a(n)=n. a(-n)=-a(n).
Series reversion of g.f. A(x) is x*C(-x)^2 where C(x) is the g.f. of A000108. - Michael Somos, Sep 04 2006
G.f. A(x) satisfies 0 = f(A(x), A(x^2)) where f(u, v) = u^2 - v - 4*u*v. - Michael Somos, Oct 03 2006
Convolution of A000012 (the all-ones sequence) with itself. - Tanya Khovanova, Jun 22 2007
a(n) = 2*a(n-1)-a(n-2); a(1)=1, a(2)=2. a(n) = 1+a(n-1). - Philippe Deléham, Nov 03 2008
a(n) = A000720(A000040(n)). - Juri-Stepan Gerasimov, Nov 29 2009
a(n+1) = Sum_{k=0..n} A101950(n,k). - Philippe Deléham, Feb 10 2012
a(n) = Sum_{d | n} phi(d) = Sum_{d | n} A000010(d). - Jaroslav Krizek, Apr 20 2012
G.f.: x * Product_{j>=0} (1+x^(2^j))^2 = x * (1+2*x+x^2) * (1+2*x^2+x^4) * (1+2*x^4+x^8) * ... = x + 2x^2 + 3x^3 + ... . - Gary W. Adamson, Jun 26 2012
a(n) = det(binomial(i+1,j), 1 <= i,j <= n). - Mircea Merca, Apr 06 2013
E.g.f.: x*E(0), where E(k) = 1 + 1/(x - x^3/(x^2 + (k+1)/E(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 03 2013
From Wolfdieter Lang, Oct 09 2013: (Start)
a(n) = Product_{k=1..n-1} 2*sin(Pi*k/n), n > 1.
a(n) = Product_{k=1..n-1} (2*sin(Pi*k/(2*n)))^2, n > 1.
These identities are used in the calculation of products of ratios of lengths of certain lines in a regular n-gon. For the first identity see the Gradstein-Ryshik reference, p. 62, 1.392 1., bringing the first factor there to the left hand side and taking the limit x -> 0 (L'Hôpital). The second line follows from the first one. Thanks to Seppo Mustonen who led me to consider n-gon lengths products. (End)
a(n) = Sum_{j=0..k} (-1)^(j-1)*j*binomial(n,j)*binomial(n-1+k-j,k-j), k>=0. - Mircea Merca, Jan 25 2014
a(n) = A052410(n)^A052409(n). - Reinhard Zumkeller, Apr 06 2014
a(n) = Sum_{k=1..n^2+2*n} 1/(sqrt(k)+sqrt(k+1)). - Pierre CAMI, Apr 25 2014
a(n) = floor(1/sin(1/n)) = floor(cot(1/(n+1))) = ceiling(cot(1/n)). - Clark Kimberling, Oct 08 2014
a(n) = floor(1/(log(n+1)-log(n))). - Thomas Ordowski, Oct 10 2014
a(k) = det(S(2,k,1)). - Ryan Stees, Dec 15 2014
a(n) = 1/(1/(n+1) + 1/(n+1)^2 + 1/(n+1)^3 + ...). - Pierre CAMI, Jan 22 2015
a(n) = Sum_{m=0..n-1} Stirling1(n-1,m)*Bell(m+1), for n >= 1. This corresponds to Bell(m+1) = Sum_{k=0..m} Stirling2(m, k)*(k+1), for m >= 0, from the fact that Stirling2*Stirling1 = identity matrix. See A048993, A048994 and A000110. - Wolfdieter Lang, Feb 03 2015
a(n) = Sum_{k=1..2n-1}(-1)^(k+1)*k*(2n-k). In addition, surprisingly, a(n) = Sum_{k=1..2n-1}(-1)^(k+1)*k^2*(2n-k)^2. - Charlie Marion, Jan 05 2016
G.f.: x/(1-x)^2 = (x * r(x) *r(x^3) * r(x^9) * r(x^27) * ...), where r(x) = (1 + x + x^2)^2 = (1 + 2x + 3x^2 + 2x^3 + x^4). - Gary W. Adamson, Jan 11 2017
a(n) = floor(1/(Pi/2-arctan(n))). - Clark Kimberling, Mar 11 2020
a(n) = Sum_{d|n} mu(n/d)*sigma(d). - Ridouane Oudra, Oct 03 2020
a(n) = Sum_{k=1..n} phi(gcd(n,k))/phi(n/gcd(n,k)). - Richard L. Ollerton, May 09 2021
a(n) = S(n-1, 2), with the Chebyshev S-polynomials A049310. - Wolfdieter Lang, Mar 09 2023
From Peter Bala, Nov 02 2024: (Start)
For positive integer m, a(n) = (1/m)* Sum_{k = 1..2*m*n-1} (-1)^(k+1) * k * (2*m*n - k) = (1/m) * Sum_{k = 1..2*m*n-1} (-1)^(k+1) * k^2 * (2*m*n - k)^2 (the case m = 1 is given above).
a(n) = Sum_{k = 0..3*n} (-1)^(n+k+1) * k * binomial(3*n+k, 2*k). (End)

Links edited by Daniel Forgues, Oct 07 2009.

A289509 Numbers k such that the gcd of the indices j for which the j-th prime prime(j) divides k is 1.

Original entry on oeis.org

2, 4, 6, 8, 10, 12, 14, 15, 16, 18, 20, 22, 24, 26, 28, 30, 32, 33, 34, 35, 36, 38, 40, 42, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 58, 60, 62, 64, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 82, 84, 85, 86, 88, 90, 92, 93, 94, 95, 96, 98, 99, 100, 102, 104

Offset: 1

Christopher J. Smyth, Jul 11 2017

nonn

Any integer k in the sequence encodes (by 'Heinz encoding' cf. A056239) a multiset of integers whose gcd is 1, namely the multiset containing r_j copies of j if k factors as Product_j prime(j)^{r_j} with gcd_j j = 1.
Clearly the sequence contains all even numbers and no odd primes or odd prime powers. It also clearly contains all numbers that are divisible by consecutive primes.
The sequence is the list of those k such that A289508(k) = 1.
It is also the list of those k such that A289506(k) = A289507(k).
Heinz numbers of integer partitions with relatively prime parts, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 13 2018

			6 is a term because 6 = p_1*p_2 and gcd(1,2) = 1.
From _Gus Wiseman_, Apr 13 2018: (Start)
Sequence of integer partitions with relatively prime parts begins:
02 : (1)
04 : (11)
06 : (21)
08 : (111)
10 : (31)
12 : (211)
14 : (41)
15 : (32)
16 : (1111)
18 : (221)
20 : (311)
22 : (51)
24 : (2111)
26 : (61)
28 : (411)
30 : (321)
32 : (11111)
33 : (52)
34 : (71)
35 : (43)
36 : (2211)
38 : (81)
40 : (3111)
(End)

Alois P. Heinz, Table of n, a(n) for n = 1..20000

Cf. A001222, A007359, A051424, A056239, A289506, A289507, A289508, A296150, A302696, A302697, A302698, A302796.

p:=1:for ind to 10000 do p:=nextprime(p);primeindex[p]:=ind;od:
out:=[]:for n from 2 to 100 do m:=[];f:=ifactors(n)[2];g:=0;
for k to nops(f) do mk:=primeindex[f[k][1]];m:=[op(m),mk];
g:=gcd(g,mk);od; if g=1 then out:=[op(out),n];fi;od:out;

Select[Range[200],GCD@@PrimePi/@FactorInteger[#][[All,1]]===1&] (* Gus Wiseman, Apr 13 2018 *)

isok(n) = my(f=factor(n)); gcd(apply(x->primepi(x), f[,1])) == 1; \\ Michel Marcus, Jul 19 2017

from sympy import gcd, primepi, primefactors
def ok(n): return gcd([primepi(p) for p in primefactors(n)]) == 1
print([n for n in range(1, 151) if ok(n)]) # Indranil Ghosh, Aug 06 2017

A289508 a(n) is the GCD of the indices j for which the j-th prime p_j divides n.

Original entry on oeis.org

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

Offset: 1

Christopher J. Smyth, Jul 11 2017

The number n = Product_j p_j can be regarded as an index for the multiset of all the j's, occurring with multiplicity corresponding to the highest power of p_j dividing n. Then a(n) is the gcd of the elements of this multiset. Compare A056239, where the same encoding for integer multisets('Heinz encoding') is used, but where A056239(n) is the sum, rather than the gcd, of the elements of the corresponding multiset (partition) of the j's. Cf. also A003963, for which A003963(n) is the product of the elements of the corresponding multiset.

a(m*n) = gcd(a(m),a(n)). - Robert Israel, Jul 19 2017

			a(n) = 1 for all even n as 2 = p_1. Also a(p_j) = j.
Further, a(703) = 4 because 703 = p_8.p_{12} and gcd(8,12) = 4.

Alois P. Heinz, Table of n, a(n) for n = 1..20000

Cf. A289506, A289507.

f:=  n -> igcd(op(map(numtheory:-pi, numtheory:-factorset(n)))):
map(f, [$1..100]); # Robert Israel, Jul 19 2017

Table[GCD @@ Map[PrimePi, FactorInteger[n][[All, 1]] ], {n, 2, 83}] (* Michael De Vlieger, Jul 19 2017 *)

a(n) = my(f=factor(n)); gcd(apply(x->primepi(x), f[,1])); \\ Michel Marcus, Jul 19 2017

from sympy import primefactors, primepi, gcd
def a(n):
    return gcd([primepi(d) for d in primefactors(n)])
print([a(n) for n in range(2, 101)]) # Indranil Ghosh, Jul 20 2017

a(n) = gcd_j j, where p_j divides n.

a(n) = A289506(n)/A289507(n).

A355737 Number of ways to choose a sequence of divisors, one of each prime index of n (with multiplicity), such that the result has no common divisor > 1.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 2, 1, 3, 4, 1, 1, 4, 1, 2, 4, 2, 1, 2, 3, 4, 7, 3, 1, 4, 1, 1, 4, 2, 6, 4, 1, 4, 6, 2, 1, 6, 1, 2, 8, 3, 1, 2, 5, 4, 4, 4, 1, 8, 4, 3, 5, 4, 1, 4, 1, 2, 10, 1, 6, 4, 1, 2, 6, 6, 1, 4, 1, 6, 8, 4, 6, 8, 1, 2, 15, 2, 1, 6, 4, 4

Offset: 1

Gus Wiseman, Jul 17 2022

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.

			The a(2) = 1 through a(18) = 4 choices:
  1  1  11  1  11  1  111  11  11  1  111  1  11  11  1111  1  111
               12          12  13     112     12  13           112
                           21                 14  21           121
                                                  23           122

Wikipedia, Coprime integers.

Dominated by A355731, firsts A355732, primes A355741, prime-powers A355742.
For weakly increasing instead of coprime we have A355735, primes A355745.
Positions of first appearances are A355738.
For strict instead of coprime we have A355739, zeros A355740.
A000005 counts divisors.
A001221 counts distinct prime factors, with sum A001414.
A001222 counts prime factors with multiplicity.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A120383 lists numbers divisible by all of their prime indices.
A289508 gives GCD of prime indices.
A289509 ranks relatively prime partitions, odd A302697, squarefree A302796.
A324850 lists numbers divisible by the product of their prime indices.
Cf. A000720, A007359, A051424, A076610, A289507, A296150, A302696, A302698, A355733, A355744, A355748.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[Select[Tuples[Divisors/@primeMS[n]],GCD@@#==1&]],{n,100}]

A289506 Write n as a product of primes p_{s_1}p_{s_2}p_{s_3}*... where p_i denotes the i-th prime; then a(n) = s_1^2 + s_2^2 + s_3^2 + ...

Original entry on oeis.org

0, 1, 4, 2, 9, 5, 16, 3, 8, 10, 25, 6, 36, 17, 13, 4, 49, 9, 64, 11, 20, 26, 81, 7, 18, 37, 12, 18, 100, 14, 121, 5, 29, 50, 25, 10, 144, 65, 40, 12, 169, 21, 196, 27, 17, 82, 225, 8, 32, 19, 53, 38, 256, 13, 34, 19, 68, 101, 289, 15, 324, 122, 24, 6

Offset: 1

Christopher J. Smyth, Jul 07 2017

When gcd_j(s_j) = 1, a(n) is the modulus of the determinant whose first row consists of the s_j, and whose remaining rows form a lattice basis for the space of integer solutions of Sum_j s_jx_j = 0. See A289507.
Compare A056239, where the same encoding for integer multisets ('Heinz encoding') is used, but where A056239(n) is the sum, rather than the sum of squares, of the elements of the corresponding multiset (partition).
See also A003963, for which A003963(n) is the product of the elements of the corresponding multiset.
See also A289507, where terms are (Sum_j s_j^2)/gcd_j(s_j) rather than Sum_j s_j^2 (this sequence).

			For n = 12 = 2^2 * 3 = p_1 * p_1 * p_2, the multiset is {1,1,2} and so a(12) = 1^2 + 1^2 + 2^2 = 6.
Also a(1) = 0 as n = 1 indexes the empty multiset.
Further a(p_k) = k^2 and a(2^r) = r.

Alois P. Heinz, Table of n, a(n) for n = 1..20000

Cf. A003963, A056239, A289507.

p:=1:for ind to 1000 do p:=nextprime(p);primeindex[p]:=ind;od: # so primeindex[p]:=k if p is the k-th prime
out:=[0]:for n from 2 to 100 do f:=ifactors(n)[2];
m:=[];for k to nops(f) do pow:=f[k];ind:=primeindex[pow[1]];for e to pow[2] do
m:=[op(m),ind];od;od;out:=[op(out),sum(m[jj]^2,jj=1..nops(m))];
od:print(out);
# second Maple program:
a:= n-> add(numtheory[pi](i[1])^2*i[2], i=ifactors(n)[2]):
seq(a(n), n=1..80);  # Alois P. Heinz, Aug 05 2017

Table[Total[FactorInteger[n] /. {p_, e_} /; p > 0 :> e PrimePi[p]^2], {n, 64}] (* Michael De Vlieger, Jul 12 2017 *)

a(n) = my(f=factor(n)); sum(k=1, #f~, primepi(f[k,1])^2*f[k,2]); \\ Michel Marcus, Jul 19 2017

For n = Product_k p_k^{r_k}, a(n) = Sum_k k^2 * r_k.

Also a(n) = Sum_j s_j^2, where the multiset of s_j's is the multiset of k's, each with multiplicity r_k.

A318585 Number of integer partitions of n whose sum of reciprocals squared is an integer.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 4, 4, 4, 4, 5, 5, 6, 6, 7, 8, 9, 9, 10, 10, 12, 12, 13, 14, 16, 16, 18, 19, 21, 23, 26, 27, 29, 30, 34, 35, 39, 43, 48, 51, 55, 57, 63, 67, 74, 78, 84, 89, 99, 103, 112, 119, 132, 139, 148, 156, 170, 182, 199

Offset: 1

Gus Wiseman, Aug 29 2018

nonn

From David A. Corneth, Sep 03 2018: (Start)
Let a valid tuple be a tuple of positive integers whose sum of reciprocals squared is an integer. Initially one only needs to consider tuples of positive integers where each element is > 1. After that some ones could be prepended to a valid tuple to find new valid tuples.
One could define a prime tuple as a valid tuple where no proper part with elements is a valid tuple. So (1) would be a prime tuple as no proper part of (1) has elements and is a valid tuple. Other examples of prime tuples are (2, 2, 2, 2) and (2, 2, 2, 3, 3, 6).
The list of distinct elements in a tuple could be whittled down by finding for each positive integer m the least sum of a prime tuple in which that integer is. For each m, that sum is at most m^3. (End)

			The a(26) = 7 integer partitions:
  (6332222222)
  (44442221111)
  (63322211111111)
  (22222222222211)
  (222222221111111111)
  (2222111111111111111111)
  (11111111111111111111111111)

Giovanni Resta, Table of n, a(n) for n = 1..130

Cf. A000041, A002865, A051908, A058360, A289506, A289507, A316854, A316856.

Cf. A316854, A318584, A318586, A318587, A318588, A318589.

Table[Length[Select[IntegerPartitions[n],IntegerQ[Total[#^(-2)]]&]],{n,30}]

a(61)-a(70) from Giovanni Resta, Sep 03 2018

A318586 Number of integer partitions of n whose sum of reciprocals squared is the reciprocal of an integer.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 3, 1, 2, 3, 3, 1, 4, 1, 3, 1, 2, 1, 5, 2, 1, 4, 5, 1, 5, 1, 6, 3, 2, 4, 8, 2, 4, 2, 6, 3, 9, 2, 4, 7, 5, 4, 11, 8, 7, 8, 9, 5, 12, 5, 16, 5, 10, 5, 25, 10, 9, 13, 18, 12, 18, 6, 11, 14, 22, 9, 24, 11, 21, 22, 25, 24, 23, 28, 32

Offset: 1

Gus Wiseman, Aug 29 2018

nonn

			The a(42) = 9 integer partitions:
  (42)
  (21,14,7)
  (18,9,9,6)
  (18,9,9,3,3)
  (20,10,4,4,4)
  (12,12,12,4,2)
  (10,5,5,5,5,5,5,2)
  (12,6,6,4,4,4,2,2,2)
  (6,6,4,4,4,4,3,3,3,3,2)

Giovanni Resta, Table of n, a(n) for n = 1..150

Cf. A000041, A051908, A058360, A289506, A289507, A316854, A316857.

Cf. A318584, A318585, A318587, A318588, A318589.

Table[Length[Select[IntegerPartitions[n],IntegerQ[1/Total[#^(-2)]]&]],{n,30}]

a(61)-a(80) from Giovanni Resta, Sep 03 2018

A318573 Numerator of the reciprocal sum of the integer partition with Heinz number n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 29 2018

The reciprocal sum of (y_1, ..., y_k) is 1/y_1 + ... + 1/y_k. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Antti Karttunen, Table of n, a(n) for n = 1..10000
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
Gus Wiseman, Sequences counting and ranking integer partitions by their reciprocal sums
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to Heinz numbers

Positions of 1's are A316857.

Cf. A051908, A056239, A058360, A112798, A289506, A289507, A296150, A316854, A316855, A316856, A318574, A325704.

Table[Sum[pr[[2]]/PrimePi[pr[[1]]],{pr,If[n==1,{},FactorInteger[n]]}],{n,100}]//Numerator

A318573(n) = { my(f=factor(n)); numerator(sum(i=1,#f~,f[i, 2]/primepi(f[i, 1]))); }; \\ Antti Karttunen, Nov 17 2019

If n = Product prime(x_i)^y_i is the prime factorization of n, then a(n) is the numerator of Sum y_i/x_i.

More terms from Antti Karttunen, Nov 17 2019

A318584 Number of integer partitions of n whose sum of reciprocals squared is 1.

Original entry on oeis.org

0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 2, 2, 0, 0, 0, 2, 0, 2, 1, 2, 2, 2, 1, 1, 2, 3, 0, 1, 1, 6, 2, 3, 2, 6, 2, 2, 3, 2, 6, 7, 2, 4, 3, 9, 4, 7, 5, 8, 8, 7, 9, 9, 11, 12, 7, 9, 11, 17, 9, 13, 12, 17, 16, 13, 15, 20, 26, 27, 18, 23

Offset: 0

Gus Wiseman, Aug 29 2018

nonn

The a(16) = 1 integer partition:
  (6,3,3,2,2,2)
The a(48) = 2 integer partitions:
  (18,9,9,3,3,2,2,2)
  (6,6,6,6,3,3,3,3,3,3,3,3)
The a(56) = 3 integer partitions:
  (12,6,6,4,4,4,4,4,4,4,2,2)
  (10,6,5,5,5,5,5,5,3,3,2,2)
  (6,6,4,4,4,4,4,4,4,4,3,3,3,3)
The a(60) = 6 integer partitions:
  (12,12,12,12,3,3,2,2,2)
  (8,8,8,8,6,4,4,4,3,3,2,2)
  (6,6,6,6,6,6,6,6,6,2,2,2)
  (12,12,12,4,3,3,3,3,3,3,2)
  (10,5,5,5,5,5,5,4,4,4,4,2,2)
  (6,4,4,4,4,4,4,4,4,4,4,4,4,3,3)

Cf. A000041, A051908, A058360, A289506, A289507, A316854, A316855.

Cf. A318585, A318586, A318587, A318588, A318589.

Table[Length[Select[IntegerPartitions[n],Total[#^(-2)]==1&]],{n,30}]

a(61)-a(100) from Alois P. Heinz, Aug 30 2018

A318587 Heinz numbers of integer partitions whose sum of reciprocals squared is 1.

Original entry on oeis.org

2, 81, 8775, 64827, 950625, 1953125, 7022925, 9055935, 21781575, 36020025, 50124555, 51883209, 57909033, 102984375, 118978125, 760816875, 816747435, 981059625, 1206902781, 1265058675, 1387132263, 2359670625, 3902169375, 4868424351, 5222768733, 5430160125

Offset: 1

Gus Wiseman, Aug 29 2018

nonn

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

			The sequence of integer partitions with Heinz numbers in this sequence begins: (1), (2222), (633222), (4444222), (66333322).

Cf. A001222, A051908, A056239, A112798, A289506, A289507, A316855, A316856, A316857.

Cf. A318584, A318585, A318586, A318588, A318589.

Select[Range[100000],Total[If[#==1,{},Cases[FactorInteger[#],{p_,k_}:>k/PrimePi[p]^2]]]==1&]

a(6)-a(26) from Alois P. Heinz, Aug 30 2018

cp's OEIS Frontend

A000027 The positive integers. Also called the natural numbers, the whole numbers or the counting numbers, but these terms are ambiguous.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Julia

Magma

Maple

Mathematica

Maxima

PARI

Perl

Python

R

Shell

Formula

Extensions

A289509 Numbers k such that the gcd of the indices j for which the j-th prime prime(j) divides k is 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

A289508 a(n) is the GCD of the indices j for which the j-th prime p_j divides n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

A355737 Number of ways to choose a sequence of divisors, one of each prime index of n (with multiplicity), such that the result has no common divisor > 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A289506 Write n as a product of primes p_{s_1}*p_{s_2}*p_{s_3}*... where p_i denotes the i-th prime; then a(n) = s_1^2 + s_2^2 + s_3^2 + ...

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A318585 Number of integer partitions of n whose sum of reciprocals squared is an integer.

Views

Author

Keywords

Comments

Examples

Links

A289506 Write n as a product of primes p_{s_1}p_{s_2}p_{s_3}*... where p_i denotes the i-th prime; then a(n) = s_1^2 + s_2^2 + s_3^2 + ...