login
A229769
Not of the form x^2 + P*y^2 for integers x > 0, y > 1, P a prime.
1
1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 14, 15, 18, 20, 23, 25, 26, 30, 35, 38, 39, 40, 42, 47, 50, 55, 58, 62, 65, 71, 74, 78, 83, 87, 90, 95, 98, 106, 110, 119, 122, 130, 138, 143, 146, 155, 158, 159, 167, 170, 182, 186, 190, 195, 203, 210, 215, 218, 222, 227, 230, 231, 250
OFFSET
1,2
COMMENTS
There are no terms <= 10^9 greater than a(3036) = 69808035, leading to the conjecture that the sequence is finite. If true, every sufficiently large number is expressible as x^2 + P*y^2 with x > 0, y > 1.
EXAMPLE
Since 17 can be expressed as 3^2 + 2 * 2^2, it is not in the sequence.
No such expression exists for 18, hence it is in the sequence.
Since 19 can be expressed as 1^2 + 2 * 3^2, it is not in the sequence.
PROG
(PARI) test(n)={local(z, x, p, y); for(x=1, sqrtint(n), z=n-x^2; p=core(z); y=core(z, 1)[2]; if(isprime(p)&&y>1, return(1))); }
for(n=1, 300, if(test(n)==0, print1(n", ")))
CROSSREFS
Supersequence of A074885 (x^2 + M*y^2, x > 0, y > 1, M > 0).
Supersequence of A212709 (not of the form p*c^2 + b^2, with p prime and c and b nonzero integers).
Sequence in context: A044919 A011871 A377357 * A130676 A039227 A257727
KEYWORD
nonn
AUTHOR
Chris Boyd, Sep 29 2013
STATUS
approved