Me funksionin GAP
MaxOnes := n -> Maximum(List(Filtered(AsList(GF(2)^[n,n]),
M->not ForAny(Tuples([1..n-1],2),
s->ForAny(Cartesian([1..n-s[1]],[1..n-s[2]]),
t->ForAll(Flat(M{[s[1],s[1]+t[1]]}
{[s[2],s[2]+t[2]]}),
IsOne)))),
N->Number(Flat(N),IsOne)));
ne kemi per shembull
gap> MaxOnes(1);
1
gap> MaxOnes(2);
3
gap> MaxOnes(3);
6
gap> MaxOnes(4);
9
dhe pergjigja e pyetjes eshte
gap> MaxOnes(43);
( ... ju lutem prisni ... ju lutem prisni ... !! )