A224770 -id:A224770 - cp's OEIS frontend

A008784 Numbers k such that sqrt(-1) mod k exists; or, numbers that are primitively represented by x^2 + y^2.

1, 2, 5, 10, 13, 17, 25, 26, 29, 34, 37, 41, 50, 53, 58, 61, 65, 73, 74, 82, 85, 89, 97, 101, 106, 109, 113, 122, 125, 130, 137, 145, 146, 149, 157, 169, 170, 173, 178, 181, 185, 193, 194, 197, 202, 205, 218, 221, 226, 229, 233, 241, 250, 257, 265, 269, 274, 277, 281, 289

Offset: 1

N. J. A. Sloane, Olivier Gérard

nonn

Numbers whose prime divisors are all congruent to 1 mod 4, with the exception of at most a single factor of 2. - Franklin T. Adams-Watters, Sep 07 2008
In appears that {a(n)} is the set of proper divisors of numbers of the form m^2+1. - Kaloyan Todorov (kaloyan.todorov(AT)gmail.com), Mar 25 2009 [This conjecture is correct. - Franklin T. Adams-Watters, Oct 07 2009]
If a(n) is a term of this sequence, then so too are all of its divisors (Euler). - Ant King, Oct 11 2010
From Richard R. Forberg, Mar 21 2016: (Start)
For a given a(n) > 2, there are 2^k solutions to sqrt(-1) mod n (for some k >= 1), and 2^(k-1) solutions primitively representing a(n) by x^2 + y^2.
Record setting values for the number of solutions (i.e., the next higher k values), occur at values for a(n) given by A006278.
A224450 and A224770 give a(n) values with exactly one and exactly two solutions, respectively, primitively representing integers as x^2 + y^2.
The 2^k different solutions for sqrt(-1) mod n can written as values for j, with j <= n, such that integers r = sqrt(n*j-1). However, the set of j values (listed from smallest to largest) transform into themselves symmetrically (i.e., largest to smallest) when the solutions are written as n-r. When the same 2^k solutions are written as r-j, it is clear that only 2^(k-1) distinct and independent solutions exist. (End)
Lucas uses the fact that there are no multiples of 3 in this sequence to prove that one cannot have an equilateral triangle on the points of a square lattice. - Michel Marcus, Apr 27 2020
For n > 1, terms are precisely the numbers such that there is at least one pair (m,k) where m + k = a(n), and m*k == 1 (mod a(n)), m > 0 and m <= k. - Torlach Rush, Oct 18 2020
A pair (s,t) such that s+t = a(n) and s*t == +1 (mod a(n)) as above is obtained from a square root of -1 (mod a(n)) for s and t = a(n)-s. - Joerg Arndt, Oct 24 2020
The Diophantine equation x^2 + y^2 = z^5 + z with gcd(x, y, z) = 1 has solutions iff z is a term of this sequence. See Gardiner reference, Olympiad links and A340129. - Bernard Schott, Jan 17 2021
Except for 1, numbers of the form a + b + 2*sqrt(a*b - 1) for positive integers a,b such that a*b-1 is a square. - Davide Rotondo, Nov 10 2024

B. C. Berndt & R. A. Rankin, Ramanujan: Letters and Commentary, see p. 176; AMS Providence RI 1995.
J. W. S. Cassels, Rational Quadratic Forms, Cambridge, 1978.
Leonard Eugene Dickson, History of the Theory Of Numbers, Volume II: Diophantine Analysis, Chelsea Publishing Company, 1992, pp.230-242.
A. Gardiner, The Mathematical Olympiad Handbook: An Introduction to Problem Solving, Oxford University Press, 1997, reprinted 2011, Problem 6 pp. 63 and 167-168 (1985).
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, Ch. 20.2-3.

T. D. Noe, Table of n, a(n) for n = 1..1000
J.-P. Allouche and F. M. Dekking, Generalized Beatty sequences and complementary triples, arXiv:1809.03424 [math.NT], 2018.
British Mathematical Olympiad, 1985 - Problem 6.
Édouard Lucas, Théorème sur la géométrie des quinconces, Nouvelles annales de mathématiques : journal des candidats aux écoles polytechnique et normale, Série 2, Tome 17 (1878), p. 129-130.
P. Cho-Ho Lam, Representation of integers using a^2+b^2-dc^2, J. Int. Seq. 18 (2015) 15.8.6, Theorems 2 and 3.
Richard J. Mathar, Construction of Bhaskara pairs, arXiv:1703.01677 [math.NT], 2017.
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references).
Index to sequences related to Olympiads.

Apart from the first term, a subsequence of A000404.

Cf. A001481, A022544, A020893, A037942, A034023, A057756, A076948, A045673, A004613, A340129, A192450 (complement).

import Data.List.Ordered (union)
a008784 n = a008784_list !! (n-1)
a008784_list = 1 : 2 : union a004613_list (map (* 2) a004613_list)
-- Reinhard Zumkeller, Oct 25 2015

with(numtheory); [seq(mroot(-1,2,p),p=1..300)];

data=Flatten[FindInstance[x^2+y^2==# && 0<=x<=# && 0<=y<=# && GCD[x,y]==1,{x,y},Integers]&/@Range[289],1]; x^2+y^2/.data//Union (* Ant King, Oct 11 2010 *)
Select[Range[289], And @@ (Mod[#, 4] == 1 & ) /@ (fi = FactorInteger[#]; If[fi[[1]] == {2, 1}, Rest[fi[[All, 1]]], fi[[All, 1]]])&] (* Jean-François Alcover, Jul 02 2012, after Franklin T. Adams-Watters *)

is(n)=if(n%2==0,if(n%4,n/=2,return(0)));n==1||vecmax(factor(n)[,1]%4)==1 \\ Charles R Greathouse IV, May 10 2012

list(lim)=my(v=List([1,2]),t); lim\=1; for(x=2,sqrtint(lim-1), t=x^2; for(y=0,min(x-1,sqrtint(lim-t)), if(gcd(x,y)==1, listput(v,t+y^2)))); Set(v) \\ Charles R Greathouse IV, Sep 06 2016

for(n=1,300,if(issquare(Mod(-1, n)),print1(n,", "))); \\ Joerg Arndt, Apr 27 2020

Checked by T. D. Noe, Apr 19 2007

A157228 Number of primitive inequivalent inclined square sublattices of square lattice of index n.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 2, 0, 0

Offset: 1

N. J. A. Sloane, Feb 25 2009

nonn

From Andrey Zabolotskiy, May 09 2018: (Start)
Also, the number of partitions of n into 2 distinct coprime squares.
All such sublattices (including non-primitive ones) are counted in A025441.
The primitive sublattices that have the same symmetries (including the orientation of the mirrors) as the parent lattice are not counted here; they are counted in A019590, and all such sublattices (including non-primitive ones) are counted in A053866.
For n > 2, equals A193138. (End)

Andrey Zabolotskiy, Table of n, a(n) for n = 1..5000
John S. Rutherford, Sublattice enumeration. IV. Equivalence classes of plane sublattices by parent Patterson symmetry and colour lattice group type, Acta Cryst. (2009). A65, 156-163. [See Table 5.]

Cf. A193138, A145393 (all sublattices of the square lattice), A025441, A019590, A053866, A157226, A157230, A157231, A000089, A304182, A224450, A224770, A281877, A024362.

a(n) = (A000089(n) - A019590(n)) / 2. - Andrey Zabolotskiy, May 09 2018
a(n) = 1 if n>2 is in A224450, a(n) = 2 if n is in A224770, a(n) is a higher power of 2 if n is in A281877 (first time reaches 8 at n = 32045). - Andrey Zabolotskiy, Sep 30 2018
a(n) = b(n) for odd n, a(n) = b(n/2) for even n, where b(n) = A024362(n). - Andrey Zabolotskiy, Jan 23 2022

New name and more terms from Andrey Zabolotskiy, May 09 2018

A224450 Numbers that are the primitive sum of two nonzero squares in exactly one way.

Original entry on oeis.org

2, 5, 10, 13, 17, 25, 26, 29, 34, 37, 41, 50, 53, 58, 61, 73, 74, 82, 89, 97, 101, 106, 109, 113, 122, 125, 137, 146, 149, 157, 169, 173, 178, 181, 193, 194, 197, 202, 218, 226, 229, 233, 241, 250, 257, 269, 274, 277, 281, 289, 293, 298, 313, 314, 317, 337

Offset: 1

Wolfdieter Lang, Apr 17 2013

nonn

If one includes 1 as the first entry then this sequence gives the numbers that are the primitive sum of two squares (square 0 allowed) in exactly one way, if neither the order of the squares nor the signs of the numbers to be squared matters.
Compare this sequence with A025284.
If 2 is omitted from this sequence then all members are primitively represented by two distinct nonzero squares in exactly one way.
The sequence A193138(n), n >= 3, gives the multiplicities of the primitive sums of two squares (automatically distinct and nonzero for n >= 3 if such a sum exists at all).
Numbers such that there is exactly one pair (m,k) where m + k = a(n), and m*k == 1 (mod a(n)), m > 0 and m <= k. - Torlach Rush, Oct 19 2020
A pair (s,t) such that s+t = a(n) and s*t == +1 (mod a(n)) as above is obtained from a square root of -1 (mod a(n)) for s and t = a(n)-s. - Joerg Arndt, Oct 24 2020

			a(1) = 2 because m = 2 is the first number with a unique doublet (a,b) in question, namely (1,1) (gcd(1,1) = 1).
This is the only case with equal entries a and b (the non-distinct case).
8 is not a member of this sequence (but of A025284) because the only representation is 2^2 +2^2 and (2,2) is not primitive. Similarly for 18, 20, ...
a(2) = 5 because 5 is the second smallest number satisfying the given requirements. 3 and 4 have no representation as sum of two nonzero squares, and the unique doublet for 5 is (1,2) (with distinct a and b).

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

Cf. A025284, A008784 (primitive sums of two squares with square 0 included), A224770 (exactly 2 ways), A193138 (multiplicities).

nn = 20; t = Sort[Select[Flatten[Table[If[GCD[a, b] == 1, a^2 + b^2, 0], {a, nn}, {b, a, nn}]], 0 < # <= nn^2 &]]; t2 = Transpose[Select[Tally[t], #[[2]] == 1 &]][[1]] (* T. D. Noe, Apr 20 2013 *)

This sequence gives the increasingly ordered numbers m which satisfy m = a^2 + b^2, with a and b integers, 0 < a <= b, gcd(a,b) = 1, and there is only one such representation, denoted by one doublet (a,b).

A281877 Numbers that are a primitive sum of two squares in more than 2 ways.

Original entry on oeis.org

1105, 1885, 2210, 2405, 2465, 2665, 3145, 3445, 3485, 3770, 3965, 4505, 4745, 4810, 4930, 5185, 5330, 5365, 5525, 5785, 5945, 6205, 6290, 6305, 6409, 6565, 6890, 6970, 7085, 7345, 7565, 7585, 7685, 7930, 8177, 8245, 8585, 8845, 8905, 9010, 9061, 9265, 9425, 9490, 9605, 9685, 9805

Offset: 1

R. J. Mathar, Feb 01 2017

nonn

"Primitive" means that x and y are coprime in the representations x^2+y^2.

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

Cf. A224450 (exactly 1 way), A224770 (exactly 2 ways), A008784, A097102.

N:= 10^4: # to get all terms <= N
V:= Vector(N):
for a from 1 to floor(sqrt(N)) do
  for b from 1 to min(a, floor(sqrt(N-a^2))) do
    if igcd(a,b) > 1 then next fi;
    r:= a^2 + b^2;
    V[r]:= V[r]+1;
od od:
select(n -> V[n] > 2, [$1..N]); # Robert Israel, Feb 07 2017

A037942 Numbers of the form x^2 + y^2 with x >= 0, y >= 0, gcd(x,y)=1, with multiplicity.

Original entry on oeis.org

1, 2, 5, 10, 13, 17, 25, 26, 29, 34, 37, 41, 50, 53, 58, 61, 65, 65, 73, 74, 82, 85, 85, 89, 97, 101, 106, 109, 113, 122, 125, 130, 130, 137, 145, 145, 146, 149, 157, 169, 170, 170, 173, 178, 181, 185, 185, 193, 194

Offset: 1

Nissim Broudo (broudo(AT)brain.math.fau.edu)

nonn

N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references).
Eric Weisstein's World of Mathematics, Postage Stamp Problem.

If repetitions are deleted, same as A008784.

Cf. A224770 (multiplicity 2), A281877 (multiplicity >= 3).

s=14; v=[]; for(y=0,s,for(x=0,y,if(gcd(x,y)==1,r=x*x+y*y; if(r<=s*s,v=concat(v,r))))); vecsort(v) \\ Craig Clapp (craig_clapp(AT)earthlink.net), Apr 30 2010

Missing second instance of 85 (6^2+7^2 = 2^2+9^2 = 85) added by Craig Clapp (craig_clapp(AT)earthlink.net), Apr 30 2010

cp's OEIS Frontend

A008784 Numbers k such that sqrt(-1) mod k exists; or, numbers that are primitively represented by x^2 + y^2.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

Extensions

A157228 Number of primitive inequivalent inclined square sublattices of square lattice of index n.

Views

Author

Keywords

Comments

Links

Crossrefs

Formula

Extensions

A224450 Numbers that are the primitive sum of two nonzero squares in exactly one way.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A281877 Numbers that are a primitive sum of two squares in more than 2 ways.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

A037942 Numbers of the form x^2 + y^2 with x >= 0, y >= 0, gcd(x,y)=1, with multiplicity.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Extensions