rsk3[p0_, m_, b_, k_] := Module[{j, r, p}, j = 0; p = p0; Label[o1]; p = seek4[p, b, m]; j = j + 1; If[j < k, Goto[o1]]; Return[p]; ] seek4[pb_, j_, k_] := Module[{m, pq, a1, b1, a2, b2, x, y}, pq = pertc[pb]; m = 1; Label[o1]; {a1, b1} = vq[pb,vv1]; {a2, b2} = vq[pq,vv1]; x = j*(a2 - a1) + (b2 - b1); y = a2 - a1 + b2 - b1; If[x < 0 && y < 0, Goto[o3]]; m = m + 1; If[m < k, Goto[o2]]; Goto[o4]; Label[o2]; pq = 2*pq + pb; pq = rat[pq, 20]; pq = N[pq, 20]; pq = pq/s4sum[pq]; Goto[o1]; Label[o4]; pq = pb; Label[o3]; Return[pq]; ] pertc[pa_] := Module[{m, p1, p2, p3, p4, x, y, z, w}, {p1, p2, p3, p4} = pa; x = R[]; y = R[]; z = R[]; w = R[]; z = R[]; m = {1 + x, 1 + y, 1 + z, 1 + w}; x = m*pa; y = s4sum[x]; z = x/y; Return[z]; ] R[] := Module[{s1, s2, s3}, s1 = Ra[]; s2 = Ra[]; s3 = Ra[]; s1 = 100000*s1 + s2; s3 = 100000*s1 + s3; s2 = s3/10^15; s1 = N[s2, 17]; Return[s1]; ] Ra[] := Random[Integer, 100000] s4sum[kk_] := Module[{r, f1, f2, f3, f4}, {f1, f2, f3, f4} = kk; r = f1 + f2 + f3 + f4; Return[r]; ] vq[pd_, vvd_] := vqc[xc[pd, vvd]] vqc[ks_] := Module[{r, g1, g2, g3, g4, g5, g6, h1, h2, h3, h4, h5, s, m}, {g1, g2, g3, g4, g5, g6} = ks; r = g1; h1 = g1 - g2; h2 = g2 - g3; h3 = g3 - g4; h4 = g4 - g5; h5 = g5 - g6; m = {h1, h2, h3, h4, h5}; s = Sort[m, Greater]; {g1, g2, g3, g4, g5} = s; m = {r, g3}; Return[m]; ] xc[pc_, vvc_] := Module[{r, j1, j2, j3, j4, j5, j6, bl7, bl8, bl9, bl10, bl11, bl12, bl13, bl14}, {j1, j2, j3, j4, j5, j6, bl7, bl8, bl9, bl10, bl11, bl12, bl13, bl14} = xb[pc, vvc]; r = {j1, j2, j3, j4, j5, j6}; Return[r]; ] xb[pb_, vvb_] := Sort[xa[pb, vvb], Greater] xa[u_, v_] := xa4[u, v] xa4[pa_, vva_] := Module[{r, w1, w2, w3, w4, w12, w13, w23, w14, w24, w34, w123, w124, w134, w234, q1, q2, q3, q4}, {q1, q2, q3, q4} = pa; {w1, w2, w3, w4, w12, w13, w23, w14, w24, w34, w123, w124, w134, w234} = vva; r = {w1 - q1, w2 - q2, w3 - q3, w4 - q4, w12 - q1 - q2, w13 - q1 - q3, w23 - q2 - q3, w14 - q1 - q4, w24 - q2 - q4, w34 - q3 - q4, w123 - q1 - q2 - q3, w124 - q1 - q2 - q4, w134 - q1 - q3 - q4, w234 - q2 - q3 - q4}; Return[r]; ] bet[pa_, pb_] := fp[xsf[pb] - xsf[pa]] fp[s_] := firpos[s] firpos[kka_] := Module[{k1, k2, k3, k4, k5, k6}, {k1, k2, k3, k4, k5, k6, k7, k8, k9, k10, k11, k12, k13, k14} = kka; If[k1 < 0, Goto[o1]]; If[k1 > 0, Goto[o2]]; If[k2 < 0, Goto[o1]]; If[k2 > 0, Goto[o2]]; If[k3 < 0, Goto[o1]]; If[k3 > 0, Goto[o2]]; If[k4 < 0, Goto[o1]]; If[k4 > 0, Goto[o2]]; If[k5 < 0, Goto[o1]]; If[k5 > 0, Goto[o2]]; If[k6 < 0, Goto[o1]]; If[k6 > 0, Goto[o2]]; If[k7 < 0, Goto[o1]]; If[k7 > 0, Goto[o2]]; If[k8 < 0, Goto[o1]]; If[k8 > 0, Goto[o2]]; If[k9 < 0, Goto[o1]]; If[k9 > 0, Goto[o2]]; If[k10 < 0, Goto[o1]]; If[k10 > 0, Goto[o2]]; If[k11 < 0, Goto[o1]]; If[k11 > 0, Goto[o2]]; If[k12 < 0, Goto[o1]]; If[k12 > 0, Goto[o2]]; If[k13 < 0, Goto[o1]]; If[k13 > 0, Goto[o2]]; If[k14 < 0, Goto[o1]]; If[k14 > 0, Goto[o2]]; Label[o1]; Return[False]; Label[o2]; Return[True]; ] xsf[px_] := xb[px, vv1] rat[x_, k_Integer] := Rationalize[x, 1/10^(k + 2)]