|
%I #16 Jan 22 2021 14:22:49
%S 10001,10010,10100,11000,20002,20020,20200,22000,30003,30030,30300,
%T 33000,40004,40040,40400,44000,50005,50050,50500,55000,60006,60060,
%U 60600,66000,70007,70070,70700,77000,80008,80080,80800,88000,90009,90090,90900
%N Numbers such that each digit is the sum of two or more other digits.
%C E. Angelini calls these "W numbers", see the link for more information.
%H E. Angelini, <a href="http://www.cetteadressecomportecinquantesignes.com/Wnumbers.htm">W and M numbers</a>, January 2012.
%H E. Angelini, <a href="/A203591/a203591.pdf">W and M numbers</a> [Cached copy, with permission]
%o (PARI) is_W(n)={ my(c=vector(10), OK(n,c) = c[2] >= n & return(1);
%o forstep( d=n-1,2,-1, forstep( m=min(c[d+1], n\d),1,-1, OK(n-m*d,c) & return(1))));
%o for( i=1,#n=Vecsmall( Str( n )), c[n[i]-47]++ ); c[1]<3 & return; c[2]==1 & return;
%o for( j=3,10, c[j] != 1 | OK( j-1, c ) | return); 1} \\ _M. F. Hasler_, Jan 03 2012
%Y Cf. A203592 for the primes in this sequence.
%K nonn,base
%O 1,1
%A _M. F. Hasler_, Jan 03 2012
|
