x6/5=c6; x7/5=c7; x8/5=c8; x9/5=c9; x10/5=c10; x11/5=c11; x12/5=c12; x13/5=c13;
@gin(c1); @gin(c2); @gin(c3); @gin(c4); @gin(c5); @gin(c6); @gin(c7); @gin(c8); @gin(c9); @gin(c10); @gin(c11); @gin(c12); @gin(c13);
(x1/2)*(x1/2)+(y1/2)*(y1/2)=325*325;
(x2/2)*(x2/2)+(y1/2+y2)*(y1/2+y2)=325*325;
(x3/2)*(x3/2)+(y1/2+y2+y3)*(y1/2+y2+y3)=325*325;
(x4/2)*(x4/2)+(y1/2+y2+y3+y4)*(y1/2+y2+y3+y4)=325*325;
(x5/2)*(x5/2)+(y1/2+y2+y3+y4+y5)*(y1/2+y2+y3+y4+y5)=325*325;
(x6/2)*(x6/2)+(y1/2+y2+y3+y4+y5+y6)*(y1/2+y2+y3+y4+y5+y6)=325*325;
(x7/2)*(x7/2)+(y1/2+y2+y3+y4+y5+y6+y7)*(y1/2+y2+y3+y4+y5+y6+y7)=325*325; (x8/2)*(x8/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8)*(y1/2+y2+y3+y4+y5+y6+y7+y8)=325*325; (x9/2)*(x9/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9)=325*325;
(x10/2)*(x10/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10)=325*325;
(x11/2)*(x11/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11)=325*325;
(x12/2)*(x12/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12)=325*325;
(x13/2)*(x13/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12+y13)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12+y13)=325*325;
h1=y1+y2+y3+y3+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12+y13; 求得结果为：
Local optimal solution found.
Objective value: 319923.9 Objective bound: 319923.9 Infeasibilities: 0.1916722E-02 Extended solver steps: 123 Total solver iterations: 24684
Variable Value Reduced Cost X1 640.0000 0.000000 Y1 113.5782 X2 620.0000 Y2 40.80702 X3 595.0000 Y3 33.24253 X4 565.0000 Y4 29.84350 X5 530.0000 Y5 27.46674 X6 495.0000 Y6 22.49003 X7 455.0000 Y7 21.45751 X8 410.0000 Y8 20.09398 X9 365.0000 Y9 16.73105 X10 315.0000 Y10 15.36493 X11 255.0000 Y11 14.65968 X12 190.0000 Y12 11.85934 X13 110.0000 Y13 9.506939 C1 128.0000 C2 124.0000 C3 119.0000 C4 113.0000 C5 106.0000 C6 99.00000 C7 91.00000 C8 82.00000 C9 73.00000 C10 63.00000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -4.402540 -11.02568 8.643625 9.237668 -28.18703 10.09900 5.968964 -18.04228 2.348676 12.55633
C11 51.00000 -7.984821 C12 38.00000 3.669596 C13 22.00000 -0.6303581 H1 443.5865 0.000000
Row Slack or Surplus Dual Price 1 319923.9 1.000000 2 10.00000 0.000000 3 20.00000 0.000000 4 25.00000 5 30.00000 6 35.00000 7 35.00000 8 40.00000 9 45.00000 10 45.00000 11 50.00000 12 60.00000 13 65.00000 14 80.00000 15 90.00000 16 0.000000 17 0.000000 18 0.000000 19 0.000000 20 0.000000 21 0.000000 22 0.000000 23 0.000000 24 0.000000 25 0.000000 26 0.000000 27 0.000000 28 0.000000 29 0.000000 30 0.000000 31 0.000000 32 0.000000 33 0.000000 34 0.000000 35 0.000000 36 0.000000 37 0.000000 38 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 4.402540 11.02568 -8.643625 -9.237668 28.18703 -10.09900 -5.968964 18.04228 -2.348676 -12.55633 7.984821 -3.669596 0.6303581 0.3521803 0.2561577 0.2292901 0.2178213 0.1860229 0.1898984 0.1938849 0.1784366 0.1859279 0.2110548
39 -0.1147564E-06 0.2174305 40 -0.1916722E-02 0.2573958 41 0.000000 0.3434148 42 0.000000 0.000000
附录三：
针对级数为14
max=x1*y1+2*x2*y2+2*x3*y3+2*x4*y4+2*x5*y5+2*x6*y6+2*x7*y7+2*x8*y8+2*x9*y9+2*x10*y10+2*x11*y11+2*x12*y12+2*x13*y13+2*x14*y14; x1<=650; x1>=x2; x2>=x3; x3>=x4; x4>=x5; x5>=x6; x6>=x7; x7>=x8; x8>=x9; x9>=x10; x10>=x11; x11>=x12; x12>=x13; x13>=x14; x14>=20; x1/5=c1; x2/5=c2; x3/5=c3; x4/5=c4; x5/5=c5; x6/5=c6; x7/5=c7; x8/5=c8; x9/5=c9; x10/5=c10; x11/5=c11; x12/5=c12; x13/5=c13; x14/5=c14; @gin(c1); @gin(c2); @gin(c3); @gin(c4); @gin(c5); @gin(c6);
@gin(c7); @gin(c8); @gin(c9); @gin(c10); @gin(c11); @gin(c12); @gin(c13); @gin(c14);
(x1/2)*(x1/2)+(y1/2)*(y1/2)=325*325;
(x2/2)*(x2/2)+(y1/2+y2)*(y1/2+y2)=325*325;
(x3/2)*(x3/2)+(y1/2+y2+y3)*(y1/2+y2+y3)=325*325;
(x4/2)*(x4/2)+(y1/2+y2+y3+y4)*(y1/2+y2+y3+y4)=325*325;
(x5/2)*(x5/2)+(y1/2+y2+y3+y4+y5)*(y1/2+y2+y3+y4+y5)=325*325;
(x6/2)*(x6/2)+(y1/2+y2+y3+y4+y5+y6)*(y1/2+y2+y3+y4+y5+y6)=325*325;
(x7/2)*(x7/2)+(y1/2+y2+y3+y4+y5+y6+y7)*(y1/2+y2+y3+y4+y5+y6+y7)=325*325; (x8/2)*(x8/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8)*(y1/2+y2+y3+y4+y5+y6+y7+y8)=325*325; (x9/2)*(x9/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9)=325*325;
(x10/2)*(x10/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10)=325*325;
(x11/2)*(x11/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11)=325*325;
(x12/2)*(x12/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12)=325*325;
(x13/2)*(x13/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12+y13)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12+y13)=325*325;
(x14/2)*(x14/2)+(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12+y13+y14)*(y1/2+y2+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12+y13+y14)=325*325;
h1=y1+y2+y3+y3+y3+y4+y5+y6+y7+y8+y9+y10+y11+y12+y13+y14;
求得结果为：
Local optimal solution found.
Objective value: 320751.3 Objective bound: 320751.3 Infeasibilities: 0.7412484E-02 Extended solver steps: 76 Total solver iterations: 22205
Variable Value Reduced Cost X1 640.0000 0.000000 Y1 113.5782 0.000000 X2 620.0000 0.000000