|

楼主 |
发表于 2021-7-10 10:24
|
显示全部楼层
本帖最后由 天山草 于 2021-7-14 23:01 编辑
可以使用常新德老师(河南永城职业学院)在论文《有重复元素的圆排列和环排列计数》中给出的公式,在使用这个公式时要注意,
就是每一种颜色珠子的数目有奇有偶,当奇数的个数多于 2 个时,公式中的 M=0,而不必按公式进行计算,如果计算就出错了。
只有当奇数的个数为0、1、2 时,才能按公式计算 M 的值。
以下是按常新德公式计算的 mathematica 程序:
- Clear["Global`*"];
- m = 3; n[1] = 4; n[2] = 4; n[3] = 1;
- k = 1;(*n[1]、n[2]、n[3] 中共有 k 个奇数 *)
- s = \!\(
- \*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(m\)]\(n[i]\)\);
- Q = (Sum[EulerPhi[d] *(s/d)!/\!\(
- \*UnderoverscriptBox[\(\[Product]\), \(i = 1\), \(m\)]\(\((
- \*FractionBox[\(n[i]\), \(d\)])\)!\)\), {d,
- Divisors[GCD[n[1], n[2], n[3]]]}])/s;(*圆排列数*)
- If[k < 3, M = (\!\(
- \*UnderoverscriptBox[\(\[Sum]\), \(i = 1\), \(m\)]\(\[LeftFloor]
- \*FractionBox[\(n[i]\), \(2\)]\[RightFloor]\)\))!/\!\(
- \*UnderoverscriptBox[\(\[Product]\), \(i = 1\), \(m\)]\(\[LeftFloor]
- \*FractionBox[\(n[i]\), \(2\)]\[RightFloor]!\)\),
- M = 0];(*圆排列中的对称排列数*)
- \[CapitalPhi] = ((Q + M)/2)(*环排列数*)
复制代码
计算结果为 38。
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有帐号?注册
x
|