哪位大侠会算这个呀?? 指点一二，谢谢

Kitaisky

题目如下，重点是程序在3分钟之内计算完毕。

求详细代码和解释，请指点一下

一个N位的十进制正整数，如果它的每个位上的数字的N次方的和等于这个数本身，则称其为花朵数。
例如：
当N=3时，153就满足条件，因为 1^3 + 5^3 + 3^3 = 153，这样的数字也被称为水仙花数（其中，“^”表示乘方，5^3表示5的3次方，也就是立方）。
当N=4时，1634满足条件，因为 1^4 + 6^4 + 3^4 + 4^4 = 1634。
当N=5时，92727满足条件。
实际上，对N的每个取值，可能有多个数字满足条件。

程序的任务是：求N=21时，所有满足条件的花朵数。注意：这个整数有21位，它的各个位数字的21次方之和正好等于这个数本身。
如果满足条件的数字不只有一个，请从小到大输出所有符合条件的数字，每个数字占一行。因为这个数字很大，请注意解法时间上的可行性。要求程序在3分钟内运行完毕

zhlong8

这完全是个数学问题好吧

yinyuemi

回复 1# Kitaisky


    C/C++版上有人问过

Kitaisky

回复 2# zhlong8


回版主的话，是数学问题，Perl的算法我也有一个，但算到第8位就两分半了；求优化:

#!/usr/bin/perl

use strict;
use warnings;

my $now = time;

my $n = shift;
my $start = 10**($n-1)+1;
my $end = 10**$n-1;

my @result;

my %cache;
$cache{$_} = ($_)**$n for (0..9);

for my $number ($start..$end) {

    my $sum = 0;

    for (split '', $number) {
        $sum += $cache{$_};
    }   

    push @result, $number if ($sum == $number);
}

print "For n = $n, the numbers are: ", "@result\n";
print "Time: ", time()-$now, "secs\n";

zhlong8

回复  zhlong8


回版主的话，是数学问题，Perl的算法我也有一个，但算到第8位就两分半了；求优化:

...
Kitaisky 发表于 2011-08-18 14:36

    用循环到宇宙毁灭也算不完，所以说是数学问题

yinyuemi

回复 4# Kitaisky


   for 循环的开始可以优化下：
  当$_**$n*$n>10**$n时，从$start=$(n-1) x n开始循环判断;

    可以加个判断结束循环
  last if ($sum >10**$n-1)

cecial

参考一下
http://www.iteye.com/topic/763764

Kitaisky

回复 7# cecial


    谢谢 yinyuemi 和 cecial...

    居然发现:

    for (small_number..large_number) {

      # 有 implementation 的上限；
      # Back to being a C-programmer
}

hitsubunnu

换 天河一号  30秒就算出来了

ssfjhh

提供个python版的

1. def narcissistic(n):
   
2.     length = len(str(n))
   
3.     sumn = sum([int(i)**length for i in list(str(n))])
   
4.     if n == sumn:
   
5.         return True
   
6.     else:
   
7.         return False
   
8. 
   
9. if __name__ == '__main__':
   
10.     for i in range(10000):
    
11.         if narcissistic(i):
    
12.             print(i)
    

复制代码