I am working with variables over GF(2) and when I list the summands of certain polynomials, I want to have a list of all distinct polynomials, but it does not work.
from sage.crypto.sbox import SBox
from sage.crypto.boolean_function import BooleanFunction
F=GF(2^4,'a');
V=F.vector_space();
V=sorted(V);
sF=sorted(F);
S= SBox(12,5,6,11,9,0,10,13,3,14,15,8,4,7,1,2); #Present SBox
f=S.interpolation_polynomial(); #univariate form
L=[S.component_function(i).algebraic_normal_form() for i in [1,2,4,8]];
L # S=(f1,f2,f3,f4)
#L=[x0 + x1*x2 + x2 + x3,
x0*x1*x2 + x0*x1*x3 + x0*x2*x3 + x1*x3 + x1 + x2*x3 + x3,
x0*x1*x3 + x0*x1 + x0*x2*x3 + x0*x3 + x1*x3 + x2 + x3 + 1,
x0*x1*x2 + x0*x1*x3 + x0*x2*x3 + x0 + x1*x2 + x1 + x3 + 1]
S=set(srange(4)); #[0,1,2,3]
C2=Combinations(S,2).list();
C1=Combinations(S,1).list();
A=list(var('a%d' % i,domain='integer') for i in srange(len(C2)+len(C1)));
R.<x0,x1,x2,x3>=GF(2)[];
X=[x0,x1,x2,x3];
def g(X):
return [X[i[0]]*X[i[1]] for i in C2]+[X[i[0]] for i in C1];
L_S=[g(L)[i] for i in [0..len(A)-1]];
L_X=[g(X)[i] for i in [0..len(A)-1]];
M0=[(L_S[i]+L_X[i]) for i in [0..len(A)-1]]
Now, this is what M0 looks like:
[x0*x1*x2 + x0*x1*x3 + x0*x3 + x2*x3 + x3,
x0*x1 + x0*x3 + x0 + x1*x2*x3 + x1*x3,
x0*x1*x2 + x0*x1 + x0*x2*x3 + x0*x2 + x0*x3 + x1*x2 + x1*x3 + x2*x3 + x2,
x0*x1*x2 + x0*x1*x3 + x0*x1 + x0*x3 + x1,
x0*x1*x3 + x0*x1 + x0*x2*x3 + x0*x3 + x1*x2,
x0*x1*x3 + x0*x1 + x0*x2 + x0 + x1*x2 + x1 + x2 + x3 + 1,
x1*x2 + x2 + x3,
x0*x1*x2 + x0*x1*x3 + x0*x2*x3 + x1*x3 + x2*x3 + x3,
x0*x1*x3 + x0*x1 + x0*x2*x3 + x0*x3 + x1*x3 + x3 + 1,
x0*x1*x2 + x0*x1*x3 + x0*x2*x3 + x0 + x1*x2 + x1 + 1]
When I flatten and use set() I get that set(flatten(M0)) equals:
{1,
1,
x0,
x0,
x0,
x0*x1,
x0*x1,
x0*x1,
x0*x1*x2,
x0*x1*x2,
x0*x1*x2,
x0*x1*x3,
x0*x1*x3,
x0*x1*x3,
x0*x1*x3,
x0*x2,
x0*x2,
x0*x2*x3,
x0*x2*x3,
x0*x2*x3,
x0*x2*x3,
x0*x3,
x0*x3,
x0*x3,
x1,
x1,
x1,
x1*x2,
x1*x2,
x1*x2,
x1*x2,
x1*x2*x3,
x1*x3,
x1*x3,
x1*x3,
x2,
x2,
x2*x3,
x2*x3,
x3,
x3,
x3}
Why are there duplicate values if I have set() operation?
