Pseudo-random F generator

PHI := proc(x, y)

local m, i, j, s, t;

    m := y;

    s := x;

    for j to 160 do

        i := BBS2(s, n);

        if m mod 2 = 0 then s := i[1]; m := 1/2*m

        else s := i[2]; m := 1/2*m - 1/2

        fi

    od;

    if m mod 2 = 0 then RETURN(i[1]) else RETURN(i[2]) fi

end

 

> PHI(x,0);

1357469249910833676661224042358551080614102081812314467127237047431286956132944238889493

 

> PHI(x,1);

117398909908260911131996112842895804188668870134335923337271692925313320578737719081841

 

> PHI(x,2);

1329411526465534016320144291726413409421572390831016521897300010679573651210026786041686

 

> PHI(x,4);

1258298423181230303566383830747046385865664693955246483692588981548888926435119030755904

 

> PHI(x,64);

1487457608820828772469699800723028681412606791663963919484258889644608112329664706894429

 

> PHI(x,65);

287067700226218067318776122090885459316550216733650028037186059295752238894619562693658

 

> PHI(x,75876532123567988798965332132);

1936993187289509274181501852942573712429893873164958009583340602403704048575994095483107