|
MATHEMATICA
|
Invariant[a_List] := Module[{i=1, j=2, xMin, xMax, aa, n, invar=0, signs, x}, xMin=Abs[a[[i]]-a[[j]]]; xMax=a[[i]]+a[[j]]; aa=Complement[a, {a[[i]], a[[j]]}]; n=Length[aa]; Do[signs=(2*IntegerDigits[k, 2, n]-1); x=aa.signs; If[x>xMin&&x<xMax, invar+=Times@@signs], {k, 0, 2^n-1}]; invar]; Table[theSet=Table[N[Log[Prime[i]]], {i, 1, n}]; Invariant[theSet], {n, 3, 23, 2}]
|