[seqfan] Re: defn of A167415
Don Reble
djr at nk.ca
Sat Jun 15 02:57:11 CEST 2019
> How is A167415 defined?
It looks like,
M such that (for all x,y except x=y=0) (x+y)^2+3xy isn't == 0 mod M.
At any rate, this Python jot produces the sequence data.
#!/usr/bin/python3
for mod in range(2,300):
halfmod = (mod // 2) + 1
for x in range(halfmod):
xsq = x * x
for y in range(max(x,1),halfmod):
sumsq = y * y + xsq
xy3 = x * y * 3
val = (sumsq + xy3) % mod
if val == 0: break
val = (sumsq - xy3) % mod
if val == 0: break
if val == 0: break
if val != 0: print(mod)
> %C Prime numbers of this sequence are congruent to {2,3} modulo 5.
BTW, that formula equals (x^2 + 5xy + y^2);
that may be why modulo 5 matters.
--
Don Reble djr at nk.ca
More information about the SeqFan
mailing list