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.

A084109 n is congruent to 1 (mod 4) and is not the sum of two squares.

This page as a plain text file.
%I A084109 #25 Jul 01 2016 22:46:24
%S A084109 21,33,57,69,77,93,105,129,133,141,161,165,177,189,201,209,213,217,
%T A084109 237,249,253,273,285,297,301,309,321,329,341,345,357,381,385,393,413,
%U A084109 417,429,437,453,465,469,473,489,497
%N A084109 n is congruent to 1 (mod 4) and is not the sum of two squares.
%C A084109 Alternatively, n is congruent to 1 (mod 4) with at least 2 distinct prime factors congruent to 3 (mod 4) in the squarefree part of n. - Comment corrected by _Jean-Christophe Hervé_, Oct 25 2015
%C A084109 Applications to the theory of optimal weighing designs and maximal determinants: An (n+1) X (n+1) conference matrix is impossible.
%C A084109 The upper bound of Ehlich/Wojtas on the determinant of a (0,1) matrix of order congruent to 1 (mod 4) cannot be achieved for n X n matrices.
%C A084109 The bound of Ehlich/Wojtas on the determinant of a (-1,1) matrix of order congruent to 2 (mod 4) cannot be achieved for (n+1) X (n+1) matrices.
%C A084109 Numbers with only odd prime factors, of which a strictly positive even number are raised to an odd power and congruent to 3 (mod 4). - _Jean-Christophe Hervé_, Oct 24 2015
%D A084109 F. J. MacWilliams and N. J. A. Sloane, The Theory of Error-Correcting Codes, Elsevier-North Holland, 1978, p. 56.
%H A084109 Jean-Christophe Hervé, <a href="/A084109/b084109.txt">Table of n, a(n) for n = 1..1000</a>
%H A084109 H. Ehlich, <a href="http://dx.doi.org/10.1007/BF01111249">Determinantenabschätzungen für binäre Matrizen</a>, Math. Z. 83 (1964) 123-132.
%H A084109 D. Raghavarao, <a href="http://dx.doi.org/10.1214/aoms/1177705664">Some aspects of weighing designs</a>, Ann. Math. Stat. 31 (1960) 878-884.
%e A084109 a(1) = 3*7 = 21, a(2) = 3*11 = 33, a(3) = 3*19 = 57, a(14) = 3^3*7 = 189.
%p A084109 N:= 1000: # to get all entries <= N
%p A084109 S:= {seq(i,i=1..N,4)} minus
%p A084109    {seq(seq(i^2+j^2, j=1..floor(sqrt(N-i^2)),2),i=0..floor(sqrt(N)),2)}:
%p A084109 sort(convert(S,list)); # _Robert Israel_, Oct 25 2015
%t A084109 a[m_] := Complement[Range[1, m, 4], Union[Flatten[Table[j^2+k^2, {j, 1, Sqrt[m], 2}, {k, 0, Sqrt[m], 2}]]]]
%o A084109 (PARI) is(n)=if(n%4!=1, return(0)); my(f=factor(n)); for(i=1, #f~, if(f[i, 1]%4==3 && f[i, 2]%2, return(1))); 0 \\ _Charles R Greathouse IV_, Jul 01 2016
%Y A084109 Cf. A000952, A003432, A003433, A022544, A001481.
%K A084109 easy,nonn
%O A084109 1,1
%A A084109 _William P. Orrick_, Jun 18 2003