A031215 -id:A031215 - cp's OEIS frontend

A155750 First differences of A031215.

4, 6, 6, 10, 8, 6, 10, 8, 10, 8, 10, 12, 6, 6, 18, 8, 12, 12, 10, 8, 12, 6, 24, 6, 10, 12, 12, 8, 10, 12, 18, 6, 20, 12, 10, 14, 10, 14, 12, 12, 12, 10, 14, 6, 16, 12, 12, 18, 20, 16, 12, 8, 16, 8, 12, 6, 22, 6, 12, 14, 10, 18, 18, 14, 10, 14, 12, 18, 22, 12, 6, 12, 18, 6, 18, 6, 24

Offset: 1

Paul Curtz, Jan 26 2009

All terms are even. Do all even numbers (A005843) appear at least once?

G. C. Greubel, Table of n, a(n) for n = 1..5000

Cf. A005843, A031215.

Cf. A001223, A031131, A110854, A155067.

Table[Prime[2*n+2] - Prime[2 n], {n, 80}] (* G. C. Greubel, Jun 05 2021 *)

[nth_prime(2*n+2) - nth_prime(2*n) for n in (1..80)] # G. C. Greubel, Jun 05 2021

a(n) = A001223(2n) + A001223(2n+1). - R. J. Mathar, Feb 27 2009

a(n) = A031131(2n). - R. J. Mathar, Feb 27 2009

Edited and extended by R. J. Mathar, Feb 27 2009

A066207 All primes that divide n are of the form prime(2k), where prime(k) is k-th prime.

Original entry on oeis.org

1, 3, 7, 9, 13, 19, 21, 27, 29, 37, 39, 43, 49, 53, 57, 61, 63, 71, 79, 81, 87, 89, 91, 101, 107, 111, 113, 117, 129, 131, 133, 139, 147, 151, 159, 163, 169, 171, 173, 181, 183, 189, 193, 199, 203, 213, 223, 229, 237, 239, 243, 247, 251, 259, 261, 263, 267, 271, 273

Offset: 1

Leroy Quet, Dec 16 2001

nonn

The partitions into even parts, encoded by their Heinz numbers. 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. Example: 63 ( = 3*3*7) is in the sequence because it is the Heinz number of the partition [2, 2, 4]. - Emeric Deutsch, May 19 2015
Numbers divisible only by primes with even indices (A031215). - Michael De Vlieger, Dec 11 2017
Numbers that are not divisible by any prime with an odd index (A031368). - Antti Karttunen, Jul 18 2020
For every positive integer m there exists a unique ordered pair of positive integers (j,k) such that m = a(j)*A066208(k). - Christopher Scussel, Jul 01 2023

			39 is included because 3 * 13 = prime(2) * prime(6) and 2 and 6 are both even.

Antti Karttunen, Table of n, a(n) for n = 1..20000 (original first 1000 terms from Harry J. Smith)

Cf. A031215, A031368, A066208, A215366.
Cf. A297002 (a permutation).
Numbers in the odd bisection of A336321.

Select[Range[273], AllTrue[PrimePi@ FactorInteger[#][[All, 1]], EvenQ] &] (* Michael De Vlieger, Dec 11 2017, range adjusted by Antti Karttunen, Jul 18 2020 *)

{ n=0; for (m=2, 10^9, f=factor(m); b=1; for(i=1, matsize(f)[1], if (primepi(f[i, 1])%2, b=0; break)); if (b, write("b066207.txt", n++, " ", m); if (n==1000, return)) ) } \\ Harry J. Smith, Feb 06 2010

isA066207(n) = (!#select(p -> (primepi(p)%2), factor(n)[,1])); \\ Antti Karttunen, Jul 18 2020

Offset changed from 0 to 1 by Harry J. Smith, Feb 06 2010
a(53)-a(58) from Harry J. Smith, Feb 06 2010
a(1) = 1 inserted (and the indexing of the rest of terms changed) by Antti Karttunen, Jul 18 2020

A005917 Rhombic dodecahedral numbers: a(n) = n^4 - (n - 1)^4.

Original entry on oeis.org

1, 15, 65, 175, 369, 671, 1105, 1695, 2465, 3439, 4641, 6095, 7825, 9855, 12209, 14911, 17985, 21455, 25345, 29679, 34481, 39775, 45585, 51935, 58849, 66351, 74465, 83215, 92625, 102719, 113521, 125055, 137345, 150415, 164289, 178991

Offset: 1

N. J. A. Sloane

Final digits of a(n), i.e., a(n) mod 10, are repeated periodically with period of length 5 {1,5,5,5,9}. There is a symmetry in this list since the sum of two numbers equally distant from the ends is equal to 10 = 1 + 9 = 5 + 5 = 2*5. Last two digits of a(n), i.e., a(n) mod 100, are repeated periodically with period of length 50. - Alexander Adamchuk, Aug 11 2006
a(n) = VarScheme(n,2) in the scheme displayed in A128195. - Peter Luschny, Feb 26 2007
If Y is a 3-subset of a 2n-set X then, for n >= 2, a(n-2) is the number of 4-subsets of X intersecting Y. - Milan Janjic, Nov 18 2007
The numbers are the constant number found in magic squares of order n, where n is an odd number, see the comment in A006003. A Magic Square of side 1 is 1; 3 is 15; 5 is 65 and so on. - David Quentin Dauthier, Nov 07 2008
Two times the area of the triangle with vertices at (0,0), ((n - 1)^2, n^2), and (n^2, (n - 1)^2). - J. M. Bergot, Jun 25 2013
Bisection of A006003. - Omar E. Pol, Sep 01 2018
Construct an array M with M(0,n) = 2*n^2 + 4*n + 1 = A056220(n+1), M(n,0) = 2*n^2 + 1 = A058331(n) and M(n,n) = 2*n*(n+1) + 1 = A001844(n). Row(n) begins with all the increasing odd numbers from A058331(n) to A001844(n) and column(n) begins with all the decreasing odd numbers from A056220(n+1) to A001844(n). The sum of the terms in row(n) plus those in column(n) minus M(n,n) equals a(n+1). The first five rows of array M are [1, 7, 17, 31, 49, ...]; [3, 5, 15, 29, 47, ...]; [9, 11, 13, 27, 45, ...]; [19, 21, 23, 25, 43, ...]; [33, 35, 37, 39, 41, ...]. - J. M. Bergot, Jul 16 2013 [This contribution was moved here from A047926 by Petros Hadjicostas, Mar 08 2021.]
For n>=2, these are the primitive sides s of squares of type 2 described in A344332. - Bernard Schott, Jun 04 2021
(a(n) + 1) / 2 = A212133(n) is the number of cells in the n-th rhombic-dodecahedral polycube. - George Sicherman, Jan 21 2024

J. H. Conway and R. K. Guy, The Book of Numbers, p. 53.
E. Deza and M. M. Deza, Figurate Numbers, World Scientific Publishing, 2012, pp. 123-124.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 1..10000
Mario Defranco and Paul E. Gunnells, Hypergraph matrix models and generating functions, arXiv:2204.11361 [math.CO], 2022.
Milan Janjic, Two Enumerative Functions
T. P. Martin, Shells of atoms, Phys. Rep., 273 (1996), 199-241, eq. (9).
Andy Nicol, Illustration of Rhombic Dodecahedral Numbers
C. J. Pita Ruiz V., Some Number Arrays Related to Pascal and Lucas Triangles, J. Int. Seq. 16 (2013) #13.5.7
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985), 4545-4558.
Eric Weisstein's World of Mathematics, Rhombic Dodecahedral Number.
Eric Weisstein's World of Mathematics, Nexus Number.
D. Zeitlin, A family of Galileo sequences, Amer. Math. Monthly 82 (1975), 819-822.
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).

(1/12)*t*(2*n^3 - 3*n^2 + n) + 2*n - 1 for t = 2, 4, 6, ... gives A049480, A005894, A063488, A001845, A063489, A005898, A063490, A057813, A063491, A005902, A063492, A063493, A063494, A063495, A063496.
Column k=3 of A047969.
Cf. A128195, A176850, A005408, A176271, A212133.
Cf. A001844, A000583, A000290.
Cf. A007588, A000292, A000332, A002817, A342354.
Cf. A031215, A008292.
Cf. A016754, A344330, A344332.

a005917 n = a005917_list !! (n-1)
a005917_list = map sum $ f 1 [1, 3 ..] where
   f x ws = us : f (x + 2) vs where (us, vs) = splitAt x ws
-- Reinhard Zumkeller, Nov 13 2014

[n^4 - (n-1)^4: n in [1..50]]; // Vincenzo Librandi, Aug 01 2011

Table[n^4-(n-1)^4,{n,40}]  (* Harvey P. Dale, Apr 01 2011 *)
#[[2]]-#[[1]]&/@Partition[Range[0,40]^4,2,1] (* More efficient than the above Mathematica program because it only has to calculate each 4th power once *) (* Harvey P. Dale, Feb 07 2015 *)
Differences[Range[0,40]^4] (* Harvey P. Dale, Aug 11 2023 *)

a(n)=n^4-(n-1)^4 \\ Charles R Greathouse IV, Jul 31 2011

A005917_list, m = [], [24, -12, 2, 1]
for _ in range(10**2):
    A005917_list.append(m[-1])
    for i in range(3):
        m[i+1] += m[i] # Chai Wah Wu, Dec 15 2015

a(n) = (2*n - 1)*(2*n^2 - 2*n + 1).
Sum_{i=1..n} a(i) = n^4 = A000583(n). First differences of A000583.
G.f.: x*(1+x)*(1+10*x+x^2)/(1-x)^4. - Simon Plouffe in his 1992 dissertation
More generally, g.f. for n^m - (n - 1)^m is Euler(m, x)/(1 - x)^m, where Euler(m, x) is Eulerian polynomial of degree m (cf. A008292). E.g.f.: x*(exp(y/(1 - x)) - exp(x*y/(1 - x)))/(exp(x*y/(1 - x))-x*exp(y/(1 - x))). - Vladeta Jovovic, May 08 2002
a(n) = sum of the next (2*n - 1) odd numbers; i.e., group the odd numbers so that the n-th group contains (2*n - 1) elements like this: (1), (3, 5, 7), (9, 11, 13, 15, 17), (19, 21, 23, 25, 27, 29, 31), ... E.g., a(3) = 65 because 9 + 11 + 13 + 15 + 17 = 65. - Xavier Acloque, Oct 11 2003
a(n) = 2*n - 1 + 12*Sum_{i = 1..n} (i - 1)^2. - Xavier Acloque, Oct 16 2003
a(n) = (4*binomial(n,2) + 1)*sqrt(8*binomial(n,2) + 1). - Paul Barry, Mar 14 2004
Binomial transform of [1, 14, 36, 24, 0, 0, 0, ...], if the offset is 0. - Gary W. Adamson, Dec 20 2007
Sum_{i=1..n-1}(a(i) + a(i+1)) = 8*Sum_{i=1..n}(i^3 + i) = 16*A002817(n-1) for n > 1. - Bruno Berselli, Mar 04 2011
a(n+1) = a(n) + 2*(6*n^2 + 1) = a(n) + A005914(n). - Vincenzo Librandi, Mar 16 2011
a(n) = -a(-n+1). a(n) = (1/6)*(A181475(n) - A181475(n-2)). - Bruno Berselli, Sep 26 2011
a(n) = A045975(2*n-1,n) = A204558(2*n-1)/(2*n - 1). - Reinhard Zumkeller, Jan 18 2012
a(n+1) = Sum_{k=0..2*n+1} (A176850(n,k) - A176850(n-1,k))*(2*k + 1), n >= 1. - L. Edson Jeffery, Nov 02 2012
a(n) = A005408(n-1) * A001844(n-1) = (2*(n - 1) + 1) * (2*(n - 1)*n + 1) = A000290(n-1)*12 + 2 + a(n-1). - Bruce J. Nicholson, May 17 2017
a(n) = A007588(n) + A007588(n-1) = A000292(2n-1) + A000292(2n-2) + A000292(2n-3) = A002817(2n-1) - A002817(2n-2). - Bruce J. Nicholson, Oct 22 2017
a(n) = A005898(n-1) + 6*A000330(n-1) (cf. Deza, Deza, 2012, p. 123, Section 2.6.2). - Felix Fröhlich, Oct 01 2018
a(n) = A300758(n-1) + A005408(n-1). - Bruce J. Nicholson, Apr 23 2020
G.f.: polylog(-4, x)*(1-x)/x. See the Simon Plouffe formula above (with expanded numerator), and the g.f. of the rows of A008292 by Vladeta Jovovic, Sep 02 2002. - Wolfdieter Lang, May 10 2021

A031368 Odd-indexed primes: a(n) = prime(2n-1).

Original entry on oeis.org

2, 5, 11, 17, 23, 31, 41, 47, 59, 67, 73, 83, 97, 103, 109, 127, 137, 149, 157, 167, 179, 191, 197, 211, 227, 233, 241, 257, 269, 277, 283, 307, 313, 331, 347, 353, 367, 379, 389, 401, 419, 431, 439, 449, 461, 467, 487, 499, 509, 523, 547, 563

Offset: 1

Jeff Burch

Appeared as a puzzle in "Stickelers", a nationally distributed feature, by Terry Stickels, Sep 28 2006. - Franklin T. Adams-Watters, Sep 28 2006
Also every second prime starting with 2. - Cino Hilliard, Dec 02 2007
Central terms of the triangle in A005145. - Reinhard Zumkeller, Aug 05 2009

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

Cf. A000040, A031215 (even-indexed primes), A005408.

First differences are A155067.

a031368 = a000040 . ((subtract 1) . (* 2))
a031368_list = map a000040 [1, 3 ..]  -- Reinhard Zumkeller, Nov 25 2012

[ NthPrime(2*n-1): n in [1..1000] ]; // Vincenzo Librandi, Apr 11 2011

A031368 := n->ithprime(2*n-1): seq(A031368(n), n=1..100);

Table[ Prime[ 2n -1], {n, 52}] (* Robert G. Wilson v, Dec 01 2013 *)

a(n) = prime(2*n-1) \\ Jianing Song, Jun 03 2021

a(n) = A219603(n) / A000040(n). - Reinhard Zumkeller, Nov 25 2012

A244990 After 1, numbers whose greatest prime factor is a prime with an even index; n such that A061395(n) is even.

Original entry on oeis.org

1, 3, 6, 7, 9, 12, 13, 14, 18, 19, 21, 24, 26, 27, 28, 29, 35, 36, 37, 38, 39, 42, 43, 48, 49, 52, 53, 54, 56, 57, 58, 61, 63, 65, 70, 71, 72, 74, 76, 78, 79, 81, 84, 86, 87, 89, 91, 95, 96, 98, 101, 104, 105, 106, 107, 108, 111, 112, 113, 114, 116, 117, 122, 126, 129, 130, 131, 133, 139, 140, 142, 143, 144

Offset: 1

Antti Karttunen, Jul 21 2014

nonn

Equally, after 1, natural numbers n such that A006530(n) is in A031215.

A122111 maps each one of these numbers to a unique term of A028260 and vice versa.

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

Complement: A244991.

Cf. A006530, A028260, A031215, A061395, A122111, A244988, A244321, A244322.

For all n, A244988(a(n)) = n.

A360006 Least positive integer whose prime indices have median n/2. a(1) = 1.

Original entry on oeis.org

1, 2, 6, 3, 14, 5, 26, 7, 38, 11, 58, 13, 74, 17, 86, 19, 106, 23, 122, 29, 142, 31, 158, 37, 178, 41, 202, 43, 214, 47, 226, 53, 262, 59, 278, 61, 302, 67, 326, 71, 346, 73, 362, 79, 386, 83, 398, 89, 446, 97, 458, 101, 478, 103, 502, 107, 526, 109, 542, 113

Offset: 1

Gus Wiseman, Jan 24 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.

The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length).

Position of first appearance of n in A360005.
The sorted version is A360007, for mean A360008.
A112798 lists prime indices, length A001222, sum A056239.
A316413 lists numbers whose prime indices have integer mean.
A325347 = partitions w/ integer median, strict A359907, complement A307683.
A326567/A326568 gives mean of prime indices.
A359893 counts partitions by median, cf. A359901, A359902.
A359908 = numbers w/ integer median of prime indices, complement A359912.
Cf. A026424, A359889, A359890, A360009.

nn=100;
prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
seq=Table[If[n==1,1,2*Median[prix[n]]],{n,nn}];
Table[Position[seq,k][[1,1]],{k,Count[Differences[Union[seq]],1]}]

Consists of 1 followed by A000040 interleaved with 2*A031215.

A066066 a(n) = prime(2n) - 2prime(n).

Original entry on oeis.org

-1, 1, 3, 5, 7, 11, 9, 15, 15, 13, 17, 15, 19, 21, 19, 25, 21, 29, 29, 31, 35, 35, 33, 45, 35, 37, 45, 49, 53, 55, 39, 49, 43, 59, 51, 57, 59, 57, 63, 63, 63, 71, 61, 71, 69, 81, 69, 57, 67, 83, 91, 91, 95, 91, 87, 87, 81, 99, 93, 97, 107, 97, 87, 97, 107, 109, 95, 95, 93, 111, 115, 109, 105, 111, 105, 115, 109, 117, 127, 123, 115

Offset: 1

Reinhard Zumkeller, Dec 01 2001

sign

a(n) = A022457(n) for n > 1.
a(n) = A031215(n)-A100484(n) = A072473(n)-A000040(n); see A179740 for primes. - Reinhard Zumkeller, Jul 25 2010
Asymptotically, a(n) ~ log(4) n, with log(4) = 2 log 2 = 1.38629436111989... = A016627. - M. F. Hasler, Oct 19 2013

Harry J. Smith, Table of n, a(n) for n = 1..1000

Table[Prime[2n]-2Prime[n],{n,100}] (* Harvey P. Dale, Aug 21 2016 *)

{ for (n = 1, 1000, a=prime(2*n) - prime(n)*2; write("b066066.txt", n, " ", a) ) } \\ Harry J. Smith, Nov 09 2009

A066066(n)=prime(2*n)-2*prime(n) \\ M. F. Hasler, Oct 19 2013

A338906 Semiprimes whose prime indices sum to an even number.

Original entry on oeis.org

4, 9, 10, 21, 22, 25, 34, 39, 46, 49, 55, 57, 62, 82, 85, 87, 91, 94, 111, 115, 118, 121, 129, 133, 134, 146, 155, 159, 166, 169, 183, 187, 194, 203, 205, 206, 213, 218, 235, 237, 247, 253, 254, 259, 267, 274, 289, 295, 298, 301, 303, 314, 321, 334, 335, 339

Offset: 1

Gus Wiseman, Nov 28 2020

nonn

A semiprime is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.

			The sequence of terms together with their prime indices begins:
      4: {1,1}      87: {2,10}    183: {2,18}    274: {1,33}
      9: {2,2}      91: {4,6}     187: {5,7}     289: {7,7}
     10: {1,3}      94: {1,15}    194: {1,25}    295: {3,17}
     21: {2,4}     111: {2,12}    203: {4,10}    298: {1,35}
     22: {1,5}     115: {3,9}     205: {3,13}    301: {4,14}
     25: {3,3}     118: {1,17}    206: {1,27}    303: {2,26}
     34: {1,7}     121: {5,5}     213: {2,20}    314: {1,37}
     39: {2,6}     129: {2,14}    218: {1,29}    321: {2,28}
     46: {1,9}     133: {4,8}     235: {3,15}    334: {1,39}
     49: {4,4}     134: {1,19}    237: {2,22}    335: {3,19}
     55: {3,5}     146: {1,21}    247: {6,8}     339: {2,30}
     57: {2,8}     155: {3,11}    253: {5,9}     341: {5,11}
     62: {1,11}    159: {2,16}    254: {1,31}    358: {1,41}
     82: {1,13}    166: {1,23}    259: {4,12}    361: {8,8}
     85: {3,7}     169: {6,6}     267: {2,24}    365: {3,21}

A031215 looks at primes instead of semiprimes.
A098350 has this as union of even-indexed antidiagonals.
A300061 looks at all numbers (not just semiprimes).
A338904 has this as union of even-indexed rows.
A338907 is the odd version.
A338908 is the squarefree case.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A056239 gives the sum of prime indices (Heinz weight).
A084126 and A084127 give the prime factors of semiprimes.
A087112 groups semiprimes by greater factor.
A289182/A115392 list the positions of odd/even terms in A001358.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.
A338911 lists products of pairs of primes both of even index.
A339114/A339115 give the least/greatest semiprime of weight n.
Cf. A000040, A001222, A024697, A037143, A112798, A300063, A319242, A320655, A332765, A338910, A339004.

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

from math import isqrt
from sympy import primepi, primerange
def A338906(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum((primepi(x//p)-a>>1) for a,p in enumerate(primerange(isqrt(x)+1),-1))
    return bisection(f,n,n) # Chai Wah Wu, Apr 03 2025

A338906 \/ A338907 = A001358.

A338908 Squarefree semiprimes whose prime indices sum to an even number.

Original entry on oeis.org

10, 21, 22, 34, 39, 46, 55, 57, 62, 82, 85, 87, 91, 94, 111, 115, 118, 129, 133, 134, 146, 155, 159, 166, 183, 187, 194, 203, 205, 206, 213, 218, 235, 237, 247, 253, 254, 259, 267, 274, 295, 298, 301, 303, 314, 321, 334, 335, 339, 341, 358, 365, 371, 377, 382

Offset: 1

Gus Wiseman, Nov 28 2020

nonn

A squarefree semiprime is a product of any two distinct prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.

			The sequence of terms together with their prime indices begins:
     10: {1,3}     115: {3,9}     213: {2,20}
     21: {2,4}     118: {1,17}    218: {1,29}
     22: {1,5}     129: {2,14}    235: {3,15}
     34: {1,7}     133: {4,8}     237: {2,22}
     39: {2,6}     134: {1,19}    247: {6,8}
     46: {1,9}     146: {1,21}    253: {5,9}
     55: {3,5}     155: {3,11}    254: {1,31}
     57: {2,8}     159: {2,16}    259: {4,12}
     62: {1,11}    166: {1,23}    267: {2,24}
     82: {1,13}    183: {2,18}    274: {1,33}
     85: {3,7}     187: {5,7}     295: {3,17}
     87: {2,10}    194: {1,25}    298: {1,35}
     91: {4,6}     203: {4,10}    301: {4,14}
     94: {1,15}    205: {3,13}    303: {2,26}
    111: {2,12}    206: {1,27}    314: {1,37}

A031215 looks at primes instead of semiprimes.
A300061 and A319241 (squarefree) look all numbers (not just semiprimes).
A338905 has this as union of even-indexed rows.
A338906 is the nonsquarefree version.
A338907 is the odd version.
A001358 lists semiprimes, with odd/even terms A046315/A100484.
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes, with odd/even terms A046388/A100484.
A024697 is the sum of semiprimes of weight n.
A025129 is the sum of squarefree semiprimes of weight n.
A056239 gives the sum of prime indices of n.
A289182/A115392 list the positions of odd/even terms in A001358.
A320656 counts factorizations into squarefree semiprimes.
A332765 gives the greatest squarefree semiprime of weight n.
A338898, A338912, and A338913 give the prime indices of semiprimes, with product A087794, sum A176504, and difference A176506.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes, with difference A338900.
A338904 groups semiprimes by weight.
A338911 lists products of pairs of primes both of even index.
A339114/A339115 give the least/greatest semiprime of weight n.
A339116 groups squarefree semiprimes by greater prime factor.
Cf. A000040, A001221, A001222, A087112, A098350, A112798, A168472, A338901, A338904, A339004, A339005.

Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&& EvenQ[Total[PrimePi/@First/@FactorInteger[#]]]&]

A316430 Heinz numbers of integer partitions whose length is equal to the GCD of all the parts.

Original entry on oeis.org

1, 2, 9, 21, 39, 57, 87, 91, 111, 125, 129, 159, 183, 203, 213, 237, 247, 267, 301, 303, 321, 325, 339, 377, 393, 417, 427, 453, 489, 519, 543, 551, 553, 559, 575, 579, 597, 669, 687, 689, 707, 717, 753, 789, 791, 813, 817, 843, 845, 879, 923, 925, 933, 951, 973

Offset: 1

Gus Wiseman, Jul 02 2018

nonn

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

2 is the only even term in the sequence. 3k is in the sequence if and only if k is in A031215. 5k is in the sequence if and only if k = pq with p and q in A031336.

			Sequence of integer partitions whose length is equal to their GCD begins: (), (1), (2,2), (4,2), (6,2), (8,2), (10,2), (6,4), (12,2), (3,3,3), (14,2), (16,2), (18,2), (10,4), (20,2), (22,2), (8,6), (24,2), (14,4), (26,2), (28,2), (6,3,3).

Subsequence of A004280.

Cf. A056239, A067538, A074761, A143773, A289508, A289509, A296150, A316413, A316431, A316432, A316433, A031215, A031336.

Select[Range[200],PrimeOmega[#]==GCD@@Cases[FactorInteger[#],{p_,k_}:>PrimePi[p]]&]

is(n,f=factor(n))=gcd(apply(primepi,f[,1]))==vecsum(f[,2]) \\ Charles R Greathouse IV, Jul 25 2024

a(n) << n log^2 n, can this be improved? - Charles R Greathouse IV, Jul 25 2024

cp's OEIS Frontend

A155750 First differences of A031215.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Sage

Formula

Extensions

A066207 All primes that divide n are of the form prime(2k), where prime(k) is k-th prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Extensions

A005917 Rhombic dodecahedral numbers: a(n) = n^4 - (n - 1)^4.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

Python

Formula

A031368 Odd-indexed primes: a(n) = prime(2n-1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Formula

A244990 After 1, numbers whose greatest prime factor is a prime with an even index; n such that A061395(n) is even.

Views

Author

Keywords

Comments

Links

Crossrefs

Formula

A360006 Least positive integer whose prime indices have median n/2. a(1) = 1.

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

Formula

A066066 a(n) = prime(2*n) - 2*prime(n).

Views

Author

Keywords

Comments

Links

Programs

Mathematica

A066066 a(n) = prime(2n) - 2prime(n).