Кстати, надо было мне сразу чуть большее расширение взять, тогда бы многие вопросы сами собой отпали:
Код:
> alias(alpha=RootOf(cyclotomic(51,z))): simplify( subs(alpha=exp(2*Pi*I/51),factors(x^16 -x^15 - 16*x^14 +16*x^13 +103*x^12 -103*x^11 -339*x^10 +339*x^9+596*x^8-596*x^7-526*x^6+526*x^5+188*x^4-188*x^3-16*x^2+16*x+1,alpha)) );
[1, [
[2*cos(1/51*Pi)+x, 1],
[x+exp(2/17*I*Pi)-exp(28/51*I*Pi)+exp(40/51*I*Pi), 1],
[1/2*I*3^(1/2)-2*I*sin(1/51*Pi)-2*I*sin(7/51*Pi)+2*I*sin(11/51*Pi)+x+1/2-exp(2/51*I*Pi)+exp(2/17*I*Pi)-exp(8/51*I*Pi)+exp(4/17*I*Pi)-exp(14/51*I*Pi)+exp(6/17*I*Pi)-exp(20/51*I*Pi)+exp(8/17*I*Pi)-exp(26/51*I*Pi)+exp(10/17*I*Pi)-exp(32/51*I*Pi)+2*I*sin(5/51*Pi)-2*exp(38/51*I*Pi), 1],
[x-exp(14/51*I*Pi)+exp(20/51*I*Pi)+exp(-16/17*I*Pi), 1],
[x+exp(2/51*I*Pi)-exp(32/51*I*Pi)+exp(12/17*I*Pi), 1],
[x+exp(8/51*I*Pi)-exp(26/51*I*Pi)+exp(14/17*I*Pi), 1],
[x-exp(10/51*I*Pi)+exp(8/17*I*Pi)+exp(-44/51*I*Pi), 1],
[2*cos(3/17*Pi)+2*cos(1/17*Pi)+x-1-exp(2/51*I*Pi)+exp(4/51*I*Pi)-exp(2/17*I*Pi)+exp(10/51*I*Pi)-exp(4/17*I*Pi)+exp(16/51*I*Pi)-exp(6/17*I*Pi)+exp(22/51*I*Pi)-exp(8/17*I*Pi)+exp(28/51*I*Pi)-exp(10/17*I*Pi)+exp(32/51*I*Pi)-exp(12/17*I*Pi)+exp(38/51*I*Pi)+exp(-46/51*I*Pi)+exp(44/51*I*Pi)+exp(-40/51*I*Pi)+exp(50/51*I*Pi), 1],
[x+exp(14/51*I*Pi)-exp(20/51*I*Pi)+exp(16/17*I*Pi), 1],
[-1/2*I*3^(1/2)+2*I*sin(1/51*Pi)+2*I*sin(7/51*Pi)-2*I*sin(11/51*Pi)+x-1/2+exp(2/51*I*Pi)-exp(4/51*I*Pi)-exp(2/17*I*Pi)+exp(8/51*I*Pi)-exp(4/17*I*Pi)+exp(14/51*I*Pi)-exp(6/17*I*Pi)+exp(20/51*I*Pi)-exp(8/17*I*Pi)+exp(26/51*I*Pi)+exp(32/51*I*Pi)-2*I*sin(5/51*Pi)+exp(38/51*I*Pi), 1],
[x+exp(4/17*I*Pi)-exp(22/51*I*Pi)+exp(46/51*I*Pi), 1],
[x-exp(8/51*I*Pi)+exp(26/51*I*Pi)+exp(-14/17*I*Pi), 1],
[2*cos(5/51*Pi)+x, 1],
[2*cos(11/51*Pi)+x, 1],
[x-exp(16/51*I*Pi)+exp(6/17*I*Pi)+exp(-50/51*I*Pi), 1],
[2*cos(7/51*Pi)+x, 1]
]]