%I #7 Aug 07 2015 03:06:02
%S 1,7,17,101,241,157,191,1019,311,277,2689,1723,1297,1087,883,809,7433,
%T 3169,1907,6709,2909,5843,3121,6871,2557,4759,3079,3373,12899,6703,
%U 5003,4073,4241,15881,8783,9173,25457,14593,7603,9007,46751,21521,36107
%N Diagonal of triangle in A082737.
%p A082738 := proc(nmax) local T,a,n,r,i,rsum,c,j ; T := [1,2,7] ; a := [1,7] ; n := 3 ; i := 1 ; while nops(a)< nmax do r := [] ; for c from 1 to n-1 do while ithprime(i) in T or ithprime(i) in r do i:= i+1 ; od ; r := [op(r),ithprime(i)] ; i:= i+1 ; od ; j := i+1 ; rsum := sum(op(k,r),k=1..nops(r)) ; while not issqr( rsum+ithprime(j)) do j := j+1 ; od ; a := [op(a),ithprime(j)] ; r := [op(r),ithprime(j)] ; T := [op(T),op(r)] ; n := n+1 ; od ; RETURN(a) ; end: a := A082738(80) : for n from 1 to nops(a) do printf("%d,",op(n,a)) ; od ; # _R. J. Mathar_, Nov 12 2006
%Y Cf. A082737, A082739, A082740.
%K nonn
%O 1,2
%A _Amarnath Murthy_, Apr 14 2003
%E More terms from _R. J. Mathar_, Nov 12 2006
|