Skip to content

Instantly share code, notes, and snippets.

@mango314
Last active August 29, 2015 14:20
Show Gist options
  • Save mango314/d1d9c1b595ca7df231b9 to your computer and use it in GitHub Desktop.
Save mango314/d1d9c1b595ca7df231b9 to your computer and use it in GitHub Desktop.
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]]]
13
34
81
113
120
13
34
81
164
298
335
336
13
34
81
182
374
650
1044
1500
1824
2048
2156
2177
2184
13
34
81
182
374
758
1510
2606
3880
4594
4850
4878
4890
4895
4896
13
34
81
182
374
744
1412
2412
3606
4930
6212
6754
6834
6839
6840
13
34
81
182
374
744
1412
2534
4268
6638
9406
11280
12028
12126
12138
12143
12144
13
34
81
182
374
756
1508
2692
4612
7622
11854
16974
21538
23638
24080
24240
24308
24340
24356
24360
13
34
81
182
374
758
1526
3044
5952
10300
16502
23440
27740
29294
29666
29738
29754
29759
29760
13
34
81
182
374
758
1526
2890
5354
9586
15944
25092
34664
42348
47142
49574
50464
50594
50614
50616
13
34
81
182
374
758
1526
3062
6130
11588
20714
34424
50478
63046
68100
68812
68878
68880
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
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment