Talk:Hofstadter Figure-Figure sequences: Difference between revisions

→‎timings for the REXX solutions: added/changed whitespace and comments, change a variable's name.
m (→‎timings for the REXX solutions: edited the source.)
(→‎timings for the REXX solutions: added/changed whitespace and comments, change a variable's name.)
Line 32:
 
==timings for the REXX solutions==
 
I normally don't including timings for the REXX solutions that I post, but when I saw the 2<sup>nd</sup> REXX example's timings, <br>
I decided to go back and include the timings here as the REXX 2<sup>nd</sup> example's timings seemed a bit high.
 
<br>I didn't expect a difference of three orders of magnitude.
<lang rexx>/*REXX pgmprogram calculates &and verifies the Hofstadter Figure-FigureFigure─Figure sequences.*/
call time 'Reset████████████████████████████████████████████████████████████████████████████████'
call time 'Reset'
parse arg x hightop bot . /*obtain anyoptional C.L.arguments specifications.from the CL*/
if x=='' then x=10;| if highx==',' then high x=1000 10 /*Not specified? Then use the defaults?default.*/
if bottop=='' | then bottop=40=',' then top=1000 /* " " " /* " " " */
lowif bot==1;'' | if x<0bot==',' then lowbot=abs(x) 40 /*only show" " " " " a single " │X│ value.*/
r.low=01; r.1=1; rr.=r.; if x<0 rr.1then low=1abs(x) /*only display a /*initialize thesingle R│X│ & RRvalue? arrays.*/
sr.=0; sr.1=2 1; rr.=r.; rr.1=1; s.=r.; s.1=2 /*initialize the R " " RR S arrayarrays. */
errs=0; $.=0
do i=low to abs(x) /*showdisplay the 1st first X values of R & S .*/
say right('R('i") =", 20) right(ffr(i), 7), /*show nice*/
right('S('i") =", 20) right(ffs(i), 7) /* R & S */
end /*i*/
if x<1 then exit /*if x ≤ 0, then we're all done.*/
do m=1 for bot; r=ffr(m); $.r=1
end /*m*/ /* [] verifycalculate the 1st 40 all 1≤#≤1kR presentvalues.*/
 
do mn=1 for top-bot; r=ffr(m); $.r=1; end /*calculate 1st 40 R values.*/ s=ffs(n)
if $.s then call ser 'duplicate number in R and S lists:' s; $.s=1
end /*n*/ /* [↑] calculate the 1st 960 S values.*/
 
do v=1 for top
@='duplicate number in R and S lists:' /* [↓] calc. 1st 960 S values.*/
do n=1 for high-bot; s=ffs(n); if \$.sv then call sErrser @ s; $.s=1; end'missing R │ S:' v
end /*v*/ /* [↑] are all 1≤ numbers ≤1k present?*/
 
/* [↓] verify all 1≤#≤1k present*/
do v=1 for high; if \$.v then call sErr 'missing R │ S:' v; end
say 'took' format(time('Elapsed'),,2) "seconds."
say
if errs==0 then say 'verification completed for all numbers from 1 ──►' hightop " [inclusive]."
else say 'verification failed with' errs "errors."
say
say 'took' format(time('Elapsed█████████████████████████████████████████████████████████████████'),,2) "seconds."
exit /*stick a fork in it, we're done.*/
/*────────────────────────────────────────────────────────────────────────────*/
/*──────────────────────────────────FFR subroutine──────────────────────*/
ffr: procedure expose r. s. rr.; parse arg n /*obtain the number parsefrom argthe narg.*/
if r.n\==0 then return r.n /*Defined? Then return the value.*/
_=ffr(n-1) + ffs(n-1) /*calculate the FFR & FFS values.*/
r.n=_; rr._=1; return _ /*assign value to R & RR; return.*/
/*────────────────────────────────────────────────────────────────────────────*/
/*──────────────────────────────────FFS subroutine──────────────────────*/
ffs: procedure expose r. s. rr.; parse arg n /*search for ¬null R│SR or S #number.*/
if s.n==0 then do k=1 for n /* [↓] 1st IF is a shortSHORT circuitCIRCUIT*/
if s.k\==0 then if r.k\==0 then iterate
call ffr k /*define R.k via FFR subroutine*/
km=k-1; _=s.km+1 /*the next S number, possibly. */
_=_+rr._; s.k=_ /*define an element of S array.*/
end end /*k*/
return s.n /*return theS.n value to the invoker*/
/*────────────────────────────────────────────────────────────────────────────*/
/*──────────────────────────────────SERR subroutine─────────────────────*/
sErrser: errs=errs+1; say '***error***!' arg(1); return</lang>
'''output''' &nbsp; when using the defaultsdefault inputs:
<pre>
R(1) = 1 S(1) = 2
Line 92 ⟶ 94:
R(9) = 56 S(9) = 13
R(10) = 69 S(10) = 14
took 0.22 seconds.
 
verification completed for all numbers from 1 ──► 1000 [inclusive].
 
took 0.22 seconds.
</pre>
The (above) example was run under Windows 7 on aan air-gap HP box (3.2GHz) using Regina REXX version 3.9.1.
<br><br>