例如:
當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次方之和正好等于這個數本身。
如果滿足條件的數字不只有一個,請從小到大輸出所有符合條件的數字,每個數字占一行。因為這個數字很大,請注意解法時間上的可行性。
優化:不管一個數是多少位,它都是由0~9這10個數字組成,所以我們可以先分別求出這10個數的各自的21次方,
? ? ? ? ? ?又因為10個9的21次方之和超過了21位數,所以從9個9,12個8開始一次往下運算即可
這個題意我是看懂了,上網百度了一下別人的代碼,但是有三處我看不懂。
#include <iostream>
#include <stdio.h>
#include<time.h>
#include <string.h>
using namespace std;
int A[10][21],count[100][21],cnt=0;
void Get()//將0-9這10個數的每個數的21次方存入A數組
{
? ? memset(A,0,sizeof(A));
? ? for(int k=0; k<10; k++) //一共有10個數
? ? {
? ? ? ? A[k][20]=1; ?//因為是乘法所以每個數的底數設為1
? ? ? ? int sum=0;
? ? ? ? for(int i=1; i<=21; i++) //21次方
? ? ? ? {
? ? ? ? ? ? for(int j=20; j>=0; j--)
? ? ? ? ? ? {
? ? ? ? ? ? ? ? sum=sum+A[k][j]*k;
? ? ? ? ? ? ? ? A[k][j]=sum%10;
? ? ? ? ? ? ? ? sum=sum/10;
? ? ? ? ? ? }
? ? ? ? }
? ? }
}
void Cal()
{
? ? int B[21]; //b[N]存從99……99到100……00之間的數 初始化為9個9后面全是8,因為10個9的21次方之和將超過21位數
? ? for(int i=0; i<9; i++) ?B[i]=9;
? ? for(int i=9; i<21; i++) B[i]=8;
? ? while(B[0]!=0)
? ? {
? ? ? ? int flag=1;
? ? ? ? int sum[21]= {0},sort[21]= {0},E[10]= {0};
? ? ? ? for(int i=0; i<21; i++) //求B中21個數的21次方之和
? ? ? ? {
? ? ? ? ? ? int temp=0;
? ? ? ? ? ? for(int j=20; j>=0; j--)
? ? ? ? ? ? {
? ? ? ? ? ? ? ? temp=temp+A[B[i]][j]+sum[j];
? ? ? ? ? ? ? ? sum[j]=temp%10;
? ? ? ? ? ? ? ? temp=temp/10;
? ? ? ? ? ? ? ? if(j==0&&temp>0)
? ? ? ? ? ? ? ? {
? ? ? ? ? ? ? ? ? ? flag=0;
? ? ? ? ? ? ? ? ? ? break;//位數超過21位舍棄
? ? ? ? ? ? ? ? }
? ? ? ? ? ? }
? ? ? ? ? ? if(flag==0)break;
? ? ? ? }
? ? ? ? if(sum[0]==0) flag=0;//位數不夠21位舍棄
? ? ? ? if(flag==1)//位數正好21位看看兩個數是不是相等
? ? ? ? {
? ? ? ? ? ? for(int i=0; i<21; i++) ? ? ? ??
? ? ? ? ? ? ? ? E[sum[i]]++; ? ? //統計sum中分別有幾個0到9 存到E[N]中
? ? ? ? ? ? for(int i=9,x=0; i>=0; i--)
? ? ? ? ? ? ? ? for(int j=0; j<E[i]; j++)
? ? ? ? ? ? ? ? ? ? sort[x++]=i;
? ? ? ? ? ? //比較二者是否相等
? ? ? ? ? ? for(int i=0; i<21; i++)
? ? ? ? ? ? {
? ? ? ? ? ? ? ? if(B[i]!=sort[i]) ?//難處1:這個B不應該跟sum比較嗎,關sort什么事
? ? ? ? ? ? ? ? {
? ? ? ? ? ? ? ? ? ? flag=0;
? ? ? ? ? ? ? ? ? ? break;
? ? ? ? ? ? ? ? }
? ? ? ? ? ? }
? ? ? ? }
? ? ? ? if(flag==1)//兩個數相等就存起來
? ? ? ? {
? ? ? ? ? ? for(int j=0; j<21; j++)
? ? ? ? ? ? ? ? count[cnt][j]=sum[j];//若這個地方換成了count[cnt][j]=B[j]那結果直接就不對,根據題意B應該跟sum是一樣的??????
? ? ? ? ? ? cnt++;
? ? ? ? }
? ? ? ? //找下一個B[N] ?(找下一個B的找法我是一點都不理解啊)
? ? ? ? int i;
? ? ? ? for(i=20; i>=0; i--)
? ? ? ? ? ? if(B[i]!=0)
? ? ? ? ? ? ? ? break;
? ? ? ? int temp=B[i];
? ? ? ? while(i<21)
? ? ? ? {
? ? ? ? ? ? B[i++]=temp-1; //保證去掉重復的數 ?并按從大到小順序 進行處理(尤其是這個地方直接不懂怎么去的重復的數)
? ? ? ? }
? ? }
? ? for(int i=cnt-1; i>=0; i--)
? ? {
? ? ? ? for(int j=0; j<21; j++)
? ? ? ? ? ? cout<<count[i][j];
? ? ? ? cout<<endl;
? ? }
}
int main()
{
? ? Get();
? ? Cal();
? ? printf("\n程序運行了%.2lf秒\n",(double)clock()/CLOCKS_PER_SEC);//只是為了測試一下程序運行的時間
? ? return 0;
}