n = 25 (m = 12): trajectories parallel to the side X_0 X_(n-1); exact arithmetic in Q(zeta_50)
j= 1  simple preclosed piece: 2 segments, length 3.9372 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 21 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 2  simple preclosed piece: 2 segments, length 7.6269 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 17 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 3  simple preclosed piece: 2 segments, length 10.8375 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 13 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 4  simple preclosed piece: 2 segments, length 13.3671 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  9 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 5  simple preclosed piece: 2 segments, length 15.0568 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  5 = -4j mod n; factor  5 = n/gcd(n,j); closed trajectory: 10 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 6  simple preclosed piece: 2 segments, length 15.8004 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  1 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 7  simple preclosed piece: 2 segments, length 15.5512 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 22 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 8  simple preclosed piece: 2 segments, length 14.3249 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 18 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j= 9  simple preclosed piece: 2 segments, length 12.1985 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 14 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=10  simple preclosed piece: 2 segments, length 9.3056 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b = 10 = -4j mod n; factor  5 = n/gcd(n,j); closed trajectory: 10 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=11  simple preclosed piece: 2 segments, length 5.8280 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  6 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
j=12  simple preclosed piece: 2 segments, length 1.9842 = lambda sin(2j pi/n)/sin(pi/n) (exact); x_2 = r^b x_0 with b =  2 = -4j mod n; factor 25 = n/gcd(n,j); closed trajectory: 50 segments   [t = 3/10, 1/7, 5/9; midpoint t = 1/2: 1 segment, half the length, same factor]
factors d_1..d_m: [25, 25, 25, 25, 5, 25, 25, 25, 25, 5, 25, 25] ; gcd = 5 ; a quotient equals n: False
DONE [0.7s]
