1,5

a(5)=6 because there are 6 essentially different ways: {1, 4, 9, 64, 49, 100, 81, 16, 25, 36}, {1, 4, 49, 64, 9, 100, 81, 16, 25, 36}, {1, 16, 81, 100, 9, 4, 49, 64, 25, 36}, {1, 16, 81, 100, 9, 64, 49, 4, 25, 36}, {1, 16, 81, 100, 49, 4, 9, 64, 25, 36} and {1, 16, 81, 100, 49, 64, 9, 4, 25, 36}

$RecursionLimit=500; try[lev_] := Module[{t, j}, If[lev>2n, (*then make sure the sum of the first and last is prime*) If[PrimeQ[soln[[1]]^2+soln[[2n]]^2]&&soln[[2]]<=soln[[2n]], (*Print[soln]; *) cnt++ ], (*else append another number to the soln list*) t=soln[[lev-1]]; For[j=1, j<=Length[s[[t]]], j++, If[ !MemberQ[soln, s[[t]][[j]]], soln[[lev]]=s[[t]][[j]]; try[lev+1]; soln[[lev]]=0]]]]; For[lst={}; n=1, n<=7, n++, s=Table[{}, {2n}]; For[i=1, i<=2n, i++, For[j=1, j<=2n, j++, If[i!=j&&PrimeQ[i^2+j^2], AppendTo[s[[i]], j]]]]; soln=Table[0, {2n}]; soln[[1]]=1; cnt=0; try[2]; AppendTo[lst, cnt]]; lst

Cf. A051252, A073451.

Sequence in context: A145223 A365979 A219948 * A085511 A187525 A187696

Adjacent sequences: A072126 A072127 A072128 * A072130 A072131 A072132

nonn,more

Santi Spadaro, Jun 25 2002

Corrected and extended by T. D. Noe, Jul 03 2002

a(13) corrected and a(14)-a(15) from Sean A. Irvine, Sep 03 2024

approved