我正在使用NTL库对Galois Field进行试验。 GF2
是整数mod 2,GF2X
是GF2
的多项式,GF2E
是GF2
的环/场扩展。
我面临的问题是我按如下方式初始化不可约多项式
GF2X irreduc;
SetCoeff(irreduc, 128, 1);
SetCoeff(irreduc, 7, 1);
SetCoeff(irreduc, 2, 1);
SetCoeff(irreduc, 1, 1);
SetCoeff(irreduc, 0, 1);
GF2E::init(irreduc);
然后我还要初始化两个多项式:
GF2X a;
SetCoeff(a, 120);
SetCoeff(a, 22);
GF2X b;
SetCoeff(b, 128);
SetCoeff(b, 51);
std::cout << "a: " << a << '\n';
std::cout << "b: " << b << '\n';
并将它们相乘:
std::cout << "\ndeg(a * b): " << deg(a * b) << '\n';
输出为deg(a * b):
248 ,它超出了2^128
的域/环,由不可约多项式定义。
我知道我可能缺少明显的东西,但是我对这个领域很陌生,所以请耐心等待。
谢谢!
答案 0 :(得分:1)
正如您已经说过的那样,GF2X
代表GF2
上的多项式,因此它们不会因初始化GF2E的多项式而减少。您需要将多项式转换为GF2E
,然后一切正常。
因此将您的最后一行更改为
std::cout << "\ndeg(a * b): " << deg(conv<GF2X>(conv<GF2E>(a) * conv<GF2E>(b))) << '\n';
输出结果
deg(a * b): 124
这种转换很难看。我不确定是否有更好的方法以及NTL的文档记录方式,很难找到适合您要执行的功能。我只找到了GF2E::degree()
,但这只能为您提供多项式不可约的阶数。当您找到正确的方法时,请告诉我。