C程序设计例解(03)

核心提示： 03. 找一个最小的自然数，使它等于不同的两组三个自然数的三次幂之和，C程序设计例解(03)，即找最小的x，使得： x=a*a*a+b*b*b+c*c*c+d*d*d+e*e*e+f*f*f 其中，则ih增1*/ { ih++; p[ih]=ih*ih*ih; /*为ih平面设定j的变化范

03. 找一个最小的自然数，使它等于不同的两组三个自然数的三次幂之和，即找最小的x，使得：
x=a*a*a+b*b*b+c*c*c+d*d*d+e*e*e+f*f*f
其中，a,b,c,d,e,f都是自然数,a<=b<=c<=d<=e<=f; [a,b,c]!=[d,e,f]
解：
利用上一问题的求解思想，上一问题在正方形平面下三角区内找解，本题在正立方体的下三角棱内找解。记i为三角棱体的平面，j为某平面的行，k为某行上的列。当前考察的下三角棱体的范围由最上平面至最下平面控制；对应每个平面的下三角区域，在每个下三角区域内当前待考查的行可由行的下界和上界控制，每个有效行上的候选列由其当前列来表示。因此有如下解法：
算法---找一个最小的自然数x，使它等于不同的两组三个自然数的三次幂之和
{
以三角棱体的顶点为最初候选者；
为最初寻找平面设定行的变化范围和列的初值；
do
{
保存上一个候选者；
if(当前候选者在最下平面)
{
寻找平面范围的最下平面向下进一层；
为新平面设定行的变化范围；
}
if(在最上平面最下角点找到候选者)
寻找平面范围的最上平面向下进一层；
else
{
if(在第一列找到候选者)
{
当前平面的行的变化上增1；
置当前平面的最高行的列为1；
}
if(在对角线上找到候选者)
当前平在的行的变化下界增1；
else
调整当前平面当前行的列号值；
}
在当前最上平面至当前最下平面范围内寻找最小值的候选者；
}while(两候选者对应的值不相等);
输出解；
}
因每个平面有行变化的下界和上界，程序分别用两个一维数组来存贮；每个平面的每行都有一个当前列，程序用一个二维数组来存贮；为避免反复计算一个整数的三次幂，另引入一个一维数组，对应第i下标位置存贮i*i*i。令当前找到的候选者为i1,j1,k表示在i1平面的第j1行的k1列找到的候选者。因候选者限制在三角棱内，i1,j1,k1满足条件：
i1>=j1>=k1
当候选者在最下平面时，则最下平面向下进一层，并为新平面设定行的变化范围和对应列号；当前最上平面的最下角点找到候选者时，最上平面向下进一层；当在第一列找到候选者时，当前平面的行的上界增，并为新的行设定初始列号；当在某行的对角线上找到候选者时，该行不应该再被考虑，当前平面的行的下界增1；其它情况，当前行的下一列将会被考虑，为该行调整当前列。在调整当前平面的行的下界和上界时，应不能超过当前平面号。为在三角棱体的当前有效平面内找最小值的候选者，先假定最上平面的最小行的当前列为下一个候选者，然后自最上平面至最下平面，每个平面自最小行至最大行，寻找最小值所在平面号、行号和列号。

程序代码如下：
#include<stdio.h>
#define N 100
void main()
{
int i,j,il,ih,i0,j0,k0,i1,j1,k1;
int jl[N],jh[N]; /*第i层平面的行的变化范围，自jl[i]至jh[i]*/
int k[N][N]; /*第i层平面中，对应行j，当前的列号值为k[i][j]*/
int p[N], min; /*p[i]=i*i*i*/
i1=1;j1=1;k1=1; /*首先只局限下三角棱体的顶点*/
il=1;ih=1; /*预置i的变化范围初值il<=i<=ih*/
jl[1]=1;jh[1]=1; /*对应i层平面的行的变化范围*/
k[il][jl[il>=1; /*第i层平面中，对应行的列的初值*/
p[1]=1;
do
{
min=p[i1]+p[j1]+p[k1];
i0=i1;j0=j1;k0=k1;
if(i1==ih) /*当前候选者在ih平面，则ih增1*/
{
ih++;
p[ih]=ih*ih*ih;
/*为ih平面设定j的变化范围和对应k值*/
jl[ih]=1;jh[ih]=1;k[ih][1]=1;
}
if(i1==il&&j1==il&&k1==il)
il++; /*在il平面最下角点找到候选者，il增1*/
else
{
if(k1==1&&jh[i1]<i1)
{ /*在第一列找到候选者，i1平面的行的上界增1*/
jh[i1]++;
k[i1][jh[i1>=1;
}
if(k1==j1&&jl[i1]<i1)
jl[i1]++; /*在对角线上找到候选者，il平面的行的下界增1*/
else
k[i1][j1]++; /*调整i1平面当前行的列号*/
}
i1=il; /*预定最上平面的最小行的当前列为下一个候选者*/
j1=jl[i1];
k1=k[i1][j1];
for(i=il;i<=ih;i++) /*寻找最小值所在平面号、行号和列号*/
{
for(j=jl[i];j<=jh[i];j++)
if(p[i]+p[j]+p[k[i][j><p[i1]+p[j1]+p[k1])

{
i1=i;j1=j;k1=k[i][j];
}
}
}while(p[i1]+p[j1]+p[k1]!=min&&ih!=N);
if(p[i1]+p[j1]+p[k1]==min)
PRintf("%4d=%2d^3+%d^3+%d^3=%2d^3+%d^3+%d^3\n",min,i0,j0,k0,i1,j1,k1);
else printf("The %d is too small.\n",N);
}

程序运行结果如下：
251 = 5^3 + 5^3 + 1^3 = 6^3 + 3^3 + 2^3