天平秤球的最佳规格化算法.docVIP

天平秤球的最佳规格化算法 mjwu mjwu1940@126.com 天平秤球问题作为一个智力锻炼过程，在网上讨论甚为热烈，经年不息。大多集中在12球问题上。叙文繁多，方法各别，有文本论述，有向量法，有表格式，更有图文并茂者。可谓五花八门，但秤球过程烦杂，实难记忆。本文将提供一种针对任意球数m的行之有效的快速查寻方法，过程规一，方便记忆；且适合于计算机拟人化查寻。 问题：m个球中含有一个次球（或轻或重），通过天平n次秤球，查寻次球并判定其轻重。 两个关于天平秤球的整数函数 给定n，通过天平n次秤球能实现查寻次球并判定其轻重的最大球数m，记为M(n), m=M(n)=(3^n-3)/2 (^表示n次方)。 如n=3,m=12;n=5,m=120;…… 给定m，要实现查寻次球并判定其轻重的天平秤球最小次数n，记为N(m), n=N(m)=【log(2m+3)/log3】 (此处【x】为进位取整函数: 即当x为整数时，【x】=x;当x不为整数时，【x】=[x]+1，即x有小数点位时进1。x的进位取整是不小于x的最小整数)。如m=39,n=4;m=40,n=5。 （证明从略，见作者的“天平秤球中的信息论”）。 二、天平秤球中的准三等分法(WKL) 天平秤球过程中需要反复对疑重球数z、疑轻球数q、待定球数d（不知其为轻或为重）分别分成为三个数放置在天平的左盘、右盘及不放盘。下面先介绍一种整数的准三等分法。 准三等分法（WKL）包括K三分法及L三分法。 1，K三分法。 在天平秤球过程中，在无先验知识条件下，一个球放在左盘，或右盘，或不放盘应是等可能的；所以，一组待处理的球应尽可能地三等分，并置于左盘、右盘与不放盘。设待处理球数p，p=pa+pb+pc，p整三等分为整数pa,pb,pc；如多余1，则pa+1;如多余2，则pa+1,pc+1。如，6=2+2+2；7=3+2+2；8=3+2+3；9=3+3+3。pa,pb,pc分别对应于左盘球数，右盘球数，不放盘球数。 分析与实验表明，在秤球分盘的过程中始终保持K三分法（pa=pc=pb）能有效降低秤球次数，达到快速查寻次球并判定其轻重之目的。（当然也可反过来，pb=pc=pa）。 2，L三分法。在天平秤球时，第一次秤球因无参考正品球，应使左盘球数=右盘球数。 设待处理球数p, p=pa+pb+pc，p整三等分为整数pa,pb,pc；如多余1，则pc+1;如多余2，pa+1,pb+1。如，6=2+2+2；7=2+2+3；8=3+3+2；9=3+3+3。 在天平秤球的实际过程中，第一次秤球时采用L三分法（左盘球数=右盘球数）；所有的其他次秤球（非第一次）时都采用K三分法。 天平秤球最佳规化算法的基本思想 对于任意个球数m（编号为1，…，m），初始疑重球数z=0、疑轻球数q=0、待定球数d=m。估判向量Y简记为（z,q,d）。初始为（0,0，m）。 分盘。除第一次秤球采用L三分法分解d=da+db+dc外，其余的各次秤球，无论是（z=q=0,d0），或者（z+q0,d=0），都采用K三分法，分解d=da+db+dc或者分解z=za+zb+zc, q=qa+qb+qc，并分别放置于天平的左盘、右盘及不放盘（注意：当秤球的左盘与右盘的球数不等时，可适当补充参考正品球） 秤球。结果：或平衡，或左重右轻，或左轻右重。 3, 更新估判向量Y，即（z,q,d）。根据原来的估判向量Y，本次秤球的分布状态P以及秤球的结果更新估判向量Y（如，在第一次称球后，平衡时d=dc,z=q=0; 左重右轻时z=da,q=db,d=0; 右重左轻时z=db,q=da,d=0）。 4， 反复上述过程，直到(z=1,q=0,d=0)或(q=1,z=0,d=0)为止，即可确定次球号ci及其轻重 。 （详见下节）。 天平秤球最佳规格化的算法流程 （一）初始化： 输入球数m； 计算秤球次数n，n=N(m)=【log(2m+3)/log3】； 随机生成一个次球号ci及随机生成其对应值-1（轻)或1（重); 设置m维球状态向量X=(x1,…,xm)，（xi=0表示i号球为正品；而次球号ci的值1为重球，-1为轻球），如本例m=8 ci=4, ci的值为-1, X=(0,0,0,-1,0,0,0,0)。 设置m维球估计向量Y=(y1,…,ym)，（yi=0表示i号球判定为正品；1为疑重球；-1为疑轻球，2表示i号球待定不知其轻重），其初值全部为2，如本例Y=(2,2,2,2,2,22,2)。 初始疑重球数z=0、疑轻球数q=0、待定球数d=m，如本例d=m=8; 秤球序数的初值cn=1; （二）第一次秤球： 采用L三分法分解d, d=da+db+dc分别放置于天平的左盘、右盘及不放盘，如8=3+3