OFFSET
1,1
COMMENTS
Prompted by a question from Shiv K. Gupta to the Number Theory mailing list.
LINKS
D. Broadhurst in reply to S. K. Gupta, Hypotenuses of Integral Right Triangles - again, NMBRTHRY list, Feb 24 2009
FORMULA
a(n) = min { A002144(k) | A145010(k) > A145010(k+1) > ... > A145010(k+n-1)}. - M. F. Hasler, Feb 26 2009
EXAMPLE
Comment from M. F. Hasler, Feb 24 2009:
The first sequence of 12 such primes is the one starting at a(12) =
931953301 = [27050, 14151]^2 ; area = 203431499448450450
931953389 = [26050, 15917]^2 ; area = 176325413694076350
931953397 = [25239, 17174]^2 ; area = 148267841956285170
931953409 = [24528, 18175]^2 ; area = 120941067830427600
931953433 = [30332, 3453 ]^2 ; area = 95111855933417940
931953437 = [23846, 19061]^2 ; area = 93319265825216970
931953469 = [30462, 2005 ]^2 ; area = 56429222392003890
931953509 = [30478, 1745 ]^2 ; area = 49241224048436490
931953569 = [30487, 1580 ]^2 ; area = 44651199683914740
931953637 = [22166, 20991]^2 ; area = 23594434443844350
931953709 = [30525 , 422 ]^2 ; area = 12000420304268550
931953733 = [21793, 21378]^2 ; area = 8346882442487610
PROG
(PARI) A144954( n, p=5, verbose=0, L=[0])={ for( i=1, n-1, while(( p=nextprime(p+2)) % 4 !=1, ); mn=sum2sqr_prime(p); L=if( L[i] > A=mn[1]*mn[2]*abs(mn[1]^2-mn[2]^2), concat( L, A), i=0; [A]) ); for( i=0, n-1, i & while( 1 != (p=precprime(p-2)) % 4, ); verbose & print( p" = " sum2sqr_prime(p) "^2 ; area = " L[n-i])); p} \\ M. F. Hasler, Feb 24 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
David Broadhurst, Feb 24 2009
STATUS
approved