精确大有限域线性代数库（例如GF（2 ^ 128）/ GF（2 ^ 256））

Question

常规

我正在寻找一个能够在大型有限域上进行精确计算的库，例如GF（2 ¹²⁸ ）/ _{2 ¹²⁸} 和GF（2 ²⁵⁶ ）/ _{2 ²⁵⁶} 。我列出了我需要的功能以及下面很酷的功能。显然，图书馆应尽可能快:-)。啊，既然我不是C ++大师（可能大部分库都是C ++），那么的示例代码会生成一个随机元素/一个常数并将其乘以它的乘法逆

必备功能

• 添加字段元素
• 字段元素的乘法
• 找到字段元素的乘法逆

很高兴有功能

• Vector / Matrix支持
• 随机元素支持

我已经看过的库可能不工作

• FFLAS/FFPACK，似乎不适用于如此大的有限域
• Givaro，似乎不适用于如此大的有限域

我已经看过可以工作的库（但我无法使用）

• NTL，我无法反转一个元素，但它应该真的有效，因为SAGE似乎在定义GF（2 ^ 256）时使用这个库，并且元素可以使用x^(-1)
进行反转
• PARI/GP，我无法在文档中找到我需要的所有内容，但是SAGE文档说它应该可以正常工作

其他说明

• 我正在编写一个Haskell程序，稍后将接口该库，因此更容易Haskell接口更好： - ）

Answer 1

NTL库似乎可以工作，使用这个（抱歉，我无法用C ++编程）代码

#include <NTL/GF2E.h>
#include <NTL/GF2EX.h>
#include <NTL/GF2X.h>
#include <NTL/GF2XFactoring.h>

NTL_CLIENT

int main()
{
    GF2X P = BuildIrred_GF2X(256);
    GF2E::init(P);

    GF2E zero = GF2E::zero();
    GF2E one;
    GF2E r = random_GF2E();
    GF2E r2 = random_GF2E();
    conv(one, 1L);
    cout << "Cardinality: " << GF2E::cardinality() << endl;
    cout << "ZERO: " << zero << " --> " << IsZero(zero) << endl;
    cout << "ONE:  " << one  << " --> " << IsOne(one)   << endl;
    cout << "1/r:  " << 1/r  << ", r * (1/r): " << (r * (1/r)) << endl;
    cout << "1/r2:  " << 1/r2  << ", r2 * (1/r2): " << (r2 * (1/r2)) << endl;
}

它似乎工作，证明（该程序的输出）：

Cardinality: 115792089237316195423570985008687907853269984665640564039457584007913129639936
ZERO: [] --> 1
ONE:  [1] --> 1
1/r:  [0 1 0 1 1 0 1 1 1 0 1 1 1 0 0 1 1 0 1 1 1 0 1 1 0 0 0 0 0 1 0 1 0 1 1 0 1 1 0 0 0 0 0 0 1 1 1 0 1 1 1 0 1 0 1 0 0 0 1 1 1 0 1 1 1 1 0 1 0 1 0 1 1 0 1 1 1 0 0 0 1 0 0 1 0 1 1 1 0 1 1 0 1 1 0 0 0 0 0 1 1 0 1 0 1 0 0 0 0 0 1 0 0 1 1 0 0 1 0 0 1 0 1 1 1 1 0 0 1 1 0 1 0 1 1 1 1 1 1 0 1 1 0 0 0 0 0 1 1 0 1 0 0 1 1 1 0 1 1 1 1 1 0 1 0 1 0 0 0 1 1 0 0 1 1 0 0 1 0 1 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 1 0 1 0 0 0 0 1 1 1 0 1 1 1 0 1 1 1 1 1 1 1 1 0 1 0 1 0 0 1 1 0 1 1 0 1 1 1 1 1 0 0 1 1 0 1 0 1 0 0 0 0 1 1 0 0 1 1 1 0 1], r * (1/r): [1]
1/r2:  [1 0 1 1 0 0 0 0 1 0 1 0 0 0 1 0 0 0 1 1 0 0 1 0 1 0 0 0 1 1 1 0 0 0 1 1 1 1 1 0 1 0 1 1 0 0 1 1 1 0 1 0 1 0 0 1 0 0 0 0 1 1 1 0 0 0 1 1 1 1 1 0 0 1 0 0 0 1 1 0 1 0 1 1 1 0 0 1 0 1 0 1 0 0 1 0 0 0 1 0 0 1 1 1 1 1 0 0 0 0 1 1 1 0 1 0 1 0 1 0 0 0 1 0 1 0 1 1 0 0 0 1 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 1 1 1 0 1 1 0 0 0 0 1 1 0 1 1 1 0 1 0 0 0 0 0 1 1 0 1 1 1 0 0 0 0 1 1 0 1 0 0 0 0 1 0 0 0 0 1 1 1 1 1 0 1 0 1 1 0 1 0 1 0 1 1 1 1 0 0 1 1 0 1 1 1 1 1 0 1 1 1 0 1 0 0 0 0 1 0 1 1 0 0 0 1 1 0 0 1 1 0 1 0 0 1 0 1 0 0 1 1], r2 * (1/r2): [1]

即使反转似乎也有效（在上面的输出示例中尽可能向右滚动）： - ）