|
Search: id:A155562
|
|
| |
|
| 0, 1, 2, 4, 8, 9, 16, 17, 18, 25, 32, 34, 36, 41, 49, 50, 64, 68, 72, 73, 81, 82, 89, 97, 98, 100, 113, 121, 128, 136, 137, 144, 146, 153, 162, 164, 169, 178, 193, 194, 196, 200, 225, 226, 233, 241, 242, 256, 257, 272, 274, 281, 288, 289, 292, 306, 313, 324, 328
(list; graph; listen)
|
|
|
OFFSET
|
1,3
|
|
|
COMMENT
|
Contains A155561 as a subsequence (obtained by restricting a,b,c,d to be nonzero). Also contains A000290 (squares) and A001105 (twice the squares) as subsequence.
Contribution from Warut Roonguthai (warut822(AT)yahoo.com), Oct 13 2009: (Start)
N is also of the form x^2 - 2y^2.
N = (p^2-q^2-2*r*s)^2+(r^2-s^2-2*p*q)^2
= (p^2+q^2-r^2-s^2)^2+2*(p*r-p*s-q*r-q*s)^2
= (p^2+q^2+r^2+s^2)^2-2*(p*r+p*s+q*r-q*s)^2
for some non-negative integers p, q, r, s. (End)
|
|
PROGRAM
|
(PARI) isA155562(n, /* use optional 2nd arg to get other analogous sequences */c=[2, 1]) = { for(i=1, #c, for(b=0, sqrtint(n\c[i]), issquare(n-c[i]*b^2) & next(2)); return); 1}
for( n=1, 500, isA155562(n) & print1(n", "))
|
|
CROSSREFS
|
Sequence in context: A035258 A115813 A048300 this_sequence A048715 A028982 A071601
Adjacent sequences: A155559 A155560 A155561 this_sequence A155563 A155564 A155565
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
M. F. Hasler (MHasler(AT)univ-ag.fr), Jan 24 2009
|
|
|
Search completed in 0.002 seconds
|
