自方幂数问题.docVIP

自方幂数问题 江跃勇 200421001010 问题的提出： 如果一个n位正整数恰好等于它的n个位上数字的n次方和，则称该数为n位自方幂数。三位自方幂数又称水仙花数；四位自方幂数又称玫瑰花数；五位自方幂数又称五角星数；六位自方幂数又称六合数。 （1）设计算法求所有水仙花数； （2）设计算法求所有玫瑰花数； （3）设计算法综合求3~6位自方幂数。 解答： 求水仙花数的算法实现: o=1;p=1;q=1; for n=100:999 o=floor(n/100);p=floor((n-100.*o)/10);q=n-100.*o-10.*p; if n==o^3+p^3+q^3 n end end n = 153 n = 370 n = 371 n = 407 （2）求玫瑰花数算法实现： r=1; o=1;p=1;q=1; for n=1000:9999 o=floor(n/1000);p=floor((n-1000.*o)/100);q=floor((n-1000.*o-100.*p)/10); r=n-1000.*o-100.*p-10.*q; if n==o^4+p^4+q^4+r^4 n end end n = 1634 n = 8208 n = 9474 （3）求3~6位自方幂数的算法实现： o=1;p=1;q=1;r=1;s=1;t=1; for n=100:999999 if n=999 o=floor(n/100);p=floor((n-100.*o)/10);q=n-100.*o-10.*p; if n==o^3+p^3+q^3 n end elseif n=9999 o=floor(n/1000);p=floor((n-1000.*o)/100);q=floor((n-1000.*o-100.*p)/10); r=n-1000.*o-100.*p-10.*q; if n==o^4+p^4+q^4+r^4 n end elseif n=99999 o=floor(n/10000);p=floor((n-10000.*o)/1000); q=floor((n-10000.*o-1000.*p)/100);r=floor((n-10000.*o-1000.*p-100.*q)/10); s=n-10000.*o-1000.*p-100.*q-10.*r; if n==o^5+p^5+q^5+r^5+s^5 n end else o=floor(n/100000);p=floor((n-100000.*o)/10000); q=floor((n-100000.*o-10000.*p)/1000); r=floor((n-100000.*o-10000.*p-1000.*q)/100); s=((n-100000.*o-10000.*p-1000.*q-100.*r)/10); t=n-100000.*o-10000.*p-1000.*q-100.*r-10.*s; if n==o^6+p^6+q^6+r^6+s^6+t^6 n end end end n = 153 n = 370 n = 371 n = 407 n = 1634 n = 8208 n = 9474 n = 54748 n = 92727 n = 93084