自然数前段的均衡样本

hujunhua

又申请了6个序列，A390082~390087，
A390082，即127#中的d(n)的奇数项。
A390083，即127#中的d(n)的偶数项
A390084，既约平衡子集的最大长度，奇数项。我提交了14项，编辑给延长到了第50项。
A390085，既约平衡子集的最大长度，偶数项。我提交了14项，与其它三个既有序列还不能区分。

前面三个序列很快就通过了，第四个是本月3号提交的，今天才通过。因为我提交的项数不足以与既有的序列相区别。这一点在编辑时系统自动有警示，但我也没办法。现在问题没解决，先给通过了。

A390086， 打算给A390084的简并，即序列发生阶跃的地方.  
A390087， 打算给A390085的简并，即序列发生阶跃的地方.  

我没打算把A390082与A390083合并申请一个新的序列，也不打算合并A390084和A390085。

hujunhua

从127#开始，我们将自然数前段平移成了对称的零和整数集，例如自然数前段{1..2n+1}平移为{-n..n}。
随之，平衡样本转化成了零和整数集。
有关问题的阶段性结果包括：
1、{-n..n}有多少个本原零和子集，0~37项见A389803
2、{-n..n}有多少个完全划分？0~28项见A389806
3、{-n..n}包含 n 但不包含 -n的本原零和集有多少个，1~37项见A390082
     A390082就是A389803的向前差分，通过计算A390082，累计即得A389803。
4、{-n..n}的本原零和子集的最大长度，0~50项见A390084
     这是为编程计算问题3服务的。如果不考虑它而直接计算A390082(n)时，需要搜索的子集长度范围就是3~n, 但是观察A390084的前50项，显示A390084(n)≈`\sqrt{6n}`, 也就是在计算计算A390082(n)时，需要搜索的子集长度不超过`\sqrt{6n}`。当n较大时，这个差别就很大了。
5、n元本原零和集最小规模，也就是问题4的逆问题。1~16项见A390086

A390084的前50项：1, 2, 2, 3, 4, 4, 5, 5, 6, 6, 6, 7, 7, 7, 8, 8, 8, 8, 9, 9, 9, 10, 10, 10, 10, 11, 11, 11, 11, 11, 12, 12, 12, 12, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 16
前0~14项是本论坛贡献的，第15~50项是OEIS的一位编辑贡献的。
忽略那些平台项，发生增长的项及相应的最长本原零和子集如下：
  a(0) = 1: {0}
  a(1) = 2: {-1, 1}
  a(3) = 3: ±{3, -2, -1}
  a(4) = 4: ±{4, -3, -2, 1}
  a(6) = 5: ±{6, -5, -4, 2, 1}
  a(8) = 6: ±{8, -7, 6, -5, -4, 2}
  a(11) = 7: ±{11, 10, -9, -8, -7, 2, 1}
  a(14) = 8: ±{14, 13, -12, -11, -10, 3, 2, 1}
  a(18) = 9: ±{18, 17, 16, -15, -14, -13, -12, 2, 1}
  a(21) = 10: ±{21, 20, 19, -18, -17, -16, -15, 3, 2, 1}
  a(25) = 11: ±{25, 24, 23, -22, -21, -20, -19, 4, 3, 2, 1}
  a(30) = 12: ±{30, 29, 28, 27, -26, -25, -24, -23, -22, 3, 2, 1}
  a(34) = 13: ±{34, 33, 32, 31, -30, -29, -28, -27, -26, 4, 3, 2, 1}
  a(39) = 14: ±{39, 38, 37, 36, -35, -34, -33, -32, -31, 5, 4, 3, 2, 1}
  a(45) = 15: ±{45, 44, 43, 42, 41, -40, -39, -38, -37, -36, -35, 4, 3, 2, 1}
  a(50) = 16: ±{50, 49, 48, 47, 46, -45, -44, -43, -42, -41, -40, 5, 4, 3, 2, 1}

我们把它取反函数，就得到A390086，1~16项如下：0, 1, 3, 4, 6, 8, 11, 14, 18, 21, 25, 30, 34, 39, 45, 50
即 给定长度的本原零和整数集，它的最大元素(按绝对值）所能达到的最小值。观察这前16项，我们得到\[A390086(n)=\lfloor\frac{n(n+3)}6\rfloor=\begin{cases}\D\frac{n(n+3)}6, \text{ if }n\equiv0\\\D\frac{n(n+3)-4}6,\text{ if }n\neq0\end{cases}\pmod{3}\]
不过，上面的公式不是来自数据拟合，而是从上面观察到的最大子集构成模式中推导出来的。
从a(11)开始，最长子集具有以下构成模式：\[\{a_k, ..., a_1, -b_{k+1}, ..., -b_1, c_m, ...,c_1\}\]
忽略{bi}的负号，整个序列按降序排列，分为正负正的三段，并且 \[\sum_{i=1}^k a_i+\sum_{i=1}^m c_i=\sum_{i=1}^{k+1} b_i, b_1>\sum_{i=1}^m c_i\]易证这种模式必是本原零和集。

hujunhua

假定它有较小的零和真子集，必定是
total(subset(a))+total(subset(c))=total(subset(b)), 即 a 段的一个子段加 c 段的一个子段等于 b 段的一个子段，
并且length(subset(b))>length(subset(a))， 即 b 段的那个子段长于 a 段的那个子段，因为 b 段的数较小嘛，得靠长度来凑。
但是，subset(b)哪怕只比 subset(a)多一个数，假定 subset(b)=s+1， subset(a)=s, 都有 total(subset(a))+total(subset(c))<total(subset(b)).
因为 total(subset(b))≥b_1+...+b_{s+1}, total(subset(a))≤a_k+...+a_{k-s+1}, b_{s+2}+...+b_{k+1}<a_{k-s}+...+a_1
由a+c=b, b_{s+2}+...+b_{k+1}<a_{k-s}+...+a_1 → a_k+...+a_{k-s+1}+c<b_1+...+b_{s+1} →  total(subset(a))+total(subset(c))<total(subset(b)).

hujunhua

为了使 a_k尽可能小，a 段, b 段, c 段都应该是连续的，b 段与 a 段宜接近， c 段的数字宜尽可小，所以可以上述模式可以更具体：\[
\{a, a-1, ..., a-k+1, -(a-k), -(a-k-1), ..., -(a-2k), c, c-1, ..., 1\}
\]由零和可得\[
(2a-k+1)k/2+c(c+1)/2=(2a-3k)(k+1)/2\tag{1}
\]给定长度 n, 所以\[2k+1+c=n\tag2\]
由以上两式解二次式最值可得\[
k=\lfloor\frac{n}3\rfloor, c=\lfloor\frac{n}3\rfloor+r-1,a=\lfloor\frac{3n(n+3)}6\rfloor
\]这里 r 是 n 除以 3 的余数。