%I #6 Aug 07 2015 02:48:10
%S 4,12,13,14,17,18,19,20,21,24,25,26,27,28,29,32,33,34,35,36,37,38,39,
%T 42,43,44,45,46,47,48,49,50,58,59,60,61,62,63,66,67,68,69,70,71,72,73,
%U 81,82,83,84,85,86,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104
%N Numbers which are not members of A074137.
%p A074137 := proc(nmax) local a,ipiv,apiv,incsum,nexta ; a := [1,2,3] ; ipiv := 2 ; while nops(a) < nmax do apiv := op(ipiv,a) ; while nops(a) < ipiv+apiv-1 do a := [op(a), op(-1,a)+1] ; od ; incsum := add(op(i,a),i=ipiv+1..nops(a)) ; nexta := apiv - ( incsum mod apiv ) ; while nexta <= op(-1,a) do nexta := nexta+apiv ; od ; a := [op(a), nexta] ; ipiv := ipiv+1 ; od ; a ; end ; a074137 := A074137(80) ; for n from 1 to op(-1,a074137) do if not n in a074137 then printf("%d,",n) ; fi ; od ; # _R. J. Mathar_, Mar 30 2007
%Y Cf. A074137.
%K nonn
%O 1,1
%A _Amarnath Murthy_, Aug 28 2002
%E More terms from _R. J. Mathar_, Mar 30 2007
|