1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
|
/*
Copyright (C) 2011 Sebastian Pancratz
Copyright (C) 2012 Lina Kulakova
This file is part of FLINT.
FLINT is free software: you can redistribute it and/or modify it under
the terms of the GNU Lesser General Public License (LGPL) as published
by the Free Software Foundation; either version 2.1 of the License, or
(at your option) any later version. See <http://www.gnu.org/licenses/>.
*/
#ifndef FMPZ_MOD_POLY_FACTOR_H
#define FMPZ_MOD_POLY_FACTOR_H
#ifdef FMPZ_MOD_POLY_FACTOR_INLINES_C
#define FMPZ_MOD_POLY_FACTOR_INLINE FLINT_DLL
#else
#define FMPZ_MOD_POLY_FACTOR_INLINE static __inline__
#endif
#undef ulong
#define ulong ulongxx /* interferes with system includes */
#include <stdio.h>
#undef ulong
#include <gmp.h>
#define ulong mp_limb_t
#include "flint.h"
#include "fmpz.h"
#ifdef __cplusplus
extern "C" {
#endif
/* Factoring ****************************************************************/
typedef struct
{
fmpz_mod_poly_struct *poly;
slong *exp;
slong num;
slong alloc;
} fmpz_mod_poly_factor_struct;
typedef fmpz_mod_poly_factor_struct fmpz_mod_poly_factor_t[1];
FLINT_DLL void fmpz_mod_poly_factor_init(fmpz_mod_poly_factor_t fac);
typedef struct
{
fmpz_mod_poly_struct * baby;
fmpz_mod_poly_struct * res;
fmpz_mod_poly_struct * H;
fmpz_mod_poly_struct * v;
fmpz_mod_poly_struct * vinv;
fmpz * tmp;
slong m;
}
fmpz_mod_poly_interval_poly_arg_t;
FLINT_DLL void fmpz_mod_poly_factor_clear(fmpz_mod_poly_factor_t fac);
FLINT_DLL void fmpz_mod_poly_factor_realloc(fmpz_mod_poly_factor_t fac, slong alloc);
FLINT_DLL void fmpz_mod_poly_factor_fit_length(fmpz_mod_poly_factor_t fac, slong len);
FLINT_DLL void fmpz_mod_poly_factor_set(fmpz_mod_poly_factor_t res, const fmpz_mod_poly_factor_t fac);
FMPZ_MOD_POLY_FACTOR_INLINE
void fmpz_mod_poly_factor_swap(fmpz_mod_poly_factor_t a, fmpz_mod_poly_factor_t b)
{
fmpz_mod_poly_factor_struct t = *a;
*a = *b;
*b = t;
}
FLINT_DLL void fmpz_mod_poly_factor_insert(fmpz_mod_poly_factor_t fac,
const fmpz_mod_poly_t poly, slong exp);
FLINT_DLL void fmpz_mod_poly_factor_print(const fmpz_mod_poly_factor_t fac);
FLINT_DLL void fmpz_mod_poly_factor_concat(fmpz_mod_poly_factor_t res,
const fmpz_mod_poly_factor_t fac);
FLINT_DLL void fmpz_mod_poly_factor_pow(fmpz_mod_poly_factor_t fac, slong exp);
FLINT_DLL int fmpz_mod_poly_is_irreducible(const fmpz_mod_poly_t f);
FLINT_DLL int fmpz_mod_poly_is_irreducible_ddf(const fmpz_mod_poly_t f);
FLINT_DLL int fmpz_mod_poly_is_irreducible_rabin(const fmpz_mod_poly_t f);
FLINT_DLL int fmpz_mod_poly_is_irreducible_rabin_f(fmpz_t fac,
const fmpz_mod_poly_t f);
FLINT_DLL int _fmpz_mod_poly_is_squarefree(const fmpz * f, slong len, const fmpz_t p);
FLINT_DLL int _fmpz_mod_poly_is_squarefree_f(fmpz_t fac,
const fmpz * f, slong len, const fmpz_t p);
FLINT_DLL int fmpz_mod_poly_is_squarefree(const fmpz_mod_poly_t f);
FLINT_DLL int fmpz_mod_poly_is_squarefree_f(fmpz_t fac, const fmpz_mod_poly_t f);
FLINT_DLL int fmpz_mod_poly_factor_equal_deg_prob(fmpz_mod_poly_t factor,
flint_rand_t state, const fmpz_mod_poly_t pol, slong d);
FLINT_DLL void fmpz_mod_poly_factor_equal_deg(fmpz_mod_poly_factor_t factors,
const fmpz_mod_poly_t pol, slong d);
FLINT_DLL void fmpz_mod_poly_factor_distinct_deg(fmpz_mod_poly_factor_t res,
const fmpz_mod_poly_t poly, slong * const *degs);
FLINT_DLL void fmpz_mod_poly_factor_squarefree(fmpz_mod_poly_factor_t res,
const fmpz_mod_poly_t f);
FLINT_DLL void fmpz_mod_poly_factor_distinct_deg_threaded(fmpz_mod_poly_factor_t res,
const fmpz_mod_poly_t poly, slong * const * degs);
FLINT_DLL void fmpz_mod_poly_factor_squarefree(fmpz_mod_poly_factor_t res,
const fmpz_mod_poly_t f);
FLINT_DLL void fmpz_mod_poly_factor(fmpz_mod_poly_factor_t res,
const fmpz_mod_poly_t f);
FLINT_DLL void fmpz_mod_poly_factor_cantor_zassenhaus(fmpz_mod_poly_factor_t res,
const fmpz_mod_poly_t f);
FLINT_DLL void fmpz_mod_poly_factor_kaltofen_shoup(fmpz_mod_poly_factor_t res,
const fmpz_mod_poly_t poly);
FLINT_DLL void fmpz_mod_poly_factor_berlekamp(fmpz_mod_poly_factor_t factors,
const fmpz_mod_poly_t f);
FLINT_DLL void _fmpz_mod_poly_interval_poly_worker(void * arg_ptr);
/* Roots *********************************************************************/
FLINT_DLL void fmpz_mod_poly_roots(fmpz_mod_poly_factor_t r,
const fmpz_mod_poly_t f, int with_multiplicity);
FLINT_DLL int fmpz_mod_poly_roots_factored(fmpz_mod_poly_factor_t r,
const fmpz_mod_poly_t f, int with_multiplicity, const fmpz_factor_t n);
/* Inlines *******************************************************************/
FLINT_DLL void fmpz_mod_poly_factor_get_fmpz_mod_poly(fmpz_mod_poly_t z, fmpz_mod_poly_factor_t fac, slong i);
#ifdef __cplusplus
}
#endif
#endif
|
