Cycle detection: Difference between revisions
m
→{{header|Wren}}: Changed to Wren S/H
Not a robot (talk | contribs) (Add BCPL) |
m (→{{header|Wren}}: Changed to Wren S/H) |
||
(14 intermediate revisions by 11 users not shown) | |||
Line 1:
{{
;Task:
Line 33:
=={{header|11l}}==
{{trans|D}}
<
print(‘Cycle length = ’len)
print(‘Start index = ’start)
V i = x0
L 1..start
i = f(i)
V cycle = [0] * len
L 0.<len
cycle[L.index] = i
i = f(i)
print(‘Cycle: ’, end' ‘’)
print(cycle)
F brent(f, x0)
Int cycle_length
V hare = x0
Line 59 ⟶ 72:
print_result(x0, f, cycle_length, cycle_start)
brent(i -> (i * i + 1) % 255, 3)</syntaxhighlight>
{{out}}
<pre>
Line 81:
=={{header|8086 Assembly}}==
<
org 100h
section .text
Line 187:
ix: db 'Index: $'
len: db 'Length: $'
nl: db 13,10,'$'</
{{out}}
<pre>3 10 101 2 5 26 167 95 101 2 5 26 167 95 101 2 5 26 167 95
Line 194:
=={{header|Ada}}==
This implementation is split across three files. The first is the specification of a generic package:
<
generic
type Element_Type is private;
Line 201:
procedure Brent(F : Brent_Function; X0 : Element_Type; Lambda : out Integer; Mu : out Integer);
end Brent;
</syntaxhighlight>
The second is the body of the generic package:
<
package body Brent is
procedure Brent (F : Brent_Function; X0 : Element_Type; Lambda : out Integer; Mu : out Integer) is
Line 233:
end Brent;
end Brent;
</syntaxhighlight>
By using a generic package, this implementation can be made to work for any unary function, not just integers. As a demonstration two instances of the test function are created and two instances of the generic package. These are produced inside the main procedure:
<
with Brent;
with Text_Io;
Line 277:
end loop;
end Main;
</syntaxhighlight>
{{out}}
<pre> 3, 10, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5
Line 288:
=={{header|APL}}==
{{works with|Dyalog APL}}
<
f←⍺⍺
lam←⊃{
Line 307:
}
(255 | 1 + ⊢×⊢) task 3</
{{out}}
<pre>3 10 101 2 5 26 167 95 101 2 5 26 167 95 101 2 5 26 167 95
Line 315:
=={{header|BCPL}}==
<
// Brent's algorithm. 'fn' is a function pointer,
Line 370:
// print the cycle
printRange(f, 3, mu, mu+lam)
$)</
{{out}}
<pre>3 10 101 2 5 26 167 95 101 2 5 26 167 95 101 2 5 26 167 95
Line 376:
Index: 2
101 2 5 26 167 95</pre>
=={{header|BQN}}==
<syntaxhighlight lang="bqn">_Brent← {F _𝕣 x0:
p←l←1
(I ← {p=l?
l↩1 ⋄ p×↩2 ⋄ 𝕩I𝕩 ;
l+↩1 ⋄ 𝕨I𝕩
}⍟≠⟜F) x0
m←0
{m+↩1 ⋄ 𝕨𝕊⍟≠○F𝕩}⟜(F⍟l) x0
l‿m‿(F⍟(m+↕l)x0)
}</syntaxhighlight>
{{out|Example use}}
<pre> (255|1+ט)_Brent 3
⟨ 6 2 ⟨ 101 2 5 26 167 95 ⟩ ⟩</pre>
=={{header|C}}==
{{trans|Modula-2}}
<
#include <stdlib.h>
Line 457 ⟶ 472:
return 0;
}</
{{out}}
<pre>[3, 10, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5]
Line 464 ⟶ 479:
Cycle = [101, 2, 5, 26, 167, 95]</pre>
=={{header|C_sharp|C#
This solution uses generics, so may find cycles of any type of data, not just integers.
<
// First file: Cycles.cs
Line 563 ⟶ 578:
}
</syntaxhighlight>
=={{header|C++}}==
<
int val;
ListNode *next;
Line 602 ⟶ 617:
}
}
}</
=={{header|CLU}}==
<syntaxhighlight lang="clu">% Find a cycle in f starting at x0 using Brent's algorithm
brent = proc [T: type] (f: proctype (T) returns (T), x0: T)
returns (int,int)
where T has equal: proctype (T,T) returns (bool)
pow: int := 1
lam: int := 1
tort: T := x0
hare: T := f(x0)
while tort ~= hare do
if pow = lam then
tort := hare
pow := pow * 2
lam := 0
end
hare := f(hare)
lam := lam + 1
end
tort := x0
hare := x0
for i: int in int$from_to(1,lam) do
hare := f(hare)
end
mu: int := 0
while tort ~= hare do
tort := f(tort)
hare := f(hare)
mu := mu + 1
end
return(lam, mu)
end brent
% Iterate over a function starting at x0 for N steps,
% starting at a given index
iterfunc = iter [T: type] (f: proctype (T) returns (T), x0: T, n, start: int)
yields (T)
for i: int in int$from_to(1, start) do
x0 := f(x0)
end
for i: int in int$from_to(1, n) do
yield(x0)
x0 := f(x0)
end
end iterfunc
% Iterated function
step_fn = proc (x: int) returns (int)
return((x*x + 1) // 255)
end step_fn
start_up = proc ()
po: stream := stream$primary_output()
% Print the first 20 items of the sequence
for i: int in iterfunc[int](step_fn, 3, 20, 0) do
stream$puts(po, int$unparse(i) || " ")
end
stream$putl(po, "")
% Find a cycle
lam, mu: int := brent[int](step_fn, 3)
% Print the length and index
stream$putl(po, "Cycle length: " || int$unparse(lam))
stream$putl(po, "Start index: " || int$unparse(mu))
% Print the cycle
stream$puts(po, "Cycle: ")
for i: int in iterfunc[int](step_fn, 3, lam, mu) do
stream$puts(po, int$unparse(i) || " ")
end
end start_up</syntaxhighlight>
{{out}}
<pre>3 10 101 2 5 26 167 95 101 2 5 26 167 95 101 2 5 26 167 95
Cycle length: 6
Start index: 2
Cycle: 101 2 5 26 167 95</pre>
=={{header|Cowgol}}==
<
typedef N is uint8;
Line 663 ⟶ 761:
print("Cycle length: "); print_i32(length as uint32); print_nl();
print("Cycle start: "); print_i32(start as uint32); print_nl();
PrintRange(x2_plus1_mod255, 3, start, length+start);</
{{out}}
<pre>3 10 101 2 5 26 167 95 101 2 5 26 167 95 101 2 5 26 167 95
Line 672 ⟶ 770:
=={{header|D}}==
{{trans|Java}}
<
import std.stdio;
Line 714 ⟶ 812:
auto iterate(int start, int function(int) f) {
return only(start).chain(generate!(() => start=f(start)));
}</
{{out}}
<pre>Cycle length: 6
Line 721 ⟶ 819:
=={{header|Elixir}}==
{{trans|Ruby}}
<
def find_cycle(x0, f) do
lambda = find_lambda(f, x0, f.(x0), 1, 1)
Line 751 ⟶ 849:
# Test the find_cycle function
{clength, cstart} = Cycle_detection.find_cycle(3, f)
IO.puts "Cycle length = #{clength}\nStart index = #{cstart}"</
{{out}}
Line 762 ⟶ 860:
=={{header|Factor}}==
This is a strict translation of the Python code from the Wikipedia article. Perhaps a more idiomatic version could be added in the future, although there is value in showing how Factor's lexical variables differ from most languages. A variable binding with <code>!</code> at the end is mutable, and subsequent uses of that name followed by <code>!</code> change the value of the variable to the value at the top of the data stack.
<
: cyclical-function ( n -- m ) dup * 1 + 255 mod ;
Line 792 ⟶ 890:
3 [ 20 [ dup , cyclical-function ] times ] { } make nip .
3 [ cyclical-function ] brent
"Cycle length: %d\nCycle start: %d\n" printf</
{{out}}
<pre>
Line 801 ⟶ 899:
=={{header|FOCAL}}==
<
01.20 S X=X0;F I=1,20;T X;D 2
01.30 D 3;T !" START",M,!,"LENGTH",L,!
Line 823 ⟶ 921:
03.80 I (T-H)3.9,3.99,3.9
03.90 S X=T;D 2;S T=X;S X=H;D 2;S H=X;S M=M+1;G 3.8
03.99 R</
{{out}}
<pre>= 3= 10= 101= 2= 5= 26= 167= 95= 101= 2= 5= 26= 167= 95= 101= 2= 5= 26= 167= 95
Line 832 ⟶ 930:
=={{header|Forth}}==
Works with gforth.
<syntaxhighlight lang="forth">
: cycle-length { x0 'f -- lambda } \ Brent's algorithm stage 1
1 1 x0 dup 'f execute
Line 869 ⟶ 967:
." The cycle is " 3 ' f(x) .cycle cr
bye
</syntaxhighlight>
{{Out}}
<pre>
Line 879 ⟶ 977:
===Brent's algorithm===
{{trans|Python}}
<
' compile with: fbc -s console
Line 950 ⟶ 1,048:
Print : Print "hit any key to end program"
Sleep
End</
{{out}}
<pre> Brent's algorithm
Line 959 ⟶ 1,057:
===Tortoise and hare. Floyd's algorithm===
{{trans|Wikipedia}}
<
' compile with: fbc -s console
Line 1,039 ⟶ 1,137:
Print : Print "hit any key to end program"
Sleep
End</
=={{header|Go}}==
Line 1,046 ⟶ 1,144:
Run it on the [https://play.golang.org/p/unOtxuwZfg go playground], or on [https://goplay.space/#S1pQZSuJci go play space].
<syntaxhighlight lang="go">
package main
Line 1,104 ⟶ 1,202:
}
fmt.Println("Cycle:", a[µ:µ+λ])
}</
{{out}}
Line 1,113 ⟶ 1,211:
=={{header|Groovy}}==
{{trans|Java}}
<
class CycleDetection {
Line 1,161 ⟶ 1,259:
println()
}
}</
{{out}}
<pre>Cycle length: 6
Line 1,172 ⟶ 1,270:
Haskellers, being able to handle conceptually infinite structures, traditionally consider totality as an important issue. The following solution is total for total inputs (modulo totality of iterated function) and allows nontermination only if input is explicitly infinite.
<
findCycle :: Eq a => [a] -> Maybe ([a], Int, Int)
Line 1,189 ⟶ 1,287:
| pow == lam = loop (2*pow) 1 y ys
| otherwise = loop pow (1+lam) x ys
in loop 1 1 x xs</
'''Examples'''
Line 1,218 ⟶ 1,316:
Brute force implementation:
<
r=. ~.@(,u@{:)^:_ y
n=. u{:r
(,(#r)-])r i. n
)</
In other words: iterate until we stop getting new values, find the repeated value, and see how it fits into the sequence.
Line 1,228 ⟶ 1,326:
Example use:
<
2 6</
(Note that 1 0 1 are the coefficients of the polynomial <code>1 + (0 * x) + (1 * x * x)</code> or, if you prefer <code>(1 * x ^ 0) + (0 * x ^ 1) + (1 * x ^ 2)</code> - it's easier and probably more efficient to just hand the coefficients to p. than it is to write out some other expression which produces the same result.)
Line 1,235 ⟶ 1,333:
=={{header|Java}}==
{{works with|Java|8}}
<
import static java.util.stream.IntStream.*;
Line 1,277 ⟶ 1,375:
.forEach(n -> System.out.printf("%s ", n));
}
}</
<pre>Cycle length: 6
Cycle: 101 2 5 26 167 95</pre>
=={{header|jq}}==
{{trans|Julia}}
{{works with|jq}}
'''Works with gojq, the Go implementation of jq'''
<syntaxhighlight lang="jq">def floyd(f; x0):
{tort: (x0|f)}
| .hare = (.tort|f)
| until( .tort == .hare;
.tort |= f
| .hare |= (f|f) )
| .mu = 0
| .tort = x0
| until( .tort == .hare;
.tort |= f
| .hare |= f
| .mu += 1)
| .lambda = 1
| .hare = (.tort|f)
| until (.tort == .hare;
.hare |= f
| .lambda += 1 )
| {lambda, mu} ;
def task(f; x0):
def skip($n; stream):
foreach stream as $s (0; .+1; select(. > $n) | $s);
floyd(f; x0)
| .,
"Cycle:",
skip(.mu; limit((.lambda + .mu); 3 | recurse(f)));
</syntaxhighlight>
'''The specific function and task'''
<syntaxhighlight lang="jq">
def f: (.*. + 1) % 255;
task(f;3)
</syntaxhighlight>
{{out}}
<pre>
{
"lambda": 6,
"mu": 2
}
Cycle:
3
10
101
2
5
26
167
95
</pre>
=={{header|Julia}}==
Line 1,287 ⟶ 1,440:
Following the Wikipedia article:
<
function floyd(f, x0)
Line 1,322 ⟶ 1,475:
x -> Iterators.take(x, λ) |>
collect
println("Cycle length: ", λ, "\nCycle start index: ", μ, "\nCycle: ", join(cycle, ", "))</
{{out}}
Line 1,330 ⟶ 1,483:
=={{header|Kotlin}}==
<
typealias IntToInt = (Int) -> Int
Line 1,378 ⟶ 1,531:
println("Start index = $mu")
println("Cycle = $cycle")
}</
{{out}}
Line 1,390 ⟶ 1,543:
=={{header|Lua}}==
Fairly direct translation of the Wikipedia code, except that the sequence is stored in a table and passed back as a third return value.
<
local tortoise, hare, mu = x0, f(x0), 0
local cycleTab, power, lam = {tortoise, hare}, 1, 1
Line 1,418 ⟶ 1,571:
print("Sequence:", table.concat(sequence, " "))
print("Cycle length:", cycleLength)
print("Start Index:", startIndex)</
{{out}}
<pre>Sequence: 3 10 101 2 5 26 167 95 101 2 5 26 167 95
Cycle length: 6
Start Index: 2</pre>
=={{header|Mathematica}}/{{header|Wolfram Language}}==
<syntaxhighlight lang="mathematica">s = {3, 10, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5,
26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101,
2, 5, 26, 167, 95, 101, 2, 5};
{transient, repeat} = FindTransientRepeat[s, 2];
Print["Starting index: ", Length[transient] + 1]
Print["Cycles: ", Floor[(Length[s] - Length[transient])/Length[repeat]]]</syntaxhighlight>
{{out}}
<pre>Starting index: 3
Cycles: 6</pre>
=={{header|Modula-2}}==
{{trans|Kotlin}}
<
FROM FormatString IMPORT FormatString;
FROM Terminal IMPORT WriteString,WriteLn,ReadChar;
Line 1,521 ⟶ 1,685:
ReadChar
END CycleDetection.</
=={{header|Nim}}==
Translation of Wikipedia Python program.
<syntaxhighlight lang="nim">import strutils, sugar
func brent(f: int -> int; x0: int): (int, int) =
# Main phase: search successive powers of two.
var
power, λ = 1
tortoise = x0
hare = f(x0)
while tortoise != hare:
if power == λ:
# Time to start a new power of two.
tortoise = hare
power *= 2
λ = 0
hare = f(hare)
inc λ
# Find the position of the first repetition of length λ.
tortoise = x0
hare = x0
for i in 0..<λ:
hare = f(hare)
# The distance between the hare and tortoise is now λ.
# Next, the hare and tortoise move at same speed until they agree.
var μ = 0
while tortoise != hare:
tortoise = f(tortoise)
hare = f(hare)
inc μ
result = (λ, μ)
when isMainModule:
func f(x: int): int = (x * x + 1) mod 255
let x0 = 3
let (λ, μ) = brent(f, x0)
echo "Cycle length: ", λ
echo "Cycle start index: ", μ
var cycle: seq[int]
var x = x0
for i in 0..<λ+μ:
if i >= μ: cycle.add x
x = f(x)
echo "Cycle: ", cycle.join(" ")</syntaxhighlight>
{{out}}
<pre>Cycle length: 6
Cycle start index: 2
Cycle: 101 2 5 26 167 95</pre>
=={{header|ooRexx}}==
<
x=3
list=x
Line 1,541 ⟶ 1,764:
End
Exit
f: Return (arg(1)**2+1)//255 </
{{out}}
<pre>
Line 1,552 ⟶ 1,775:
=={{header|Perl}}==
{{trans|Raku}}
<
sub cyclical_function { ($_[0] * $_[0] + 1) % 255 }
Line 1,600 ⟶ 1,823:
print "Cycle length $l\n";
print "Cycle start index $s\n";
print show_range($s,$s+$l-1) . "\n";</
{{out}}
<pre>3 10 101 2 5 26 167 95 101 2 5 26 167 95 101 2 5 26 167 95
Line 1,609 ⟶ 1,832:
=={{header|Phix}}==
Translation of the Wikipedia code, but using the more descriptive len and pos, instead of lambda and mu, and adding a limit.
<!--<syntaxhighlight lang="phix">(phixonline)-->
<span style="color: #008080;">function</span> <span style="color: #000000;">f</span><span style="color: #0000FF;">(</span><span style="color: #004080;">integer</span> <span style="color: #000000;">x</span><span style="color: #0000FF;">)</span>
<span style="color: #008080;">return</span> <span style="color: #7060A8;">mod</span><span style="color: #0000FF;">(</span><span style="color: #000000;">x</span><span style="color: #0000FF;">*</span><span style="color: #000000;">x</span><span style="color: #0000FF;">+</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #000000;">255</span><span style="color: #0000FF;">)</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">function</span>
<span style="color: #008080;">function</span> <span style="color: #000000;">brent</span><span style="color: #0000FF;">(</span><span style="color: #004080;">integer</span> <span style="color: #000000;">x0</span><span style="color: #0000FF;">)</span>
<span style="color: #004080;">integer</span> <span style="color: #000000;">pow2</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">1</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">len</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">1</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">pos</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">1</span>
<span style="color: #004080;">integer</span> <span style="color: #000000;">tortoise</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">x0</span><span style="color: #0000FF;">,</span>
<span style="color: #000000;">hare</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">f</span><span style="color: #0000FF;">(</span><span style="color: #000000;">x0</span><span style="color: #0000FF;">)</span>
<span style="color: #004080;">sequence</span> <span style="color: #000000;">s</span> <span style="color: #0000FF;">=</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">tortoise</span><span style="color: #0000FF;">,</span><span style="color: #000000;">hare</span><span style="color: #0000FF;">}</span> <span style="color: #000080;font-style:italic;">-- (kept for output only)
-- main phase: search successive powers of two</span>
<span style="color: #008080;">while</span> <span style="color: #000000;">tortoise</span><span style="color: #0000FF;">!=</span><span style="color: #000000;">hare</span> <span style="color: #008080;">do</span>
<span style="color: #008080;">if</span> <span style="color: #000000;">pow2</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">len</span> <span style="color: #008080;">then</span>
<span style="color: #000000;">tortoise</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">hare</span>
<span style="color: #000000;">pow2</span> <span style="color: #0000FF;">*=</span> <span style="color: #000000;">2</span>
<span style="color: #008080;">if</span> <span style="color: #000000;">pow2</span><span style="color: #0000FF;">></span><span style="color: #000000;">16</span> <span style="color: #008080;">then</span>
<span style="color: #008080;">return</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">s</span><span style="color: #0000FF;">,</span><span style="color: #000000;">0</span><span style="color: #0000FF;">,</span><span style="color: #000000;">0</span><span style="color: #0000FF;">}</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
<span style="color: #000000;">len</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">0</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">if</span>
<span style="color: #000000;">hare</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">f</span><span style="color: #0000FF;">(</span><span style="color: #000000;">hare</span><span style="color: #0000FF;">)</span>
<span style="color: #000000;">s</span> <span style="color: #0000FF;">&=</span> <span style="color: #000000;">hare</span>
<span style="color: #000000;">len</span> <span style="color: #0000FF;">+=</span> <span style="color: #000000;">1</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">while</span>
<span style="color: #000080;font-style:italic;">-- Find the position of the first repetition of length len</span>
<span style="color: #000000;">tortoise</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">x0</span>
<span style="color: #000000;">hare</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">x0</span>
<span style="color: #008080;">for</span> <span style="color: #000000;">i</span><span style="color: #0000FF;">=</span><span style="color: #000000;">1</span> <span style="color: #008080;">to</span> <span style="color: #000000;">len</span> <span style="color: #008080;">do</span>
<span style="color: #000000;">hare</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">f</span><span style="color: #0000FF;">(</span><span style="color: #000000;">hare</span><span style="color: #0000FF;">)</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">for</span>
<span style="color: #000080;font-style:italic;">-- The distance between the hare and tortoise is now len.
-- Next, the hare and tortoise move at same speed until they agree</span>
<span style="color: #008080;">while</span> <span style="color: #000000;">tortoise</span><span style="color: #0000FF;"><></span><span style="color: #000000;">hare</span> <span style="color: #008080;">do</span>
<span style="color: #000000;">tortoise</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">f</span><span style="color: #0000FF;">(</span><span style="color: #000000;">tortoise</span><span style="color: #0000FF;">)</span>
<span style="color: #000000;">hare</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">f</span><span style="color: #0000FF;">(</span><span style="color: #000000;">hare</span><span style="color: #0000FF;">)</span>
<span style="color: #000000;">pos</span> <span style="color: #0000FF;">+=</span> <span style="color: #000000;">1</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">while</span>
<span style="color: #008080;">return</span> <span style="color: #0000FF;">{</span><span style="color: #000000;">s</span><span style="color: #0000FF;">,</span><span style="color: #000000;">len</span><span style="color: #0000FF;">,</span><span style="color: #000000;">pos</span><span style="color: #0000FF;">}</span>
<span style="color: #008080;">end</span> <span style="color: #008080;">function</span>
<span style="color: #004080;">sequence</span> <span style="color: #000000;">s</span>
<span style="color: #004080;">integer</span> <span style="color: #000000;">len</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">pos</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">x0</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">3</span>
<span style="color: #0000FF;">{</span><span style="color: #000000;">s</span><span style="color: #0000FF;">,</span><span style="color: #000000;">len</span><span style="color: #0000FF;">,</span><span style="color: #000000;">pos</span><span style="color: #0000FF;">}</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">brent</span><span style="color: #0000FF;">(</span><span style="color: #000000;">x0</span><span style="color: #0000FF;">)</span>
<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">" Brent's algorithm\n"</span><span style="color: #0000FF;">)</span>
<span style="color: #0000FF;">?</span><span style="color: #000000;">s</span>
<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">" Cycle starts at position: %d (1-based)\n"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">pos</span><span style="color: #0000FF;">})</span>
<span style="color: #7060A8;">printf</span><span style="color: #0000FF;">(</span><span style="color: #000000;">1</span><span style="color: #0000FF;">,</span><span style="color: #008000;">" The length of the Cycle = %d\n"</span><span style="color: #0000FF;">,{</span><span style="color: #000000;">len</span><span style="color: #0000FF;">})</span>
<span style="color: #0000FF;">?</span><span style="color: #000000;">s</span><span style="color: #0000FF;">[</span><span style="color: #000000;">pos</span><span style="color: #0000FF;">..</span><span style="color: #000000;">pos</span><span style="color: #0000FF;">+</span><span style="color: #000000;">len</span><span style="color: #0000FF;">-</span><span style="color: #000000;">1</span><span style="color: #0000FF;">]</span>
<!--</syntaxhighlight>-->
{{out}}
<pre>
Line 1,670 ⟶ 1,895:
=={{header|PL/M}}==
<
/* CP/M CALLS */
BDOS: PROCEDURE (FN, ARG); DECLARE FN BYTE, ARG ADDRESS; GO TO 5; END BDOS;
Line 1,767 ⟶ 1,992:
CALL PRINT$RANGE(.F$ADDR, 3, MU, MU+LAM);
CALL EXIT;
EOF</
{{out}}
<pre>3 10 101 2 5 26 167 95 101 2 5 26 167 95 101 2 5 26 167 95
Line 1,776 ⟶ 2,001:
===Procedural===
Function from the Wikipedia article:
<
def brent(f, x0):
Line 1,818 ⟶ 2,043:
print("Cycle length: %d" % lam)
print("Cycle start index: %d" % mu)
print("Cycle: %s" % list(itertools.islice(iterate(f, x0), mu, mu+lam)))</
{{out}}
<pre>
Line 1,827 ⟶ 2,052:
A modified version of the above where the first stage is restructured for clarity:
<
def brent_length(f, x0):
Line 1,872 ⟶ 2,097:
print("Cycle length: %d" % lam)
print("Cycle start index: %d" % mu)
print("Cycle: %s" % list(itertools.islice(iterate(f, x0), mu, mu+lam)))</
{{out}}
<pre>Cycle length: 6
Line 1,882 ⟶ 2,107:
{{Trans|Haskell}}
{{Works with|Python|3.7}}
<
from itertools import (chain, cycle, islice)
Line 2,115 ⟶ 2,340:
# MAIN ---
if __name__ == '__main__':
main()</
{{Out}}
<pre>First cycle detected, if any:
Line 2,129 ⟶ 2,354:
Recursive ''until'':
<
def until(p):
'''The result of repeatedly applying f until p holds.
Line 2,135 ⟶ 2,360:
def go(f, x):
return x if p(x) else go(f, f(x))
return lambda f: lambda x: go(f, x)</
''cycleLength'' refactored in terms of ''until'':
<
def cycleLength(xs):
'''Just the length of the first cycle found,
Line 2,164 ⟶ 2,389:
) if ys else Nothing()
else:
return Nothing()</
Iterative reimplementation of ''until'':
<
def until_(p):
'''The result of repeatedly applying f until p holds.
Line 2,176 ⟶ 2,401:
v = f(v)
return v
return lambda f: lambda x: go(f, x)</
Line 2,182 ⟶ 2,407:
The Python no longer falls out of the tree at the sight of an ouroboros, and we can happily search for cycles in lists of several thousand items:
{{Works with|Python|3.7}}
<
from itertools import (chain, cycle, islice)
Line 2,437 ⟶ 2,662:
# MAIN ---
if __name__ == '__main__':
main()</
{{Out}}
<pre>First cycle detected, if any:
Line 2,445 ⟶ 2,670:
[1..10000] -> No cycle found
f(x) = (x*x + 1) modulo 255 -> [101,2,5,26,167,95] (from:2, length:6)</pre>
=={{header|Quackery}}==
<syntaxhighlight lang="quackery"> [ stack ] is fun ( --> s )
[ stack ] is pow ( --> s )
[ stack ] is len ( --> s )
[ fun put
1 pow put
1 len put
dup fun share do
[ 2dup != while
len share
pow share = if
[ nip dup
pow share
pow tally
0 len replace ]
fun share do
1 len tally
again ]
2drop
fun release
pow release
len take ] is cyclelen ( n x --> n )
[ 0 temp put
dip [ fun put dup ]
times [ fun share do ]
[ 2dup != while
fun share do
dip [ fun share do ]
1 temp tally
again ]
2drop
fun release
temp take ] is cyclepos ( n x n --> n )
[ 2dup cyclelen
dup dip cyclepos ] is brent ( n x --> n n )
[ 2dup
20 times
[ over echo sp
tuck do swap ]
cr cr
2drop
brent
say "cycle length is "
echo cr
say "cycle starts at "
echo ] is task ( n x --> )
3 ' [ 2 ** 1+ 255 mod ] task</syntaxhighlight>
{{out}}
<pre>3 10 101 2 5 26 167 95 101 2 5 26 167 95 101 2 5 26 167 95
cycle length is 6
cycle starts at 2
</pre>
=={{header|Racket}}==
I feel a bit bad about overloading λ, but it’s in the spirit of the algorithm.
<
#lang racket/base
Line 2,488 ⟶ 2,776:
(let-values (([µ λ] (brent f 3)))
(printf "Cycle length = ~a~%Start Index = ~a~%" µ λ)))
</syntaxhighlight>
{{out}}
Line 2,502 ⟶ 2,790:
Pretty much a line for line translation of the Python code on the Wikipedia page.
<syntaxhighlight lang="raku"
my ( $l, $s ) = brent( &cyclical-function, 3 );
Line 2,535 ⟶ 2,823:
}
return $λ, $μ;
}</
{{out}}
<pre>3, 10, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, ...
Line 2,544 ⟶ 2,832:
=={{header|REXX}}==
===Brent's algorithm===
<
init= 3; $= init
do until length($)>79; $= $ f( word($, words($) ) )
Line 2,577 ⟶ 2,865:
return # mu
/*──────────────────────────────────────────────────────────────────────────────────────*/
f: return ( arg(1) **2 + 1) // 255 /*this defines/executes the function F.*/</
{{out|output|text= when using the defaults:}}
<pre>
Line 2,588 ⟶ 2,876:
===sequential search algorithm===
<
x= 3; $= x /*initial couple of variables*/
do until cycle\==0; x= f(x) /*calculate another number. */
Line 2,601 ⟶ 2,889:
exit /*stick a fork in it, we're all done. */
/*──────────────────────────────────────────────────────────────────────────────────────*/
f: return ( arg(1) **2 + 1) // 255 /*this defines/executes the function F.*/</
{{out|output|:}}
<pre>
Line 2,613 ⟶ 2,901:
===hash table algorithm===
This REXX version is a lot faster (than the sequential search algorithm) if the ''cycle length'' and/or ''start index'' is large.
<
!.= .; !.x= 1; x= 3; $= x /*assign initial value. */
do #=1+words($); x= f(x); $= $ x /*add the number to list.*/
Line 2,626 ⟶ 2,914:
exit /*stick a fork in it, we're all done. */
/*──────────────────────────────────────────────────────────────────────────────────────*/
f: return ( arg(1) **2 + 1) // 255 /*this defines/executes the function F.*/</
{{out|output|text= is identical to the 2<sup>nd</sup> REXX version.}} <br><br>
Line 2,639 ⟶ 2,927:
<br>test the hash table algorithm. A divisor which is <big> ''two raised to the 49<sup>th</sup> power'' </big> was chosen; it
<br>generates a cyclic sequence that contains over 1.5 million numbers.
<
parse arg power . /*obtain optional args from C.L. */
if power=='' | power="," then power=8 /*Not specified? Use the default*/
Line 2,663 ⟶ 2,951:
exit /*stick a fork in it, we're all done. */
/*──────────────────────────────────────────────────────────────────────────────────────*/
f: return ( arg(1) **2 + 1) // 255 /*this defines/executes the function F.*/</
{{out|output|text= when the input (power of two) used is: <tt> 49 </tt>}}
<pre>
Line 2,679 ⟶ 2,967:
There is more information in the "hash table"<br>
and f has no "side effect".
<
x=3; list=x; p.=0; p.x=1
Do q=2 By 1
Line 2,696 ⟶ 2,984:
Exit
/*-------------------------------------------------------------------*/
f: Return (arg(1)**2+1)// 255; /*define the function F*/</
=={{header|RPL}}==
Translation of the Brent algorithm given in Wikipedia.
{{works with|HP|48}}
≪ 1 1 0 → f x0 power lam mu
≪ x0 DUP f EVAL <span style="color:grey">@ Main phase: search successive powers of two</span>
'''WHILE''' DUP2 ≠ '''REPEAT'''
'''IF''' power lam == '''THEN''' <span style="color:grey">@ time to start a new power of two?</span>
SWAP DROP DUP
2 'power' STO*
0 'lam' STO
'''END'''
f EVAL
1 'lam' STO+
'''END'''
DROP2 x0 DUP <span style="color:grey">@ Find the position of the first repetition of length λ</span>
0 lam 1 - '''START'''
f EVAL '''NEXT''' <span style="color:grey">@ distance between the hare and tortoise is now λ</span>
'''WHILE''' DUP2 ≠ '''REPEAT''' <span style="color:grey">@ the hare and tortoise move at same speed until they agree</span>
f EVAL SWAP
f EVAL SWAP
1 'mu' STO+
'''END'''
DROP2 lam mu
≫ ≫ '<span style="color:blue">CYCLEB</span>' STO
≪ SQ 1 + 255 MOD ≫ 0 <span style="color:blue">CYCLEB</span>
{{out}}
<pre>
2: 6
1: 2
</pre>
=={{header|Ruby}}==
{{works with|ruby|2.0}}
<
# Purpose:
# Find the cycle length and start position of a numerical seried using Brent's cycle algorithm.
Line 2,758 ⟶ 3,078:
# Test the findCycle function
clength, cstart = findCycle(3) { |x| f(x) }
puts "Cycle length = #{clength}\nStart index = #{cstart}"</
{{out}}
Line 2,770 ⟶ 3,090:
=== Procedural ===
{{Out}}Best seen in running your browser either by [https://scalafiddle.io/sf/6O7WjnO/0 ScalaFiddle (ES aka JavaScript, non JVM)] or [https://scastie.scala-lang.org/kPCg0fxOQQCZPkOnmMR0Kg Scastie (remote JVM)].
<
def brent(f: Int => Int, x0: Int): (Int, Int) = {
Line 2,814 ⟶ 3,134:
println(s"Cycle = ${cycle.force}")
}</
=== Functional ===
<syntaxhighlight lang="scala">
import scala.annotation.tailrec
Line 2,880 ⟶ 3,200:
}
</syntaxhighlight>
{{out}}
Line 2,891 ⟶ 3,211:
=={{header|Sidef}}==
{{trans|Raku}}
<
var power = 1
var λ = 1
Line 2,934 ⟶ 3,254:
say "Cycle length #{l}.";
say "Cycle start index #{s}."
say [seq[s .. (s + l - 1)]]</
{{out}}
<pre>3, 10, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, 101, 2, 5, 26, 167, 95, ...
Line 2,943 ⟶ 3,263:
=={{header|Visual Basic .NET}}==
{{trans|C#}}
<
Function FindCycle(Of T As IEquatable(Of T))(x0 As T, yielder As Func(Of T, T)) As Tuple(Of Integer, Integer)
Line 3,001 ⟶ 3,321:
End Sub
End Module</
{{out}}
<pre>3,10,101,2,5,26,167,95,101,2,5,26,167,95,101,2,5,26,167,95,101,2,5,26,167,95,101,2,5,26,167,95,101,2,5,26,167,95,101,2,5
Line 3,009 ⟶ 3,329:
=={{header|Wren}}==
Working from the code in the Wikipedia article:
<
var lam = 1
var power = 1
Line 3,048 ⟶ 3,368:
System.print("Cycle length = %(lam)")
System.print("Start index = %(mu)")
System.print("Cycle = %(seq[mu...mu+lam])")</
{{out}}
Line 3,060 ⟶ 3,380:
=={{header|zkl}}==
Algorithm from the Wikipedia
<
# main phase: search successive powers of two
power:=lam:=1;
Line 3,082 ⟶ 3,402:
}
return(lam,mu);
}</
<
{{out}}
<pre>
|