|
|
A335767
|
|
Numbers m such that m = S_af(m) = af(d_1)+af(d_2)+...+af(d_k) where d_1 d_2 ... d_n is the decimal expansion of m and af(m) = m!-(m-1)!+(m-2)!+...1! (alternating factorial) with S_af(0) = 1.
|
|
0
|
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Largest k such that S_af(k) > k is 1599999. That's why there are only four numbers such that S_af(m) = m. Proved by computer calculations.
If m has eight or more digits then S_af(m) < m. Proved directly.
|
|
LINKS
|
|
|
EXAMPLE
|
For m = 643, S_af(643) = af(6)+af(4)+af(3) = 619+19+5 = 643.
|
|
MATHEMATICA
|
af[0] = af[1] = 1; af[n_] := af[n] = n! - af[n - 1]; Select[Range[400000], Total[af /@ IntegerDigits[#]] == # &] (* Amiram Eldar, Jun 24 2020 *)
|
|
PROG
|
(C++)
#include <iostream>
using namespace std;
int S_af(int n) { const int af[]={1, 1, 1, 5, 19, 101, 619, 4421, 35899, 326981}; int s=0; while (n) {s+=af[n%10]; n/=10; } return s; }
int main() {int n=1; while(n<=1599999){if(n==S_af(n)){cout<<n<<endl; }n++; }}
|
|
CROSSREFS
|
Cf. A005165 (alternating factorial).
|
|
KEYWORD
|
nonn,base,fini,full
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|