cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-6 of 6 results.

A075889 Primes in A075888, as they appear.

Original entry on oeis.org

3, 2, 5, 3, 7, 13, 5, 17, 13, 7, 23, 43, 17, 37, 43, 67, 23, 127, 137, 47, 103, 167, 127, 193, 223, 163, 167, 283, 103, 107, 257, 137, 293, 313, 487, 337, 563, 613, 617, 643, 647, 433, 773, 523, 283, 313, 1033, 347, 373, 757, 1187, 397, 1193, 797, 1277, 443
Offset: 1

Views

Author

Zak Seidov, Oct 17 2002

Keywords

Comments

Previous name was: Prime differences of successive primes squared divided by 24, (prime(n+1)^2-prime(n)^2)/24.
For n>=3, prime(n+1)^2-prime(n)^2 is always divisible by 24 and for many n's (prime(n+1)^2-prime(n)^2)/24 is prime.
Subsequence of primes of A075888. - Michel Marcus, Oct 03 2013

Examples

			a(1)=3 because 3 is the first prime value obtained for (prime(n+1)^2-prime(n)^2)/24 and n=5; next prime value is a(2)=2 and corresponds to n=6.
		

Programs

  • PARI
    lista(nn) = {pr = primes(nn); for (n = 4, #pr, if (isprime(q = (pr[n]^2 - pr[n-1]^2)/24), print1(q, ", ")););} \\ Michel Marcus, Oct 03 2013

Extensions

New name from Michel Marcus, Oct 05 2013

A026741 a(n) = n if n odd, n/2 if n even.

Original entry on oeis.org

0, 1, 1, 3, 2, 5, 3, 7, 4, 9, 5, 11, 6, 13, 7, 15, 8, 17, 9, 19, 10, 21, 11, 23, 12, 25, 13, 27, 14, 29, 15, 31, 16, 33, 17, 35, 18, 37, 19, 39, 20, 41, 21, 43, 22, 45, 23, 47, 24, 49, 25, 51, 26, 53, 27, 55, 28, 57, 29, 59, 30, 61, 31, 63, 32, 65, 33, 67, 34, 69, 35, 71, 36, 73, 37, 75, 38
Offset: 0

Views

Author

J. Carl Bellinger (carlb(AT)ctron.com)

Keywords

Comments

a(n) is the size of largest conjugacy class in D_2n, the dihedral group with 2n elements. - Sharon Sela (sharonsela(AT)hotmail.com), May 14 2002
a(n+1) is the composition length of the n-th symmetric power of the natural representation of a finite subgroup of SL(2,C) of type D_4 (quaternion group). - Paul Boddington, Oct 23 2003
For n > 1, a(n) is the greatest common divisor of all permutations of {0, 1, ..., n} treated as base n + 1 integers. - David Scambler, Nov 08 2006 (see the Mathematics Stack Exchange link below).
From Dimitrios Choussos (choussos(AT)yahoo.de), May 11 2009: (Start)
Sequence A075888 and the above sequence are fitting together.
First 2 entries of this sequence have to be taken out.
In some cases two three or more sequenced entries of this sequence have to be added together to get the next entry of A075888.
Example: Sequences begin with 1, 3, 2, 5, 3, 7, 4, 9 (4 + 9 = 13, the next entry in A075888).
But it works out well up to primes around 50000 (haven't tested higher ones).
As A075888 gives a very regular graph. There seems to be a regularity in the primes. (End)
Starting with 1 = triangle A115359 * [1, 2, 3, ...]. - Gary W. Adamson, Nov 27 2009
From Gary W. Adamson, Dec 11 2009: (Start)
Let M be an infinite lower triangular matrix with (1, 1, 1, 0, 0, 0, ...) in every column, shifted down twice. This sequence starting with 1 = M * (1, 2, 3, ...)
M =
1;
1, 0;
1, 1, 0;
0, 1, 0, 0;
0, 1, 1, 0, 0;
0, 0, 1, 0, 0, 0;
0, 0, 1, 1, 0, 0, 0;
...
A026741 = M * (1, 2, 3, ...); but A002487 = lim_{n->infinity} M^n, a left-shifted vector considered as a sequence. (End)
A particular case of sequence for which a(n+3) = (a(n+2) * a(n+1)+q)/a(n) for every n > n0. Here n0 = 1 and q = -1. - Richard Choulet, Mar 01 2010
For n >= 2, a(n+1) is the smallest m such that s_n(2*m*(n-1))/(n-1) is even, where s_b(c) is the sum of digits of c in base b. - Vladimir Shevelev, May 02 2011
A001477 and A005408 interleaved. - Omar E. Pol, Aug 22 2011
Numerator of n/((n-1)*(n-2)). - Michael B. Porter, Mar 18 2012
Number of odd terms of n-th row in the triangles A162610 and A209297. - Reinhard Zumkeller, Jan 19 2013
For n >= 3, a(n) is the periodic of integer of spiral length ratio of spiral that have (n-1) circle centers. See illustration in links. - Kival Ngaokrajang, Dec 28 2013
This is the sequence of Lehmer numbers u_n(sqrt(R), Q) with the parameters R = 4 and Q = 1. It is a strong divisibility sequence, that is, gcd(a(n), a(m)) = a(gcd(n, m)) for all natural numbers n and m. Cf. A005013 and A108412. - Peter Bala, Apr 18 2014
The sequence of convergents of the 2-periodic continued fraction [0; 1, -4, 1, -4, ...] = 1/(1 - 1/(4 - 1/(1 - 1/(4 - ...)))) = 2 begins [0/1, 1/1, 4/3, 3/2, 8/5, 5/3, 12/7, ...]. The present sequence is the sequence of denominators; the sequence of numerators of the continued fraction convergents [0, 1, 4, 3, 8, 5, 12, ...] is A022998, also a strong divisibility sequence. - Peter Bala, May 19 2014
For n >= 3, (a(n-2)/a(n))*Pi = vertex angle of a regular n-gon. See illustration in links. - Kival Ngaokrajang, Jul 17 2014
For n > 1, the numerator of the harmonic mean of the first n triangular numbers. - Colin Barker, Nov 13 2014
The difference sequence is a permutation of the integers. - Clark Kimberling, Apr 19 2015
From Timothy Hopper, Feb 26 2017: (Start)
Given the function a(n, p) = n/p if n mod p = 0, else n, then a possible formula is: a(n, p) = n*(1 + (p-1)*((n^(p-1)) mod p))/p, p prime, (n^(p-1)) mod p = 1, n not divisible by p. (Fermat's Little Theorem). Examples: p = 2; a(n), p = 3; A051176(n), p = 5; A060791(n), p = 7; A106608(n).
Conjecture: lcm(n, p) = p*a(n, p), gcd(n, p) = n/a(n, p). (End)
Let r(n) = (a(n+1) + 1)/a(n+1) if n mod 2 = 1, a(n+1)/(a(n+1) + 2) otherwise; then lim_{k->oo} 2^(k+2) * Product_{n=0..k} r(n)^(k-n) = Pi. - Dimitris Valianatos, Mar 22 2021
Number of integers k from 1 to n such that gcd(n,k) is odd. - Amiram Eldar, May 18 2025

Examples

			G.f. = x + x^2 + 3*x^3 + 2*x^4 + 5*x^5 + 3*x^6 + 7*x^7 + 4*x^8 + ...
		

References

  • David Wells, Prime Numbers: The Most Mysterious Figures in Math. Hoboken, New Jersey: John Wiley & Sons (2005), p. 53.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, 2nd Ed. Penguin (1997), p. 79.

Crossrefs

Signed version is in A030640. Partial sums give A001318.
Cf. A051176, A060819, A060791, A060789 for n / gcd(n, k) with k = 3..6. See also A106608 thru A106612 (k = 7 thru 11), A051724 (k = 12), A106614 thru A106621 (k = 13 thru 20).
Cf. A013942.
Cf. A227042 (first column). Cf. A005013 and A108412.

Programs

  • Haskell
    import Data.List (transpose)
    a026741 n = a026741_list !! n
    a026741_list = concat $ transpose [[0..], [1,3..]]
    -- Reinhard Zumkeller, Dec 12 2011
    
  • Magma
    [2*n/(3+(-1)^n): n in [0..70]]; // Vincenzo Librandi, Aug 14 2011
    
  • Maple
    A026741 := proc(n) if type(n,'odd') then n; else n/2; end if; end proc: seq(A026741(n), n=0..76); # R. J. Mathar, Jan 22 2011
  • Mathematica
    Numerator[Abs[Table[Det[DiagonalMatrix[Table[1/i^2 - 1, {i, 1, n - 1}]] + 1], {n, 20}]]] (* Alexander Adamchuk, Jun 02 2006 *)
    halfMax = 40; Riffle[Range[0, halfMax], Range[1, 2halfMax + 1, 2]] (* Harvey P. Dale, Mar 27 2011 *)
    a[ n_] := Numerator[n / 2]; (* Michael Somos, Jan 20 2017 *)
    Array[If[EvenQ[#],#/2,#]&,80,0] (* Harvey P. Dale, Jul 08 2023 *)
  • PARI
    a(n) = numerator(n/2) \\ Rick L. Shepherd, Sep 12 2007
    
  • Python
    def A026741(n): return n if n % 2 else n//2 # Chai Wah Wu, Apr 02 2021
  • Sage
    [lcm(n, 2) / 2 for n in range(77)] # Zerinvary Lajos, Jun 07 2009
    

Formula

G.f.: x*(1 + x + x^2)/(1-x^2)^2. - Len Smiley, Apr 30 2001
a(n) = 2*a(n-2) - a*(n-4) for n >= 4.
a(n) = n * 2^((n mod 2) - 1). - Reinhard Zumkeller, Oct 16 2001
a(n) = 2*n/(3 + (-1)^n). - Benoit Cloitre, Mar 24 2002
Multiplicative with a(2^e) = 2^(e-1) and a(p^e) = p^e, p > 2. - Vladeta Jovovic, Apr 05 2002
a(n) = n / gcd(n, 2). a(n)/A045896(n) = n/((n+1)*(n+2)).
For n > 0, a(n) = denominator of Sum_{i=1..n-1} 2/(i*(i+1)), numerator=A022998. - Reinhard Zumkeller, Apr 21 2012, Jul 25 2002 [thanks to Phil Carmody who noticed an error]
For n > 1, a(n) = GCD of the n-th and (n-1)-th triangular numbers (A000217). - Ross La Haye, Sep 13 2003
Euler transform of finite sequence [1, 2, -1]. - Michael Somos, Jun 15 2005
G.f.: x * (1 - x^3) / ((1 - x) * (1 - x^2)^2) = Sum_{k>0} k * (x^k - x^(2*k)). - Michael Somos, Jun 15 2005
a(n+3) + a(n+2) = 3 + a(n+1) + a(n). a(n+3) * a(n) = - 1 + a(n+2) * a(n+1). a(n) = -a(-n) for all n in Z. - Michael Somos, Jun 15 2005
For n > 1, a(n) is the numerator of the average of 1, 2, ..., n - 1; i.e., numerator of A000217(n-1)/(n-1), with corresponding denominators [1, 2, 1, 2, ...] (A000034). - Rick L. Shepherd, Jun 05 2006
Equals A126988 * (1, -1, 0, 0, 0, ...). - Gary W. Adamson, Apr 17 2007
For n >= 1, a(n) = gcd(n,A000217(n)). - Rick L. Shepherd, Sep 12 2007
a(n) = numerator(n/(2*n-2)) for n >= 2; A022998(n-1) = denominator(n/(2*n-2)) for n >= 2. - Johannes W. Meijer, Jun 18 2009
a(n) = A167192(n+2, 2). - Reinhard Zumkeller, Oct 30 2009
a(n) = A106619(n) * A109012(n). - Paul Curtz, Apr 04 2011
From R. J. Mathar, Apr 18 2011: (Start)
a(n) = A109043(n)/2.
Dirichlet g.f.: zeta(s-1)*(1 - 1/2^s). (End)
a(n) = A001318(n) - A001318(n-1) for n > 0. - Jonathan Sondow, Jan 28 2013
a((2*n+1)*2^p - 1) = 2^p - 1 + n*A151821(p+1), p >= 0 and n >= 0. - Johannes W. Meijer, Feb 03 2013
a(n+1) = denominator(H(n, 1)), n >= 0, with H(n, 1) = 2*n/(n+1) the harmonic mean of n and 1. a(n+1) = A227042(n, 1). See the formula a(n) = n/gcd(n, 2) given above. - Wolfdieter Lang, Jul 04 2013
a(n) = numerator(n/2). - Wesley Ivan Hurt, Oct 02 2013
a(n) = numerator(1 - 2/(n+2)), n >= 0; a(n) = denominator(1 - 2/n), n >= 1. - Kival Ngaokrajang, Jul 17 2014
a(n) = Sum_{i = floor(n/2)..floor((n+1)/2)} i. - Wesley Ivan Hurt, Apr 27 2016
Euler transform of length 3 sequence [1, 2, -1]. - Michael Somos, Jan 20 2017
G.f.: x / (1 - x / (1 - 2*x / (1 + 7*x / (2 - 9*x / (7 - 4*x / (3 - 7*x / (2 + 3*x))))))). - Michael Somos, Jan 20 2017
From Peter Bala, Mar 24 2019: (Start)
a(n) = Sum_{d|n, n/d odd} phi(d), where phi(n) is the Euler totient function A000010.
O.g.f.: Sum_{n >= 1} phi(n)*x^n/(1 - x^(2*n)). (End)
a(n) = A256095(2*n,n). - Alois P. Heinz, Jan 21 2020
E.g.f.: x*(2*cosh(x) + sinh(x))/2. - Stefano Spezia, Apr 28 2023
From Ctibor O. Zizka, Oct 05 2023: (Start)
For k >= 0, a(k) = gcd(k + 1, k*(k + 1)/2).
If (k mod 4) = 0 or 2 then a(k) = (k + 1).
If (k mod 4) = 1 or 3 then a(k) = (k + 1)/2. (End)
Sum_{n=1..oo} 1/a(n)^2 = 7*Pi^2/24. - Stefano Spezia, Dec 02 2023
a(n)*a(n+1) = A000217(n). - Rémy Sigrist, Mar 19 2025

Extensions

Better description from Jud McCranie
Edited by Ralf Stephan, Jun 04 2003

A024702 a(n) = (prime(n)^2 - 1)/24.

Original entry on oeis.org

1, 2, 5, 7, 12, 15, 22, 35, 40, 57, 70, 77, 92, 117, 145, 155, 187, 210, 222, 260, 287, 330, 392, 425, 442, 477, 495, 532, 672, 715, 782, 805, 925, 950, 1027, 1107, 1162, 1247, 1335, 1365, 1520, 1552, 1617, 1650, 1855, 2072, 2147, 2185, 2262, 2380, 2420, 2625, 2752, 2882, 3015
Offset: 3

Views

Author

Clark Kimberling, Dec 11 1999

Keywords

Comments

Note that p^2 - 1 is always divisible by 24 since p == 1 or 2 (mod 3), so p^2 == 1 (mod 3) and p == 1, 3, 5, or 7 (mod 8) so p^2 == 1 (mod 8). - Michael B. Porter, Sep 02 2016
For n > 3 and m > 1, a(n) = A000330(m)/(2*m + 1), where 2*m + 1 = prime(n). For example, for m = 8, 2*m + 1 = 17 = prime(7), A000330(8) = 204, 204/17 = 12 = a(7). - Richard R. Forberg, Aug 20 2013
For primes => 5, a(n) == 0 or 2 (mod 5). - Richard R. Forberg, Aug 28 2013
The only primes in this sequence are 2, 5 and 7 (checked up to n = 10^7). The set of prime factors, however, appears to include all primes. - Richard R. Forberg, Feb 28 2015
Subsequence of generalized pentagonal numbers (cf. A001318): a(n) = k_n*(3*k_n - 1)/2, for k_n in {1, -1, 2, -2, 3, -3, 4, 5, -5, -6, 7, -7, 8, 9, 10, -10, ...} = A024699(n-2)*((A000040(n) mod 6) - 3)/2, n >= 3. - Daniel Forgues, Aug 02 2016
The only primes in this sequence are indeed 2, 5 and 7. For a prime p >= 5, if both p + 1 and p - 1 contains a prime factor > 3, then (p^2 - 1)/24 = (p + 1)*(p - 1)/24 contains at least 2 prime factors, so at least one of p + 1 and p - 1 is 3-smooth. Let's call it s. Also, If (p^2 - 1)/24 is a prime, then A001222(p^2-1) = 5. Since A001222(p+1) and A001222(p-1) are both at least 2, A001222(s) <= 5 - 2 = 3. From these we can see the only possible cases are p = 7, 11 and 13. - Jianing Song, Dec 28 2018

Examples

			For n = 6, the 6th prime is 13, so a(6) = (13^2 - 1)/24 = 168/24 = 7.
		

Crossrefs

Subsequence of generalized pentagonal numbers A001318.
Cf. A075888.

Programs

Formula

a(n) = (A000040(n)^2 - 1)/24 = (A001248(n) - 1)/24. - Omar E. Pol, Dec 07 2011
a(n) = A005097(n-1)*A006254(n-1)/6. - Bruno Berselli, Dec 08 2011
a(n) = A084920(n)/24. - R. J. Mathar, Aug 23 2013
a(n) = A127922(n)/A000040(n) for n >= 3. - César Aguilera, Nov 01 2019

A283532 Primes p such that (q^2 - p^2) / 24 is prime, where q is the next prime after p.

Original entry on oeis.org

7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 67, 83, 101, 109, 127, 131, 137, 251, 271, 281, 307, 331, 379, 383, 443, 487, 499, 563, 617, 641, 769, 821, 877, 937, 971, 1009, 1123, 1223, 1231, 1283, 1291, 1297, 1543, 1567, 1697, 1877, 2063, 2081, 2237, 2269, 2371, 2381, 2383, 2389, 2551, 2657, 2659, 2801, 2851, 2857
Offset: 1

Views

Author

Thomas Ordowski and Altug Alkan, Mar 10 2017

Keywords

Comments

This sequence is union of primes of the form:
6t-1 such that 6t+1 and t are both prime,
6t-1 such that 6t+5 and 3t+1 are both prime and 6t+1 is composite,
6t+1 such that 6t+5 and 2t+1 are both prime,
6t+1 such that 6t+7 and 3t+2 are both prime and 6t+5 is composite.

Examples

			7 is a term since 11 is the next prime and (11^2 - 7^2)/24 = 3 is prime.
		

Crossrefs

A060213 is a subsequence.
Cf. A075888.

Programs

  • Maple
    N:= 10000: # to get all terms <= N
    Primes:= select(isprime, [seq(i,i=3..N,2)]):
    f:= proc(p,q)
      local r;
      r:= (q^2-p^2)/24;
      if r::integer and isprime(r) then p fi
    end proc:
    seq(f(Primes[i],Primes[i+1]),i=1..nops(Primes)-1); # Robert Israel, Mar 10 2017
  • Mathematica
    Select[Prime@ Range@ 415, PrimeQ[(NextPrime[#]^2 - #^2)/24] &] (* Michael De Vlieger, Mar 13 2017 *)
  • PARI
    is(n) = n>3 && isprime(n) && isprime((nextprime(n+1)^2-n^2)/24);

A282445 For n>=5, a(n) is the smallest m>=3 such that odd part of ((prime(n)^2 - prime(m)^2)/3) is prime, or a(n)=0 if there is no such m

Original entry on oeis.org

4, 3, 3, 3, 4, 3, 4, 3, 4, 3, 7, 3, 12, 6, 8, 4, 13, 7, 8, 4, 11, 3, 20, 5, 6, 22, 11, 23, 13, 16, 14, 9, 10, 10, 24, 29, 6, 40, 31, 0, 3, 4, 40, 11, 32, 45, 13, 7, 30, 3, 53, 20, 6, 30, 35, 27, 54, 26, 0, 63, 46, 57, 16, 67, 67, 38, 0, 39, 52, 5, 61, 75, 3
Offset: 5

Views

Author

Vladimir Shevelev, Feb 15 2017

Keywords

Comments

a(n) = 0 for n: 44, 63, 71, 80, 89, 95, 97, 108, 118, 122, 132, 141, 150, etc. Robert G. Wilson v, Feb 15 2017

Examples

			Let n=9, prime(9)=23. If m=3, then odd part of (23^2 - 5^2)/24 is 21, while if m=4, then odd part of (23^2 - 7^2)/24 is 5 which is prime. So a(9)=4.
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Block[{m = 3, p = Prime[n]^2}, While[q = (p - Prime[m]^2)/3; m < n && ! PrimeQ[q/2^IntegerExponent[q, 2]], m++]; If[m < n, m, 0]]; Array[f, 73, 5] (* Robert G. Wilson v, Feb 15 2017 *)

Extensions

More terms from Peter J. C. Moses, Feb 15 2017

A286679 a(n) = (2*prime(n)^2 + 1)/3.

Original entry on oeis.org

17, 33, 81, 113, 193, 241, 353, 561, 641, 913, 1121, 1233, 1473, 1873, 2321, 2481, 2993, 3361, 3553, 4161, 4593, 5281, 6273, 6801, 7073, 7633, 7921, 8513, 10753, 11441, 12513, 12881, 14801, 15201, 16433, 17713, 18593, 19953, 21361, 21841
Offset: 3

Views

Author

Dimitris Valianatos, May 12 2017

Keywords

Comments

For primes p other than 3, p == 1 or 2 (mod 3) and p^2 == 1 (mod 3). Thus 2*p^2 + 1 is a multiple of 3.

Crossrefs

Programs

  • Magma
    [(2*NthPrime(n)^2+1)/3: n in [3..50]]; // Vincenzo Librandi, May 15 2017
  • Mathematica
    (2Prime[Range[3, 50]]^2 + 1)/3 (* Alonso del Arte, May 12 2017 *)
  • PARI
    {
    forprime(n=5,300,
             print1((2*n^2+1)/3", ")
            )
    }
    

Formula

Product_{n >= 3} (3*a(n) + 1) / (3*a(n) - 1) = (26/25) * (50/49) * (122/121) * ... = 54/(5*Pi^2) = 1.0942687833372479315938982026650585002 (constant).
a(3) = 17; a(n + 1) = a(n) + 16 * A075888(n-2) for n > 3.
Numbers of the form 16k + 1 for some k. In particular, k belongs to A001318, excluding those for which sqrt(24 * A001318(k) + 1) are composites.
Showing 1-6 of 6 results.