例 x^(23^50)+1=(x+1)*Product(Sum((-x)^(i*m^j),i=0..22),j=0..49)

> m:=9:
ithprime_m:=ithprime(m);
n:=50;
power:=ithprime_m^n;

[Maple Math]

[Maple Math]

[Maple Math]

> LHS:=x^(ithprime(m)^n)+1;
RHS:=(x+1)*Product(Sum((-x)^(i*ithprime(m)^j),i=0..ithprime(m)-1),j=0..n-1);

[Maple Math]

[Maple Math]

> LHS_RHS=normal(expand(value(LHS-RHS)));

[Maple Math]

>