9# hujunhua
我是借助软件找出来的规律，
当时只是考虑大于1的奇数，
所以没能准确的提供表达式。。。
:L :L
借助软件，利用非常充分的数据验证了一下，得到三元，四元，五元，六元的均衡样本个数的递推公式分别如下：
(注释：K为奇数时，偶数n的样本个数为0，所以下面的3，5元反映的是去掉所有的偶数项的0之后的新序列)
三元: an=2an-1-2an-3+an-4
初始值{0, 1, 2, 5}
四元：an=2an-1-an-3-an-4+2an-6-an-7
初始值{0, 0, 0, 1, 1, 3, 5}
五元：an=2an-1-an-3-2an-5+2an-6+an-8-2an-10+an-11
初始值{0, 0, 1, 3, 12, 32, 73, 141, 252, 414, 649}
六元：an=an-1+2an-2-an-3-an-4-an-5-an-6+2an-8+2an-9-an-11-an-12-an-13-an-14+2an-15+an-16-an-17
初始值{0, 0, 0, 0, 0, 1, 1, 4, 8, 18, 32, 58, 94, 151, 227, 338, 480}
比如，六元均衡样本，用mathematica进行迭代：
LinearRecurrence[{1, 2, -1, -1, -1, -1, 0, 2, 2, 0, -1, -1, -1, -1, 2,1, -1}, {0, 0, 0, 0, 0, 1, 1, 4, 8, 18, 32, 58, 94, 151, 227, 338, 480}, 100]
{0, 0, 0, 0, 0, 1, 1, 4, 8, 18, 32, 58, 94, 151, 227, 338, 480, 676, 920, 1242, 1636, 2137, 2739, 3486, 4370, 5444, 6698, 8196, 9926, 11963, 14293, 17002, 20076, 23612, 27594, 32134, 37212, 42955, 49341, 56512, 64444, 73294, 83036, 93844, 105690, 118765, 133037, 148718, 165772, 184430, 204654, 226694, 250510, 276373, 304239, 334402, 366814, 401792, 439284, 479632, 522780, 569093, 618513, 671430, 727782, 787986, 851974, 920192, 992568, 1069575, 1151137, 1237756, 1329352, 1426456, 1528984, 1637498, 1751908, 1872809, 2000105, 2134424, 2275666, 2424490, 2580792, 2745266, 2917802, 3099129,3289131, 3488574, 3697336, 3916220, 4145098, 4384810, 4635224, 4897217, 5170651, 5456444, 5754450, 6065628, 6389826, 6728044}
发主帖前，偶只粗略算了一下3元均衡样本数，当时也是去掉了偶零项，得到通项公式Ceiling[(n-1)2/2]。为了消除取整函数，形成了递推公式a+a=n2-n+1。
根据wayne所给的递推公式，相应的特征多项式是：
k=3时，x^4 - 2 x^3 + 2 x - 1=(x-1)^3(x+1)=(x-1)^2(x^2-1)
k=4时，x^7 - 2 x^6 + x^4 + x^3 - 2 x + 1=(x-1)^4(x+1)(x^2+x+1)=(x-1)^2(x^2-1)(x^3-1)
k=5时，x^11 - 2 x^10 + x^8 + 2 x^6 - 2 x^5 - x^3 + 2 x - 1=(x-1)^5(x+1)^2(x^2+1)(x^2+x+1)=(x-1)^2(x^2-1)(x^3-1)(x^4-1)
这些特征多项式的不可约因式都是分圆多项式，并且刚好可凑成1~(k-1)阶分圆方程之积，仅有因式(x-1)多了一重。多么美妙的规律啊，我满怀期望地分解k=6的特征多项式，
1 - x - 2 x^2 + x^3 + x^4 + x^5 + x^6 - 2 x^8 - 2 x^9 + x^11 + x^12 + x^13 + x^14 - 2 x^15 - x^16 + x^17
=(x-1 )^6 (1 + x)^3 (1 + x^2) (1 + x + x^2) (1 + x + x^2 + x^3 + x^4)
=(x-1)(x^2-1)^2(x^3-1)(x^4-1)(x^5-1)
结果这回多了一重的是(x^2-1).虽然仍然只有分圆多项式，仍然能够凑成1~(k-1)阶分圆方程之积，但是重数的一点点不同仍然不啻一瓢冷水当头浇下。
wayne, 检查一下你的6元公式，能否将特征多项式降到16次。
能否再接再厉，给出7元、8元、9元、…特征方程来。
呵呵，一刷新才发现 大大捷足先登了啊
我再整理整理
六元的我感觉已经比较充分了。
有一点疑问，奇数元的我们都忽略了偶零项，这会不会就是不和谐的根源？
不对，把偶数项加进去的话就更乱了
偶零项加进去的话，公式恐怕会很怪。暂时不加为好。
既然k为奇数时忽略了偶零项，那么可以考虑k为偶数时忽略n偶项。
也就是，可以暂限于讨论Range(2n-1)的均衡样本数。
Range(2n-1)的总体平均值为n，我发主帖时，第一版其实就是Range(2n-1)，后来为了一般性改成了Range(n).
14# hujunhua
k=7时,公式会不会是(x-1)(x^2-1)(x^3-1)^2(x^4-1)(x^5-1)(x^6-1)
也许。即使有出入，也不会很大，我坚信不可约因式都是分圆多项式这一点不会错，虽然现在还没进入证明阶段。我期望这些分圆多项式可以凑成1~(k-1)阶分圆方程之积，并且重数有简明的分布规律。