login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A155560 Intersection of A000404 and A092572: N = a^2 + b^2 = c^2 + 3d^2 with a,b,c,d>0 16

%I #7 Jan 11 2015 10:51:07

%S 13,37,52,61,73,97,100,109,117,148,157,169,181,193,208,229,241,244,

%T 277,292,313,325,333,337,349,373,388,397,400,409,421,433,436,457,468,

%U 481,541,549,577,592,601,613,628,637,657,661,673,676,709,724,733,757,769

%N Intersection of A000404 and A092572: N = a^2 + b^2 = c^2 + 3d^2 with a,b,c,d>0

%C Nonsquare terms of A155563. - _Joerg Arndt_, Jan 11 2015

%e a(1)=13 is the least number that can be written as A+B and C+3D where A,B,C,D are positive squares (namely 13 = 2^2 + 3^2 = 1^2 + 3*2^2).

%e a(2)=37 is the second smallest number which figures in A000404 and in A092572 as well.

%o (PARI) isA155560(n /* omit optional 2nd arg for the present sequence */, c=[3,1]) = { for(i=1,#c,for(b=1,sqrtint((n-1)\c[i]),issquare(n-c[i]*b^2)&next(2));return);1}

%o for( n=1,10^3, isA155560(n) & print1(n","))

%o (PARI) is(n)=!issquare(n) && #bnfisintnorm(bnfinit(z^2+z+1), n) && #bnfisintnorm(bnfinit(z^2+1), n);

%o select(n->is(n), vector(1500,j,j)) \\ _Joerg Arndt_, Jan 11 2015

%K nonn

%O 1,1

%A _M. F. Hasler_, Jan 24 2009

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 05:48 EDT 2024. Contains 371265 sequences. (Running on oeis4.)