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.

A165918 a(n) = number of squarefree quadratic non-residues of n.

This page as a plain text file.
%I A165918 #6 Jul 23 2013 08:17:42
%S A165918 0,0,1,2,2,2,3,5,4,3,4,7,5,5,7,10,7,8,6,11,10,8,10,15,10,8,11,15,10,
%T A165918 12,12,18,15,11,15,21,13,13,19,25,16,19,14,25,24,17,19,29,19,18,21,28,
%U A165918 17,22,25,33,25,18,19,35,21,22,32,34,28,29,20,36,33,31,29,44,26,23
%N A165918 a(n) = number of squarefree quadratic non-residues of n.
%H A165918 C. H. Gribble, <a href="/A165918/b165918.txt">Table of n, a(n) for n=1,...,1000</a>.
%o A165918 (PARI) squares(n) = {local(s = Set()); for (j=1, n, s = setunion(s, Set(j^2 % n));); return (s);}
%o A165918 qnonr(n) = {local(s = Set()); sq = squares(n);for (j=0, n-1, if (length(setintersect(Set(j), sq))==0, s = setunion(s, Set(j)));); return (s);}
%o A165918 a(n) = {s = Set(); qnr = qnonr(n); for (j=1, #qnr, if (issquarefree(eval(qnr[j])), s = setunion(s, Set(qnr[j])));); return (#s);}  \\ _Michel Marcus_, Jul 23 2013
%Y A165918 Cf. A165916.
%K A165918 nonn
%O A165918 1,4
%A A165918 _Christopher Hunt Gribble_, Sep 30 2009
%E A165918 a(1) added by _Christopher Hunt Gribble_, Oct 05 2009