A006254 -id:A006254 - cp's OEIS frontend

A305901 Filter sequence for all such sequences b, for which b(A006254(k)) = constant for all k >= 3.

1, 2, 3, 4, 5, 4, 4, 6, 4, 4, 7, 4, 8, 9, 4, 4, 10, 11, 4, 12, 4, 4, 13, 4, 14, 15, 4, 16, 17, 4, 4, 18, 19, 4, 20, 4, 4, 21, 22, 4, 23, 4, 24, 25, 4, 26, 27, 28, 4, 29, 4, 4, 30, 4, 4, 31, 4, 32, 33, 34, 35, 36, 37, 4, 38, 4, 39, 40, 4, 4, 41, 42, 43, 44, 4, 4, 45, 46, 4, 47, 48, 4, 49, 4, 50, 51, 4, 52, 53, 4, 4, 54, 55, 56, 57, 4, 4, 58, 4, 4, 59, 60, 61

Offset: 1

Antti Karttunen, Jun 14 2018

nonn

Restricted growth sequence transform of A305900(A064216(n)).
For all i, j:
  a(i) = a(j) => A278223(i) = A278223(j).
  a(i) = a(j) => A253786(i) = A253786(j).

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

Cf. A006254, A064216, A305900.

Cf. also A305902.

up_to = 1000;
partialsums(f,up_to) = { my(v = vector(up_to), s=0); for(i=1,up_to,s += f(i); v[i] = s); (v); }
v_partsums = partialsums(x -> isprime(x+x-1), up_to);
A305901(n) = if(n<=3,n,if(isprime(n+n-1),4,3+n-v_partsums[n]));

For n <= 3, a(n) = n, and for n >= 4, a(n) = 4 if 2n-1 is a prime (for all n in A006254[3..] = 4, 6, 7, 9, 10, 12, 15, ...), and for all other n (numbers n such that 2n-1 is composite), a(n) = running count from 5 onward.

A081385 Incorrect version of A006254 (or A111333).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 16, 17, 18, 19, 20, 26, 27, 28, 29, 30, 36, 37, 38, 39, 40, 46, 47, 48, 49, 50, 51, 52, 53, 54, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65

Offset: 1

Michael Joseph Halm, Apr 20 2003

dead

Name was: Numbers just more than half-prime.

			a(11) = 16 because [(31)/2 + 1/2]=16.

M. J. Halm, More Sequences, Mpossibilities 83, 2003.

M. J. Halm, Sequences

Cf. A000040.

a(n) = ceil(prime(n)/2); \\ Michel Marcus, Aug 13 2013

a(n) = [p(n)/2 + 1/2].

A001248 Squares of primes.

Original entry on oeis.org

4, 9, 25, 49, 121, 169, 289, 361, 529, 841, 961, 1369, 1681, 1849, 2209, 2809, 3481, 3721, 4489, 5041, 5329, 6241, 6889, 7921, 9409, 10201, 10609, 11449, 11881, 12769, 16129, 17161, 18769, 19321, 22201, 22801, 24649, 26569, 27889, 29929, 32041, 32761, 36481

Offset: 1

N. J. A. Sloane

Also 4, together with numbers n such that Sum_{d|n}(-1)^d = -A048272(n) = -3. - Benoit Cloitre, Apr 14 2002
Also, all solutions to the equation sigma(x) + phi(x) = 2x + 1. - Farideh Firoozbakht, Feb 02 2005
Unique numbers having 3 divisors (1, their square root, themselves). - Alexandre Wajnberg, Jan 15 2006
Smallest (or first) new number deleted at the n-th step in an Eratosthenes sieve. - Lekraj Beedassy, Aug 17 2006
Subsequence of semiprimes A001358. - Lekraj Beedassy, Sep 06 2006
Integers having only 1 factor other than 1 and the number itself. Every number in the sequence is a multiple of 1 factor other than 1 and the number itself. 4 : 2 is the only factor other than 1 and 4; 9 : 3 is the only factor other than 1 and 9; and so on. - Rachit Agrawal (rachit_agrawal(AT)daiict.ac.in), Oct 23 2007
The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime. - Omar E. Pol, May 06 2008
There are 2 Abelian groups of order p^2 (C_p^2 and C_p x C_p) and no non-Abelian group. - Franz Vrabec, Sep 11 2008
Also numbers n such that phi(n) = n - sqrt(n). - Michel Lagneau, May 25 2012
For n > 1, n is the sum of numbers from A006254(n-1) to A168565(n-1). - Vicente Izquierdo Gomez, Dec 01 2012
A078898(a(n)) = 2. - Reinhard Zumkeller, Apr 06 2015
Let r(n) = (a(n) - 1)/(a(n) + 1); then Product_{n>=1} r(n) = (3/5) * (4/5) * (12/13) * (24/25) * (60/61) * ... = 2/5. - Dimitris Valianatos, Feb 26 2019
Numbers k such that A051709(k) = 1. - Jianing Song, Jun 27 2021

Ray Chandler, Table of n, a(n) for n = 1..10000 (first 5000 terms from N. J. A. Sloane)
Nicholas John Bizzell-Browning, LIE scales: Composing with scales of linear intervallic expansion, Ph. D. Thesis, Brunel Univ. (UK, 2024). See p. 143.
R. P. Boas and N. J. A. Sloane, Correspondence, 1974
Marius Coman, On the special relation between the numbers of the form 505+ 1008k and the squares of primes, 2015.
Brady Haran and Matt Parker, Squaring Primes, Numberphile video (2018).
Eric Weisstein's World of Mathematics, Prime Power.
OEIS Wiki, Index entries for number of divisors
Index to sequences related to prime signature

Cf. A000005, A000040, A001358, A001694, A002033, A005408, A005722, A006254.
Cf. A008864, A033273, A049001, A048272, A051709, A059956, A060800, A062799.
Cf. A078898, A082020, A085548, A087112, A095874, A168565, A192134, A211110.
Cf. A258599.
Subsequence of A000430, A001358, and of A190641.
Cf. A024450 (partial sums), A069482 (first differences).

a001248 n = a001248_list !! (n-1)
a001248_list = map (^ 2) a000040_list -- Reinhard Zumkeller, Sep 23 2011

[p^2: p in PrimesUpTo(300)]; // Vincenzo Librandi, Mar 27 2014

A001248:=n->ithprime(n)^2; seq(A001248(k), k=1..50); # Wesley Ivan Hurt, Oct 11 2013

Prime[Range[30]]^2 (* Zak Seidov, Dec 07 2011 *)
Select[Range[40000], DivisorSigma[0, #] == 3 &] (* Carlos Eduardo Olivieri, Jun 01 2015 *)

forprime(p=2,1e3,print1(p^2", ")) \\ Charles R Greathouse IV, Jun 10 2011

A001248(n)=prime(n)^2  \\ M. F. Hasler, Sep 16 2012

from sympy import prime
def A001248(n): return prime(n)**2 # Chai Wah Wu, Aug 09 2024

[n^2 for n in prime_range(1,301)] # G. C. Greubel, May 02 2024

n such that A062799(n) = 2. - Benoit Cloitre, Apr 06 2002
A000005(a(n)^(k-1)) = A005408(k) for all k>0. - Reinhard Zumkeller, Mar 04 2007
a(n) = A000040(n)^(3-1)=A000040(n)^2, where 3 is the number of divisors of a(n). - Omar E. Pol, May 06 2008
A000005(a(n)) = 3 or A002033(a(n)) = 2. - Juri-Stepan Gerasimov, Oct 10 2009
A033273(a(n)) = 3. - Juri-Stepan Gerasimov, Dec 07 2009
For n > 2: (a(n) + 17) mod 12 = 6. - Reinhard Zumkeller, May 12 2010
A192134(A095874(a(n))) = A005722(n) + 1. - Reinhard Zumkeller, Jun 26 2011
For n > 2: a(n) = 1 (mod 24). - Zak Seidov, Dec 07 2011
A211110(a(n)) = 2. - Reinhard Zumkeller, Apr 02 2012
a(n) = A087112(n,n). - Reinhard Zumkeller, Nov 25 2012
a(n) = prime(n)^2. - Jon E. Schoenfield, Mar 29 2015
Product_{n>=1} a(n)/(a(n)-1) = Pi^2/6. - Daniel Suteu, Feb 06 2017
Sum_{n>=1} 1/a(n) = P(2) = 0.4522474200... (A085548). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(2)/zeta(4) = 15/Pi^2 (A082020).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(2) = 6/Pi^2 (A059956). (End)

A005097 (Odd primes - 1)/2.

Original entry on oeis.org

1, 2, 3, 5, 6, 8, 9, 11, 14, 15, 18, 20, 21, 23, 26, 29, 30, 33, 35, 36, 39, 41, 44, 48, 50, 51, 53, 54, 56, 63, 65, 68, 69, 74, 75, 78, 81, 83, 86, 89, 90, 95, 96, 98, 99, 105, 111, 113, 114, 116, 119, 120, 125, 128, 131, 134, 135, 138, 140, 141, 146, 153, 155, 156

Offset: 1

N. J. A. Sloane

Or, numbers k such that 2k+1 is prime.
Also numbers not of the form 2xy + x + y. - Jose Brox (tautocrona(AT)terra.es), Dec 29 2005
This sequence arises if you factor the product of a large number of the first odd numbers into the form 3^n(3)5^n(5)7^n(7)11^n(11)... Then n(3)/n(5) = 2, n(3)/n(7) = 3, n(3)/n(11) = 5, ... . - Andrzej Staruszkiewicz (astar(AT)th.if.uj.edu.pl), May 31 2007
Kohen shows: A king invites n couples to sit around a round table with 2n+1 seats. For each couple, the king decides a prescribed distance d between 1 and n which the two spouses have to be seated from each other (distance d means that they are separated by exactly d-1 chairs). We will show that there is a solution for every choice of the distances if and only if 2n+1 is a prime number [i.e., iff n is in A005097], using a theorem known as Combinatorial Nullstellensatz. - Jonathan Vos Post, Jun 14 2010
Starting from 6, positions at which new primes are seen for Goldbach partitions. E.g., 31 is first seen at 34 from 31+3, so position = 1 + (34-6)/2 = 15. - Bill McEachen, Jul 05 2010
Perfect error-correcting Lee codes of word length n over Z: it is conjectured that these always exist when 2n+1 is a prime, as mentioned in Horak. - Jonathan Vos Post, Sep 19 2011
Also solutions to:  A000010(2*n+1) = n * A000005(2*n+1). - Enrique Pérez Herrero, Jun 07 2012
A193773(a(n)) = 1. - Reinhard Zumkeller, Jan 02 2013
I conjecture that the set of pairwise sums of terms of this sequence (A005097) is the set of integers greater than 1, i.e.: 1+1=2, 1+2=3, ..., 5+5=10, ... (This is equivalent to Goldbach's conjecture: every even integer greater than or equal to 6 can be expressed as the sum of two odd primes.) - Lear Young, May 20 2014
See conjecture and comments from Richard R. Forberg, in Links section below, on the relationship of this sequence to rules on values of c that allow both p^q+c and p^q-c to be prime, for an infinite number of primes p. - Richard R. Forberg, Jul 13 2016
The sequence represents the minimum number Ng of gears which are needed to draw a complete graph of order p using a Spirograph(R), where p is an odd prime. The resulting graph consists of Ng hypotrochoids whose respective nodes coincide. If the teethed ring has a circumference p then Ng = (p-1)/2. Examples: A complete graph of order three can be drawn with a Spirograph(R) using a ring with 3n teeth and one gear with n teeth. n is an arbitrary number, only related to the geometry of the gears. A complete graph of order 5 can be drawn using a ring with diameter 5 and 2 gears with diameters 1 and 2 respectively. A complete graph of order 7 can be drawn using a ring with diameter 7 and 3 gears with diameters 1, 2 and 3 respectively. - Bob Andriesse, Mar 31 2017

T. D. Noe, Table of n, a(n) for n = 1..1000
Richard R. Forberg, Comments on A005097
Peter Horak and Bader F. AlBdaiwi, Diameter Perfect Lee Codes, arXiv:1109.3475 [cs.IT], 2011-2012.
Daniel Kohen and Ivan Sadofschi, A New Approach on the Seating Couples Problem, arXiv:1006.2571 [math.CO], 2010.
Dhananjay P. Mehendale, On Hamilton Decompositions, arXiv:0806.0251 [math.GM], 2008.
Eric Weisstein's World of Mathematics, Legendre Symbol

Complement of A047845. Cf. A000040, A006005, A006093.
A130290 is an essentially identical sequence.
Cf. A005384 (subsequence of primes), A266400 (their indices in this sequence).
Numbers n such that 2n+k is prime: this seq(k=1), A067076 (k=3), A089038 (k=5), A105760 (k=7), A155722 (k=9), A101448 (k=11), A153081 (k=13), A089559 (k=15), A173059 (k=17), A153143 (k=19).
Numbers n such that 2n-k is prime: A006254 (k=1), A098090 (k=3), A089253 (k=5), A089192 (k=7), A097069 (k=9), A097338 (k=11), A097363 (k=13), A097480 (k=15), A098605 (k=17), A097932 (k=19).
Cf. also A266409, A294507.

a005097 = (`div` 2) . a065091 -- Reinhard Zumkeller, Jan 02 2013

[n: n in [1..160] |IsPrime(2*n+1)]; // Vincenzo Librandi, Feb 16 2015

with(numtheory): p:=n-> ithprime(n):seq((p(n+1)^2-1)/(2*sigma(p(n+1))), n= 1..64) # Gary Detlefs, May 02 2012

Table[p=Prime[n];(p-1)/2, {n, 2, 22}] (* Vladimir Joseph Stephan Orlovsky, Apr 29 2008 *)
(Prime[Range[2,70]]-1)/2 (* Harvey P. Dale, Jul 11 2020 *)

forprime(p=3,1e4,print1(p>>1", ")) \\ Charles R Greathouse IV, Jun 16 2011

from sympy import prime
def A005097(n): return prime(n+1)//2 # Chai Wah Wu, Jun 04 2022

a(n) = A006093(n)/2 = A000010(A000040(n+1))/2.
a(n) = (prime(n+1)^2-1)/(2*sigma(prime(n+1))) = (A000040(n+1)^2-1)/(2*A000203(A000040(n+1))). - Gary Detlefs, May 02 2012
a(n) = (A065091(n) - 1) / 2. - Reinhard Zumkeller, Jan 02 2013
a(n) ~ n*log(n)/2. - Ilya Gutkovskiy, Jul 11 2016
a(n) = A294507(n) (mod prime(n+1)). - Jonathan Sondow, Nov 04 2017
a(n) = A130290(n+1). - Chai Wah Wu, Jun 04 2022

A005382 Primes p such that 2p-1 is also prime.

Original entry on oeis.org

2, 3, 7, 19, 31, 37, 79, 97, 139, 157, 199, 211, 229, 271, 307, 331, 337, 367, 379, 439, 499, 547, 577, 601, 607, 619, 661, 691, 727, 811, 829, 877, 937, 967, 997, 1009, 1069, 1171, 1237, 1279, 1297, 1399, 1429, 1459, 1531, 1609, 1627, 1657, 1759, 1867, 2011

Offset: 1

N. J. A. Sloane

Sequence gives values of p such Sum_{i=1..p} gcd(p,i) = A018804(p) is prime. - Benoit Cloitre, Jan 25 2002
Let q = 2n-1. For these n (and q), the sum of two cyclotomic polynomials can be written as a product of cyclotomic polynomials and as a cyclotomic polynomial in x^2: Phi(q,x) + Phi(2q,x) = 2 Phi(n,x) Phi(2n,x) = 2 Phi(n,x^2). - T. D. Noe, Nov 04 2003
Primes in A006254. - Zak Seidov, Mar 26 2013
If a(n) is in A168421 then A005383(n) is a twin prime with a Ramanujan prime, A005383(n) - 2. If this sequence has an infinite number of terms in A168421, then the twin prime conjecture can be proved. - John W. Nicholson, Dec 05 2013
Records subsequence of A023509 (n >= 2). - David James Sycamore, May 05 2025

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

T. D. Noe, Table of n, a(n) for n = 1..10000
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972, p. 870.
R. P. Boas and N. J. A. Sloane, Correspondence, 1974
Ajeet Kumar, Subhamoy Maitra, and Chandra Sekhar Mukherjee, On approximate real mutually unbiased bases in square dimension, Cryptography and Communications (2020) Vol. 13, 321-329.
Sagar Mandal, Divisibility and Sequence Properties of sigma+ and phi+, arXiv:2508.11660 [math.GM], 2025. See p. 8.
Marius Tărnăuceanu, Arithmetic progressions in finite groups, arXiv:2003.10060 [math.GR], 2020.
Wikipedia, Cunningham chain

Cf. A005383, A005384 (2p+1), A057326, A057327, A057328, A057329, A057330, A005603, A063908 (2p-3), A063909 (2p-5), A023204 (2p+3), A000384, A001358.
Cf. A010051, A000040, A053685 (subsequence), A006254.
Cf. A023509.

a005382 n = a005382_list !! (n-1)
a005382_list = filter
   ((== 1) . a010051 . (subtract 1) . (* 2)) a000040_list
-- Reinhard Zumkeller, Oct 03 2012

[n: n in [0..1000] | IsPrime(n) and IsPrime(2*n-1)]; // Vincenzo Librandi, Nov 18 2010

f := proc(Q) local t1,i,j; t1 := []; for i from 1 to 500 do j := ithprime(i); if isprime(2*j-Q) then t1 := [op(t1),j]; fi; od: t1; end; f(1);
# second Maple program:
q:= p-> andmap(isprime, [p, 2*p-1]):
select(q, [$2..2500])[];  # Alois P. Heinz, Dec 16 2024

Select[Prime[Range[300]], PrimeQ[2#-1]&]

select(p->isprime(2*p-1),primes(500)) \\ Charles R Greathouse IV, Apr 26 2012

forprime(n=2, 10^3, if(ispseudoprime(2*n-1), print1(n, ", "))) \\ Felix Fröhlich, Jun 15 2014

a(n) = A129521(n) / A005383(n). - Reinhard Zumkeller, Apr 19 2007

a(n) = (A005383(n) + 1)/2. - Zak Seidov, Nov 04 2010

A067076 Numbers k such that 2*k + 3 is a prime.

Original entry on oeis.org

0, 1, 2, 4, 5, 7, 8, 10, 13, 14, 17, 19, 20, 22, 25, 28, 29, 32, 34, 35, 38, 40, 43, 47, 49, 50, 52, 53, 55, 62, 64, 67, 68, 73, 74, 77, 80, 82, 85, 88, 89, 94, 95, 97, 98, 104, 110, 112, 113, 115, 118, 119, 124, 127, 130, 133, 134, 137, 139, 140, 145, 152, 154, 155

Offset: 1

David Williams, Aug 17 2002

The following sequences (allowing offset of first term) all appear to have the same parity: A034953, triangular numbers with prime indices; A054269, length of period of continued fraction for sqrt(p), p prime; A082749, difference between the sum of next prime(n) natural numbers and the sum of next n primes; A006254, numbers n such that 2n-1 is prime; A067076, 2n+3 is a prime. - Jeremy Gardiner, Sep 10 2004
n is in the sequence iff none of the numbers (n-3k)/(2k+1), 1 <= k <= (n-1)/5, is positive integer. - Vladimir Shevelev, May 31 2009
Zeta(s) = Sum_{n>=1} 1/n^s = 1/1 - 2^(-s) * Product_{p=prime=(2*A067076)+3} 1/(1 - (2*A067076+3)^(-s)). - Eric Desbiaux, Dec 15 2009
This sequence is a subsequence of A047949. - Jason Kimberley, Aug 30 2012

Harry J. Smith, Table of n, a(n) for n = 1..1000
Mutsumi Suzuki, Vincenzo Librandi's method for sequential primes (Librandi's description in Italian).

Cf. A086801, A047949.
Numbers n such that 2n+k is prime: A005097 (k=1), this seq(k=3), A089038 (k=5), A105760 (k=7), A155722 (k=9), A101448 (k=11), A153081 (k=13), A089559 (k=15), A173059 (k=17), A153143 (k=19). - Jason Kimberley, Sep 07 2012
Numbers n such that 2n-k is prime: A006254 (k=1), A098090 (k=3), A089253 (k=5), A089192 (k=7), A097069 (k=9), A097338 (k=11), A097363 (k=13), A097480 (k=15), A098605 (k=17), A097932 (k=19).

Filtered([0..200], k-> IsPrime(2*k+3) ); # G. C. Greubel, May 21 2019

[n: n in [0..200]| IsPrime(2*n+3)]; // Vincenzo Librandi, Feb 23 2012

select(t -> isprime(2*t+3), [$0..1000]); # Robert Israel, Feb 19 2015

(Prime[Range[100]+1]-3)/2  (* Vladimir Joseph Stephan Orlovsky, Sep 08 2008, modified by G. C. Greubel, May 21 2019 *)
Select[Range[0,200],PrimeQ[2#+3]&] (* Harvey P. Dale, Jun 10 2014 *)

[k | k<-[0..99], isprime(2*k+3)] \\ for illustration

A067076(n) = (prime(n+1)-3)/2 \\ M. F. Hasler, Feb 14 2024

[n for n in (0..200) if is_prime(2*n+3) ] # G. C. Greubel, May 21 2019

a(n) = A006254(n) - 2 = A086801(n+1)/2. [Corrected by M. F. Hasler, Feb 14 2024]
a(n) = A089253(n) - 4. - Giovanni Teofilatto, Dec 14 2003
Conjecture: a(n) = A008507(n) + n - 1 = A005097(n) - 1 = A102781(n+1) - 1. - R. J. Mathar, Jul 07 2009
a(n) = A179893(n) - A000040(n). - Odimar Fabeny, Aug 24 2010

Offset changed from 0 to 1 in 2008: some formulas here and elsewhere may need to be corrected.

A034953 Triangular numbers (A000217) with prime indices.

Original entry on oeis.org

3, 6, 15, 28, 66, 91, 153, 190, 276, 435, 496, 703, 861, 946, 1128, 1431, 1770, 1891, 2278, 2556, 2701, 3160, 3486, 4005, 4753, 5151, 5356, 5778, 5995, 6441, 8128, 8646, 9453, 9730, 11175, 11476, 12403, 13366, 14028, 15051, 16110, 16471, 18336, 18721, 19503

Offset: 1

Patrick De Geest, Oct 15 1998

The following sequences (allowing offset of first term) all appear to have the same parity: A034953, triangular numbers with prime indices; A054269, length of period of continued fraction for sqrt(p), p prime; A082749, difference between the sum of next prime(n) natural numbers and the sum of next n primes; A006254, numbers n such that 2n-1 is prime; A067076, 2n+3 is a prime. - Jeremy Gardiner, Sep 10 2004
Given a rectangular prism with sides 1, p, p^2 for p = n-th prime (n > 1), the area of the six sides divided by the volume gives a remainder which is 4*a(n). - J. M. Bergot, Sep 12 2011
The infinite sum over the reciprocals is given by 2*A179119. - Wolfdieter Lang, Jul 10 2019

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Triangular Number.

Cf. A000217, A034954, A034955, A011756, A179119, A195678, A307868.

Cf. A054269, A067076, A082749.

a034953 n = a034953_list !! (n-1)
a034953_list = map a000217 a000040_list
-- Reinhard Zumkeller, Sep 23 2011

a:= n-> (p-> p*(p+1)/2)(ithprime(n)):
seq(a(n), n=1..65);  # Alois P. Heinz, Apr 20 2022

t[n_] := n(n + 1)/2; Table[t[Prime[n]], {n, 44}] (* Robert G. Wilson v, Aug 12 2004 *)
(#(# + 1))/2&/@Prime[Range[50]] (* Harvey P. Dale, Feb 27 2012 *)
With[{nn=200},Pick[Accumulate[Range[nn]],Table[If[PrimeQ[n],1,0],{n,nn}],1]] (* Harvey P. Dale, Mar 05 2023 *)

forprime(p=2,1e3,print1(binomial(p+1,2)", ")) \\ Charles R Greathouse IV, Jul 19 2011

apply(n->binomial(n+1,2),primes(100)) \\ Charles R Greathouse IV, Jun 04 2013

a(n) = A000217(A000040(n)). - Omar E. Pol, Jul 27 2009
a(n) = Sum_{k=1..prime(n)} k. - Wesley Ivan Hurt, Apr 27 2021
Product_{n>=1} (1 - 1/a(n)) = A307868. - Amiram Eldar, Nov 07 2022

A040040 Average of twin prime pairs (A014574), divided by 2. Equivalently, 2a(n)-1 and 2a(n)+1 are primes.

Original entry on oeis.org

2, 3, 6, 9, 15, 21, 30, 36, 51, 54, 69, 75, 90, 96, 99, 114, 120, 135, 141, 156, 174, 210, 216, 231, 261, 285, 300, 309, 321, 330, 405, 411, 414, 429, 441, 510, 516, 525, 531, 546, 576, 615, 639, 645, 651, 660, 714, 726, 741, 744, 804, 810, 834, 849, 861, 894

Offset: 1

N. J. A. Sloane

Intersection of A005097 and A006254. - Zak Seidov, Mar 18 2005
The only possible pairs for 2a(n)+-1 are prime/prime (this sequence), not prime/not prime (A104278), prime/notprime (A104279) and not prime/prime (A104280), ... this sequence + A104280 + A104279 + A104278 = the odd numbers.
These numbers are never k mod (2k+1) or (k+1) mod (2k+1) with 2k+1 < a(n). - Jon Perry, Sep 04 2012
Excluding the first term, all remaining terms have digital root 3, 6 or 9. - J. W. Helkenberg, Jul 24 2013
Positive numbers x such that the difference between x^2 and adjacent squares are prime (both x^2-(x-1)^2 and (x+1)^2-x^2 are prime). - Doug Bell, Aug 21 2015

T. D. Noe, Table of n, a(n) for n = 1..10001

Cf. A001359, A006512, A014574, A054735, A111046, A045753 (even terms halved), A002822 (terms divided by 3).
Cf. A221310.
Cf. A260689, A264526.

a040040 = flip div 2 . a014574  -- Reinhard Zumkeller, Nov 17 2015

P := select(isprime,[$1..1789]): map(p->(p+1)/2, select(p->member(p+2,P),P)); # Peter Luschny, Mar 03 2011

Select[Range[900], And @@ PrimeQ[{-1, 1} + 2# ] &] (* Ray Chandler, Oct 12 2005 *)

p=2; forprime(b=3, 1e4, if(b-p==2, print1((p+1)/2", ")); p=b) \\ Altug Alkan, Nov 10 2015

a(n) = A014574(n)/2 = A054735(n+1)/4 = A111046(n+1)/8.
For n > 1, a(n) = 3*A002822(n-1). - Jason Kimberley, Nov 06 2015
A260689(a(n),1) = A264526(a(n)) = 1. - Reinhard Zumkeller, Nov 17 2015
From Michael G. Kaarhus, Aug 19 2022: (Start)
a(n) = (A001359(n) + 1)/2.
a(n) = (A006512(n) - 1)/2.
For n > 1, a(n) = A167379(n-1) * 3/2. (End)

More terms from Cino Hilliard, Oct 21 2002
Title corrected by Daniel Forgues, Jun 01 2009
Edited by Daniel Forgues, Jun 21 2009
Comment corrected by Daniel Forgues, Jul 12 2009

A061286 Smallest integer for which the number of divisors is the n-th prime.

Original entry on oeis.org

2, 4, 16, 64, 1024, 4096, 65536, 262144, 4194304, 268435456, 1073741824, 68719476736, 1099511627776, 4398046511104, 70368744177664, 4503599627370496, 288230376151711744, 1152921504606846976

Offset: 1

Labos Elemer, May 22 2001

Seems to be the same as "Even numbers with prime number of divisors" - Jason Earls, Jul 04 2001

Except for the first term, smallest number == 1 (mod prime(n)) having n divisors (by Fermat's little theorem). - Amarnath Murthy and Meenakshi Srikanth (menakan_s(AT)yahoo.com), Jun 20 2003

Karl V. Keller, Jr., Table of n, a(n) for n = 1..460
Index to divisibility sequences

Cf. A000005, A005179, A003680, A061283, A061286, A006093, A005097, A006254, A119523, A196202.

Table[2^(p-1),{p,Table[Prime[n],{n,1,18}]}] (* Geoffrey Critzer, May 26 2013 *)

forstep(n=2,100000000,2,x=numdiv(n); if(isprime(x),print(n)))

a(n)=2^(prime(n)-1) \\ Charles R Greathouse IV, Apr 08 2012

from sympy import isprime, divisor_count as tau
[2] + [2**(2*n) for n in range(1, 33) if isprime(tau(2**(2*n)))] # Karl V. Keller, Jr., Jul 10 2020

a(n) = 2^(prime(n)-1) = 2^A006093(n).
a(n) = A005179(prime(n)). - R. J. Mathar, Aug 09 2019
Sum_{n>=1} 1/a(n) = A119523. - Amiram Eldar, Aug 11 2020

A089038 Nonnegative numbers k such that 2k+5 is prime.

Original entry on oeis.org

0, 1, 3, 4, 6, 7, 9, 12, 13, 16, 18, 19, 21, 24, 27, 28, 31, 33, 34, 37, 39, 42, 46, 48, 49, 51, 52, 54, 61, 63, 66, 67, 72, 73, 76, 79, 81, 84, 87, 88, 93, 94, 96, 97, 103, 109, 111, 112, 114, 117, 118, 123, 126, 129, 132, 133, 136, 138, 139, 144, 151, 153, 154, 156, 163

Offset: 1

Ray Chandler, Nov 29 2003

Vincenzo Librandi, Table of n, a(n) for n = 1..1000

Cf. A086304.
Numbers n such that 2n+k is prime: A005097 (k=1), A067076 (k=3), this seq (k=5), A105760 (k=7), A155722 (k=9), A101448 (k=11), A153081 (k=13), A089559 (k=15), A173059 (k=17), A153143 (k=19).
Numbers n such that 2n-k is prime: A006254 (k=1), A098090 (k=3), A089253 (k=5), A089192 (k=7), A097069 (k=9), A097338 (k=11), A097363 (k=13), A097480 (k=15), A098605 (k=17), A097932 (k=19).

Filtered([0..200], k-> IsPrime(2*k+5) ); # G. C. Greubel, May 21 2019

[n: n in [0..200]| IsPrime(2*n+5)]; // Vincenzo Librandi, Nov 17 2010

(Prime[Range[3,100]]-5)/2 (* Vladimir Joseph Stephan Orlovsky, Feb 08 2010 *)
Select[Range[0, 200], PrimeQ[2*# + 5] &] (* Vincenzo Librandi, Oct 16 2012 *)

is(n)=isprime(2*n+5) \\ Charles R Greathouse IV, Feb 16 2017

[n for n in (0..200) if is_prime(2*n+5) ] # G. C. Greubel, May 21 2019

a(n) = (A086304(n-1) + 1)/2, n > 1.

Removed wrong comment by Ralf Stephan, Nov 17 2010

Definition clarified by Zak Seidov, Jul 11 2014

cp's OEIS Frontend

A305901 Filter sequence for all such sequences b, for which b(A006254(k)) = constant for all k >= 3.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A081385 Incorrect version of A006254 (or A111333).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

PARI

Formula

A001248 Squares of primes.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

SageMath

Formula

A005097 (Odd primes - 1)/2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

A005382 Primes p such that 2p-1 is also prime.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Formula

A067076 Numbers k such that 2*k + 3 is a prime.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

A040040 Average of twin prime pairs (A014574), divided by 2. Equivalently, 2a(n)-1 and 2a(n)+1 are primes.