|
|
A125180
|
|
a(n) = 2*a(n-1) + prime(n) - prime(n-1), a(1)=2, where prime(n) denotes the n-th prime.
|
|
3
|
|
|
2, 5, 12, 26, 56, 114, 232, 466, 936, 1878, 3758, 7522, 15048, 30098, 60200, 120406, 240818, 481638, 963282, 1926568, 3853138, 7706282, 15412568, 30825142, 61650292, 123300588, 246601178, 493202360, 986404722, 1972809448, 3945618910
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Limit_{n->oo} a(n)/a(n-1) = 2.
|
|
LINKS
|
|
|
FORMULA
|
a(n) = prime(n) + Sum_{j=1..n-1} 2^(n-j-1)*prime(j), where prime(k) denotes the k-th prime.
a(1) = 2 and a(n) = prime(n) + Sum_{i=1..n-1} a(i) for n > 1. - Alexandre Herrera, Dec 10 2023
|
|
EXAMPLE
|
a(4)=26 because 4*prime(1)+2*prime(2)+prime(3)+prime(4) = 8+6+5+7 = 26.
|
|
MAPLE
|
a[1]:=2: for n from 2 to 35 do a[n]:=2*a[n-1]+ithprime(n)-ithprime(n-1) od: seq(a[n], n=1..35);
|
|
MATHEMATICA
|
a[1] = 2; a[n_] := 2*a[n - 1] + Prime[n] - Prime[n - 1]; Table[a[n], {n, 1, 31}] (* James C. McMahon, Dec 10 2023 *)
|
|
PROG
|
(Magma) [n eq 1 select 2 else 2*Self(n-1)+NthPrime(n)-NthPrime(n-1):n in [1..31]]; // Marius A. Burtea, Oct 17 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|