%I #12 Dec 26 2023 12:11:34
%S 5,41,61,109,137,149,157,229,269,281,317,389,397,401,461,509,557,601,
%T 617,653,701,769,773,797,829,857,877,937,941,977,997,1009,1013,1061,
%U 1093,1097,1153,1181,1213,1217,1229,1249,1277,1289,1321,1409,1453,1489
%N Primes of the form x^2+y^2 such that L(x)*L(y) = -1, where L is the Liouville lambda-function A008836.
%e 5 = 1 + 2^2 and L(1)*L(2)= (1) *(-1) = -1.
%e 41 = 4^2 + 5^2 and L(4)*L(5)= (1)*(-1) = -1
%p isA174054 := proc(n)
%p local x,y ;
%p if not isprime(n) then
%p return false;
%p end if;
%p for x from 1 do
%p if x^2 > n then
%p return false;
%p end if;
%p if issqr(n-x^2) then
%p y := sqrt(n-x^2) ;
%p if A008836(x) * A008836(y) = -1 then
%p return true;
%p end if;
%p end if;
%p end do:
%p end proc:
%p for n from 1 to 2000 do
%p if isA174054(n) then
%p printf("%d,\n",n) ;
%p end if;
%p end do: # _R. J. Mathar_, Jul 09 2012
%t lambdaQ[{x_, y_}] := LiouvilleLambda[x]*LiouvilleLambda[y] == -1; Select[ Prime /@ Range[300], Or @@ lambdaQ /@ PowersRepresentations[#, 2, 2] &] (* _Jean-François Alcover_, Jul 30 2013 *)
%Y Cf. A174050, A002313.
%K nonn
%O 1,1
%A _Michel Lagneau_, Mar 06 2010
|