求21位數的水仙花數

/*
 * 求21位數的水仙花數
 */
import java.math.BigInteger;

public class Demo01_BigInteger {
    // 求 每個 i 的 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){
        if(cur==9){
            nn[9] = 21 - use;
            ji_suan(pw,nn);
            return;
        }
        // 對當前位置所有可能進行枚舉
        for(int i=0;i<21-use;i++){
            nn[cur] = i;
            f(pw,nn,cur+1,use+i);
        }
    }
    public static void main(String[] args){
        long startTime = System.currentTimeMillis();    // 程序開始時間
        BigInteger pw[] = new BigInteger[10];
        for(int i=0;i<pw.length;i++){
            pw[i] = p(i);
        }
        int nn[] = new int[10];
        f(pw,nn,0,0);
        System.out.println("OK");
        long endTime = System.currentTimeMillis();  // 程序結束時間
        System.out.println((endTime-startTime)/1000f+"秒");  // 運行總時間
    }
}

運行結果:

128468643043731391252
449177399146038697307