def coset_unrank(self, rank, af=False):
"""
unrank using Schreier-Sims representation
coset_unrank is the inverse operation of coset_rank
if 0 <= rank < order; otherwise it returns None.
"""
u = self.coset_repr()
if rank < 0 or rank >= self.order():
return None
un = self._coset_repr_n
base = self._base
m = len(u)
nb = len(base)
assert nb == len(un)
v = [0] * m
for i in range(nb - 1, -1, -1):
j = base[i]
rank, c = divmod(rank, un[i])
v[j] = c
a = [u[i][v[i]] for i in range(m)]
h = perm_af_muln(*a)
if af:
return h
else:
return _new_from_array_form(h)
def is_normal(self, gr):
"""
test if G=self is a normal subgroup of gr
G is normal in gr if
for each g2 in G, g1 in gr, g = g1*g2*g1**-1 belongs to G
It is sufficient to check this for each g1 in gr.generator and
g2 g2 in G.generator
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([1, 2, 0])
>>> b = Permutation([1, 0, 2])
>>> G = PermutationGroup([a, b])
>>> G1 = PermutationGroup([a, Permutation([2, 0, 1])])
>>> G1.is_normal(G)
True
"""
gens2 = [p.array_form for p in self.generators]
gens1 = [p.array_form for p in gr.generators]
for g1 in gens1:
for g2 in gens2:
p = perm_af_muln(g1, g2, perm_af_invert(g1))
if not self.coset_decomposition(p):
return False
return True
def coset_unrank(self, rank, af=False):
"""
unrank using Schreier-Sims representation
coset_unrank is the inverse operation of coset_rank
if 0 <= rank < order; otherwise it returns None.
"""
u = self.coset_repr()
if rank < 0 or rank >= self.order():
return None
un = self._coset_repr_n
base = self._base
m = len(u)
nb = len(base)
assert nb == len(un)
v = [0]*m
for i in range(nb-1, -1,-1):
j = base[i]
rank, c = divmod(rank, un[i])
v[j] = c
a = [u[i][v[i]] for i in range(m)]
h = perm_af_muln(*a)
if af:
return h
else:
return _new_from_array_form(h)
def is_normal(self, gr):
"""
test if G=self is a normal subgroup of gr
G is normal in gr if
for each g2 in G, g1 in gr, g = g1*g2*g1**-1 belongs to G
It is sufficient to check this for each g1 in gr.generator and
g2 g2 in G.generator
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([1, 2, 0])
>>> b = Permutation([1, 0, 2])
>>> G = PermutationGroup([a, b])
>>> G1 = PermutationGroup([a, Permutation([2, 0, 1])])
>>> G1.is_normal(G)
True
"""
gens2 = [p.array_form for p in self.generators]
gens1 = [p.array_form for p in gr.generators]
for g1 in gens1:
for g2 in gens2:
p = perm_af_muln(g1, g2, perm_af_invert(g1))
if not self.coset_decomposition(p):
return False
return True
def coset_rank(self, g):
"""
rank using Schreier-Sims representation
The coset rank of `g` is the ordering number in which
it appears in the lexicographic listing according to the
coset decomposition, see coset_decomposition;
the ordering is the same as in G.generate(method='coset').
If `g` does not belong to the group it returns None
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([[0, 1, 3, 7, 6, 4], [2, 5]])
>>> b = Permutation([[0, 1, 3, 2], [4, 5, 7, 6]])
>>> G = PermutationGroup([a, b])
>>> c = Permutation([[0, 1, 2, 3, 4], [5, 6, 7]])
>>> G.coset_rank(c)
>>> c = Permutation([[0, 6], [1, 7], [2, 4], [3, 5]])
>>> G.coset_rank(c)
40
>>> G.coset_unrank(40, af=True)
[6, 7, 4, 5, 2, 3, 0, 1]
"""
u = self.coset_repr()
if isinstance(g, Permutation):
g = g.array_form
g1 = g
m = len(u)
a = []
un = self._coset_repr_n
n = self.degree
rank = 0
base = [1]
for i in un[m:0:-1]:
base.append(base[-1]*i)
base.reverse()
a1 = [0]*m
i1 = -1
for i in self._base:
i1 += 1
x = g1[i]
for j, h in enumerate(u[i]):
if h[i] == x:
a.append(h)
a1[i] = j
rank += j*base[i1]
p2 = perm_af_invert(h)
g1 = perm_af_mul(p2, g1)
break
else:
return None
if perm_af_muln(*a) == g:
return rank
return None
def coset_rank(self, g):
"""
rank using Schreier-Sims representation
The coset rank of `g` is the ordering number in which
it appears in the lexicographic listing according to the
coset decomposition, see coset_decomposition;
the ordering is the same as in G.generate(method='coset').
If `g` does not belong to the group it returns None
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([[0, 1, 3, 7, 6, 4], [2, 5]])
>>> b = Permutation([[0, 1, 3, 2], [4, 5, 7, 6]])
>>> G = PermutationGroup([a, b])
>>> c = Permutation([[0, 1, 2, 3, 4], [5, 6, 7]])
>>> G.coset_rank(c)
>>> c = Permutation([[0, 6], [1, 7], [2, 4], [3, 5]])
>>> G.coset_rank(c)
40
>>> G.coset_unrank(40, af=True)
[6, 7, 4, 5, 2, 3, 0, 1]
"""
u = self.coset_repr()
if isinstance(g, Permutation):
g = g.array_form
g1 = g
m = len(u)
a = []
un = self._coset_repr_n
n = self.degree
rank = 0
base = [1]
for i in un[m:0:-1]:
base.append(base[-1] * i)
base.reverse()
a1 = [0] * m
i1 = -1
for i in self._base:
i1 += 1
x = g1[i]
for j, h in enumerate(u[i]):
if h[i] == x:
a.append(h)
a1[i] = j
rank += j * base[i1]
p2 = perm_af_invert(h)
g1 = perm_af_mul(p2, g1)
break
else:
return None
if perm_af_muln(*a) == g:
return rank
return None
def test_muln():
n = 6
m = 8
a = [Permutation.unrank_nonlex(n, i).array_form for i in range(m)]
h = range(n)
for i in range(m):
h = perm_af_mul(h, a[i])
h2 = perm_af_muln(*a[:i+1])
assert h == h2
def test_muln():
n = 6
m = 8
a = [Permutation.unrank_nonlex(n, i).array_form for i in range(m)]
h = range(n)
for i in range(m):
h = perm_af_mul(h, a[i])
h2 = perm_af_muln(*a[:i + 1])
assert h == h2
def stabilizer(self, alpha):
"""
return the stabilizer subgroup leaving alpha fixed.
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([0, 2, 1, 3])
>>> b = Permutation([2, 0, 1, 3])
>>> g = PermutationGroup([a, b])
>>> g2 = g.stabilizer(2)
>>> g2
PermutationGroup([Permutation([1, 0, 2, 3])])
"""
if alpha == 0:
self.schreier_sims()
gens = self._stabilizers_gens
gens = [_new_from_array_form(p) for p in gens]
return PermutationGroup(gens)
# h[alpha] = 0; h[0] = alpha
n = self.degree
h = range(n)
h[0] = alpha
h[alpha] = 0
genv = [p.array_form for p in self.generators]
# conjugate the group to breng alpha to 0
gens = [_new_from_array_form(perm_af_muln(h, p, h)) for p in genv]
G = PermutationGroup(gens)
G.schreier_sims()
# stabilizers for 0
gens1 = G._stabilizers_gens
# conjugate the group
gens2 = [_new_from_array_form(perm_af_muln(h, p, h)) for p in gens1]
return PermutationGroup(gens2)
def stabilizer(self, alpha):
"""
return the stabilizer subgroup leaving alpha fixed.
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([0, 2, 1, 3])
>>> b = Permutation([2, 0, 1, 3])
>>> g = PermutationGroup([a, b])
>>> g2 = g.stabilizer(2)
>>> g2
PermutationGroup([Permutation([1, 0, 2, 3])])
"""
if alpha == 0:
self.schreier_sims()
gens = self._stabilizers_gens
gens = [_new_from_array_form(p) for p in gens]
return PermutationGroup(gens)
# h[alpha] = 0; h[0] = alpha
n = self.degree
h = range(n)
h[0] = alpha
h[alpha] = 0
genv = [p.array_form for p in self.generators]
# conjugate the group to breng alpha to 0
gens = [_new_from_array_form(perm_af_muln(h, p, h)) for p in genv]
G = PermutationGroup(gens)
G.schreier_sims()
# stabilizers for 0
gens1 = G._stabilizers_gens
# conjugate the group
gens2 = [_new_from_array_form(perm_af_muln(h, p, h)) for p in gens1]
return PermutationGroup(gens2)
def coset_decomposition(self, g):
"""
Decompose `g` as h_0*...*h_{len(u)}
The Schreier-Sims coset representation u of `G`
gives a univoque decomposition of an element `g`
as h_0*...*h_{len(u)}, where h_i belongs to u[i]
Output: [h_0, .., h_{len(u)}] if `g` belongs to `G`
False otherwise
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([[0, 1, 3, 7, 6, 4], [2, 5]])
>>> b = Permutation([[0, 1, 3, 2], [4, 5, 7, 6]])
>>> G = PermutationGroup([a, b])
>>> c = Permutation([[0, 1, 2, 3, 4], [5, 6, 7]])
>>> G.coset_decomposition(c)
False
>>> c = Permutation([[0, 6], [1, 7], [2, 4], [3, 5]])
>>> G.coset_decomposition(c)
[[6, 4, 2, 0, 7, 5, 3, 1], [0, 4, 1, 5, 2, 6, 3, 7], [0, 1, 2, 3, 4, 5, 6, 7]]
>>> G.has_element(c)
True
"""
u = self.coset_repr()
if isinstance(g, Permutation):
g = g.array_form
g1 = g
n = len(u)
a = []
for i in range(n):
x = g1[i]
for h in u[i]:
if h[i] == x:
a.append(h)
p2 = perm_af_invert(h)
g1 = perm_af_mul(p2, g1)
break
else:
return False
if perm_af_muln(*a) == g:
return a
return False
def insert(self, g, alpha):
"""
insert permutation `g` in stabilizer chain at point alpha
"""
n = len(g)
if not g == self.idn:
vertex = self.vertex
jg = self.jg
i = _smallest_change(g, alpha)
ig = g[i]
nn = vertex[i].index_neighbor[ig]
if nn >= 0: # if ig is already neighbor of i
jginn = jg[vertex[i].perm[nn]]
if g != jginn:
# cycle consisting of two edges;
# replace jginn by g and insert h = g**-1*jginn
g1 = perm_af_invert(g)
h = perm_af_mul(g1, jginn)
jg[ vertex[i].perm[nn] ] = g
self.insert(h, alpha)
else: # new edge
self.insert_edge(g, i, ig)
self.cycle = [i]
if self.find_cycle(i, i, ig, -1):
cycle = self.cycle
cycle.append(cycle[0])
# find the smallest point (vertex) of the cycle
minn = min(cycle)
cmin = cycle.index(minn)
# now walk around the cycle starting from the smallest
# point, and multiply around the cycle to obtain h
# satisfying h[cmin] = cmin
ap = []
for c in range(cmin, len(cycle)-1) + range(cmin):
i = cycle[c]
j = cycle[c+1]
nn = vertex[i].index_neighbor[j]
p = jg[ vertex[i].perm[nn] ]
if i > j:
p = perm_af_invert(p)
ap.append(p)
ap.reverse()
h = perm_af_muln(*ap)
self.remove_edge(cycle[cmin], cycle[cmin + 1])
self.insert(h, alpha)
def coset_decomposition(self, g):
"""
Decompose `g` as h_0*...*h_{len(u)}
The Schreier-Sims coset representation u of `G`
gives a univoque decomposition of an element `g`
as h_0*...*h_{len(u)}, where h_i belongs to u[i]
Output: [h_0, .., h_{len(u)}] if `g` belongs to `G`
False otherwise
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([[0, 1, 3, 7, 6, 4], [2, 5]])
>>> b = Permutation([[0, 1, 3, 2], [4, 5, 7, 6]])
>>> G = PermutationGroup([a, b])
>>> c = Permutation([[0, 1, 2, 3, 4], [5, 6, 7]])
>>> G.coset_decomposition(c)
False
>>> c = Permutation([[0, 6], [1, 7], [2, 4], [3, 5]])
>>> G.coset_decomposition(c)
[[6, 4, 2, 0, 7, 5, 3, 1], [0, 4, 1, 5, 2, 6, 3, 7], [0, 1, 2, 3, 4, 5, 6, 7]]
>>> G.has_element(c)
True
"""
u = self.coset_repr()
if isinstance(g, Permutation):
g = g.array_form
g1 = g
n = len(u)
a = []
for i in range(n):
x = g1[i]
for h in u[i]:
if h[i] == x:
a.append(h)
p2 = perm_af_invert(h)
g1 = perm_af_mul(p2, g1)
break
else:
return False
if perm_af_muln(*a) == g:
return a
return False
def insert(self, g, alpha):
"""
insert permutation `g` in stabilizer chain at point alpha
"""
n = len(g)
if not g == self.idn:
vertex = self.vertex
jg = self.jg
i = _smallest_change(g, alpha)
ig = g[i]
nn = vertex[i].index_neighbor[ig]
if nn >= 0: # if ig is already neighbor of i
jginn = jg[vertex[i].perm[nn]]
if g != jginn:
# cycle consisting of two edges;
# replace jginn by g and insert h = g**-1*jginn
g1 = perm_af_invert(g)
h = perm_af_mul(g1, jginn)
jg[vertex[i].perm[nn]] = g
self.insert(h, alpha)
else: # new edge
self.insert_edge(g, i, ig)
self.cycle = [i]
if self.find_cycle(i, i, ig, -1):
cycle = self.cycle
cycle.append(cycle[0])
# find the smallest point (vertex) of the cycle
minn = min(cycle)
cmin = cycle.index(minn)
# now walk around the cycle starting from the smallest
# point, and multiply around the cycle to obtain h
# satisfying h[cmin] = cmin
ap = []
for c in range(cmin, len(cycle) - 1) + range(cmin):
i = cycle[c]
j = cycle[c + 1]
nn = vertex[i].index_neighbor[j]
p = jg[vertex[i].perm[nn]]
if i > j:
p = perm_af_invert(p)
ap.append(p)
ap.reverse()
h = perm_af_muln(*ap)
self.remove_edge(cycle[cmin], cycle[cmin + 1])
self.insert(h, alpha)
def test_coset_decomposition():
a = Permutation([2,0,1,3,4,5])
b = Permutation([2,1,3,4,5,0])
g = PermutationGroup([a, b])
assert g.order() == 360
rep = g.coset_repr()
d = Permutation([1,0,2,3,4,5])
assert not g.coset_decomposition(d.array_form)
assert not g.has_element(d)
c = Permutation([1,0,2,3,5,4])
v = g.coset_decomposition(c)
assert perm_af_muln(*v) == [1,0,2,3,5,4]
assert g.has_element(c)
a = Permutation([0,2,1])
g = PermutationGroup([a])
c = Permutation([2,1,0])
assert not g.coset_decomposition(c)
assert g.coset_rank(c) == None
def test_coset_decomposition():
a = Permutation([2,0,1,3,4,5])
b = Permutation([2,1,3,4,5,0])
g = PermutationGroup([a, b])
assert g.order() == 360
rep = g.coset_repr()
d = Permutation([1,0,2,3,4,5])
assert not g.coset_decomposition(d.array_form)
assert not g.has_element(d)
c = Permutation([1,0,2,3,5,4])
v = g.coset_decomposition(c)
assert perm_af_muln(*v) == [1,0,2,3,5,4]
assert g.has_element(c)
a = Permutation([0,2,1])
g = PermutationGroup([a])
c = Permutation([2,1,0])
assert not g.coset_decomposition(c)
assert g.coset_rank(c) == None
def normal_closure(self, gens):
"""
normal closure in self of a list gens2 of permutations
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([1, 2, 0])
>>> b = Permutation([1, 0, 2])
>>> G = PermutationGroup([a, b])
>>> G.order()
6
>>> G1 = G.normal_closure([a])
>>> list(G1.generate(af=True))
[[0, 1, 2], [1, 2, 0], [2, 0, 1]]
"""
G2 = PermutationGroup(gens)
if G2.is_normal(self):
return G2
gens1 = [p.array_form for p in self.generators]
b = 0
while not b:
gens2 = [p.array_form for p in G2.generators]
b = 1
for g1 in gens1:
if not b:
break
for g2 in gens2:
p = perm_af_muln(g1, g2, perm_af_invert(g1))
p = Permutation(p)
if not G2.has_element(p):
gens2 = G2.generators + [p]
G2 = PermutationGroup(gens2)
b = 0
break
return G2
def normal_closure(self, gens):
"""
normal closure in self of a list gens2 of permutations
Examples
========
>>> from sympy.combinatorics.permutations import Permutation
>>> from sympy.combinatorics.perm_groups import PermutationGroup
>>> a = Permutation([1, 2, 0])
>>> b = Permutation([1, 0, 2])
>>> G = PermutationGroup([a, b])
>>> G.order()
6
>>> G1 = G.normal_closure([a])
>>> list(G1.generate(af=True))
[[0, 1, 2], [1, 2, 0], [2, 0, 1]]
"""
G2 = PermutationGroup(gens)
if G2.is_normal(self):
return G2
gens1 = [p.array_form for p in self.generators]
b = 0
while not b:
gens2 = [p.array_form for p in G2.generators]
b = 1
for g1 in gens1:
if not b:
break
for g2 in gens2:
p = perm_af_muln(g1, g2, perm_af_invert(g1))
p = Permutation(p)
if not G2.has_element(p):
gens2 = G2.generators + [p]
G2 = PermutationGroup(gens2)
b = 0
break
return G2
