问题提出:
设x,y为非负整数,试计算集合M={2^x*3^y|x>=0,y>=0}的元素小于指定整数n的个数,并求这些元素从小到大排序的第m项。
以下给出案例求解的三种设计:
1.穷举求解
设计要点:
集合元素由2的幂与3的幂及其乘积组成,设元素从小到大排序的第k项为f(k)。显然,f(1)=1,f(2)=2。
设置a从3开始递增至n循环,对每一个a(赋值给j),逐次试用2试商,然后逐次试用3试商。
试商后若j>1,说明原a有2,3以外的因数,不属于该序列。
若j=1,说明原a只有2,3的因数,为序列第k项赋值。
由于实施从小到大穷举赋值,所得项无疑是从小到大的序列。
当达到指定的n,退出循环,输出指定项f(m)。
代码如下:
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![ExpandedBlockStart.gif](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
2 void main()
3 { int k,m;
4 long a,j,n,f[ 1000 ];
5 printf( " 计算小于n的项数,请指定n: " );
6 scanf( " %ld " , & n);
7 printf( " 输出序列的第m项,请指定m: " );
8 scanf( " %d " , & m);
9 f[ 1 ] = 1 ;f[ 2 ] = 2 ;k = 2 ;
10 for (a = 3 ;a <= n;a ++ )
11 { j = a;
12 while (j % 2 == 0 ) j = j / 2 ; // 反复用2试商
13 while (j % 3 == 0 ) j = j / 3 ; // 反复用3试商
14 if (j == 1 )
15 { k ++ ;f[k] = a;}
16 }
17 printf( " 幂序列中小于%ld的项数为:%d\n " ,n,k);
18 if (m <= k)
19 printf( " 从小到大排序的第%d项为:%ld\n " ,m,f[m]);
20 else
21 printf( " 所输入的m大于序列的项数!\n " );
22 }
23
2. 递推排序求解
设计要点:
为探索x+y=i时各项与x+y=i−1时各项之间的递推规律,剖析x+y的前若干项情形:
x+y=0时,元素为1(初始条件);
x+y=1时,元素为2*1=2,3*1=3,共2项;
x+y=2时,序列有2*2=4,2*3=6,3*3=9,共3项;
x+y=3时,序列有2*2*2=8,2*2*3=12,2*3*3=18,3*3*3=27,共4项;
……
可归纳出以下递推关系:
x+y=i时,序列共i+1项,其中前i项是x+y=i−1时的所有i项分别乘2所得;最后一项为x+y=i−1时的最后一项乘3所得(即t=3^i)。
注意,对x+y=i−1的所有i项分别乘2,设为f[h]*2,必须检测是否小于n而大于0。同样,对t也必须检测是否小于n而大于0。只有小于n且大于0时才能赋值。
这里要指出,最后若干行可能不是完整的,即可能只有前若干项能递推出新项。为此设置变量u: 当一行有递推项时u=1;否则u=0。只有当u=0时停止,否则会影响序列的项数。
代码如下:
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![ExpandedBlockStart.gif](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
2 void main()
3 { int i,j,h,k,m,u,c[ 100 ];
4 double d,n,t,f[ 1000 ];
5 printf( " 计算小于n的项数,请指定n: " );
6 scanf( " %lf " , & n);
7 printf( " 输出序列的第m项,请指定m: " );
8 scanf( " %d " , & m);
9 k = 1 ;t = 1.0 ; i = 1 ;
10 c[ 0 ] = 1 ; f[ 1 ] = 1.0 ;
11 while ( 1 )
12 { u = 0 ;
13 for (j = 0 ;j <= i - 1 ;j ++ )
14 { h = c[i - 1 ] + j;
15 if (f[h] * 2 < n && f[h] > 0 ) // 第i行各项为前一行各项乘2
16 { k ++ ;f[k] = f[h] * 2 ;u = 1 ;
17 if (j == 0 ) c[i] = k; // 该行的第1项的项数值赋给c(i)
18 }
19 else break ;
20 }
21 t = t * 3 ; // 最后一项为3的幂
22 if (t < n && t > 0 )
23 { k ++ ;f[k] = t; } // 用t给f[k]赋值
24 if (u == 0 ) break ;
25 i ++ ;
26 }
27 for (i = 1 ;i < k;i ++ ) // 逐项比较排序
28 for (j = i + 1 ;j <= k;j ++ )
29 if (f[i] > f[j])
30 { d = f[i];f[i] = f[j];f[j] = d;}
31 printf( " 幂序列中小于%f的项数为:%d\n " ,n,k);
32 if (m <= k)
33 printf( " 从小到大排序的第%d项为:%.0f\n " ,m,f[m]);
34 else
35 printf( " 所输入的m大于序列的项数!\n " );
36 }
37
3.递推结合比较赋值求解
设计要点:
从u=1, f(u)=1开始, 在已求得f(u)的基础上,可递推地求出f(u+1):求出各大于f(u)的最小数, 取其中最小者即为f(u+1)。
递推结合比较赋值设置永真外循环,实施乘2的内循环。
首先,从p=0开始,若q[p]<=f[u],则递推得一个3的幂即q[p]=3^p,并赋给最小值标志量h。
然后转入内循环i(0——p-1)中,若q[i]<=f[u],则q[i]乘2,即q[i]=2*q[i]。
然后q[i]与h比较,即2^j*3^i(i<p)与3^p比较,取较小者为h。
若h≤n,则h赋值给序列新的项,用u标记的项数。
若h>n,表明递推结合比较赋值完成,退出外循环,输出序列的项数u与序列中指定的项f[m]后结束。
代码如下:
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![ExpandedBlockStart.gif](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
2 void main()
3 { int i,m,u,p;
4 double n,h,f[ 1000 ],q[ 100 ];
5 printf( " 计算小于n的项数,请指定n: " );
6 scanf( " %lf " , & n);
7 printf( " 输出序列的第m项,请指定m: " );
8 scanf( " %d " , & m);
9 u = 1 ; f[u] = 1.0 ;
10 p = 0 ; q[p] = 1.0 ;
11 while ( 1 )
12 { if (q[p] <= f[u])
13 { p ++ ; q[p] = 3 * q[p - 1 ]; }
14 h = q[p]; // 递推3的幂,q[p]=3^p
15 for (i = 0 ;i < p;i ++ )
16 { if (q[i] <= f[u]) q[i] *= 2 ; // 幂积q[i]=2^j*3^i,j=1,2,3,…
17 if (q[i] < h) h = q[i];
18 }
19 if (h > n) break ;
20 u ++ ;f[u] = h;
21 }
22 printf( " 幂序列中小于%.0f的项数为:%d\n " ,n,u);
23 if (m <= u)
24 printf( " 从小到大排序的第%d项为:%.0f\n " ,m,f[m]);
25 else
26 printf( " 所输入的m大于序列的项数!\n " );
27 }
28