|
Search: id:A137395
|
|
|
| A137395 |
|
a(1)=1. a(n) = a(n-1) + (sum of the distinct primes that are <= n and don't divide a(n-1)). |
|
+0
1
|
|
| 1, 3, 5, 10, 13, 23, 40, 50, 60, 67, 95, 118, 157, 198, 223, 264, 306, 342, 395, 467, 544, 602, 693, 772, 870, 960, 1050, 1133, 1251, 1377, 1517, 1677, 1821, 1978, 2113, 2273, 2470, 2628, 2820, 3007, 3214, 3450, 3698, 3934, 4206, 4482, 4805, 5097, 5422, 5748
(list; graph; listen)
|
|
|
OFFSET
|
1,2
|
|
|
LINKS
|
Leroy Quet, Home Page (listed in lieu of email address)
|
|
EXAMPLE
|
The primes <= 8 are 2,3,5,7. Of these, only 3 and 7 don't divide a(7)=40. So a(8) = a(7) + 3 + 7 = 50.
|
|
MAPLE
|
A137395 := proc(n) option remember ; local aprev, i, a, p; if n = 1 then RETURN(1) ; fi; aprev := A137395(n-1) ; a := aprev ; for i from 1 do p := ithprime(i) ; if p > n then break; fi ; if aprev mod p <> 0 then a := a+p ; fi ; od: a ; end: seq(A137395(n), n=1..100) ; - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), May 23 2008
|
|
CROSSREFS
|
Adjacent sequences: A137392 A137393 A137394 this_sequence A137396 A137397 A137398
Sequence in context: A034746 A080931 A031878 this_sequence A001767 A048214 A001841
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Leroy Quet Apr 10 2008
|
|
EXTENSIONS
|
More terms from R. J. Mathar (mathar(AT)strw.leidenuniv.nl), May 23 2008
|
|
|
Search completed in 0.002 seconds
|
