| version 1.2, 2000/08/21 08:31:41 |
version 1.4, 2001/03/08 07:49:13 |
|
|
| * shall be made on your publication or presentation in any form of the |
* shall be made on your publication or presentation in any form of the |
| * results obtained by use of the SOFTWARE. |
* results obtained by use of the SOFTWARE. |
| * (4) In the event that you modify the SOFTWARE, you shall notify FLL by |
* (4) In the event that you modify the SOFTWARE, you shall notify FLL by |
| * e-mail at risa-admin@flab.fujitsu.co.jp of the detailed specification |
* e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification |
| * for such modification or the source code of the modified part of the |
* for such modification or the source code of the modified part of the |
| * SOFTWARE. |
* SOFTWARE. |
| * |
* |
|
|
| * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
| * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
| * |
* |
| * $OpenXM: OpenXM_contrib2/asir2000/lib/fff,v 1.1.1.1 1999/12/03 07:39:11 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/lib/fff,v 1.3 2000/08/22 05:04:21 noro Exp $ |
| */ |
*/ |
| /* |
/* |
| fff : Univariate factorizer over a finite field. |
fff : Univariate factorizer over a finite field. |
|
|
| F2 = ugcd(F,F1); |
F2 = ugcd(F,F1); |
| /* FLAT = H/gcd(H,H') : square free part of H */ |
/* FLAT = H/gcd(H,H') : square free part of H */ |
| FLAT = sdiv(F,F2); |
FLAT = sdiv(F,F2); |
| |
FLAT /= LCOEF(FLAT); |
| I = 0; |
I = 0; |
| /* square free factorization of H */ |
/* square free factorization of H */ |
| while ( deg(FLAT,V) ) { |
while ( deg(FLAT,V) ) { |
|
|
| FLAT1 = simp_ff(1); |
FLAT1 = simp_ff(1); |
| else |
else |
| FLAT1 = ugcd(F,FLAT); |
FLAT1 = ugcd(F,FLAT); |
| |
FLAT1 /= LCOEF(FLAT1); |
| G = sdiv(FLAT,FLAT1); |
G = sdiv(FLAT,FLAT1); |
| FLAT = FLAT1; |
FLAT = FLAT1; |
| L = cons([G,I],L); |
L = cons([G,I],L); |