//---[CLO, Using ..., p. 445], Reed-Solomon [8,4,8-4+1]
//--- si possono correggere 2 errori
//--- c = x^7 + 2*x^5 + x^2 + 2x +1
//--- y = x^7 + a*x^5 + (a-1)*x^2 + 2*x + 1 
//--- posizioni di informazione: x^7,x^6,x^5,x^4
//--- posizioni di controllo: x^3,x^2,x^1,1
ring s = (3,a),x,(lp,c); minpoly = a^2+a+2;
s;
//   characteristic : 3
//   1 parameter    : a 
//   minpoly        : (a2+a-1)
//   number of vars : 1
//        block   1 : ordering lp
//                  : names    x
//        block   2 : ordering c
 poly w = x^7 + 2*x^5;
proc codifica (poly w)
{
poly g = (x-a)*(x-a^2)*(x-a^3)*(x-a^4);
return(w - reduce(w,g));
};
c = codifica(w);
poly c = codifica(w);
c;
x7-x5+x2-x+1
//---decodifica di y
poly y = x^7 + a*x^5 + (a+2)*x^2 + 2*x + 1;
proc decodifica (poly y)
{
poly syndrome = 
subst(y,x,a)+subst(y,x,a^2)*x+subst(y,x,a^3)*x^2+subst(y,x,a^4)*x^3;
vector v1 = [x^4,0]; // 2*t = 4, d = 5
vector v2 = [syndrome,1];
module K = v1, v2;
option(redSB);
module G = std(K); //base ridotta di Groebner per il modulo K
if (G[1]<G[2]) 
 {vector dG = G[1]/subst(G[1][2],x,0); 
  poly omega = dG[1];
  poly lambda = dG[2];} //lambda = polinomio delle posizioni degli errori
else
 {vector dG = G[2]/subst(G[2][2],x,0);
  poly omega = dG[1];
  poly lambda = dG[2];};
//---dividere per il termine costante di lambda
int i;
list I;
//---I = insieme delle posizioni degli errori  
for (i=1; i<=8; i++) 
{if (subst(lambda,x,a^i) == 0) {I=insert(I,8-i);}};
poly e; //calcolo del polinomio errore e con la formula di Forney
for (i=1; i<=size(I); i++)
{poly chi(i) = subst(lambda/(1-a^I[i]*x),x,a^(-I[i]));
 e = e + subst(omega,x,a^(-I[i]))/(a^I[i]*chi(i))*x^I[i];
};
return(y - e);
};
c;
x7-x5+x2-x+1
y;
x7+(a)*x5+(a-1)*x2-x+1
decodifica(y);
x7-x5+x2-x+1
monitor("");