Quelle polyfinf.gi Sprache: unbekannt

#############################################################################
##
## This file is part of GAP, a system for computational discrete algebra.
## This file's authors include Frank Celler, Alexander Hulpke.
##
## Copyright of GAP belongs to its developers, whose names are too numerous
## to list here. Please refer to the COPYRIGHT file for details.
##
## SPDX-License-Identifier: GPL-2.0-or-later
##
## This file contains functions for polynomials over finite fields
##

#############################################################################
##
#F FactorsCommonDegreePol( <R>, <f>, <d> ) . . . . . . . . . . . . . factors
##
## <f> must be a square free product of irreducible factors of degree <d>
## and leading coefficient 1. <R> must be a polynomial ring over a finite
## field of size p^k.
##
InstallGlobalFunction(FactorsCommonDegreePol,function( R, f, d )
local c, ind, br, g, h, k, i,dou;

 # if <f> has degree 0, return f
 dou:=DegreeOfLaurentPolynomial(f);
 if dou<d then
 return [];
 # if <f> has degree <d>, return irreducible <f>
 elif dou=d then
 return [f];
 fi;

 c := CoefficientsOfLaurentPolynomial(f);
 ind := IndeterminateNumberOfLaurentPolynomial(f);
 br := CoefficientsRing(R);

 # if <f> has a trivial constant term signal an error
 if c[2] <> 0 then
 Error("<f> must have a non-trivial constant term");
 fi;

 # choose a random polynomial <g> of degree less than 2*<d>
 repeat
 g := RandomPol(br,2*d-1,ind);
 until DegreeOfLaurentPolynomial(g)<>DEGREE_ZERO_LAURPOL;

 # if p = 2 take <g> + <g>^2 + <g>^(2^2) + ... + <g>^(2^(k*<d>-1))
 if Characteristic(br) = 2 then
 g := CoefficientsOfLaurentPolynomial(g);
 h := ShiftedCoeffs(g[1],g[2]);
 k := ShiftedCoeffs(c[1],c[2]);
 g := g[1];
 for i in [1..DegreeOverPrimeField(br)*d-1] do
 g := ProductCoeffs(g,g);
 ReduceCoeffs(g,k);
 ShrinkRowVector(g);
 AddCoeffs(h,g);
 od;
 h := LaurentPolynomialByCoefficients(
 FamilyObj(h[1]), h, 0, ind );

 # if p > 2 take <g> ^ ((p ^ (k*<d>) - 1) / 2) - 1
 else
 h:=PowerMod(g,(Characteristic(br)^(DegreeOverPrimeField(br)*d)-1)/2,f)
 - One(br);
 fi;

 # gcd of <f> and <h> is with probability > 1/2 a proper factor
 g := GcdOp(f,h);
 return Concatenation(
 FactorsCommonDegreePol( R, Quotient(R,f,g), d ),
 FactorsCommonDegreePol( R, g, d ) );
end);

#############################################################################
##
#M FactorsSquarefree( <R>, <f>, <opt> ) . . . . . . . . . . . . . . factors
##
## <f> must be square free and must have leading coefficient 1. <R> must be
## a polynomial ring over a finite field of size q.
##
InstallMethod( FactorsSquarefree,"univariate polynomial over finite field",
 true, [ IsFiniteFieldPolynomialRing, IsUnivariatePolynomial, IsRecord ],0,
function( R, f, opt )
local br, ind, c, facs, deg, px, cyc, gcd,d,powc,fc,fam;

 br := CoefficientsRing(R);
 ind := IndeterminateNumberOfLaurentPolynomial(f);
 c := CoefficientsOfLaurentPolynomial(f);

 # if <f> has a trivial constant term signal an error
 if c[2] <> 0 then
 Error("<f> must have a non-trivial constant term");
 fi;

 # <facs> will contain factorisation
 facs := [];

 # in the following <pow> = x ^ (q ^ (<deg>+1))
 deg := 0;
 #px := LaurentPolynomialByExtRepNC(
 # FamilyObj(f), [One(br)],1, ind );
 # pow := px;
 if IsFinite(br) and IsField(br) and Size(br)>MAXSIZE_GF_INTERNAL then
 px:=Immutable([Zero(br),-One(br)]);
 else
 px:=ImmutableVector(br, [Zero(br),-One(br)]);
 fi;
 powc:=-px;
 fc:=CoefficientsOfLaurentPolynomial(f)[1];
 fam:=FamilyObj(One(br));

 # while <f> could still have two irreducible factors
 while 2*(deg+1) <= DegreeOfLaurentPolynomial(f) do

 #pow := PowerMod(pow,Size(br),f);
 powc:=PowerModCoeffs(powc,Length(powc),Size(br),fc,Length(fc));
 # next degree and next cyclotomic polynomial x^(q^(<deg>+1))-x
 deg := deg + 1;
 if not IsBound(opt.onlydegs) or deg in opt.onlydegs then
 #cyc := pow - px;
 cyc:=ShallowCopy(powc);
 AddCoeffs(cyc,px);
 cyc:=LaurentPolynomialByCoefficients(fam,cyc,0,ind);
 # compute the gcd of <f> and <cyc>
 gcd := GcdOp( f, cyc );

 # split the gcd with 'FactorsCommonDegree'
 d:=DegreeOfLaurentPolynomial(gcd);
 if 0<d and d>=deg then
 Info(InfoPoly,3,"Factor Common Deg.",deg );
 Append(facs,FactorsCommonDegreePol(R,gcd,deg));
 f := Quotient(f,gcd);
 fi;
 fi;
 od;

 # if necessary add irreducible <f> to the list of factors
 if 0 < DegreeOfLaurentPolynomial(f) then
 Add(facs,f);
 fi;

 # return the factorisation
 return facs;

end );

#############################################################################
##
#F RootsRepresentativeFFPol( <R>, <f>, <n> )
##
InstallGlobalFunction(RootsRepresentativeFFPol,function( R, f, n )
local r, br, nu, ind, p, d, z, v, o, i, e;

 r := [];
 br := CoefficientsRing(R);
 nu := Zero(br);
 ind := IndeterminateNumberOfLaurentPolynomial(f);

 p := Characteristic(br);
 d := DegreeOverPrimeField(br);
 z := PrimitiveRoot(br);
 f := CoefficientsOfLaurentPolynomial(f);
 v := f[2];
 f := f[1];
 o := p^d-1;
 for i in [0..(Length(f)-1)/n] do
 e := f[i*n+1];
 if e = nu then
 r[i+1] := nu;
 else
 r[i+1] := z ^ ((LogFFE(e,z) / n) mod o);
 fi;
 od;
 return LaurentPolynomialByCoefficients(
 FamilyObj(nu), r, v/n, ind );

end);

#############################################################################
##
#M Factors( <R>, <f> ) . . . . . . . . . . . . . . factors of <f>
##
InstallMethod( Factors, "polynomial over a finite field",
 IsCollsElms, [ IsFiniteFieldPolynomialRing, IsUnivariatePolynomial ],0,

function(R,f)
local cr, opt, irf, i, ind, v, l, g, k, d,
 facs, h, q, char, r,
 gc, hc, fam, val;

 # parse the arguments
 cr := CoefficientsRing(R);

 opt:=ValueOption("factoroptions");
 PushOptions(rec(factoroptions:=rec())); # options do not hold for
 # subsequent factorizations
 if opt=fail then
 opt:=rec();
 fi;

 # check if we already know a factorisation
 irf := IrrFacsPol(f);
 i := PositionProperty( irf, i -> i[1] = cr );
 if i <> fail then
 PopOptions();
 return ShallowCopy(irf[i][2]);
 fi;

 # handle the trivial cases
 ind := IndeterminateNumberOfLaurentPolynomial(f);
 v := CoefficientsOfLaurentPolynomial(f);
 #fam := FamilyObj(v[1][1]);

 if DegreeOfLaurentPolynomial(f) < 2
 or DegreeOfLaurentPolynomial(f)=DEGREE_ZERO_LAURPOL then
 Add( irf, [cr,[f]] );
 PopOptions();
 return [f];

 elif Length(v[1]) = 1 then
 l:= ListWithIdenticalEntries( v[2],
 IndeterminateOfUnivariateRationalFunction( f ) );
 l[1] := l[1]*v[1][1];
 Add( irf, [cr,l] );
 PopOptions();
 return l;
 fi;

 # make the polynomial normed, remember the leading coefficient for later
 l:=LeadingCoefficient(f);
 #g := StandardAssociate(R,f);
 #l := Quotient(f,g);

 v := CoefficientsOfLaurentPolynomial(f);
 if v[2]=0 then
 k:=1/l*f;
 else
 k:=LaurentPolynomialByExtRepNC( FamilyObj(f), 1/l*v[1],0, ind );
 fi;
 v := v[2];

 # compute the derivative
 d := Derivative(k);

 # if the derivative is nonzero then $k / Gcd(k,d)$ is squarefree
 if not IsZero(d) then

 # compute the gcd of <k> and the derivative <d>
 g := GcdOp( k, d );
 if DegreeOfLaurentPolynomial(g)>0 then

 # factor the squarefree quotient and the remainder
 facs := FactorsSquarefree( R, Quotient(k,g), opt );
 else
 facs := FactorsSquarefree( R, k, opt );
 fi;

 if not (IsBound(opt.onlydegs) or IsBound(opt.stopdegs)) then
 # tell the factors they are factors
 for h in facs do
 StoreFactorsPol(cr,h,[h]);
 od;
 fi;

 if DegreeOfLaurentPolynomial(g)>0 then
 # Above w computed a square free factorization of k/g = k/k' in facs.
 # Now determine how often each factor h in facs divides g, by repeatedly
 # dividing g by h.
 #
 # The code for this is basically equivalent to the commented out code
 # snippet below; however, it avoids converting coefficient lists to
 # polynomials and back again during each loop iteration. This has a
 # significant performance effect if the polynomial is divided by a
 # larger power of a given divisor polynomial.
 #
 # for h in ShallowCopy(facs) do
 # q := Quotient( g, h );
 # while q <> fail do
 # Add( facs, h );
 # g := q;
 # q := Quotient( g, h );
 # od;
 # od;

 # convert g to coefficients list
 fam := FamilyObj( g );
 gc := CoefficientsOfLaurentPolynomial( g );
 if gc[2] > 0 then
 gc := ShiftedCoeffs(gc[1],gc[2]);
 else
 # the call to ShallowCopy below is necessary to ensure gc is mutable,
 # so that QUOTREM_LAURPOLS_LISTS can modify it
 gc := ShallowCopy(gc[1]);
 fi;

 # perform repeated divisions by the factors in facs
 for h in ShallowCopy(facs) do
 # convert h to coefficients list
 hc := CoefficientsOfLaurentPolynomial(h);
 if hc[2] > 0 then
 hc := ShiftedCoeffs(hc[1], hc[2]);
 else
 # calling ShallowCopy here is not necessary, but results in a
 # considerable speed boost
 hc := ShallowCopy(hc[1]);
 fi;

 # divide g by h as long as there is no remainder
 while true do
 # perform the actual division; since QUOTREM_LAURPOLS_LISTS modifies
 # its first argument, we pass a copy of gc in; this is necessary to
 # correctly handle the final iteration, in which division by h fails
 q := QUOTREM_LAURPOLS_LISTS( ShallowCopy(gc), hc );
 if not IsZero( q[2] ) then break; fi;
 Add( facs, h );
 gc := q[1];
 od;
 od;

 # convert coefficients list back into a polynomial
 val := RemoveOuterCoeffs( gc, fam!.zeroCoefficient );
 g := LaurentPolynomialByExtRepNC( fam, gc, val, ind );

 fi;
 if 0=DegreeOfLaurentPolynomial(g) then
 if not IsOne(g) then
 facs[1]:=facs[1]*g;
 fi;
 else
#T how shall this ever happen?
 Append( facs, Factors(R,g:factoroptions:=opt) );
 fi;

 # otherwise <k> is the -th power of another polynomial <r>
 else

 # compute the -th root of <f>
 char := Characteristic(cr);
 r := RootsRepresentativeFFPol( R, k, char );

 # factor this polynomial
 h := Factors( R, r: factoroptions:=opt );

 # each factor appears times in <k>
 facs := [];
 for i in [ 1 .. char ] do
 Append( facs, h );
 od;

 fi;

 # Sort the factorization
 Sort(facs);
 if v>0 then
 ind := IndeterminateOfUnivariateRationalFunction(f);
 facs:=Concatenation(List( [ 1 .. v ], x -> ind ),facs );
 fi;

 # return the factorization and store it
 if not IsOne(l) then
 facs[1] := facs[1]*l;
 fi;
 if not (IsBound(opt.onlydegs) or IsBound(opt.stopdegs)) then
 StoreFactorsPol(cr,f,facs);
 fi;
 Assert(2,Product(facs)=f);
 PopOptions();
 return facs;

end);

#############################################################################
##
#F ProductPP( <l>, <r> ) . . . . . . . . . . . . product of two prime powers
##
BindGlobal("ProductPP",function( l, r )
 local res, p1, p2, ps, p, i, n;

 if IsEmpty(l) then
 return r;
 elif IsEmpty(r) then
 return l;
 fi;
 res := [];
 p1 := l{ 2 * [ 1 .. Length( l ) / 2 ] - 1 };
 p2 := r{ 2 * [ 1 .. Length( r ) / 2 ] - 1 };
 ps := Set( Union( p1, p2 ) );
 for p in ps do
 n := 0;
 Add( res, p );
 i := Position( p1, p );
 if i <> fail then
 n := l[ 2*i ];
 fi;
 i := Position( p2, p );
 if i <> fail then
 n := n + r[ 2*i ];
 fi;
 Add( res, n );
 od;
 return res;

end);

#############################################################################
##
#F LcmPP( <l>, <r> ) . . . . . . . . . . . . lcm of prime powers <l> and <r>
##
BindGlobal("LcmPP",function( l, r )
 local res, p1, p2, ps, p, i, n;

 if l = [] then
 return r;
 elif r = [] then
 return l;
 fi;
 res := [];
 p1 := l{ 2 * [ 1 .. Length( l ) / 2 ] - 1 };
 p2 := r{ 2 * [ 1 .. Length( r ) / 2 ] - 1 };
 ps := Set( Union( p1, p2 ) );
 for p in ps do
 n := 0;
 Add( res, p );
 i := Position( p1, p );
 if i <> false then
 n := l[ 2*i ];
 fi;
 i := Position( p2, p );
 if i <> false and n < r[ 2*i ] then
 n := r[ 2*i ];
 fi;
 Add( res, n );
 od;
 return res;

end);

#############################################################################
##
#F FFPPowerModCheck(<g>, <pp>, <f> ) . . . . . . . . . . . . . . local
##
BindGlobal("FFPPowerModCheck",function( g, pp, f )
local qq, i;

 qq := [];
 for i in [ 1 .. Length(pp)/2 ] do
 Add( qq, pp[2*i-1] );
 Add( qq, pp[2*i] );
 g := PowerMod( g, pp[2*i-1] ^ pp[2*i], f );
 if DegreeOfLaurentPolynomial(g) = 0 then
 return [ g, qq ];
 fi;
 od;
 return [ g, qq ];

end);

#############################################################################
##
#F OrderKnownDividendList( <l>, <pp> ) . . . . . . . . . . . . . . . . local
##
## Computes an integer n such that OnSets( <l>, n ) contains only one
## element e. <pp> must be a list of prime powers of an integer d such that
## n divides d. The functions returns the integer n and the element e.
##
InstallGlobalFunction(OrderKnownDividendList,function( l, pp )
local pp1, # first half of <pp>
 pp2, # second half of <pp>
 a, # half exponent of first prime power
 k, # power of <l>
 o, o1, # computed order of <k>
 i; # loop

 # if <pp> contains no element return order 1
 if Length(pp) = 0 then
 return [ 1, l[1] ];

 # if <l> contains only one element return order 1
 elif Length(l) = 1 then
 return [ 1, l[1] ];

 # if the dividend is a prime return
 elif Length(pp) = 2 and pp[2] = 1 then
 return [ pp[1], l[1]^pp[1] ];

 # if the dividend is a prime power divide and conquer
 elif Length(pp) = 2 then
 pp := ShallowCopy(pp);
 a := QuoInt( pp[2], 2 );
 k := OnSets( l, pp[1]^a );

 # if <k> is trivial try smaller dividend
 if Length(k) = 1 then
 pp[2] := a;
 return OrderKnownDividendList( l, pp );

 # otherwise try to find order of <h>
 else
 pp[2] := pp[2] - a;
 o := OrderKnownDividendList( k, pp );
 return [ pp[1]^a*o[1], o[2] ];
 fi;

 # split different primes into two parts
 else
 a := 2 * QuoInt( Length(pp), 4 );
 pp1 := pp{[ 1 .. a ]};
 pp2 := pp{[ a+1 .. Length(pp) ]};

 # compute the order of <l>^<pp1>
 k := l;
 for i in [ 1 .. Length(pp2)/2 ] do
 k := OnSets( k, pp2[2*i-1]^pp2[2*i] );
 od;
 o1 := OrderKnownDividendList( k, pp1 );

 # compute the order of <l>^<o1> and return
 o := OrderKnownDividendList( OnSets( l, o1[1] ), pp2 );
 return [ o1[1]*o[1], o[2] ];
 fi;

end);

#############################################################################
##
#F FFPOrderKnownDividend( <R>, <g>, <f>, <pp> ) . . . . . . . . . . . local
##
## Computes an integer n such that <g>^n = const mod <f> where <g> and <f>
## are polynomials in <R> and <pp> is list of prime powers of an integer d
## such that n divides d. The functions returns the integer n and the
## element const.
##
InstallGlobalFunction(FFPOrderKnownDividend,function ( R, g, f, pp )
local l, a, h, n1, pp1, pp2, k, o, q;

 #Info( InfoPoly, 3, "FFPOrderKnownDividend started with:" );
 #Info( InfoPoly, 3, " <g> = ", g );
 #Info( InfoPoly, 3, " <f> = ", f );
 #Info( InfoPoly, 3, " <pp> = ", pp );

 # if <g> is constant return order 1
 if 0 = DegreeOfLaurentPolynomial(g) then
 #Info( InfoPoly, 3, " <g> is constant" );
 l := CoefficientsOfUnivariatePolynomial(g);
 l := [ 1, l[1] ];
 #Info( InfoPoly, 3, "FFPOrderKnownDividend returns ", l );
 return l;

 # if the dividend is a prime, we must compute g^pp[1] to get the constant
 elif Length(pp) = 2 and pp[2] = 1 then
 k := PowerMod( g, pp[1], f );
 l := CoefficientsOfUnivariatePolynomial(k);
 l := [ pp[1], l[1] ];
 #Info( InfoPoly, 3, "FFPOrderKnownDividend returns ", l );
 return l;

 # if the dividend is a prime power find the necessary power
 elif Length(pp) = 2 then
 #Info( InfoPoly, 3, "prime power, divide and conquer" );
 pp := ShallowCopy( pp );
 a := QuoInt( pp[2], 2 );
 q := pp[1] ^ a;
 h := PowerMod( g, q, f );

 # if <h> is constant try again with smaller dividend
 if 0 = DegreeOfLaurentPolynomial(h) then
 pp[2] := a;
 o := FFPOrderKnownDividend( R, g, f, pp );
 else
 pp[2] := pp[2] - a;
 l := FFPOrderKnownDividend( R, h, f, pp );
 o := [ q*l[1], l[2] ];
 fi;
 #Info( InfoPoly, 3, "FFPOrderKnownDividend returns ", o );
 return o;

 # split different primes.
 else

 # divide primes
 #Info( InfoPoly, 3, " ", Length(pp)/2, " different primes" );
 n1 := QuoInt( Length(pp), 4 );
 pp1 := pp{[ n1*2+1 .. Length(pp) ]};
 pp2 := pp{[ 1 .. n1*2 ]};
 #Info( InfoPoly, 3, " <pp1> = ", pp1 );
 #Info( InfoPoly, 3, " <pp2> = ", pp2 );

 # raise <g> to the power <pp2>
 k := FFPPowerModCheck( g, pp2, f );
 pp2 := k[2];
 k := k[1];

 # compute order for <pp1>
 o := FFPOrderKnownDividend( R, k, f, pp1 );

 # compute order for <pp2>
 k := PowerMod( g, o[1], f );
 l := FFPOrderKnownDividend( R, k, f, pp2 );
 o := [ o[1]*l[1], l[2] ];
 #Info( InfoPoly, 3, "FFPOrderKnownDividend returns ", o );
 return o;
 fi;

end);

#############################################################################
##
#F FFPUpperBoundOrder( <R>, <f> ) . . . . . . . . . . . . . . . . . . local
##
## Computes the irreducible factors f_i of a polynomial <f> over a field
## with p^n elements. It returns a list l of quadruples (f_i,a_i,pp_i,pb_i) such
## that the p-part of x mod f_i is p^a_i and the p'-part divides d_i for
## which the prime powers pp_i and not-yet-prime powers pb_i (in case
## factorization fails) are given.
##
BindGlobal("FFPUpperBoundOrder",function( R, f )
local fs, F, L, phi, B, i, d, pp, a, deg,t,pb;

 # factorize <f> into irreducible factors
 fs := Collected( Factors( R, f ) );

 # get the field over which the polynomials are written
 F := CoefficientsRing(R);

 # <phi>(m) gives ( minpol of 1^(1/m) )( F.char )
 # cache values
 if not IsBound(F!.FFPUBOVAL) then
 F!.FFPUBOVAL:=[ [PrimePowersInt( Characteristic(F)-1 ),[]] ];
 fi;

 L:=F!.FFPUBOVAL;
 phi := function( m )
 local x, pp, a, good,bad, d, i,primes;
 if not IsBound( L[m] ) then
 bad:=[];
 x := Characteristic(F)^m-1;
 primes:=PrimeDivisors(x); # use the Cunningham tables, and then store the prime
 # factors such that they end up cached below. In fact, since we
 # only divide off some known primes, this factorization really
 # shouldn't be harder than the one below.
 for d in Difference( DivisorsInt( m ), [m] ) do
 pp := phi( d );
 if Length(pp[2])>0 then
 bad:=ProductPP(pp[2],bad);
 fi;
 pp:=pp[1]; # nothing bad can happen here as d is small
 for i in [ 1 .. Length(pp)/2 ] do
 x := x / pp[2*i-1]^pp[2*i];
 od;
 od;
 a := PrimePowersInt( x:quiet );
 good:=[];
 for i in [1,3..Length(a)-1] do
 if a[i] in primes # we assume that the factorization above really gave
 # prime factors.
 or IsPrimeInt(a[i]) then
 Add(good,a[i]);
 Add(good,a[i+1]);
 else
 Add(bad,a[i]);
 Add(bad,a[i+1]);
 fi;
 od;
 good:=[good,bad];
 if Length(good[1])<Length(a) then
 Info(InfoWarning,1,"disregarded nonfactorable",bad);
 else
 L[m]:=good;
 fi;
 else
 good:=L[m];
 fi;
 return good;
 end;

 # compute a_i and pp_i
 B := [];
 for i in [ 1 .. Length(fs) ] do

 # p-part is p^Roof(Log_p(e_i)) where e_i is the multiplicity of f_i
 a := 0;
 if fs[i][2] > 1 then
 a := 1+LogInt(fs[i][2]-1,Characteristic(F));
 fi;

 # p'-part: (p^n)^d_i-1/(p^n-1) where d_i is the degree of f_i
 d := DegreeOfLaurentPolynomial(fs[i][1]);
 pp := [];
 pb:=[];
 deg := DegreeOverPrimeField(F);
 for f in DivisorsInt( d*deg ) do
 if deg mod f <> 0 then
 t:=phi(f);
 pp := ProductPP( t[1], pp );
 pb := ProductPP( t[2], pb );
 fi;
 od;

 # add <a> and <pp> to 
 Add( B, [fs[i][1],a,pp,pb] );
 od;

 # OK, that's it
 return B;

end);

#############################################################################
##
#M ProjectiveOrder( <f> ) . . . . . . . . . . . . . projective order of <f>
##
## Return an integer n and a finite field element const such that x^n =
## const mod <f>.
##
InstallOtherMethod( ProjectiveOrder,
 "divide and conquer for univariate polynomials", true,
 [ IsUnivariatePolynomial ],0,
function( f )
local v, R, U, x, O, n, g, q, o, bas;

 # <f> must not be divisible by x.
 v := CoefficientsOfLaurentPolynomial(f);
 if 0 < v[2] then
 Error( "<f> must have a non zero constant term" );
 fi;

 # if degree is zero, return
 if 0 = DegreeOfLaurentPolynomial(f) then
 return [ 1, v[1][1] ];
 fi;

 # use 'UpperBoundOrder' to split <f> into irreducibles
 #R := DefaultRing(f);
 R:=PolynomialRing(
 DefaultField(CoefficientsOfUnivariateLaurentPolynomial(f)[1]),
 [IndeterminateNumberOfLaurentPolynomial(f)]);
 U := FFPUpperBoundOrder( R, f );

 # run through the irrducibles and compute their order
 x := IndeterminateOfUnivariateRationalFunction(f);
 O := [];
 n := 1;
 for g in U do
 if Length(g[3])=0 and Length(g[4])>0 then
 # in this case `FFPOrderKnownDividend' might run in an infinite
 # recursion.
 Error("cannot compute order due to limits in the integer factorization!");
 fi;
 #o := FFPOrderKnownDividend(R,EuclideanRemainder(R,x,g[1]),g[1],g[3]);
 bas:=QuotRemLaurpols(x,g[1],2);
 o := FFPOrderKnownDividend(R,bas,g[1],g[3]);
 if Length(g[4])>0 then
 q:=DegreeOfLaurentPolynomial(PowerMod(bas,o[1],g[1]));
 if not(q=0 or q=DEGREE_ZERO_LAURPOL) then
 # in fact x^o[1] is not congruent to a constant -- we really need the
 # primes.
 Error("cannot compute order due to limits in the integer factorization!");
 fi;
 fi;
 q := Characteristic(CoefficientsRing(R))^g[2];
 n := LcmInt( n, o[1]*q );
 Add( O, [ o[1]*q, o[2]^q ] );
 od;

 # try to get the same constant in each block
 U := [];
 q := Size( CoefficientsRing(R) ) - 1;
 for g in O do
 AddSet( U, g[2]^((n/g[1]) mod q) );
 od;

 # return the order <n> times the order of 
 o := OrderKnownDividendList( U, PrimePowersInt(q) );
 return [ n*o[1], o[2] ];

end );

RedispatchOnCondition(ProjectiveOrder,true,
 [IsRationalFunction],[IsUnivariatePolynomial],0);

#############################################################################
##
#M SplittingField( <f> )
##
InstallMethod( SplittingField,"finite field polynomials",true,
 [ IsUnivariatePolynomial ],0,
function( f )
local c,b,l;
 if Characteristic(f)=0 then
 TryNextMethod();
 fi;
 c:=CoefficientsOfUnivariatePolynomial(f);
 if Length(c)=0 then
 return GF(Characteristic(f));
 fi;
 b:=DefaultField(c);
 l:=List(Factors(PolynomialRing(b),f),DegreeOfLaurentPolynomial);
 l:=Filtered(l,i->i>0); # lcm otherwise returns 0
 l:=Lcm(l);
 return GF(Characteristic(f)^(l*DegreeOverPrimeField(b)));
end);

Quelle polyfinf.gi Sprache: unbekannt

[ Dauer der Verarbeitung: 0.49 Sekunden (vorverarbeitet) ]