# Perform polynomial factorization?

Posted 10 months ago
444 Views
|
0 Replies
|
0 Total Likes
|
 Hello! I have this function: Clear[chernC]; chernC[0] = 1; chernC[1] = Subscript[ch, 1]; chernC[k_Integer /; k > 0] := chernC[k] = Simplify[(1/k)* Total[Table[(-1)^(j + 1)*Factorial[j]*Subscript[ch, j]* chernC[k - j], {j, 1, k}]]]; which when called chernC[2] gives an output like this: $\frac{1}{2}(ch_1^2-2ch_2)$, where $ch$ is the variable of the polynomials. I am trying to write it in a nice way (factor out the common terms, write all the number in prime factorization form etc.) the way I do here for another function: http://community.wolfram.com/groups/-/m/t/1305165My code that tries to do this is this: primeFactorForm[n_] := If[Length@# == 1, First@#, CenterDot @@ #] &[ Superscript @@@ FactorInteger[n]] string = With[{order = "Lexicographic"}, StringJoin[ Riffle[Table[ poly = MonomialList[chernC[i], Table[Subscript[ch, k], {k, 1, 9}], order]; gcd = GCD @@ poly /. Rational[n_, d_]*c_ :> d; ToString[ Inactive[Set][Subscript[ch, i], DisplayForm@ RowBox[{1/primeFactorForm[gcd], "(", RowBox[(List @@ Distribute[gcd*poly]) /. {Times[Rational[n_, d_], e__] :> RowBox[{If[n/d > 0, "+", "-"], primeFactorForm[Abs@n]/primeFactorForm[Abs@d], e}], Times[n_Integer, e__] :> RowBox[{If[n > 0, "+", "-"], primeFactorForm[Abs@n], e}]}], ")"}]], TeXForm], {i, 0, 9}], "\\\\"]]] CopyToClipboard[string] For the previous function it works great, but now it just gives weird output, for example the gcd function appears in the output explicitly, without actually computing the numerical gcd value. Can someone tell me how to fix this? Thank you!