Line data Source code
1 : /* Copyright (C) 2003 The PARI group.
2 :
3 : This file is part of the PARI/GP package.
4 :
5 : PARI/GP is free software; you can redistribute it and/or modify it under the
6 : terms of the GNU General Public License as published by the Free Software
7 : Foundation; either version 2 of the License, or (at your option) any later
8 : version. It is distributed in the hope that it will be useful, but WITHOUT
9 : ANY WARRANTY WHATSOEVER.
10 :
11 : Check the License for details. You should have received a copy of it, along
12 : with the package; see the file 'COPYING'. If not, write to the Free Software
13 : Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. */
14 :
15 : #include "pari.h"
16 : #include "paripriv.h"
17 :
18 : static long
19 1967 : groupelts_sumorders(GEN S)
20 : {
21 1967 : long i, s = 0;
22 43603 : for(i=1; i < lg(S); i++) s += perm_orderu(gel(S,i));
23 1967 : return s;
24 : }
25 :
26 : static long
27 700 : vecgroup_sumorders(GEN L)
28 : {
29 700 : long i, s = 0;
30 32501 : for(i=1; i<lg(L); i++) s += group_order(gel(L,i));
31 700 : return s;
32 : }
33 :
34 : static int
35 14 : indexgroupsubgroup(GEN L, long order, const long *good, const long *bad)
36 : {
37 : long i;
38 826 : for(i=1; i<lg(L); i++)
39 : {
40 812 : GEN G=gel(L,i);
41 : long idx;
42 : const long *p;
43 812 : if (group_order(G)!=order) continue;
44 210 : idx = group_ident(G,NULL);
45 420 : for(p=good; *p; p++)
46 210 : if (*p==idx) return 1;
47 210 : for(p=bad; *p; p++)
48 0 : if (*p==idx) return 0;
49 : }
50 14 : return 0;
51 : }
52 :
53 : static int
54 133 : indexgroupcentre(GEN G, GEN Z, const long *good, const long *bad)
55 : {
56 : long i;
57 357 : for(i=1;i<lg(Z);i++)
58 : {
59 357 : GEN z=gel(Z,i);
60 357 : if (perm_orderu(z)==2)
61 : {
62 182 : pari_sp btop=avma;
63 182 : GEN H = cyclicgroup(z,2);
64 182 : GEN C = group_quotient(G,H);
65 182 : GEN Q = quotient_group(C,G);
66 : const long *p;
67 182 : long idx=group_ident(Q,NULL);
68 182 : set_avma(btop);
69 413 : for(p=good;*p;p++)
70 252 : if (*p==idx) return 1;
71 259 : for(p=bad;*p;p++)
72 210 : if (*p==idx) return 0;
73 : }
74 : }
75 0 : return 0;
76 : }
77 :
78 : static GEN
79 91 : vecgroup_idxlist(GEN L, long order)
80 : {
81 91 : pari_sp av = avma;
82 91 : long i, j, l = lg(L);
83 91 : GEN V = cgetg(l, t_VECSMALL);
84 8001 : for (i = j = 1; i < l; i++)
85 7910 : if (group_order(gel(L,i)) == order) V[j++] = group_ident(gel(L,i),NULL);
86 91 : setlg(V, j); return gerepileuptoleaf(av, vecsmall_uniq(V));
87 : }
88 :
89 : /* Group identification code.
90 : * The group data are taken from GAP 4 and the SMALL GROUPS LIBRARY by Hans
91 : * Ulrich Besche, Bettina Eick and Eamonn O'Brien.
92 : */
93 :
94 : static long
95 1267 : group_ident_i(GEN G, GEN S)
96 : {
97 1267 : long n = S? lg(S)-1: group_order(G);
98 : long s;
99 : GEN F,p,e;
100 1267 : if (n==1) return 1;
101 1267 : if (!S) S = group_elts(G,group_domain(G));
102 1267 : s = groupelts_sumorders(S);/*This is used as a hash value*/
103 1267 : F = factoru(n);
104 1267 : p = gel(F,1);
105 1267 : e = gel(F,2);
106 1267 : switch(lg(p)-1)
107 : {
108 980 : case 1:/*prime power*/
109 980 : switch (e[1])
110 : {
111 98 : case 1: /* p */
112 98 : return 1;
113 35 : case 2: /* p^2 */
114 35 : return (s == 1 - p[1] + n*p[1])? 2: 1; /* pxp || p^2 */
115 231 : case 3: /* p^3 */
116 : {
117 231 : GEN H = group_abelianSNF(G, S);
118 231 : if (H) /*G is abelian*/
119 : {
120 161 : long l = lg(H)-1;
121 161 : return (l==3)?5: l; /*pxpxp || p^3 or p^2xp*/
122 : } /*G is not abelian*/
123 70 : if (p[1] == 2)
124 63 : return (s == 19)? 3: 4; /*D8 || Q8*/
125 : else
126 : {
127 7 : long q = p[1]*p[1];
128 7 : q *= q;
129 7 : return (s == q - p[1] + 1)?3 :4; /* pxp:p || p^2:p */
130 : }
131 : }
132 : }
133 616 : break;
134 273 : case 2:
135 273 : switch(e[1]+e[2])
136 : {
137 49 : case 2: /*pq*/
138 49 : return (p[2]%p[1]!=1)?1:1+group_isabelian(G); /*pq || pq or p:q*/
139 77 : case 3:
140 77 : if (p[1]==2 && e[1]==2) /* 4p */
141 : {
142 63 : long q = p[2], q2 = q*q, pmq2 = (q%4 == 1 || q==3);
143 63 : if (s==3+5*q+3*q2) return 1; /* 2p.2 */
144 63 : if (s==11-11*q+11*q2) return 2; /* 4p */
145 63 : if (s==3+q+3*q2) return 3+pmq2; /* D4p */
146 63 : if (s==7-7*q+7*q2) return 4+pmq2; /* 2px2 */
147 63 : return 3; /*A4 or p:4 */
148 : }
149 14 : else if (p[1]==2 && e[1]==1) /*2p^2*/
150 : {
151 7 : long q = p[2], q2 = q*q, q3 = q*q2, q4 = q*q3;
152 7 : if (s==1-q+3*q2-q3+q4) return 1; /* p^2:2 */
153 7 : if (s==3-3*q+3*q2-3*q3+3*q4) return 2; /* 2p^2 */
154 7 : if (s==1+q-2*q2+3*q3) return 3; /* D2pxp */
155 7 : if (s==1-q+2*q2+q3) return 4; /* p:2+p:2 */
156 7 : return 5; /* 2pxp */
157 : }
158 7 : else if (p[1]==3 && e[1]==2) /*9p, p>3*/
159 : {
160 7 : long q= p[2], q2 = q*q, p3 = (q%3 == 1), p9 = (q%9 == 1);
161 7 : if (s==7+47*q+7*q2) return 1; /* 3p.3 */
162 7 : if (s==61-61*q+61*q2) return 1+p3; /* 9p */
163 7 : if (s==1+59*q+q2) return 3; /* p:9 */
164 7 : if (s==7+11*q+7*q2) return 3+p9; /* p:3x3 */
165 7 : return 2+2*p3+p9; /* 3^2xp */
166 : }
167 0 : break;
168 : }
169 161 : case 3:
170 161 : switch(e[1]+e[2]+e[3])
171 : {
172 14 : case 3: /*pqr*/
173 14 : if (p[1]==2) /*2qr*/
174 : {
175 14 : long q = p[2],q2 = q*q, qc = 1-q+q2, qd = 1+q+q2;
176 14 : long r = p[3],r2 = r*r, rc = 1-r+r2, rd = 1+r+r2;
177 14 : long pmq = (r%q==1)? 2: 0;
178 14 : if (pmq && s==3*r*(q2-q)+rd) return 1; /* r:2q */
179 14 : if (pmq && s==3*(r2+(q2-q-1)*r+1)) return 2; /* r:qx2 */
180 7 : if (s==qd*rc) return 1+pmq; /* D2qxr */
181 7 : if (s==rd*qc) return 2+pmq; /* D2rxq */
182 7 : if (s==3*qc*rc) return 4+pmq; /* 2qr */
183 0 : return 3+pmq; /* q:2+r:2 */
184 : }
185 0 : break;
186 : }
187 763 : }
188 : {
189 763 : const long tab[]={
190 : 24, 173, 301, 99, 125, 113, 101, 97, 85, 161, 133, 189, 67, 87, 73, 105, -1,
191 : 36, 255, 671, 265, 219, 427, 243, 147, 275, 115, 127, 121, 159, 111, 175,
192 : -1,
193 : 40, 391, 903, 263, 311, 291, 271, 227, 207, 483, 399, 567, 163, 187, 315,
194 : -1,
195 : 56, 697, 1849, 585, 557, 529, 413, 385, 989, 817, 1161, 351, 357, 645, -1,
196 : 60, 945, 721, 561, 1617, 211, 497, 337, 373, 651, 581, 693, 501, 1029, -1,
197 : 75, 3647, 271, 847, -1,
198 : 84, 647, 935, 1935, 1295, 1071, 3311, 451, 699, 595, 1333, 469, 1099, 1419,
199 : 987, 2107, -1,
200 : 88, 1573, 4773, 1397, 1353, 1309, 953, 909, 2553, 2109, 2997, 865, 1665, -1,
201 : 90, 1659, 1891, 1371, 3843, 775, 1407, 735, 903, 615, 1575, -1,
202 : 104, 2143, 6751, 991, 1935, 1883, 1831, 1307, 1255, 3611, 2983, 4239, 731,
203 : 1203, 2355, -1,
204 : 105, 1785, 6321, -1,
205 : 126, 1533, 2037, 3397, 3477, 2749, 7869, 777, 937, 721, 1281, 1425, 2881,
206 : 1369, 1849, 1201, 3225, -1,
207 : 132, 4995, 2947, 2595, 8547, 1663, 3441, 2639, 3663, 2463, 5439, -1,
208 : 136, 3547, 11739, 1499, 3275, 3207, 3139, 2183, 2115, 6279, 5187, 7371,
209 : 987, 1159, 2047, 4095, -1,
210 : 140, 4429, 3885, 2989, 9933, 3053, 1613, 1697, 3297, 3569, 2849, 6321, -1,
211 : -1};
212 : long i;
213 : const long *t;
214 11277 : for(t=tab;*t!=-1;t++)
215 : {
216 10577 : if (t[0]==n)
217 : {
218 686 : for(i=1; t[i] != -1; i++)
219 686 : if (t[i]==s) return i;
220 0 : return -1;
221 : }
222 145138 : while (*t>=0) t++;
223 : }
224 : }
225 : {
226 700 : const long tab[]={
227 : 16, 1710031, 550075, 70099, 70075, 870059, 110059, 30079, 30071, 30063,
228 : 470131, 70155, 70107, 110115, 310307, -1,
229 : 32, 6830063, 150323, 1830171, /*230171*/0, 230227, 30283, 30251, 30203,
230 : /*70267*/0, 70219, 110227, /*230171*/0, /*70187*/0, /*70187*/0, 110155,
231 : 3430123, 430123, 30191, 30175, 30159, 1110387, 150563, 150387, 230323, 230411,
232 : 230299, 70619, 70467, 70355, /*70379*/0, /*70379*/0, /*70267*/0, 70291, 70555,
233 : 70331, 1750291, 230291, 430275, 70411, 70363, 70315, 110331, 30371, 30323,
234 : 950819, 150995, 150643, 230691, 30747, 30635, 632451, -1,
235 : 48, 430156, 11970124, 10219, 430324, 110324, 30396, 30444, 30348, 230300,
236 : 110300, 230396, /*70396*/0, /*70396*/0, 70540, 30412, 30364, 30380, 30332,
237 : 70492, 3850300, 490396, 490300, 6090236, 770236, 210316, 210284, 210252, 30299,
238 : 30371, 30363, 110299, 70311, 110271, 70588, 230732, 70876, 110636, 30868,
239 : 30628, 30596, 30644, 150684, 70828, 3290524, 490620, 490428, 770460, 30627,
240 : 70571, 10643, 151740, 2171228, -1,
241 : 54, 10256, 16410120, 70292, 610232, 10373, 10292, 10616, 70517, 5610228,
242 : 210309, 210228, 250448, 70616, 11696, 2370552, -1,
243 : /*64*/
244 : 72, 110307, 26230195, 210272, 30523, 110631, 30739, 70575, 30683,
245 : 14030351, 1830403, 1830299, 770346, 110586, 10750330, 10620, 210420, 71223,
246 : 9150663, 30426, 30858, 30978, 30954, 31074, 30762, 210452, 210538, 770634,
247 : 210730, 490626, 210722, 31018, 111234, 31450, 71106, 31322, 5750594, 750682,
248 : 750506, 10346, 10826, 10490, 70620, 11124, 10716, 30786, 31746, 210636, 491202,
249 : 72402, 3751122, -1,
250 : 80, 430250, 35910186, 110282, 430530, 110530, 30650, 30730, 30570,
251 : 230482, 110482, 230642, /*70642*/0, /*70642*/0, 70882, 30666, 30586, 30618,
252 : 30538, 70786, 11550450, 1470594, 1470450, 18270354, 2310354, 630474, 630426,
253 : 630378, 110562, 30562, 110722, 30722, 70546, 30546, 30946, 70962, 231202,
254 : 71442, 111042, 31426, 31026, 30978, 31058, 151106, 71346, 9870786, 1470930,
255 : 1470642, 2310690, 10467, 71266, 152866, 6511842, -1,
256 : 81,49210121, 6730319, 250427, 250319, 16450238, 610238, 70454, 70346,
257 : 70400, 70319, 5650670, 250994, 250670, 610589, 2412452, -1,
258 : /*96*/
259 : 100, 30393, 57310217, 10481, 30693, 36470341, 630408, 30968, 13310392,
260 : 210416, 10576, 11256, 10856, 11096, 630648, 31768, 8470616, -1,
261 : 108, 30552, 60170280, 610359, 30984, 38290440, 210660, 1830540, 30849,
262 : 30660, 31308, 211137, 20570532, 770721, 770532, 70769, 11636, 11813, 610647,
263 : 70647, 250674, 70674, 70890, 211092, 1830852, 31389, 31092, 32388, 211965,
264 : 13090836, 491133, 490836, 751080, 211416, 33576, 8691288, 70904, 11048, 71720,
265 : 13688, 12344, 251538, 751608, 212280, 36600, 5532024,-1,
266 : 112, 430344, 73530248, 430736, 110736, 30904, 31016, 30792, 230664,
267 : 110664, 230888, 0/*70888*/, 0/*70888*/, 71224, 30920, 30808, 30856, 30744,
268 : 71080, 23650600, 3010792, 3010600, 37410472, 4730472, 1290632, 1290568,
269 : 1290504, 71336, 231672, 72008, 111448, 31984, 31424, 31360, 31472, 151528,
270 : 71864, 20211048, 3011240, 3010856, 4730920, 30643, 153992, 13332456,-1,
271 : 120, 2310456, 770488, 110648, 63210360, 30763, 210552, 30712, 31256,
272 : 31240, 31336, 31400, 31496, 31480, 31096, 630498, 210808, 770968, 211128,
273 : 490904, 211064, 630744, 2310888, 631032, 1470840, 630984, 31128, 111368, 31608,
274 : 71224, 31464, 33810648, 4410744, 4410552, 11211, 31263, 11416, 210858, 11290,
275 : 11090, 211032, 31352, 32600, 630738, 491864, 1471704, 72664, 22051224, -1,
276 : 135, 114870240, 39270456, 1470618, 1470456, 16591104 , -1,
277 : -1};
278 : long i;
279 : const long *t;
280 700 : GEN Z = groupelts_center(S), L = group_subgroups(G);
281 700 : long scenter = groupelts_sumorders(Z), svecgroup = vecgroup_sumorders(L);
282 700 : long u = svecgroup+10000*scenter; /*This is used as a hash value*/
283 :
284 2268 : for(t=tab; *t!=-1; t++)
285 : {
286 2191 : if (t[0]==n)
287 : {
288 12894 : for(i=1; t[i] != -1; i++)
289 12971 : if (t[i]==u) return i;
290 105 : switch(n)
291 : {
292 98 : case 32:
293 : switch(u)
294 : {
295 42 : case 230171:
296 : {
297 42 : const long good[]={2,0}, bad[]={4,5,0};
298 42 : return indexgroupcentre(G,Z,good,bad)? 4: 12;
299 : }
300 7 : case 70267:
301 7 : if (s==135) return 9;
302 7 : return 32;
303 35 : case 70187:
304 : {
305 35 : const long good[]={8,0}, bad[]={7,9,0};
306 35 : return indexgroupcentre(G,Z,good,bad)? 13: 14;
307 : }
308 14 : case 70379:
309 : {
310 14 : const long good[]={4,0},bad[]={0};
311 14 : return indexgroupsubgroup(L,8,good,bad)? 31: 30;
312 : }
313 : }
314 0 : break;
315 7 : case 48: case 80:
316 : {
317 7 : const long good[]={5,8,0},bad[]={6,7,0};
318 7 : return indexgroupcentre(G,Z,good,bad)? 12: 13;
319 : }
320 0 : case 112:
321 : {
322 0 : const long good[]={7,4,0},bad[]={5,6,0};
323 0 : return indexgroupcentre(G,Z,good,bad)? 11: 12;
324 : }
325 : }
326 0 : return -1;
327 : }
328 47530 : while (*t!=-1) t++;
329 : }
330 : {
331 77 : const long tab[]={ 64, 1270001, /*4270003*/0, /*4270003*/0, 8190072,
332 : 6430073, 6670445, 5550446, 8190278, 7070279, 6350071, 5230072, 8110154,
333 : /*5870155*/0, /*5870155*/0, /*4270042*/0, /*4270042*/0, 7710246, 7390277,
334 : 6750037, 8032377, 8670266, 6750397, 11390022, 7710267, 7070277, /*3630046*/0,
335 : /*3630046*/0, 3630057, 4830196, 4830207, 4671808, 9070697, 6670700, 8750094,
336 : 6990091, 6350111, 5870115, 6671599, 5711601, 5551702, 5871512, 6351709,
337 : 5391711, /*3630046*/0, 3630047, 4111467, /*4430156*/0, /*4430156*/0, 3790166,
338 : /*2510026*/0, /*2510026*/0, 4470028, 4150300, 3830030, 13470021, 20350065,
339 : 10910041, 16514365, /*12190433*/0, 34110271, /*16514475*/0, 15230465,
340 : /*10910433*/0, 9630041, 13470233, /*16514475*/0, 20834696, /*10910433*/0,
341 : 13954343, /*12190433*/0, 19553542, 13474377, 25790466, 15870467, 18914476,
342 : 14110477, /*14590443*/0, 13310443, 11550043, /*14590443*/0, 10270043, 8990002,
343 : 8990546, 8990646, 8993647, 8356896, 13310905, 13310915, 12039018, 16990866,
344 : 12670888, 15071116, 11551217, 12038218, 16031739, 12512740, 12353138, 12993048,
345 : 15391849, 11872850, 12993551, 12353851, 8991446, 8991447, 8354830, 9951566,
346 : 9951666, 8674709, 9317039, 8031897, 8030187, 7713641, 7714641, 7074743,
347 : 9236585, 9236415, 9236586, 10990821, 9879066, 8751833, 9873399, 8751766,
348 : 10990754, 8593054, 8593087, 6830446, 6833546, 17472434, 13953445, 14432313,
349 : 16352544, 12833555, 13313423, 15635143, 13234877, 13874853, 12755287, 17870919,
350 : 14190949, 13075209, 11955310, 10835253, 9715354, 11312124, 10193135, 11074927,
351 : 12197529, 9957664, 11074970, 12196539, 9956674, 9958907, 10439497, 9479551,
352 : 9554015, 8433958, 9553915, 8434058, 8918081, 7798421, 10110856, 10110756,
353 : 9476648, 8991757, 8991857, 8356682, 10994275, 8750435, 8590474, 9230510,
354 : 10355254, 8115355, 13556790, 15790679, 11310691, 12275539, 14035061, 11795172,
355 : 8750465, 7474472, 8750475, 8114920, 6110196, 6111806, 5951808, 10191819,
356 : 9238364, 8271841, 8590736, 9390959, 8432079, 25470255, 41792701, 25470275,
357 : 20355344, 27233751, 18593673, 19717567, 23394762, 17312707, 19717633, 46115277,
358 : 31557088, 22917189, 24677288, 24039835, 24676366, 16032362, 17793529, 17153269,
359 : 38432486, 21153763, 23393635, 16037076, 27710971, 27074338, 20995174, 23396204,
360 : 20193482, 17157790, 19550231, 14751475, 17153832, 19070249, 16038080, 33391110,
361 : 25875097, 22197835, 22195018, 21070221, 24590112, 18999456, 15952565, 18356361,
362 : 17237769, 18359003, 15951169, 14832955, 16110295, 18350268, 21392354, 22030301,
363 : 18353365, 15955257, 13550032, 18430405, 18434015, 17150260, 17154128, 27234036,
364 : 23710639, 20194057, 21157900, 24198470, 20679613, 21158141, 22435065, 21318520,
365 : 20197076, 67390501, 83715011, 51070497, 54590283, 58915129, 50275230, 52035340,
366 : 263870051, -1,
367 : -1};
368 : long i;
369 : const long *t;
370 77 : GEN V=vecgroup_idxlist(L,32);
371 77 : long idxlist=vecsmall_pack(V,10,9967);
372 77 : long w=10000*svecgroup+idxlist; /*This is used as a hash value*/
373 84 : for(t=tab; *t!=-1; t++)
374 : {
375 77 : if (t[0]==n)
376 : {
377 18760 : for(i=1; t[i] != -1; i++)
378 18760 : if (t[i]==w) return i;
379 70 : switch(n)
380 : {
381 70 : case 64:
382 : switch(w)
383 : {
384 7 : case 4270003:
385 7 : return (scenter==439)? 2: 3;
386 7 : case 5870155:
387 : {
388 7 : const long good[]={8,9,0},bad[]={7,0};
389 7 : return indexgroupcentre(G,Z,good,bad)? 13: 14;
390 : }
391 7 : case 4270042:
392 : {
393 7 : const long good[]={13,0},bad[]={14,0};
394 7 : return indexgroupcentre(G,Z,good,bad)? 15: 16;
395 : }
396 7 : case 4430156:
397 : {
398 7 : const long good[]={18,20,0},bad[]={19,0};
399 7 : return indexgroupcentre(G,Z,good,bad)? 47: 48;
400 : }
401 7 : case 2510026:
402 7 : return scenter==1367? 50: 51;
403 7 : case 12190433:
404 7 : return scenter==47? 59: 70;
405 7 : case 16514475:
406 : {
407 7 : const long good[]={22,24,28,0},bad[]={23,25,27,30,0};
408 7 : return indexgroupcentre(G,Z,good,bad)? 61: 66;
409 : }
410 7 : case 10910433:
411 : {
412 7 : const long good[]={23,31,0},bad[]={25,26,29,30,33,0};
413 7 : return indexgroupcentre(G,Z,good,bad)? 63: 68;
414 : }
415 7 : case 14590443:
416 : {
417 7 : const long good[]={28,33,0},bad[]={30,34,0};
418 7 : return indexgroupcentre(G,Z,good,bad)? 77: 80;
419 : }
420 7 : case 3630046:
421 : {
422 7 : const long good[]={3,0},bad[]={12,16,0};
423 7 : if (scenter==695) return 26;
424 7 : return indexgroupcentre(G,Z,good,bad)? 27: 44;
425 : }
426 : }
427 0 : break;
428 : }
429 0 : return -1;
430 : }
431 1883 : while (*t!=-1) t++;
432 : }
433 : }
434 : {
435 7 : const long tab[]={ 96, 316010002, 252010002, 707020000, 676160124,
436 : 676170124, 87180027, 988190278, 892200028, 876030110, 876040110, 876120110,
437 : 215111237, 503062153, 972141052, 972131052, 455091156, 167101154, 620160033,
438 : 620170033, 908031033, 908121033, 908041033, 199101564, 7130153, 7140153,
439 : 812150123, 247051165, 487091566, 812151024, 391071276, 103081274, 247111377,
440 : 988180195, 892190205, 924190197, 828200207, 103050134, 679020464, 295091075,
441 : 199070145, 391060235, 103101076, 7080146, 135111157, 295020044, 684030223,
442 : 684040223, 908050274, 135060255, 7070285, 812080286, 71092475, 876102476,
443 : 908112287, 684120223, 748130243, 748140243, 620150254, 492160043, 492170043,
444 : 764180045, 700190477, 636200047, 963110003, 779050031, 935100032, 799110033,
445 : 819210003, 791260032, 246270050, 723330003, 987340003, 651360031, 623380033,
446 : 647263930, 839351534, 455321350, 178211335, 791244031, 322256848, 189340236,
447 : 130316409, 599331360, 743244548, 935295937, 551333907, 189222029, 274255883,
448 : 525275609, 82306043, 610289391, 82315641, 82307025, 647262487, 839353950,
449 : 0/*455322385*/, 0/*455322385*/, 178233588, 791245051, 322256982, 130307015,
450 : 658286619, 983297004, 983297037, 599333858, 631365165, 631376165, 535386399,
451 : 66408676, 354390966, 871428265, 775411064, 631376407, 535386309, 114430028,
452 : 823441008, 314398920, 74437993, 831420054, 42405827, 90439425, 799440830,
453 : 847426805, 767410275, 815440314, 863429143, 487360134, 487371044, 66211564,
454 : 66231664, 871295713, 66231764, 679242567, 125228724, 210253894, 18306803,
455 : 546288536, 162390938, 919437378, 871401018, 162255761, 967304398, 967313318,
456 : 413274329, 498283470, 498288163, 29345108, 967401139, 727449579, 679411219,
457 : 775352619, 583261276, 919295225, 66312839, 423381047, 2437470, 759424560,
458 : 711448550, 770222372, 109272382, 551210244, 258222592, 551230264, 295242430,
459 : 647254411, 199262266, 482272602, 871283751, 423293752, 519303751, 519312662,
460 : 71320222, 167332232, 226340245, 327352266, 167360274, 167372584, 103382587,
461 : 647392595, 455406162, 263412616, 327428742, 487438955, 295440098, 358290331,
462 : 622253358, 886280358, 322410312, 754400322, 394443122, 282440313, 354423123,
463 : 522430323, 726220349, 990273529, 470450359, 742460359, 522470032, 198470031,
464 : 282480353, 290490033, 274500353, 414470000, 614490000, 605473864, 664459790,
465 : 723464091, 893482714, 675465704, 845486215, 184493728, 653478045, 941489155,
466 : 605501588, 925482982, 264492577, 589502601, 312450472, 371466994, 285450492,
467 : 989464902, 578470486, 770489139, 994490497, 546500507, 604460529, 65270050,
468 : 684510049, 468510050, 134510562, 831510052, -1
469 : -1};
470 : long i;
471 : const long *t;
472 7 : GEN V=vecgroup_idxlist(L,48);
473 7 : long idx48=vecsmall_pack(V,10,9967);
474 7 : long idx32=vecgroup_idxlist(L,32)[1];
475 7 : long w=1000000*(svecgroup%997)+10000*idx32+idx48;
476 : /*This is used as a hash value*/
477 7 : for(t=tab; *t!=-1; t++)
478 : {
479 7 : if (t[0]==n)
480 : {
481 1127 : for(i=1; t[i] != -1; i++)
482 1127 : if (t[i]==w) return i;
483 0 : switch(n)
484 : {
485 0 : case 96:
486 : switch(w)
487 : {
488 0 : case 455322385:
489 : {
490 0 : const long good[]={37,40,0},bad[]={34,41,0};
491 0 : return indexgroupcentre(G,Z,good,bad)? 96: 97;
492 : }
493 : }
494 0 : break;
495 : }
496 0 : return -1;
497 : }
498 0 : while (*t!=-1) t++;
499 : }
500 : }
501 : }
502 0 : return 0;
503 : }
504 :
505 : long
506 1267 : group_ident(GEN G, GEN S)
507 : {
508 1267 : pari_sp av = avma;
509 1267 : long idx = group_ident_i(G, S);
510 1267 : if (idx < 0) pari_err_TYPE("group_ident [not a group]", G);
511 1267 : if (!idx) pari_err_IMPL("galoisidentify for groups of order > 127");
512 1267 : return gc_long(av, idx);
513 : }
514 :
515 : long
516 0 : group_ident_trans(GEN G, GEN S)
517 : {
518 0 : const long tab[]={
519 : 4, 1, 2, -1,
520 : 6, 2, 1, -1,
521 : 8, 1, 2, 4, 5, 3, -1,
522 : 9, 1, 2, -1,
523 : 10, 2, 1, -1,
524 : 12, 5, 1, 4, 3, 2, -1,
525 : 14, 2, 1, -1,
526 : 15, 1, -1,
527 : 16, 1, 4, 10, 8, 5, 6, 13, 12, 14, 2, 9, 7, 11, 3, -1,
528 : 18, 5, 1, 3, 4, 2, -1,
529 : 20, 2, 1, 5, 4, 3, -1,
530 : 21, 2, 1, -1,
531 : 22, 2, 1, -1,
532 : 24, 8, 1, 7, 5, 12, 13, 6, 14, 2, 15, 4, 10, 9, 11, 3, -1,
533 : 25, 1, 2, -1,
534 : 26, 2, 1, -1,
535 : 27, 1, 2, 3, 5, 4, -1,
536 : 28, 3, 1, 4, 2, -1,
537 : 30, 2, 4, 3, 1, -1,
538 : -1};
539 0 : long n = group_order(G), s;
540 : const long *t;
541 0 : if (n == 1) return 1;
542 : /* N.B. known up to 32 (Cannon-Holt) */
543 0 : if (n > 30) pari_err_IMPL("group_ident_trans [n > 30]");
544 0 : if (uisprime(n)) return 1;
545 0 : s = group_ident(G,S);
546 0 : for(t=tab;*t>=0;t++)
547 : {
548 0 : if (t[0]==n) return t[s];
549 0 : while (*t>=0) t++;
550 : }
551 : return 0; /*LCOV_EXCL_LINE*/
552 : }
|