3元序列为A000982
4元序列为A001973
6元序列为A001977
5元序列未见收录 k=4可以如下分析。
记4个数中中间两个数相邻的排列方案为b(n),中间两个数相隔1的排列数目为c(n),总排列数目为a(n)
那么a(n+1)=a(n)+b(n+1)+c(n+1)
b(n+1)-b(n)和c(n+1)-c(n)都很容易计算
这种方法应该可以扩充到更加大的k,只是这时数列b和c需要更加复杂的分析 14# hujunhua
要想继续走下去就很难了，这个需要在算法上作进一步的榨取。。。
我目前用的方法是chyanog 的改进版，即先生成整数n的所有K元不重复划分，然后再统计符合条件的划分的个数。
不知大大能否提供更优秀的算法来？ 设自然数n分拆成K个互不相同的正整数之和，有P(n,K)种情况，那么
P(n,1)=1
P(n,2)=Floor
P(n,3)=If]
自然数分拆为互异的整数之和的个数有
f(n)=P(n,1)+P(n,2)+P(n,3)+...+P(n,floor(2*n))
而f(n)又可以通过展开式 (1+x)(1+x2)...(1+xn)的xn系数得到 EvenQ是什么东西? 我曾以为这是业余爱好者容易做出成果的问题，看来是过于乐观，对问题的难度估计不足，特别是对计算量估计不足。
第1）问就这么难搞，后面3问岂不是更没戏呀？我最想进攻的是第4）问，但现在还没有着手。按顺序一步步来的做法是否必要还不得而知，有时间了准备试跳一下，直接攻第4)问。
我计算的7元序列前24项，记录一下：0,0,0,1,4,24,94,289,734,1656,3370,6375,11322,19138,30982, 48417,73316,108108,155646,219489,303748,413442, 554256,733005
计算起来确实慢，但改进空间很大，因为我是最 P 的计算者。 我算到了37个
{0, 0, 0, 1, 4, 24, 94, 289, 734, 1656, 3370, 6375, 11322, 19138, 30982, 48417, 73316, 108108, 155646, 219489, 303748, 413442, 554256, 733005, 957332, 1236222, 1579666, 1999265, 2507780, 3119876, 3851588, 4721127, 5748298, 6955424, 8366614, 10008857, 11911188}
再大的话就超过了软件默认设置的计算内存。。。
但要推出一个递推公式，这些数据还不够充分 25# qianyb
Q是 Question的缩写，意为 "ask a question"
所以，EvenQ 意思就是 "Is it Even?"-----这是偶数吗
这是Mathematica函数的特点，:) 原来如此,没用过Mathematica软件
27#wayne, 37项，多多足矣，咋算的嘛？膜拜一下。
用合情猜想的模式来推导递推公式，37项足够了。且记Y_k(x)=\prod_{i=1}^{k-1}(x^i-1)，k元均衡样本数序列的特征多项式记为S_k(x)。我们的合情猜想是：Y_k(x)|S_k(x)Y7(x)=1-x-x2 +x5+2x7 -x9-x10-x11-x12+2x14+x16-x19-x20+x21
构造bn=an-a1+n-a2+n+a5+n+2a7+n-a9+ n- a10+n-a11+n-a12+n+2a14+n+ a16+n-a19+n-a20+n+a21+n
bn满足的特征多项式就是剩余的因式。
可得到的前16项bn={3364, 3364,…, 3364},bn无疑是常数序列，特征多项式为x-1, 所以S7(x)=(x-1)Y7(x).
这就太好了，又回到了最初的轨道。也许k=奇数时，都是Sk(x)=(x-1)Yk(x).