OFFSET
1,1
COMMENTS
Alternate definition: Lists of P satisfying P = 2a^2 + 2b^2 = c^2 + d^2 (by parallelogram law) and a + b > max(c,d) (by triangle inequality) and a!=b (removes rhombus) and c!=d (removes rectangle) and gcd(a,b,c,d)=1.
EXAMPLE
a...b...c...d...P=2a^2+2b^2=c^2+d^2
4...7...7...9...130
6...7...7...11..170
5...10..9...13..250
7...9...8...14..260
8...9...11..13..290
7...11..12..14..340
8...11..9...17..370
6...13..11..17..410
10..11..9...19..442
9...13..10..20..500
MATHEMATICA
PMax=2000;
Do[Sqrt[P/2-a^2]//If[IntegerQ[#]&&GCD[a, #, c[[1]], c[[2]]]==1, P//Sow]&, {P, 2, PMax, 2}, {c, DeleteCases[PowersRepresentations[P, 2, 2], {x_, x_}]}, {a, (c[[2]]-c[[1]])/2+1, Sqrt[P]/2-1//Ceiling}]//Reap//Last//Last
CROSSREFS
KEYWORD
nonn
AUTHOR
Albert Lau, Oct 17 2014
STATUS
approved