|
|
A008057
|
|
Smallest sum of an addition chain for 2n+1.
|
|
1
|
|
|
0, 5, 10, 16, 20, 27, 31, 35, 40, 47, 51, 56, 60, 65, 74, 78, 80, 86, 92, 96, 102, 106, 110, 120, 121, 125, 134, 137, 142, 148, 153, 156, 160, 167, 171, 182, 184, 185, 192, 201, 200, 206, 210, 219, 227, 231, 233
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
EXAMPLE
|
The smallest chain for 5 is 2, 3, 5 with sum a(2) = 2+3+5 = 10.
the smallest chain for 7 is 2, 3, 4, 7 with sum a(3) = 2+3+4+7 = 16.
|
|
PROG
|
(PARI) step(V)=my(U=List(), v); for(i=1, #V, v=V[i]; for(i=1, #v, for(j=i, #v, if(v[i]+v[j]>v[#v], listput(U, concat(v, v[i]+v[j])))))); vecsort(Vec(U), , 8)
sm(v)=sum(i=2, #v, v[i])
a(n)=if(n<2, return(5*n)); n=2*n+1; my(V=[[1, 2]], U, t); while(#(U=select(v->v[#v]==n, V))==0, V=select(v->v[#v]<=n, step(V))); t=vecmin(apply(sm, U)); while(#V, V=step(select(v->sm(v)<t-n, V)); for(i=1, #V, my(v=V[i]); if(v[#v]==n, t=min(sm(v), t)))); t \\ Charles R Greathouse IV, Jul 17 2013
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,nice,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|