expansion in the group SL(2,Fp) with A = [[[1,1],[0,1]],[[1,-1],[0,1]],[[1,0],[1,1]],[[1,0],[-1,1]]]

05.txt

13
34
81
113
120

07.txt

13
34
81
164
298
335
336

13.txt

13
34
81
182
374
650
1044
1500
1824
2048
2156
2177
2184

17.txt

13
34
81
182
374
758
1510
2606
3880
4594
4850
4878
4890
4895
4896

19.txt

13
34
81
182
374
744
1412
2412
3606
4930
6212
6754
6834
6839
6840

23.txt

13
34
81
182
374
744
1412
2534
4268
6638
9406
11280
12028
12126
12138
12143
12144

29.txt

13
34
81
182
374
756
1508
2692
4612
7622
11854
16974
21538
23638
24080
24240
24308
24340
24356
24360

31.txt

13
34
81
182
374
758
1526
3044
5952
10300
16502
23440
27740
29294
29666
29738
29754
29759
29760

37.txt

13
34
81
182
374
758
1526
2890
5354
9586
15944
25092
34664
42348
47142
49574
50464
50594
50614
50616

41.txt

13
34
81
182
374
758
1526
3062
6130
11588
20714
34424
50478
63046
68100
68812
68878
68880

expand.py

A = [[[1,1],[0,1]],[[1,-1],[0,1]],[[1,0],[1,1]],[[1,0],[-1,1]]]
p = 41
B = A
for t in range(35):
    B = [dot(b,a) for b in B for a in A]
    B = [eval(y) for y in set([str(x) for x in B])]
    print len(B)
    if(len(B) == (p**2-1)*p):
        break