布尔环上的分支Gr(o)bner基算法.pdfVIP

系统科学与数学 3，S譬8。Sci．Math．s残S． 29(9)(2009，9)，1266—1277 布尔环上的分支Gr6bner基算法 孙 瑶 王定康 (孛黧辩学院数学与鬃统辩掌辑究翳数学爨城亿重悫实验鬻，魏察100190) 铺黉 众掰霹熬Gr6bner基在徭多领域都有着十分重要的应雳。近些年寒Gr5bner蒸算法 有了很大的改进，其中最著名的是№ugbre提出的F4和F5算法．这两个算法具有很高的 效窀缝透鬻需要瀵耗大囊的凑襻。鉴予忿，将给礴一今帮零环上基予zdd数糖结梅静分支 GrSbner熬算法，该算法不仅町以大大降低对内存的消耗，还能有效的控制矩阵规模，从而 撼巍算法的整体效率．详缨阐述劳证明了算泼躬基本理论，分绍该分支算法的数据结构及分 支策略．最后通过实验数据可以发现，在很多例子中耽算法都要优于Magma中的F4算法． 关键词 分支Gr6bner基，带尔耀，zdd数据结构。 MR(2000)主戆分类罨12Y05 1 引 言 多项式系统求懈是计算机代数中的一个裁本问题，通过求解多项式系统我们可以解7央很 多领域申的实际闻题，比螺：力学，税械学，密码学等。在掰有求勰多项式系统的方法之中， GrSbner基和特征列方法都是我们熟知的最有效的方法．在特征列方两，高小山的布尔环上 的特征戮算法在褥淡流密鹦等方面取褥了疆好酶效果【卜鞭，露嚣雩也给我们的分支GrSbner基 算法以很大启发．相对于特征列算法，近些年来Gr6bner基算法有了很大的改进． GrSbner基算法的改进可以分为3拿方露，兵体来说魏是：数据结构，选择策略和减少 的终止饿弼并在1979年的文章中提高了三条改迸标准澜，其中阿条是减少冗余计算，另外 一条是选择策略。 雳合冲去除冗余计算的算法由Mora等人于1993年提出燃，从l毳；大大减少了不必要的计算， 继承了Larzard的线性囊二慧法，遴过成熟的短阵处理技术大大降低了约化多项式荔挚溃耗泌时 Mora等人提爨酶台冲标臻上的改进，霹操佟健更强，效率更高． 予F5算法对予多顼式的约纯顺序有特殊要求，菜矜程度上限铡了线健{弋数技术的使用，使 +隧家自然鼹学基食(NSFF02资助． 牧稿鞲期：2009-07-30． 9期 孙瑶等：布尔环上的分支Gr5bner基算法 1267 矩阵的处理效率受到一些影响，因而除了Faug色re的实现外，尚无法获得高效的F5算法．然 而从算法角度讲，F4算法固然拥有很高的效率，但其对内存的消耗却是惊人的．例如，用 算过程中所生成的多项式，因而将消耗掉大量的存储空间以及数据读取时间．此外，随着次 数上升，F4算法生成的矩阵规模将会迅速增大，从而导致算法整体的复杂度大幅度增加， 特别的，当基域的特征为0时，多项式的系数会疾速膨胀，这也使得矩阵变得更难处理． 鉴于此，我们的算法主要是在节省内存空间以及控制矩阵规模两方面进行改进．首先， 在算法的理论方面，我们修改并证明了F5算法所提出的两条准则，降低了算法对输入多项 式顺序的依赖；其次，针对布尔环的特性，我们采用了zdd数据结构，使得多项式可以以一 种更节省的方式进行存储，从而大大降低了算法对内存的消耗；最后，通过引入分支，我们 可以有效的控制矩阵规模，提高算法的整体效率． decision 算法所采用的zdd数据结构(zero-suppressedbinary 最先使用zdd数据结构来计算布尔环中的GrSbner基，并在纯字典序下取得了较好的结果 【15】，然而，该程序对于计算全次逆字典序的GrSbner基效果不甚理想．考虑到对于同～理 想，计算全次逆字典序的GrSbner基复杂度更低，因而，我们在实现算法时，对多项式的zdd 表达形式进行了适当修改，从而使得我们的算法更适于计算全次逆字典序的GrSbner基．此 外，前面提到的高小山提出的布尔环上的特征列算法也利用了zdd数据结构． 本文的主要内容如下：第2节介绍我们GrSbner基算法的基本理论；第3节介绍算法中 多项式的zdd表达形式；第4节介绍分支算法的实现以及分支算法的分支策略；第5节是有 关算法实现的一些细节问题；第6节给出一些实验数据；第7节总结