A001479 -id:A001479 - cp's OEIS frontend

A272200 Bisection of primes congruent to 1 modulo 3 (A002476), depending on the corresponding A001479 entry being congruent to 1 modulo 3 or not. Here the first case.

13, 19, 43, 61, 97, 103, 109, 127, 157, 163, 181, 193, 241, 277, 283, 331, 349, 373, 379, 409, 433, 463, 487, 499, 523, 601, 607, 619, 631, 661, 673, 691, 727, 733, 757, 769, 787, 811, 859, 883, 937, 967, 991

Offset: 1

Wolfdieter Lang, Apr 28 2016

The other primes congruent to 1 modulo 3 are given in A272201.
Each prime == 1 (mod 3) has a unique representation A002476(m) = A(m)^2 + 3*B(m)^2 with positive A(m) = A001479(m+1) and B(m) = A001480(m+1), m >= 1 (see also A001479). The present sequence gives such primes corresponding to A(m+1) == 1 (mod 3). The ones corresponding to A(m+1) not == 1 (mod 3) (the complement) are given in A272201.
This bisection of the primes from A002476 is needed in the formula for the coefficients of the q-expansion (q = exp(2*Pi*i*z), Im(z) > 0) of the modular cusp form (eta(6*z))^4|A000727%20which%20gives%20the%20coefficients%20of%20the%20q-expansion%20of%20F(q)%20=%20Eta64(q%5E(1/6))/q%5E(1/6)%20=%20(Product">{z=z(q)} = Eta64(q) with Dedekind's eta function. See A000727 which gives the coefficients of the q-expansion of F(q) = Eta64(q^(1/6))/q^(1/6) = (Product{m>=0} (1 - q^m))^4. The coefficients F(q) = Sum_{n>=0} f(6*n+1)*q^n are given in the Finch link on p. 5, using multiplicativity. For primes congruent to 1 modulo 6 the formula involves x_p and y_p which are the present A and B for prime p == 1 (mod 3).
See also the p-defects of the elliptic curve y^2 = x^3 + 1, related to (eta(6*z))^4, given in A272198 with another (equivalent) way to find the coefficients of the Eta64(q) expansion, hence those of F(q).

Robert Israel, Table of n, a(n) for n = 1..10000
Steven R. Finch, Powers of Euler's q-Series, arXiv:math/0701251 [math.NT], 2007.

Cf. A000727, A001479, A002476, A001480, A272198, A272201 (complement relative to A002476).

filter:= proc(n) local S,x,y;
    if not isprime(n) then return false fi;
    S:= remove(hastype,[isolve(x^2+3*y^2=n)],negative);
    subs(S[1],x) mod 3 = 1
end proc:
select(filter, [seq(i,i=7..1000,6)]); # Robert Israel, Apr 29 2019

filterQ[n_] := Module[{S, x, y}, If[!PrimeQ[n], Return[False]]; S = Solve[x > 0 && y > 0 && x^2 + 3 y^2 == n, Integers]; Mod[x /. S[[1]], 3] == 1];
Select[Range[7, 1000, 6], filterQ] (* Jean-François Alcover, Apr 21 2020, after Robert Israel *)

This sequence collects the 1 (mod 3) primes p(m) = A002476(m) = A(m)^2 + 3*B(m)^2 with positive A(m) == 1 (mod 3), for m >= 1. A(m) = A001479(m+1).

A272201 Bisection of primes congruent to 1 modulo 3 (A002476), depending on the corresponding A001479 entry being congruent to 1 modulo 3 or not. Here the second case.

Original entry on oeis.org

7, 31, 37, 67, 73, 79, 139, 151, 199, 211, 223, 229, 271, 307, 313, 337, 367, 397, 421, 439, 457, 541, 547, 571, 577, 613, 643, 709, 739, 751, 823, 829, 853, 877, 907, 919, 997

Offset: 1

Wolfdieter Lang, Apr 28 2016

The other primes congruent to 1 modulo 3 are given in A272200, where also more details are given.

Each prime == 1 (mod 3) has a unique representation A002476(m) = A(m)^2 + 3*B(m)^2 with positive A(m) = A001479(m+1) and B(m) = A001480(m+1), m >= 1 (see also A001479). The present sequence gives these primes corresponding to A(m+1) not congruent 1 modulo 3. The ones corresponding to A(m+1) == 1 (mod 3) (the complement) are given in A272200.

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A000727, A001479, A002476, A001480, A272198, A272200 (complement relative to A002476).

filter:= proc(n) local S,x,y;
    if not isprime(n) then return false fi;
    S:= remove(hastype,[isolve(x^2+3*y^2=n)],negative);
    subs(S[1],x) mod 3 <> 1
end proc:
select(filter, [seq(i,i=7..1000,6)]); # Robert Israel, Apr 29 2019

filterQ[n_] := Module[{S, x, y}, If[!PrimeQ[n], Return[False]]; S = Solve[x > 0 && y > 0 && x^2 + 3 y^2 == n, Integers]; Mod[x /. S[[1]], 3] != 1];
Select[Range[7, 1000, 6], filterQ] (* Jean-François Alcover, Apr 21 2020, after Robert Israel *)

This sequence collects the 1 (mod 3) primes p(m) = A002476(m) = A(m)^2 + 3*B(m)^2 with positive A(m) not == 1 (mod 3), for m >= 1. A(m) = A001479(m+1).

A272205 A bisection of the primes congruent to 1 modulo 3 (A002476). This is the part depending on the corresponding A001479 entry being congruent to 4 or 5 modulo 6.

Original entry on oeis.org

19, 37, 43, 73, 103, 127, 163, 229, 283, 313, 331, 337, 379, 397, 421, 457, 463, 487, 499, 523, 541, 577, 607, 613, 619, 631, 691, 709, 727, 787, 811, 829, 853, 859, 877, 883, 967, 991, 997

Offset: 1

Wolfdieter Lang, May 05 2016

The other part of this bisection appears in A272204.
Each prime == 1 (mod 3) has a unique representation A002476(m) = A(m)^2 + 3*B(m)^2 with positive A(m) = A001479(m+1) and B(m) = A001480(m+1), m >= 1. The present sequence gives such primes corresponding to A(m) == 4, 5 (mod 6). The ones corresponding to A(m) == 1, 2 (mod 6) (the complement) are given in A272205.
The corresponding A001479 entries are 4, 5, 4, 5, 10, 10, 4, 11, 16, 11, 16, 17, 4, 17, 11, 5, 10, 22, 16, 4, 23, 23, 10, 5, 16, 22, 4, 11, 22, 28, 28, 23, 29, 28, 17, 4, 10, 22, 5, ...
See  A272204 for a comment on the relevance of this bisection in connection with the signs of the q-expansion coefficients of the modular cusp form eta^{12}(12*z) / (eta^4(6*z)*eta^4(24*z)).

Cf. A001479, A001480, A002476, A047239, A187076, A272203, A272204 (complement relative to A002476).

This sequence collects the 1 (mod 3) primes p(m) = A002476(m) = A(m)^2 + 3*B(m)^2 with positive A(m) == 4, 5 (mod 6), for m >= 1. A(m) = A001479(m+1).

A272204 A bisection of the primes congruent to 1 modulo 3 (A002476). This is the part depending on the corresponding A001479 entry being congruent to 1 or 2 modulo 6.

Original entry on oeis.org

7, 13, 31, 61, 67, 79, 97, 109, 139, 151, 157, 181, 193, 199, 211, 223, 241, 271, 277, 307, 349, 367, 373, 409, 433, 439, 547, 571, 601, 643, 661, 673, 733, 739, 751, 757, 769, 823, 907, 919, 937

Offset: 1

Wolfdieter Lang, May 05 2016

The other part of this bisection appears in A272205.
Each prime == 1 (mod 3) has a unique representation A002476(m) = A(m)^2 + 3*B(m)^2 with positive A(m) = A001479(m+1) and B(m) = A001480(m+1), m >= 1. The present sequence gives all such primes corresponding to A(m) == 1, 2 (mod 6). The ones corresponding to A(m) not == 1, 2 (mod 6) (the complement), that is == 4, 5 (mod 6), are given in A272205.
The corresponding A001479 entries are 2, 1, 2, 7, 8, 2, 7, 1, 8, 2, 7, 13, 1, 14, 8, 14, 7, 14, 13, 8, 7, 2, 19, 19, 1, 14, 20, 8, 13, 20, 19, 25, 25, 8, 26, 13, 1, 26, 20, 26, 13, ...
This bisection of the 1 (mod 3) primes A002476 is needed to determine the sign in the formula for the coefficients of the q-expansion (q = exp(2*Pi*i*z), Im(z) > 0) of the modular weight 2 cusp form
  eta^{12}(12*z) / (eta^4(6*z)*eta^4(24*z)) |A187076%20which%20gives%20the%20coefficients%20of%20the%20q-expansion%20of%20F(q)%20=%20Eta(q%5E%7B1/6%7D)%20/%20q%5E%7B1/6%7D%20=%20Product">{z=z(q)} =: Eta(q) with Dedekind's eta function. See A187076 which gives the coefficients of the q-expansion of F(q) = Eta(q^{1/6}) / q^{1/6} = Product{m>=0} (1 - q^(2*m))^{12} / ((1 - q^m)*(1 - q^(4*m)))^4. The q-expansion coefficients (called b(n)) of the modular cusp form are given there using multiplicativity. Note that there x can also be negative, whereas here A is positive.

Cf. A001479, A001480, A002476, A047239, A187076, A272203, A272205 (complement relative to A002476).

This sequence collects the 1 (mod 3) primes p(m) = A002476(m) = A(m)^2 + 3*B(m)^2 with positive A(m) == 1, 2 (mod 6), for m >= 1. A(m) = A001479(m+1).

A002476 Primes of the form 6m + 1.

Original entry on oeis.org

7, 13, 19, 31, 37, 43, 61, 67, 73, 79, 97, 103, 109, 127, 139, 151, 157, 163, 181, 193, 199, 211, 223, 229, 241, 271, 277, 283, 307, 313, 331, 337, 349, 367, 373, 379, 397, 409, 421, 433, 439, 457, 463, 487, 499, 523, 541, 547, 571, 577, 601, 607, 613, 619

Offset: 1

N. J. A. Sloane

Equivalently, primes of the form 3m + 1.
Rational primes that decompose in the field Q(sqrt(-3)). - N. J. A. Sloane, Dec 25 2017
Primes p dividing Sum_{k=0..p} binomial(2k, k) - 3 = A006134(p) - 3. - Benoit Cloitre, Feb 08 2003
Primes p such that tau(p) == 2 (mod 3) where tau(x) is the Ramanujan tau function (cf. A000594). - Benoit Cloitre, May 04 2003
Primes of the form x^2 + xy - 2y^2 = (x+2y)(x-y). - N. J. A. Sloane, May 31 2014
Primes of the form x^2 - xy + 7y^2 with x and y nonnegative. - T. D. Noe, May 07 2005
Primes p such that p^2 divides Sum_{m=1..2(p-1)} Sum_{k=1..m} (2k)!/(k!)^2. - Alexander Adamchuk, Jul 04 2006
A006512 larger than 5 (Greater of twin primes) is a subsequence of this. - Jonathan Vos Post, Sep 03 2006
A039701(A049084(a(n))) = A134323(A049084(a(n))) = 1. - Reinhard Zumkeller, Oct 21 2007
Also primes p such that the arithmetic mean of divisors of p^2 is an integer: sigma_1(p^2)/sigma_0(p^2) = C. (A000203(p^2)/A000005(p^2) = C). - Ctibor O. Zizka, Sep 15 2008
Fermat knew that these numbers can also be expressed as x^2 + 3y^2 and are therefore not prime in Z[omega], where omega is a complex cubic root of unity. - Alonso del Arte, Dec 07 2012
Primes of the form x^2 + xy + y^2 with x < y and nonnegative. Also see A007645 which also applies when x=y, adding an initial 3. - Richard R. Forberg, Apr 11 2016
For any term p in this sequence, let k = (p^2 - 1)/6; then A016921(k) = p^2. - Sergey Pavlov, Dec 16 2016; corrected Dec 18 2016
For the decomposition p=x^2+3*y^2, x(n) = A001479(n+1) and y(n) = A001480(n+1). - R. J. Mathar, Apr 16 2024

			Since 6 * 1 + 1 = 7 and 7 is prime, 7 is in the sequence. (Also 7 = 2^2 + 3 * 1^2 = (2 + sqrt(-3))(2 - sqrt(-3)).)
Since 6 * 2 + 1 = 13 and 13 is prime, 13 is in the sequence.
17 is prime but it is of the form 6m - 1 rather than 6m + 1, and is therefore not in the sequence.

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.
David A. Cox, Primes of the Form x^2 + ny^2. New York: Wiley (1989): 8.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 261.

Ray Chandler, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
C. Banderier, Calcul de (-3/p)
Barry Brent, Finite Field Models of Polynomials Interpolating Fourier Coefficients of Modular Functions for Hecke Groups, Integers (2024) Vol. 24, Art. No. A18. See p. 13.
F. S. Carey, On some cases of the Solutions of the Congruence z^p^(n-1)=1, mod p, Proceedings of the London Mathematical Society, Volume s1-33, Issue 1, November 1900, Pages 294-312.
A. Granville and G. Martin, Prime number races, arXiv:math/0408319 [math.NT], 2004.
Jorma K. Merikoski, Pentti Haukkanen, and Timo Tossavainen, The congruence x^n = -a^n (mod m): Solvability and related OEIS sequences, Notes. Num. Theor. Disc. Math. (2024) Vol. 30, No. 3, 516-529. See p. 526.
K. G. Reuschle, Tafeln complexer Primzahlen, Königl. Akademie der Wissenschaften, Berlin, 1875, p. 1.
Neville Robbins, On the Infinitude of Primes of the Form 3k+1, Fib. Q., 43,1 (2005), 29-30.
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
Index to sequences related to decomposition of primes in quadratic fields

Cf. A045331, A242660.
For values of m see A024899. Primes of form 3n - 1 give A003627.
These are the primes arising in A024892, A024899, A034936.
A091178 gives prime index.
Cf. A006512, A007528.
Subsequence of A016921 and of A050931.
Cf. A004611 (multiplicative closure).

Filtered(List([0..110],k->6*k+1),n-> IsPrime(n)); # Muniru A Asiru, Mar 11 2019

a002476 n = a002476_list !! (n-1)
a002476_list = filter ((== 1) . (`mod` 6)) a000040_list
-- Reinhard Zumkeller, Jan 15 2013

(#~ 1&p:) >: 6 * i.1000 NB. Stephen Makdisi, May 01 2018

[n: n in [1..700 by 6] | IsPrime(n)]; // Vincenzo Librandi, Apr 05 2011

a := [ ]: for n from 1 to 400 do if isprime(6*n+1) then a := [ op(a), n ]; fi; od: A002476 := n->a[n];

Select[6*Range[100] + 1, PrimeQ[ # ] &] (* Stefan Steinerberger, Apr 06 2006 *)

select(p->p%3==1,primes(100)) \\ Charles R Greathouse IV, Oct 31 2012

From R. J. Mathar, Apr 03 2011: (Start)
Sum_{n >= 1} 1/a(n)^2 = A175644.
Sum_{n >= 1} 1/a(n)^3 = A175645. (End)
a(n) = 6*A024899(n) + 1. - Zak Seidov, Aug 31 2016
From Vaclav Kotesovec, May 02 2020: (Start)
Product_{k>=1} (1 - 1/a(k)^2) = 1/A175646.
Product_{k>=1} (1 + 1/a(k)^2) = A334481.
Product_{k>=1} (1 - 1/a(k)^3) = A334478.
Product_{k>=1} (1 + 1/a(k)^3) = A334477. (End)
Legendre symbol (-3, a(n)) = +1 and (-3, A007528(n)) = -1, for n >= 1. For prime 3 one sets (-3, 3) = 0. - Wolfdieter Lang, Mar 03 2021

A007645 Generalized cuban primes: primes of the form x^2 + xy + y^2; or primes of the form x^2 + 3*y^2; or primes == 0 or 1 (mod 3).

Original entry on oeis.org

3, 7, 13, 19, 31, 37, 43, 61, 67, 73, 79, 97, 103, 109, 127, 139, 151, 157, 163, 181, 193, 199, 211, 223, 229, 241, 271, 277, 283, 307, 313, 331, 337, 349, 367, 373, 379, 397, 409, 421, 433, 439, 457, 463, 487, 499, 523, 541, 547, 571, 577, 601, 607, 613

Offset: 1

N. J. A. Sloane, Mira Bernstein and Robert G. Wilson v

Also, odd primes p such that -3 is a square mod p. - N. J. A. Sloane, Dec 25 2017
Equivalently, primes of the form p = (x^3 - y^3)/(x - y). If x=y+1 we get the cuban primes A002407, which is therefore a subsequence.
These are not to be confused with the Eisenstein primes, which are the primes in the ring of integers Z[w], where w = (-1+sqrt(-3))/2. The present sequence gives the rational primes which are also Eisenstein primes. - N. J. A. Sloane, Feb 06 2008
Also primes of the form x^2+3y^2 and, except for 3, x^2+xy+7y^2. See A140633. - T. D. Noe, May 19 2008
Conjecture: this sequence is Union(A002383,A162471). - Daniel Tisdale, Jul 04 2009
Primes p such that antiharmonic mean B(p) of the numbers k < p such that gcd(k, p) = 1 is not integer, where B(p) = A053818(p) / A023896(p) = A175505(p) / A175506(p) = (2p - 1) / 3. Primes p such that A175506(p) > 1. Subsequence of A179872. Union a(n) + A179891 = A179872. Example: a(6) = 37 because B(37) = A053818(37) / A023896(37) = A175505(37) / A175506(37) = 16206 / 666 = 73 / 3 (not integer). Cf. A179871, A179872, A179873, A179874, A179875, A179876, A179877, A179878, A179879, A179880, A179882, A179883, A179884, A179885, A179886, A179887, A179890, A179891, A003627, A034934. - Jaroslav Krizek, Aug 01 2010
Subsequence of Loeschian numbers, cf. A003136 and A024614; A088534(a(n)) > 0. - Reinhard Zumkeller, Oct 30 2011
Primes such that there exist a unique x, y, with 1 < x <= y < p, x + y == 1 (mod p) and x * y == 1 (mod p). - Jon Perry, Feb 02 2014
The prime factors of A002061. - Richard R. Forberg, Dec 10 2014
This sequence gives the primes p which solve s^2 == -3 (mod 4*p) (see Buell, Proposition 4.1., p. 50, for Delta = -3). p = 2 is not a solution. x^2 == -3 (mod 4) has solutions for all odd x. x^2 == -3 (mod  p) has for odd primes p, not 3, the solutions of Legendre(-3|p) = +1 which are p == {1, 7} (mod 12). For p = 3 the representative solution is x = 0. Hence the solution of s^2 == -3 (mod 4*p) are the odd primes p = 3 and p == {1, 7} (mod 12) (or the primes p = 0, 1 (mod 3)). - Wolfdieter Lang, May 22 2021

D. A. Buell, Binary Quadratic Forms. Springer-Verlag, NY, 1989, p. 50.
Conway, J. H. and Guy, R. K. The Book of Numbers. New York: Springer-Verlag, pp. 220-223, 1996.
David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989, p. 7.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Wagon, S. "Eisenstein Primes." Section 9.8 in Mathematica in Action. New York: W. H. Freeman, pp. 319-323, 1991.

T. D. Noe, Table of n, a(n) for n = 1..1000
U. P. Nair, Elementary results on the binary quadratic form a^2+ab+b^2, arXiv:math/0408107 [math.NT], 2004.
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references)
Eric Weisstein's World of Mathematics, Eisenstein Integer.

Subsequence of A003136.
Subsequences include A002407, A002648, and A201477.
Apart from initial term, same as A045331.
Cf. A001479, A001480 (x and y such that a(n) = x^2 + 3y^2).
Cf. A000040, A003627.
Primes in A003136 and A034017.

a007645 n = a007645_list !! (n-1)
a007645_list = filter ((== 1) . a010051) $ tail a003136_list
-- Reinhard Zumkeller, Jul 11 2013, Oct 30 2011

select(isprime,[3, seq(6*k+1, k=1..1000)]); # Robert Israel, Dec 12 2014

Join[{3},Select[Prime[Range[150]],Mod[#,3]==1&]] (* Harvey P. Dale, Aug 21 2021 *)

forprime(p=2,1e3,if(p%3<2,print1(p", "))) \\ Charles R Greathouse IV, Jun 16 2011

p == 0 or 1 (mod 3).
{3} UNION A002476. - R. J. Mathar, Oct 28 2008
A007645 UNION A003627 = A000040. - Juri-Stepan Gerasimov, Jan 28 2010

Entry revised by N. J. A. Sloane, Jan 29 2013

A001480 Let p = A007645(n) be the n-th generalized cuban prime and write p = x^2 + 3*y^2; a(n) = y.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 3, 2, 1, 4, 5, 4, 1, 6, 3, 5, 7, 6, 7, 2, 8, 1, 7, 3, 6, 8, 5, 6, 3, 9, 8, 5, 4, 10, 11, 2, 11, 6, 4, 10, 12, 9, 12, 11, 1, 9, 13, 2, 7, 13, 4, 12, 13, 14, 11, 7, 9, 10, 4, 15, 14, 9, 6, 15, 5, 14, 16, 1, 3, 7, 10, 2, 5, 14, 17, 13, 9, 16, 17

Offset: 1

N. J. A. Sloane

a(n) = A000196((A007645(n) - A000290(A001479(n))) / 3). - Reinhard Zumkeller, Jul 11 2013

A. J. C. Cunningham, Quadratic Partitions. Hodgson, London, 1904, p. 1.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
B. van der Pol and P. Speziali, The primes in k(rho). Nederl. Akad. Wetensch. Proc. Ser. A. {54} = Indagationes Math. 13, (1951). 9-15 (1 plate).

T. D. Noe, Table of n, a(n) for n = 1..1000
A. J. C. Cunningham, Quadratic Partitions, Hodgson, London, 1904 [Annotated scans of selected pages]
S. R. Finch, Powers of Euler's q-Series, (arXiv:math.NT/0701251).
B. van der Pol and P. Speziali, The primes in k(rho) (annotated and scanned copy)

Cf. A001479, A007645.

a001480 n = a000196 $ (`div` 3) $ (a007645 n) - (a001479 n) ^ 2
-- Reinhard Zumkeller, Jul 11 2013

nmax = 63; nextCuban[p_] := If[p1 = NextPrime[p]; Mod[p1, 3] > 1, nextCuban[p1], p1]; cubanPrimes = NestList[ nextCuban, 3, nmax ]; f[p_] := y /. ToRules[ Reduce[x > 0 && y > 0 && p == x^2 + 3*y^2, {x, y}, Integers]]; a[1] = 1; a[n_] := f[cubanPrimes[[n]]]; Table[ a[n] , {n, 1, nmax}] (* Jean-François Alcover, Oct 19 2011 *)

do(lim)=my(v=List(), q=Qfb(1,0,3)); forprime(p=2,lim, if(p%3==2,next); listput(v, qfbsolve(q,p)[2])); Vec(v) \\ Charles R Greathouse IV, Feb 07 2017

Definition revised by N. J. A. Sloane, Jan 29 2013

A259694 a(n) = Sum(k^6sigma(k)sigma(n-k),k=1..n-1).

Original entry on oeis.org

0, 1, 195, 3496, 38195, 192780, 977386, 3216320, 9860049, 27321870, 65803045, 144005856, 308944925, 635774072, 1112550390, 2153146880, 3618341556, 6391671525, 9949570455, 16725562160, 24691972080, 40979569092, 56807498030, 89450231040, 120404165825

Offset: 1

N. J. A. Sloane, Jul 03 2015

nonn

This was formerly A001479.

Colin Barker, Table of n, a(n) for n = 1..1000
J. Touchard, On prime numbers and perfect numbers, Scripta Math., 129 (1953), 35-39. [Annotated scanned copy]

Cf. A000385, A000441, A000477, A000499, A259692, A259693, A259695, A259696.

S:=(n,e)->add(k^e*sigma(k)*sigma(n-k),k=1..n-1); f:=e->[seq(S(n,e),n=1..30)]; f(6);

a(n) = sum(k=1, n-1, k^6*sigma(k)*sigma(n-k)) \\ Colin Barker, Jul 16 2015

A272198 The p-defect p - N(p) of the congruence y^2 == x^3 + 1 (mod p) for primes p, where N(p) is the number of solutions given by A272197(n).

Original entry on oeis.org

0, 0, 0, -4, 0, 2, 0, 8, 0, 0, -4, -10, 0, 8, 0, 0, 0, 14, -16, 0, -10, -4, 0, 0, 14, 0, 20, 0, 2, 0, 20, 0, 0, -16, 0, -4, 14, 8, 0, 0, 0, 26, 0, 2, 0, -28, -16, -28, 0, -22, 0, 0, 14, 0, 0, 0, 0, -28, 26, 0

Offset: 1

Wolfdieter Lang, May 02 2016

This sequence for an elliptic curve (of the Bachet-Mordell type) is discussed in the Silverman reference. In Exercise 45.5, in the table on p. 405, the p-defects are called a_p, and are shown for primes 2 to 113.
The modularity pattern series is the expansion of the 51st modular cusp form of weight 2 and level N=36, given in the table I of the Martin reference, i.e., eta^4(6*z) in powers of q = exp(2*Pi*i*z), with Im(z) > 0. Here eta is the Dedekind function. See A000727 for the expansion in powers of q^6 (after deleting a factor q^(1/6)). Note that also for the possibly bad prime 2 and the bad prime 3 this expansion gives the correct numbers 0 (the discriminant of this elliptic curve is -3^3).
See also the comment on the Martin-Ono reference in A272197 which implies that eta^4(6*z) provides the modularity sequence for this elliptic curve.
For primes p == 0 and 2 (mod 3) (A045309) a(p) = 0. The proof runs along the same line as the one given in the Silverstein reference on pp. 400 - 402 for 17 replaced by 1. From the expansion of the known modularity function eta^4(6*z) follows that only the coefficients for powers q^n with n == 1 (mod 6) are nonzero, and therefore all a(p) for primes p == 0 and 2 (mod 3) have to vanish.
If prime(n) == 1 (mod 3) = A002476(m) (for a unique m = m(n)) then prime(n) = A(m)^2 + 3*B(m)^2 with A(m) = A001479(m+1) and B(m) = A001480(m+1), m >= 1. In this case (4*prime(n) - a(n)^2)/12, seems to be a square, q(m)^2. In fact is seems that (the positive) q(m) = B(m). This is true at least for the first 80 primes 1 (mod 3), i.e. for such primes <= 997. (In the Silverman reference, in hint c) for Exercise 4.5, on p. 405, a more complicated way is suggested: 4*p is decomposed there non-uniquely instead of p uniquely.) If this conjecture is true then a(n) = 2*(+/-sqrt(prime(n) - 3*B(m)^2)) = +- 2*A(m) for prime(n) = A002476(m). This leads to a bisection of the primes 1 (mod 3) into two types: type I if the + sign applies, and type II for the - sign. Primes of type I are given in A272200: 13, 19, 43, 61, 97, ... and those of type II in A272201: 7, 31, 37, 67, 73, ...

			a(1) = 2 - A272197(1) = 0, and 2 == 2(mod 3).
a(4) = 7 - A272197(4) = 7 - 11 = -4, and 7 = A002476(1) = 2^2 + 3*1^2, 2 = A001479(1+1), 7 = A272201(1), hence a(4) = -2*2 = -4.
a(6) = 13 - A272197(6) = 13 - 11 = 2, and 13 = A002476(2) = 1^2 + 3*2^2; 1 = A001479(2+1), 13 = A272200(1), hence a(6) = +2*1 = +2.

J. H. Silverman, A Friendly Introduction to Number Theory, 3rd ed., Pearson Education, Inc, 2006, Exercise 45.5, p. 405, Exercise 47.2, p. 415, and pp. 400 - 402 (4th ed., Pearson 2014, Exercise 5, p. 371, Exercise 2, p. 385, and pp. 366 - 368).

Seiichi Manyama, Table of n, a(n) for n = 1..10000
Y. Martin, Multiplicative eta-quotients, Trans. Amer. Math. Soc. 348 (1996), no. 12, 4825-4856, see page 4852 Table I.
Yves Martin and Ken Ono, Eta-Quotients and Elliptic Curves, Proc. Amer. Math. Soc. 125, No 11 (1997), 3169-3176.

Cf. A000040, A002476, A001479, A001480, A045309, A272197, A272200, A272201.

a(n) = prime(n) - N(prime(n)), n = 1, where N(prime(n)) = A272197(n), the number of solutions of the congruence y^2 == x^3 + 1 (mod prime(n)).
a(n) = 0 for prime(n) == 0, 2 (mod 3) (see A045309).
The above given conjecture for primes 1 (mod 3) is true because Mordell proved the Ramanujan conjecture on the expansion coefficients of eta^4(6*z), and with the present a(n) the result of Ramanujan follows. See the references and a comment on A000727.
a(n) = +2*A001479(m+1) if prime(n) == A002476(m) (m is unique) is a prime of A272200 (type I).
a(n) = -2*A001479(m+1) if prime(n) == A002476(m) is from A272201 (type II).
See a comment above for the bisection of the primes 1 (mod 3) into type I and II.

A272203 P-defects p - N(p) of the congruence y^2 == x^3 - 1 (mod p) for primes p, where N(p) is the number of solutions given by A272202(n).

Original entry on oeis.org

0, 0, 0, 4, 0, 2, 0, -8, 0, 0, 4, -10, 0, -8, 0, 0, 0, 14, 16, 0, -10, 4, 0, 0, 14, 0, -20, 0, 2, 0, -20, 0, 0, 16, 0, 4, 14, -8, 0, 0, 0, 26, 0, 2, 0, 28, 16, 28, 0, -22, 0, 0, 14, 0, 0, 0, 0, 28, 26, 0, -32, 0, 16, 0, -22, 0, -32, -34, 0, 14, 0, 0, 4, 38, -8, 0, 0, -34, 0, 38, 0, -22, 0, 2, 28, 0, 0, -10, 0, -20, 0, 0, -44, 0, -32, 0, 0, 0, -8, -46, 40, 0, 0, 0, 16, -46

Offset: 1

Wolfdieter Lang, May 05 2016

The analysis of this elliptic curve runs along the same lines as in A000727, A272197 and A272198, and it is inspired by the Silverman reference where the curve y^2 = x^3 + 1 modulo primes is treated.
The series showing the modularity pattern is the expansion of the 67th modular cusp form of weight 2 and level N=144, given in the table I of the Martin reference, i.e., eta^{12}(12*z)/( eta^4(6*z)*eta^4(24*z)), symbolically 12^{12} 6^(-4) 24^{-4}, in powers of q = exp(2*Pi*i*z), with Im(z) > 0. Here eta is the Dedekind function. See A187076 for the expansion in powers of q^6 (after deleting a factor q^(1/6)). Note that also for the possibly bad prime 2 and the bad prime 3 this expansion gives the correct numbers 0 (the discriminant of this elliptic curve is -3^3).
See also the comment on the Martin-Ono reference in A272202 which implies that 12^{12} 6^(-4) 24^{-4} provides the modularity sequence for this elliptic curve.
If prime(n) == 1 (mod 3) = A002476(m) (for a unique m = m(n)) then prime(n) = A(m)^2 + 3*B(m)^2 with A(m) = A001479(m+1) and B(m) = A001480(m+1), m >= 1. In this case (4*prime(n) - a(n)^2)/12, seems to be a square, q(m)^2. In fact is seems that (the positive) q(m) = B(m). If this conjecture is true then a(n) = 2*(+-sqrt(prime(n) - 3*B(m)^2)) = +- 2*A(m) for prime(n) = A002476(m). This leads to a bisection of the primes 1 (mod 3) into two types: type I if the + sign applies, and type II for the - sign. Primes of type I are given in A272204: 7,13,31,61,67, ... and those of type II in A272205: 19,37,43,73,103, ...

			a(1) = 2 - A272202(1) = 0, and 2 == 2 (mod 3).
a(4) = 7 - A272202(4) = 7 - 3 = +4, and 7 = A002476(1) = 2^2 + 3*1^2, 2 = A001479(1+1), 7 = A272204(1), hence a(4) = +2*2 = +4.
a(8) = 19 - A272202(8) = 19 - 27 = -8, and 19 = A002476(3) = 4^2 + 3*1^2; 4=A001479(3+1), 19 = A272205(1), hence a(8) = 2*(-4) = -8.

J. H. Silverman, A Friendly Introduction to Number Theory, 3rd ed., Pearson Education, Inc, 2006, Exercise 45.5, p. 405, Exercise 47.2, p. 415, and pp. 400 - 402 (4th ed., Pearson 2014, Exercise 5, p. 371, Exercise 2, p. 385, and pp. 366 - 368).

Seiichi Manyama, Table of n, a(n) for n = 1..10000
Y. Martin, Multiplicative eta-quotients, Trans. Amer. Math. Soc. 348 (1996), no. 12, 4825-4856, see page 4852 Table I.
Yves Martin and Ken Ono, Eta-Quotients and Elliptic Curves, Proc. Amer. Math. Soc. 125, No 11 (1997), 3169-3176.

Cf. A000040, A187076, A272202, A272204, A272205.

a(n) = prime(n) - N(prime(n)), n = 1, where N(prime(n)) = A272202(n), the number of solutions of the congruence y^2 == x^3 - 1 (mod prime(n)).
a(n) = 0 for prime(n) == 0, 2 (mod 3) (see A045309).
The above given conjecture for primes 1 (mod 3) is expected to be true by analogy to the case A272198 where only the signs differ.
a(n) = +2*A001479(m+1) if prime(n) == A002476(m) (m is unique) is a prime of A272204 (type I).
a(n) = -2*A001479(m+1) if prime(n) == A002476(m) is from A272205 (type II).
See a comment above for this bisection of the primes 1 (mod 3) into type I and II.

cp's OEIS Frontend

A272200 Bisection of primes congruent to 1 modulo 3 (A002476), depending on the corresponding A001479 entry being congruent to 1 modulo 3 or not. Here the first case.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A272201 Bisection of primes congruent to 1 modulo 3 (A002476), depending on the corresponding A001479 entry being congruent to 1 modulo 3 or not. Here the second case.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A272205 A bisection of the primes congruent to 1 modulo 3 (A002476). This is the part depending on the corresponding A001479 entry being congruent to 4 or 5 modulo 6.

Views

Author

Keywords

Comments

Crossrefs

Formula

A272204 A bisection of the primes congruent to 1 modulo 3 (A002476). This is the part depending on the corresponding A001479 entry being congruent to 1 or 2 modulo 6.

Views

Author

Keywords

Comments

Crossrefs

Formula

A002476 Primes of the form 6m + 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

J

Magma

Maple

Mathematica

PARI

Formula

A007645 Generalized cuban primes: primes of the form x^2 + xy + y^2; or primes of the form x^2 + 3*y^2; or primes == 0 or 1 (mod 3).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A001480 Let p = A007645(n) be the n-th generalized cuban prime and write p = x^2 + 3*y^2; a(n) = y.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

A259694 a(n) = Sum(k^6sigma(k)sigma(n-k),k=1..n-1).