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);