# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a155709 Showing 1-1 of 1 %I A155709 #4 Jul 14 2012 11:32:23 %S A155709 22,33,73,88,97,118,121,132,150,166,177,193,198,214,225,241,249,262, %T A155709 292,294,297,313,321,337,352,358,388,393,409,433,438,441,454,457,472, %U A155709 484,502,528,537,550,577,582,600,601,649,657,664,673,681,694,708,726,753 %N A155709 Intersection of A154777 and A155716: N = a^2 + 2b^2 = c^2 + 6d^2 for some positive integers a,b,c,d. %C A155709 Subsequence of A155569 (where a,b,c,d may be zero). %o A155709 (PARI) isA155709(n,/* optional 2nd arg allows us to get other sequences */c=[6,2]) = { for(i=1,#c, for(b=1,sqrtint((n-1)\c[i]), issquare(n-c[i]*b^2) & next(2)); return);1} %o A155709 for( n=1,999, isA155709(n) & print1(n",")) %Y A155709 Cf. A000404, A154777, A092572, A097268, A154778, A155716, ... %K A155709 easy,nonn %O A155709 1,1 %A A155709 _M. F. Hasler_, Jan 25 2009 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE