OFFSET
1,3
COMMENTS
L(a/n) stands for generalized Legendre symbol, with value = 1 only if a is a quadratic residue of n.
FORMULA
a(n) = qrs2bincode(n)
MAPLE
A055094 := proc(n)
local i, z;
z := 0;
for i from 1 to n-1 do
z := z*2;
if (1 = numtheory[quadres](i, n)) then
z := z + 1;
fi;
od;
return z;
end proc:
MATHEMATICA
a[n_] := With[{rr = Table[Mod[k^2, n], {k, 1, n - 1}] // Union}, Boole[ MemberQ[rr, #]]& /@ Range[n - 1]] // FromDigits[#, 2]&; Array[a, 40] (* Jean-François Alcover, Mar 05 2016*)
PROG
(PARI) {a(n)=sum(k=1, n-1, 2^(k-1)*(0<sum(i=1, n-1, i^2%n==n-k)))} /* Michael Somos, Oct 14 2006 */
(Sage)
def A055094(n) :
Q = quadratic_residues(n)
z = 0
for i in (1..n-1) :
z = z*2
if i in Q : z += 1
return z
[A055094(n) for n in (1..33)] # Peter Luschny, Aug 08 2012
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 04 2000
STATUS
approved