?

題目:

找到滿足條件的數(shù)組
給定函數(shù)d(n)=n+n的各位之和，n為正整數(shù)，如d(78)=78+7+8=93。這樣這個函數(shù)可以看成一個生成器，如93可以看成由78生成。
定義數(shù)A：數(shù)A找不到一個數(shù)B可以由d(B)=A，即A不能由其他數(shù)生成。現(xiàn)在要寫程序，找出1至10000里的所有符合數(shù)A定義的數(shù)。
回答：
申請一個長度為10000的bool數(shù)組，每個元素代表對應(yīng)的值是否可以有其它數(shù)生成。開始時將數(shù)組中的值都初始化為false。
由于大于10000的數(shù)的生成數(shù)必定大于10000，所以我們只需遍歷1到10000中的數(shù)，計算生成數(shù)，并將bool數(shù)組中對應(yīng)的值設(shè)置為true，表示這個數(shù)可以有其它數(shù)生成。
最后bool數(shù)組中值為false的位置對應(yīng)的整數(shù)就是不能由其它數(shù)生成的。

?

---------------------------------------------------------

我嘗試用解方程的思路來看一下:

先把問題的范圍縮小到: 找出1至1000里的所有符合數(shù)A定義的數(shù), (為了便于表達(dá), 定義數(shù)A構(gòu)成一個集合SA, SA的補(bǔ)集為RA)

?

對于小于1000的正整數(shù)n(n的各位用ni表示), 有: n =? n0 + n1 * 10 + n2 * 10^2 =?∑ni * 10^i ;? i∈{0~2}, ni∈{0~9}

?

考察函數(shù) y?≡ f(n, x) = n2 * x^2 + n1 * x + n0; (x = 10, 因?yàn)槭鞘M(jìn)制數(shù))

其中, 向量n=(n2, n1, n0), ni∈{0~9}, i∈{0~2};

上述問題里的d(n) = n + ∑ni = f(n, x)+ ∑ni; (x = 10)

對于1至1000里的正整數(shù)N, 如果存在另一個正整數(shù)n, 使得d(n) = N, 那么N∈RA; 反之, 如果找不到另一個正整數(shù)n, 使得d(n) = N, 那么N?SA. .

考察方程(1): d(n) - N = 0; (即: n2 * x^2 + n1 * x + n0 + (n2 + n1 +n0) - N = 0) ..............(1)

如果把方程(1)看做是關(guān)于x的方程, 那么問題變?yōu)? 當(dāng)ni和N取哪些正整數(shù)時, 方程(1)至少有一個值為10的根.

?

?根據(jù)二次方程(A*x^2 + B*x +C = 0)求根公式:

x0 = (-b + sqrt(Δ)/(2a),? x1 = (-b - sqrt(Δ)/(2a);

Δ = B^2 - 4AC;

A = n2;

B = n1;

C = n0 + (n2 + n1 +n0) - N;

當(dāng)x = 10時, 得到N, ni的關(guān)系式:

N = 101*n2 + 11*n1 + 2*n0; ...........................(2)

這與1~1000時找N的方法是一致的.

?

?

現(xiàn)在考察問題范圍擴(kuò)大至1~10000, 這時就要涉及三次方程求根的方法了. 看看能得到什么結(jié)果

參考盛金公式:

一元三次方程aX³+bX²+cX+d=0，（a,b,c,d∈R，且a≠0） 重根判別式

總判別式Δ=B²－4AC

盛金判別法:

1.當(dāng)A=B=0時，方程有一個三重實(shí)根。 2,當(dāng)Δ=B²－4AC>0時，方程有一個實(shí)根和一對共軛虛根。 3,當(dāng)Δ=B²－4AC=0時，方程有三個實(shí)根，其中有一個二重根。 4,當(dāng)Δ=B²－4AC<0時，方程有三個不相等的實(shí)根。

?

考察函數(shù) y?≡ f(n, x) = n3 * x^3 + n2 * x^2 + n1 * x + n0; (x = 10, 因?yàn)槭鞘M(jìn)制數(shù))

考察方程(3): d(n) - N = 0; (即: n3 * x^3 + n2 * x^2 + n1 * x + n0 + (n3 + n2 + n1 +n0) - N = 0) ..............(3)

由盛金公式得:

a = n3;?????? b = n2; ? ? ?? c = n1;??????? d = n0 + (n3 + n2 + n1 +n0) - N = (n3 +n2 +n1 +2*n0 -N);

A = b^2 - 3ac = n2^2 - 3*n3*n1;

B = bc - 9ad = n2*n1 - 9* n3*(n3 + n2 + n1 +2*n0 - N);

C = c^2 - 3bd = n1^2 - 3* n2* (n3 + n2 + n1 +2*n0 - N);

Δ = B^2 - 4AC;

盛金判別法:

1.當(dāng)A=B=0時，方程有一個三重實(shí)根。

?

令x1 = x2 = x3 = 10, 得:

?

?

30a +b = 0;......(1.1)
10b +c = 0;......(1.2)
10c +3d = 0;......(1.3)
同時, A=0, B=0;得:
b^2 - 3ac = 0;......(1.4)
bc - 9ad = 0;......(1.5)
連立方程組得:
N = 2*n0 + 1271*n3;
考慮到a, b,c都是正整數(shù), 實(shí)際上可從(1.1), (1.2), (1.3) 得出a=b=c=0, 這和a≠0的假設(shè)矛盾, 所以這種情況下沒有我們想要的數(shù).

3.當(dāng)Δ=B²－4AC=0時，方程有三個實(shí)根，其中有一個二重根。 ......(3.1) ......(3.2) 其中 。 由(3.1)得:N = n2^3/(9*n3^2) - 4/9 * n2*n1/n3 + n3 +n2 - 7/3*n1 +2*n0 +10/9* n2^2/n3; 由(3.2)得: N = -(n3 + n2 + n1 + 2n0) + 20 *n2^2/(9*n3) - 20/3*n1 + 1/9* n2*n1/n3; (推導(dǎo)到這里已經(jīng)變得很惡心了啊.) 注意到, 如果令每一個加和項(xiàng)都為整數(shù)的話, 須滿足: n1 = 3 * i; n2 = 3* n3 *j; (i, j為正整數(shù)), 同時還要滿足N>0, 即使這樣, 得到的也只是答案的一個子集而以. 2.當(dāng)Δ=B²－4AC>0時，方程有一個實(shí)根和一對共軛虛根。 4.當(dāng)Δ=B²－4AC<0時，方程有三個不相等的實(shí)根。 至于2和4的解里面牽扯到sqrt, arccos等函數(shù), 在這些情況下要考慮整數(shù)解的話.........我還是決定放棄了:)

?哎, 結(jié)果不就是個 N = 1001*n3 + 101*n2 + 11*n1 + 2*n0嘛; 看來用解方程的方法還是繞遠(yuǎn)了啊,

?

轉(zhuǎn)載于:https://www.cnblogs.com/yaoyansi/p/4113002.html

總結(jié)

以上是生活随笔為你收集整理的一道算法题的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。