|
|
A195333
|
|
Numbers n that can be expressed as the sum of the arithmetic derivatives of k consecutive numbers starting from n for some k.
|
|
2
|
|
|
1, 2, 4, 6, 25, 27, 33, 42, 221, 274, 581, 1957, 3125, 11406, 47058, 823543, 1535573, 5056941, 19246541, 19571621, 36861842, 50330577, 2590282517, 45546909393
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
A051674 is a subsequence of this sequence.
|
|
LINKS
|
|
|
FORMULA
|
n = Sum_{j=1..k} (n+j-1)', for some k >= 1.
|
|
EXAMPLE
|
k=1: n=27 -> 27 = 27'.
k=2: n=33 -> 33 = 33' + 34' = 14 + 19.
k=3: n=1957 -> 1957 = 1957' + 1958' + 1959' = 122 + 1179 + 656.
|
|
MAPLE
|
with(numtheory);
local b, c, n, p;
for n from 1 to i do c:=0; b:=-1;
while c<n do b:=b+1; c:=c+(n+b)*add(op(2, p)/op(1, p), p=ifactors(n+b)[2]); od;
if n=c then print(n); fi; od; end:
|
|
MATHEMATICA
|
dn[0] = 0; dn[1] = 0; dn[n_?Negative] := -dn[-n]; dn[n_] := Module[{f = Transpose[FactorInteger[n]]}, If[PrimeQ[n], 1, Plus @@ (n*f[[2]]/f[[1]])]]; t = {}; Do[k = n; sm = dn[n]; While[sm < n, k++; sm = sm + dn[k]]; If[sm == n, AppendTo[t, n]], {n, 100000}]; t (* T. D. Noe, Jan 04 2013 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|