21位的花朵数

import java.math.BigInteger;
/*
 * 一个N位的十进制正整数，如果它的每个位上的数字的N次方的和等于这个数本身，则称其为花朵数。
例如：
当N=21时，1521就满足条件，因为 1^21 + 5^21 + 21^21 = 1521，这样的数字也被称为水仙花数（其中，“^”表示乘方，5^21表示5的21次方，也就是立方）。
当N=4时，16214满足条件，因为 1^4 + 6^4 + 21^4 + 4^4 = 16214。
当N=5时，92727满足条件。
实际上，对N的每个取值，可能有多个数字满足条件。
程序的任务是：求N=21时，所有满足条件的花朵数。注意：这个整数有21位，它的各个位数字的21次方之和正好等于这个数本身。
如果满足条件的数字不只有一个，请从小到大输出所有符合条件的数字，每个数字占一行。因为这个数字很大，请注意解法时间上的可行性。
*/
public class a_21位的花朵数 {
	
	      public static BigInteger p(int i){  
	            BigInteger base = BigInteger.valueOf(i);  
	            return base.pow(21);  
	        }  
	        public static void ji_suan(BigInteger[] pw,int[] nn){  
	            BigInteger sum = BigInteger.ZERO;  
	            for(int i=0;i<10;i++){  
	                sum = sum.add(  
	                        pw[i].multiply(BigInteger.valueOf(nn[i])));  
	            }  
	            String s = "" + sum;  
	            if(s.length()!=21) return;  
	            // 确定各数字出现的多少次  
	            int[] nn2 = new int[10];  
	            for(int i=0;i<21;i++){  
	            
	                nn2[s.charAt(i)-'0']++;  
	            }  
	            for(int i=0;i<10;i++){  
	                if(nn[i]!=nn2[i]) return;  
	            }  
	            System.out.println(s);  
	        }  
	        public static void f(BigInteger[] pw, int[] nn, int cur, int use){    //cur 当前处理的是nn[]的第几位
	        	                                                                //use 21位中的位数已经用掉了多少
 	            if(cur==9){  
	                nn[9] = 21 - use;  
	                ji_suan(pw,nn);  
	                return;  
	            }  
	            // 对当前位置所有可能进行枚举 //nn[0]  nn[1]  nn[2]  nn[3]  nn[4]  nn[5]  nn[6]  nn[7]  nn[8]  nn[9] 
	            for(int i=0;i<21-use;i++){  // 0        0     0      0      0      0      0      0      0     21
	                nn[cur] = i;            // 0        0     0      0      0      0      0      0      1     20
	                                        // 0        0     0      0      0      0      0      0      2     19 
	                f(pw,nn,cur+1,use+i);    //nn[]所有元素加起来=21  第一行代表21位数都是9
	            }  
	        }  
	        public static void main(String[] args){  
	            long startTime = System.currentTimeMillis();    // 程序开始时间  
	            BigInteger pw[] = new BigInteger[10];  
	            for(int i=0;i<pw.length;i++){        //将0到9的21次方按顺序保存在PW[]数组里
	                pw[i] = p(i);  
	            }  
	            int nn[] = new int[10];       //记录这个21位数中各个数字出现了几次
	            f(pw,nn,0,0);  
	            System.out.println("OK");  
	            long endTime = System.currentTimeMillis();  // 程序结束时间  
	            System.out.println((endTime-startTime)/1000f+"秒");  // 运行总时间  
	        }  
	    }