[go: up one dir, main page]

login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Search: a006254 -id:a006254
     Sort: relevance | references | number | modified | created      Format: long | short | data
Filter sequence for all such sequences b, for which b(A006254(k)) = constant for all k >= 3.
+20
6
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,2
COMMENTS
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).
LINKS
FORMULA
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.
PROG
(PARI)
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]));
CROSSREFS
Cf. also A305902.
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 14 2018
STATUS
approved
Incorrect version of A006254 (or A111333).
+20
0
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,2
COMMENTS
Name was: Numbers just more than half-prime.
REFERENCES
M. J. Halm, More Sequences, Mpossibilities 83, 2003.
FORMULA
a(n) = [p(n)/2 + 1/2].
EXAMPLE
a(11) = 16 because [(31)/2 + 1/2]=16.
PROG
(PARI) a(n) = ceil(prime(n)/2); \\ Michel Marcus, Aug 13 2013
CROSSREFS
Cf. A000040.
KEYWORD
dead
AUTHOR
Michael Joseph Halm, Apr 20 2003
STATUS
approved
Squares of primes.
+10
599
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,1
COMMENTS
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
LINKS
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
Brady Haran and Matt Parker, Squaring Primes, Numberphile video (2018).
Eric Weisstein's World of Mathematics, Prime Power.
FORMULA
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)
MAPLE
A001248:=n->ithprime(n)^2; seq(A001248(k), k=1..50); # Wesley Ivan Hurt, Oct 11 2013
MATHEMATICA
Prime[Range[30]]^2 (* Zak Seidov, Dec 07 2011 *)
Select[Range[40000], DivisorSigma[0, #] == 3 &] (* Carlos Eduardo Olivieri, Jun 01 2015 *)
PROG
(PARI) forprime(p=2, 1e3, print1(p^2", ")) \\ Charles R Greathouse IV, Jun 10 2011
(PARI) A001248(n)=prime(n)^2 \\ M. F. Hasler, Sep 16 2012
(Haskell)
a001248 n = a001248_list !! (n-1)
a001248_list = map (^ 2) a000040_list -- Reinhard Zumkeller, Sep 23 2011
(Magma) [p^2: p in PrimesUpTo(300)]; // Vincenzo Librandi, Mar 27 2014
(SageMath) [n^2 for n in prime_range(1, 301)] # G. C. Greubel, May 02 2024
(Python)
from sympy import prime
def A001248(n): return prime(n)**2 # Chai Wah Wu, Aug 09 2024
CROSSREFS
KEYWORD
nonn,easy
STATUS
approved
(Odd primes - 1)/2.
+10
142
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,2
COMMENTS
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
LINKS
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
FORMULA
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
MAPLE
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
MATHEMATICA
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 *)
PROG
(PARI) forprime(p=3, 1e4, print1(p>>1", ")) \\ Charles R Greathouse IV, Jun 16 2011
(Haskell)
a005097 = (`div` 2) . a065091 -- Reinhard Zumkeller, Jan 02 2013
(Magma) [n: n in [1..160] |IsPrime(2*n+1)]; // Vincenzo Librandi, Feb 16 2015
(Python)
from sympy import prime
def A005097(n): return prime(n+1)//2 # Chai Wah Wu, Jun 04 2022
CROSSREFS
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.
KEYWORD
nonn,easy
STATUS
approved
Primes p such that 2p-1 is also prime.
(Formerly M0849)
+10
122
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,1
COMMENTS
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
REFERENCES
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
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 & 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.
Marius Tărnăuceanu, Arithmetic progressions in finite groups, arXiv:2003.10060 [math.GR], 2020.
Wikipedia, Cunningham chain
FORMULA
a(n) = A129521(n) / A005383(n). - Reinhard Zumkeller, Apr 19 2007
a(n) = (A005383(n) + 1)/2. - Zak Seidov, Nov 04 2010
MAPLE
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);
MATHEMATICA
Select[Prime[Range[300]], PrimeQ[2#-1]&]
PROG
(Magma) [n: n in [0..1000] | IsPrime(n) and IsPrime(2*n-1)] // Vincenzo Librandi, Nov 18 2010
(PARI) select(p->isprime(2*p-1), primes(500)) \\ Charles R Greathouse IV, Apr 26 2012
(Haskell)
a005382 n = a005382_list !! (n-1)
a005382_list = filter
((== 1) . a010051 . (subtract 1) . (* 2)) a000040_list
-- Reinhard Zumkeller, Oct 03 2012
(PARI) forprime(n=2, 10^3, if(ispseudoprime(2*n-1), print1(n, ", "))) \\ Felix Fröhlich, Jun 15 2014
CROSSREFS
Cf. A010051, A000040, A053685 (subsequence), A006254.
KEYWORD
nonn,easy
STATUS
approved
Numbers k such that 2*k + 3 is a prime.
+10
65
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,3
COMMENTS
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
LINKS
Mutsumi Suzuki, Vincenzo Librandi's method for sequential primes (Librandi's description in Italian).
FORMULA
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
MAPLE
select(t -> isprime(2*t+3), [$0..1000]); # Robert Israel, Feb 19 2015
MATHEMATICA
(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 *)
PROG
(PARI) [k | k<-[0..99], isprime(2*k+3)] \\ for illustration
(PARI) A067076(n) = (prime(n+1)-3)/2 \\ M. F. Hasler, Feb 14 2024
(Magma)[n: n in [0..200]| IsPrime(2*n+3)]; // Vincenzo Librandi, Feb 23 2012
(Sage) [n for n in (0..200) if is_prime(2*n+3) ] # G. C. Greubel, May 21 2019
(GAP) Filtered([0..200], k-> IsPrime(2*k+3) ) # G. C. Greubel, May 21 2019
CROSSREFS
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).
KEYWORD
nonn,easy
AUTHOR
David Williams, Aug 17 2002
EXTENSIONS
Offset changed from 0 to 1 in 2008: some formulas here and elsewhere may need to be corrected.
STATUS
approved
Triangular numbers (A000217) with prime indices.
+10
50
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,1
COMMENTS
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
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Triangular Number.
FORMULA
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
MAPLE
a:= n-> (p-> p*(p+1)/2)(ithprime(n)):
seq(a(n), n=1..65); # Alois P. Heinz, Apr 20 2022
MATHEMATICA
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 *)
PROG
(PARI) forprime(p=2, 1e3, print1(binomial(p+1, 2)", ")) \\ Charles R Greathouse IV, Jul 19 2011
(PARI) apply(n->binomial(n+1, 2), primes(100)) \\ Charles R Greathouse IV, Jun 04 2013
(Haskell)
a034953 n = a034953_list !! (n-1)
a034953_list = map a000217 a000040_list
-- Reinhard Zumkeller, Sep 23 2011
KEYWORD
nonn,easy
AUTHOR
Patrick De Geest, Oct 15 1998
STATUS
approved
Average of twin prime pairs (A014574), divided by 2. Equivalently, 2*a(n)-1 and 2*a(n)+1 are primes.
+10
43
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,1
COMMENTS
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
FORMULA
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)
MAPLE
P := select(isprime, [$1..1789]): map(p->(p+1)/2, select(p->member(p+2, P), P)); # Peter Luschny, Mar 03 2011
MATHEMATICA
Select[Range[900], And @@ PrimeQ[{-1, 1} + 2# ] &] (* Ray Chandler, Oct 12 2005 *)
PROG
(PARI) p=2; forprime(b=3, 1e4, if(b-p==2, print1((p+1)/2", ")); p=b) \\ Altug Alkan, Nov 10 2015
(Haskell)
a040040 = flip div 2 . a014574 -- Reinhard Zumkeller, Nov 17 2015
CROSSREFS
Cf. A001359, A006512, A014574, A054735, A111046, A045753 (even terms halved), A002822 (terms divided by 3).
Cf. A221310.
KEYWORD
nonn,easy
EXTENSIONS
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
STATUS
approved
Nonnegative numbers k such that 2k+5 is prime.
+10
34
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,3
LINKS
FORMULA
a(n) = (A086304(n-1) + 1)/2, n > 1.
MATHEMATICA
(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 *)
PROG
(Magma) [n: n in [0..200]| IsPrime(2*n+5)] // Vincenzo Librandi, Nov 17 2010
(PARI) is(n)=isprime(2*n+5) \\ Charles R Greathouse IV, Feb 16 2017
(Sage) [n for n in (0..200) if is_prime(2*n+5) ] # G. C. Greubel, May 21 2019
(GAP) Filtered([0..200], k-> IsPrime(2*k+5) ) # G. C. Greubel, May 21 2019
CROSSREFS
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).
KEYWORD
easy,nonn
AUTHOR
Ray Chandler, Nov 29 2003
EXTENSIONS
Removed wrong comment by Ralf Stephan, Nov 17 2010
Definition clarified by Zak Seidov, Jul 11 2014
STATUS
approved
Numbers k such that 2k-3 is prime.
+10
34
3, 4, 5, 7, 8, 10, 11, 13, 16, 17, 20, 22, 23, 25, 28, 31, 32, 35, 37, 38, 41, 43, 46, 50, 52, 53, 55, 56, 58, 65, 67, 70, 71, 76, 77, 80, 83, 85, 88, 91, 92, 97, 98, 100, 101, 107, 113, 115, 116, 118, 121, 122, 127, 130, 133, 136, 137, 140, 142, 143, 148, 155, 157, 158
OFFSET
1,1
COMMENTS
Supersequence of A063908.
Left edge of the triangle in A065305. - Reinhard Zumkeller, Jan 30 2012
LINKS
FORMULA
Half of p + 3, where p is a prime greater than 2.
A122845(a(n), 3) = 3; a(n) = A113935(n+1)/2. - Reinhard Zumkeller, Sep 14 2006
MATHEMATICA
(Prime[Range[2, 100]]+3)/2 (* Vladimir Joseph Stephan Orlovsky, Feb 08 2010 *)
Select[Range[200], PrimeQ[2#-3]&] (* Harvey P. Dale, Mar 05 2022 *)
PROG
(Magma) [ n: n in [1..200] | IsPrime(2*n-3) ]; // Vincenzo Librandi, Dec 26 2010
(PARI) is(n)=isprime(2*n-3) \\ Charles R Greathouse IV, Feb 17 2017
(Sage) [n for n in (1..200) if is_prime(2*n-3) ] # G. C. Greubel, May 21 2019
(GAP) Filtered([1..200], k-> IsPrime(2*k-3) ) # G. C. Greubel, May 21 2019
CROSSREFS
Numbers n such that 2n+k is prime: A005097 (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), this sequence (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).
KEYWORD
easy,nonn
AUTHOR
Douglas Winston (douglas.winston(AT)srupc.com), Sep 14 2004
STATUS
approved

Search completed in 0.037 seconds