本帖最后由 王守恩 于 2023-3-6 09:43 编辑
(公式1): Table[FindInstance[{(x^2 + x + y^2 + y)/(m n - 2 m - n + 4)
== (2 m)/(3 + Cos[n \[Pi]]), y >= x > 0}, {x, y}, Integers, 1], {n, 3, 9}, {m, 1, 9}]
{{{}, {}, {{x -> 2, y -> 2}}, {}, {}, {{x -> 3, y -> 5}}, {}, {{x -> 5, y -> 6}}, {}}, {{},
{{x -> 1, y -> 1}}, {}, {}, {}, {{x -> 2, y -> 5}}, {}, {}, {}}, {{}, {}, {{x -> 3, y -> 3}},
{{x -> 1, y -> 6}}, {}, {{x -> 3, y -> 9}}, {{x -> 5, y -> 10}}, {{x -> 1, y -> 13}},
{}}, {{}, {}, {}, {}, {}, {}, {}, {{x -> 5, y -> 9}}, {}}, {{}, {{x -> 1, y -> 3}}, {{x -> 2, y -> 5}},
{{x -> 3, y -> 7}}, {{x -> 4, y -> 9}}, {{x -> 2, y -> 12}}, {{x -> 6, y -> 13}},
{{x -> 7, y -> 15}}, {{x -> 8, y -> 17}}}, {{}, {{x -> 1, y -> 2}}, {}, {{x -> 4, y -> 4}},
{}, {{x -> 2, y -> 9}}, {}, {{x -> 4, y -> 12}}, {}}, {{}, {{x -> 2, y -> 3}},
{{x -> 2, y -> 6}}, {{x -> 1, y -> 9}}, {}, {{x -> 3, y -> 14}}, {{x -> 1, y -> 17}}, {}, {}}}
(公式2): Ceiling[FractionalPart[(Table[Dimensions[FindInstance[
{(x^2 + x + y^2 + y)/(m n - 2 m - n + 4) == (2 m)/(3 + Cos[n \[Pi]]),
y >= x > 0}, {x, y}, Integers, 1]] // Total, {n, 3, 9}, {m, 1, 9}] + 1)^-1]]
{{0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 0},(3 边形数)
{0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},(4 边形数)
{0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1},(5 边形数)
{0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0},(6 边形数)
{0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},(7 边形数)
{0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0},(8 边形数)
{0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0}}(9 边形数)
(公式2)只是在(公式1)的外面套个框,(公式2)太臭了! 能不能让(公式2)漂亮一点。
|