8 static const char *prog
;
10 __attribute__((noreturn
))
11 static void barf(const char *msg
, ...)
16 fprintf(stderr
, "%s: ", prog
);
17 vfprintf(stderr
, msg
, ap
);
23 int main(int argc
, char *argv
[])
29 if (argc
!= 2) { fprintf(stderr
, "usage: %s N\n", prog
); exit(2); }
30 pari_init(16ul*1024ul*1024ul, 0);
31 n
= gp_read_str(argv
[1]);
32 if (typ(n
) != t_INT
|| signe(n
) <= 0) barf("expected a positive integer");
33 sd_factor_proven("1", d_SILENT
);
35 assert(typ(ff
) == t_MAT
);
37 pp
= gel(ff
, 1); ee
= gel(ff
, 2);
38 assert(typ(pp
) == t_COL
);
39 assert(typ(ee
) == t_COL
);
40 nf
= lg(pp
); assert(nf
== lg(ee
));
44 pari_printf("%Ps %Ps\n", *pp
, *ee
);