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.

A007528 Primes of the form 6k-1.

Original entry on oeis.org

5, 11, 17, 23, 29, 41, 47, 53, 59, 71, 83, 89, 101, 107, 113, 131, 137, 149, 167, 173, 179, 191, 197, 227, 233, 239, 251, 257, 263, 269, 281, 293, 311, 317, 347, 353, 359, 383, 389, 401, 419, 431, 443, 449, 461, 467, 479, 491, 503, 509, 521, 557, 563, 569, 587
Offset: 1

Views

Author

Keywords

Comments

For values of k see A024898.
Also primes p such that p^q - 2 is not prime where q is an odd prime. These numbers cannot be prime because the binomial p^q = (6k-1)^q expands to 6h-1 some h. Then p^q-2 = 6h-1-2 is divisible by 3 thus not prime. - Cino Hilliard, Nov 12 2008
a(n) = A211890(3,n-1) for n <= 4. - Reinhard Zumkeller, Jul 13 2012
There exists a polygonal number P_s(3) = 3s - 3 = a(n) + 1. These are the only primes p with P_s(k) = p + 1, s >= 3, k >= 3, since P_s(k) - 1 is composite for k > 3. - Ralf Steiner, May 17 2018
From Bernard Schott, Feb 14 2019: (Start)
A theorem due to Andrzej Mąkowski: every integer greater than 161 is the sum of distinct primes of the form 6k-1. Examples: 162 = 5 + 11 + 17 + 23 + 47 + 59; 163 = 17 + 23 + 29 + 41 + 53. (See Sierpiński and David Wells.)
{2,3} Union A002476 Union {this sequence} = A000040.
Except for 2 and 3, all Sophie Germain primes are of the form 6k-1.
Except for 3, all the lesser of twin primes are also of the form 6k-1.
Dirichlet's theorem on arithmetic progressions states that this sequence is infinite. (End)
For all elements of this sequence p=6*k-1, there are no (x,y) positive integers such that k=6*x*y-x+y. - Pedro Caceres, Apr 06 2019

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 870.
  • A. Mąkowski, Partitions into unequal primes, Bull. Acad. Polon. Sci. Sér. Sci. Math. Astr. Phys. 8 (1960), 125-126.
  • Wacław Sierpiński, Elementary Theory of Numbers, p. 144, Warsaw, 1964.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers, Penguin Books, Revised edition, 1997, p. 127.

Crossrefs

Intersection of A016969 and A000040.
Prime sequences A# (k,r) of the form k*n+r with 0 <= r <= k-1 (i.e., primes == r (mod k), or primes p with p mod k = r) and gcd(r,k)=1: A000040 (1,0), A065091 (2,1), A002476 (3,1), A003627 (3,2), A002144 (4,1), A002145 (4,3), A030430 (5,1), A045380 (5,2), A030431 (5,3), A030433 (5,4), A002476 (6,1), this sequence (6,5), A140444 (7,1), A045392 (7,2), A045437 (7,3), A045471 (7,4), A045458 (7,5), A045473 (7,6), A007519 (8,1), A007520 (8,3), A007521 (8,5), A007522 (8,7), A061237 (9,1), A061238 (9,2), A061239 (9,4), A061240 (9,5), A061241 (9,7), A061242 (9,8), A030430 (10,1), A030431 (10,3), A030432 (10,7), A030433 (10,9), A141849 (11,1), A090187 (11,2), A141850 (11,3), A141851 (11,4), A141852 (11,5), A141853 (11,6), A141854 (11,7), A141855 (11,8), A141856 (11,9), A141857 (11,10), A068228 (12,1), A040117 (12,5), A068229 (12,7), A068231 (12,11).
Cf. A034694 (smallest prime == 1 (mod n)).
Cf. A038700 (smallest prime == n-1 (mod n)).
Cf. A038026 (largest possible value of smallest prime == r (mod n)).
Cf. A001359 (lesser of twin primes), A005384 (Sophie Germain primes).

Programs

  • GAP
    Filtered(List([1..100],n->6*n-1),IsPrime); # Muniru A Asiru, May 19 2018
  • Haskell
    a007528 n = a007528_list !! (n-1)
    a007528_list = [x | k <- [0..], let x = 6 * k + 5, a010051' x == 1]
    -- Reinhard Zumkeller, Jul 13 2012
    
  • Maple
    select(isprime,[seq(6*n-1,n=1..100)]); # Muniru A Asiru, May 19 2018
  • Mathematica
    Select[6 Range[100]-1,PrimeQ]  (* Harvey P. Dale, Feb 14 2011 *)
  • PARI
    forprime(p=2, 1e3, if(p%6==5, print1(p, ", "))) \\ Charles R Greathouse IV, Jul 15 2011
    
  • PARI
    forprimestep(p=5,1000,6, print1(p", ")) \\ Charles R Greathouse IV, Mar 03 2025
    

Formula

A003627 \ {2}. - R. J. Mathar, Oct 28 2008
Conjecture: Product_{n >= 1} ((a(n) - 1) / (a(n) + 1)) * ((A002476(n) + 1) / (A002476(n) - 1)) = 3/4. - Dimitris Valianatos, Feb 11 2020
From Vaclav Kotesovec, May 02 2020: (Start)
Product_{k>=1} (1 - 1/a(k)^2) = 9*A175646/Pi^2 = 1/1.060548293.... =4/(3*A333240).
Product_{k>=1} (1 + 1/a(k)^2) = A334482.
Product_{k>=1} (1 - 1/a(k)^3) = A334480.
Product_{k>=1} (1 + 1/a(k)^3) = A334479. (End)
Legendre symbol (-3, a(n)) = -1 and (-3, A002476(n)) = +1, for n >= 1. For prime 3 one sets (-3, 3) = 0. - Wolfdieter Lang, Mar 03 2021

A045469 Primes congruent to {1, 4} mod 7.

Original entry on oeis.org

11, 29, 43, 53, 67, 71, 109, 113, 127, 137, 151, 179, 193, 197, 211, 239, 263, 277, 281, 337, 347, 379, 389, 421, 431, 449, 463, 487, 491, 547, 557, 571, 599, 613, 617, 631, 641, 659, 673, 683, 701, 739, 743
Offset: 1

Views

Author

Keywords

Comments

The same as A045461 except for the initial 7 (the only prime congruent to 0 mod 7). - Alonso del Arte, Feb 18 2012

Crossrefs

Programs

  • Magma
    [ p: p in PrimesUpTo(1000) | p mod 7 in {1, 4} ]; // Vincenzo Librandi, Aug 06 2012
  • Mathematica
    Select[Prime[Range[120]], MemberQ[{1, 4}, Mod[#, 7]] &] (* Alonso del Arte, based on Orlovsky's program for A042992, Feb 18 2012 *)

A141849 Primes congruent to 1 mod 11.

Original entry on oeis.org

23, 67, 89, 199, 331, 353, 397, 419, 463, 617, 661, 683, 727, 859, 881, 947, 991, 1013, 1123, 1277, 1321, 1409, 1453, 1607, 1783, 1871, 2003, 2069, 2113, 2179, 2267, 2311, 2333, 2377, 2399, 2531, 2663, 2707, 2729, 2861, 2927, 2971, 3037, 3169, 3191, 3257
Offset: 1

Views

Author

N. J. A. Sloane, Jul 11 2008

Keywords

Comments

Conjecture: Also primes p such that ((x+1)^11-1)/x has 10 distinct irreducible factors of degree 1 over GF(p). - Federico Provvedi, Apr 17 2018
Primes congruent to 1 mod 22. - Chai Wah Wu, Apr 28 2025

Crossrefs

Prime sequences A# (k,r) of the form k*n+r with 0 <= r <= k-1 (i.e., primes == r (mod k), or primes p with p mod k = r) and gcd(r,k)=1: A000040 (1,0), A065091 (2,1), A002476 (3,1), A003627 (3,2), A002144 (4,1), A002145 (4,3), A030430 (5,1), A045380 (5,2), A030431 (5,3), A030433 (5,4), A002476 (6,1), A007528 (6,5), A140444 (7,1), A045392 (7,2), A045437 (7,3), A045471 (7,4), A045458 (7,5), A045473 (7,6), A007519 (8,1), A007520 (8,3), A007521 (8,5), A007522 (8,7), A061237 (9,1), A061238 (9,2), A061239 (9,4), A061240 (9,5), A061241 (9,7), A061242 (9,8), A030430 (10,1), A030431 (10,3), A030432 (10,7), A030433 (10,9), this sequence (11,1), A090187 (11,2), A141850 (11,3), A141851 (11,4), A141852 (11,5), A141853 (11,6), A141854 (11,7), A141855 (11,8), A141856 (11,9), A141857 (11,10), A068228 (12,1), A040117 (12,5), A068229 (12,7), A068231 (12,11).
Cf. A034694 (smallest prime == 1 (mod n)).
Cf. A038700 (smallest prime == n-1 (mod n)).
Cf. A038026 (largest possible value of smallest prime == r (mod n)).

Programs

Formula

a(n) ~ 10n log n. - Charles R Greathouse IV, Jul 02 2016

A319101 Number of solutions to x^7 == 1 (mod n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 7
Offset: 1

Views

Author

Jianing Song, Sep 10 2018

Keywords

Comments

All terms are powers of 7. Those n such that a(n) > 1 are in A066502.

Examples

			Solutions to x^7 == 1 (mod 29): x == 1, 7, 16, 20, 23, 24, 25 (mod 29).
Solutions to x^7 == 1 (mod 43): x == 1, 4, 11, 16, 21, 35, 41 (mod 43).
Solutions to x^7 == 1 (mod 49): x == 1, 8, 15, 22, 29, 36, 43 (mod 49) (x == 1 (mod 7)).
		

Crossrefs

Number of solutions to x^k == 1 (mod n): A060594 (k=2), A060839 (k=3), A073103 (k=4), A319099 (k=5), A319100 (k=6), this sequence (k=7), A247257 (k=8).
Mobius transform gives A307382.

Programs

  • Mathematica
    f[p_, e_] := If[Mod[p, 7] == 1, 7, 1]; f[7, 1] = 1; f[7, e_] := 7; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Aug 10 2023 *)
  • PARI
    a(n)=my(Z=znstar(n)[2]); prod(i=1, #Z, gcd(7, Z[i]))

Formula

Multiplicative with a(7) = 1, a(7^e) = 7 if e >= 2; for other primes p, a(p^e) = 7 if p == 1 (mod 7), a(p^e) = 1 otherwise.
If the multiplicative group of integers modulo n is isomorphic to C_{k_1} x C_{k_2} x ... x C_{k_m}, where k_i divides k_j for i < j; then a(n) = Product_{i=1..m} gcd(7, k_i).
a(n) = A000010(n)/A293484(n). - Jianing Song, Nov 10 2019

A263770 Smallest prime q such that (prime(n)^2 + q*prime(n))/(prime(n) + 1) is an integer.

Original entry on oeis.org

7, 5, 7, 17, 13, 29, 19, 41, 73, 31, 97, 191, 43, 89, 97, 109, 61, 311, 137, 73, 149, 241, 337, 181, 197, 103, 313, 109, 331, 229, 257, 397, 139, 281, 151, 457, 317, 821, 337, 349, 181, 547, 193, 389, 199, 401, 1061, 449, 229, 461, 937, 241, 727, 757, 1033, 1321, 271, 1361, 557
Offset: 1

Views

Author

Juri-Stepan Gerasimov, Oct 25 2015

Keywords

Comments

Least prime q such that q == 1 (mod prime(n) + 1).

Crossrefs

Programs

  • Mathematica
    Table[q = 2; While[! IntegerQ[(Prime[n]^2 + q Prime@ n)/(Prime@ n + 1)], q = NextPrime@ q]; q, {n, 59}] (* Michael De Vlieger, Oct 26 2015 *)
  • PARI
    a(n) = {p = prime(n); q = 2; while ((p^2 + p*q) % (p + 1), q = nextprime(q+1)); q;} \\ Michel Marcus, Oct 26 2015

Formula

5 is in this sequence because (prime(2)^2 + 5*prime(2))/(prime(2) + 1) = 6 and 5 is prime.

A291747 Nonprimes of the form 7*k + 1.

Original entry on oeis.org

1, 8, 15, 22, 36, 50, 57, 64, 78, 85, 92, 99, 106, 120, 134, 141, 148, 155, 162, 169, 176, 183, 190, 204, 218, 225, 232, 246, 253, 260, 267, 274, 288, 295, 302, 309, 316, 323, 330, 344, 351, 358, 365, 372, 386, 393, 400, 407, 414, 428, 435, 442, 456, 470, 477, 484, 498
Offset: 1

Views

Author

Vincenzo Librandi, Aug 31 2017

Keywords

Comments

Subsequence of A018252.

Crossrefs

Programs

  • Magma
    [n: n in [1..500 by 7] | not IsPrime(n)];
  • Mathematica
    DeleteCases[7 Range[0, 200] + 1, _?PrimeQ]
Showing 1-6 of 6 results.