给定一个内角的等周多边形问题

王守恩

已知 n 边形一内角为90度,  周长为 n,  求它的最大面积。

lihpb00

显然为等腰直角三角形时面积最大，设直角边长为a，则有
2a+√2a=3→a=3(2-√2)/2
所以最大面积即a^2/2=9(3/2-√2)

lihpb00

显然就是正方形，面积最大值为1

hujunhua

应该是一个内接于圆的n边形。直角三角形占了圆的一半，斜边是圆的直径，另一半圆弧被剩下的n-3个点平分。
所以，是半个正方形与半个正(n-2)边形合成的图形。
比如五边形应该下图这样的：


本来想用canvas画图的，结果在发现本版块的附加选项中，HTML代码选项是灰色的，勾不上。

aimisiyou

hujunhua 发表于 2024-8-9 18:37
应该是一个内接于圆的n边形。直角三角形占了圆的一半，斜边是圆的直径，另一半圆弧被剩下的n-3个点平分。
...

是猜想还是证明？

王守恩

单位边长正n边形的面积等于$n/(4*tan(π/n))$,
[0.433012702, 1.000000000, 1.720477400, 2.598076212, 3.633912443, 4.828427124, 6.181824193, 7.694208842, 9.365639904, 11.19615242],
四舍五入后得到A134030——2024年3月11日
0, 1, 2, 3, 4, 5, 6, 8, 9, 11, 13, 15, 18, 20, 23, 26, 28, 32, 35, 38, 42, 46, 49, 54, 58, 62, 67, 71, 76, 81, 86, 92, 97, 103, 109, 115, 121, 127, 134, 140, 147, 154, 161, 168,
176, 183, 191, 199, 207, 215, 223, 232, 240, 249, 258, 267, 277, 286, 296, 306, 316, 326, 336, 346, 357, 368, 379, 390, 401, 412, 424, 436, 447, 459, 472, 484, 496, 509,
522, 535, 548, 561, 575, 588, 602, 616, 630, 644, 659, 673, 688, 703, 718, 733, 748, 764, 780, 796, 812, 828, 844, 860, 877, 894, 911, 928, 945, 963, 980, 998, 1016, ...}

1. Table[Round[Cot[Pi/n]*n/4], {n, 3, 112}]

复制代码

我们也可以仿此得到一个圆整序列提交上去。

mathe

容易看出除去直角三角形部分必须是等腰直角三角形，余下部分(n-2)条边必须相等而且内接一个圆，但是最优解直角不一定内接圆。

于是如图，必然有
\(\begin{cases}a=2R\sin(t)\\ \sqrt{2}d=2R|\sin((n-2)t)|\\(n-2)a+2d=n\end{cases}\)
而且显然要求\(0\lt 2(n-2)t \le 2\pi\)
所以我们可以得出
\(R=\frac{n}{2(n-2)\sin(t)+2\sqrt{2}\sin((n-2)t)}\)
而面积为
\(S=\frac{d^2}2+\frac{R^2}2((n-2)\sin(2t)-\sin(2(n-2)t))\)
比如n=7时，t的最优弧度应该略小于0.4, 最大面积大于3.5

不过面积表达式对于pari/gp有点太复杂，直接求导它算不出来，所以现在我还没有找到准确最值。
而对于n=10时，图会变成如下

可以看出对于n充分大时，应该取d=0,也就是使用正n-2边形.
比如n=6时我取了一个看上去比较靠近最值的配置，得到

(上图圆半径为0.949079053149547)

王守恩

谢谢 mathe！7楼错啦！！！我再试试(心里还是没底)。
已知三边形一内角为90度,  周长为3,  三边形面积最大值=0.386038969。
已知四边形一内角为90度,  周长为4,  四边形面积最大值=1.000000000。
已知五边形一内角为90度,  周长为5,  五边形面积最大值=1.701738567。
已知六边形一内角为90度,  周长为6,  六边形面积最大值=2.536906888。
已知七边形一内角为90度,  周长为7,  七边形面积最大值=3.514780784。
已知八边形一内角为90度,  周长为8,  八边形面积最大值=4.638479869。
已知九边形一内角为90度,  周长为9,  九边形面积最大值=5.909348892。
......
{0.386038969, 1.000000000, 1.701738567, 2.536906888, 3.514780784, 4.638479869, 5.909348892, 7.328062981, 8.894998347, 10.61038125, 12.47435594, 14.48701873, 16.64843639}
四舍五入后得到这样一串数。
0, 1, 2, 3, 4, 5, 6, 7, 9, 11, 12, 14, 17, 19, 21, 24, 27, 30, 33, 36, 39, 43, 46, 50, 54, 58, 63, 67, 71, 76, 81, 86, 91, 96, 102, 107, 113, 119, 125, 131, 138, 144, 151, 158,
164, 172, 179, 186, 194, 201, 209, 217, 225, 233, 242, 250, 259, 268, 277, 286, 296, 305, 315, 324, 334, 344, 355, 365, 375, 386, 397, 408, 419, 430, 442, 453, 465, 477,
489, 501, 513, 525, 538, 551, 564, 577, 590, 603, 617, 630, 644, 658, 672, 686, 701, 715, 730, 745, 760, 775, 790, 806, 821, 837, 853, 869, 885, 901, 918, 934, 951, ......}

1. Table[(Sqrt[8]Sin[3(n-2)Pi/(4(n-1))]Sin[(2n-5)Pi/(4(n-1))]+(n-2)Sin[(2n-5)Pi/(2(n-1))])n^2/(8(Sqrt[2]Sin[3(n-2)Pi/(4(n-1))]+(n-2)Cos[(2n-5)Pi/(4(n-1))])^2),{n,5,113}]

复制代码

特别地,  6楼: 9楼。
3/4 + 1/Pi=1.068309886183791

王守恩

1. Table[NMaximize[{((Sqrt[8]Sin[b]Sin[Pi/4+b]+(n-2)Sin[2 a])R^2)/2,Sqrt[8]R Sin[b]+2(n-2) R Sin[a]==n,(n-2)a+b==Pi,R>0,Pi>a>0,Pi>b>0},{R,a,b}],{n,3,9}]

复制代码

n=3, {0.386039, {R -> 0.936776, a -> 0.725159,  b -> 2.41643}},
n=4, {1.000000, {R -> 0.707107, a -> 0.785398,  b -> 1.5708}},
n=5, {1.70174,  {R -> 0.81867,   a -> 0.589039,  b -> 1.37448}},
n=6, {2.53691,  {R -> 0.949083, a -> 0.471232,  b -> 1.25667}},
n=7, {3.51478,  {R -> 1.08696,  a -> 0.392699,  b -> 1.1781}},
n=8, {4.63848,  {R -> 1.22856, a -> 0.336599,  b -> 1.122}},
n=9, {5.90935,  {R -> 1.37228, a -> 0.294524, b -> 1.07992}}}

1. NMaximize[{(2 Sqrt[2] R Sin[b]*R*Sin[\[Pi]/4 + b] + 5*R^2 Sin[2 a])/2,2*Sqrt[2] R Sin[b] + 5*2 R Sin[a] == 7, 5 a + b == \[Pi], R > 0, \[Pi]/2 > a > 0, \[Pi]/2 > b > 0}, {R, a, b}]

复制代码

{3.51478, {R -> 1.08696, a -> 0.392698, b -> 1.1781}}

1. NMaximize[{(2 Sqrt[2] R Sin[b]*R*Sin[\[Pi]/4 + b] + 5*R^2 Sin[2 a])/2,2*Sqrt[2] R Sin[b] + 5*2 R Sin[a] == 9, 5 a + b == \[Pi], R > 0, \[Pi]/2 > a > 0, \[Pi]/2 > b > 0}, {R, a, b}]

复制代码

{5.81015, {R -> 1.39752, a -> 0.392702, b -> 1.17808}}