伽罗华域上(Galois Field)上的四则运算-笔记

原网页：https://blog.csdn.net/shelldon/article/details/54729687

---

单位元：$a\ast e=e\ast a=a$

本原多项式：域中不可约多项式是不能够进行因子分解的多项式。

1、多项式加减法：$(x^2+1)+(x+1)=x^2+x+1$
2、多项式乘法：$(x^2+1)\ast (x+1)=x^2+2x+1$
3、多项式除法：$\left(x^3-12x^2-42\right)÷\left(x-3\right)=x^2-9x-27\cdots \cdots -123$

---

$GF\left(2^w\right)$多项式计算，多项式系数只能取$0,1$；
$GF\left(3^w\right)$多项式计算，多项式系数只能取$0,1,2$；
在$GF\left(2^w\right)$多项式加法中，合并阶数相同的同类项时为异或运算，减法等同于加法。

---

$GF\left(p\right)$中的$p$一定要为素数，加法$\left(a+b\right)modp$，乘法$\left(a\ast b\right)modp$。(素数$=1\ast$本身，不可被$1$与本身之外的整数整除)
$p$为素数时能保证任意一个$a$在$GF\left(2^w\right)$中有$b$，使$a\ast b=e($不含$0)$。

---

$GF\left(2^w\right)$上的加法和乘法运算，不再是一般的加法和乘法，而是多项式运算。

---

伽罗华域的元素可以通过该域上的的本原多项式生成，通过本原多项式得到的域，其加法单位元都是0，乘法单位元是1。

伽罗华域上的运算：

加法和减法：异或运算

乘法和除法：$mod\left(p\right)$

当为0时，小于多项式，不需要取模运算
当为1时，一定小于等于，不需要取模运算

---

$GF\left(2^w\right)$是一个有限域，就是元素个数是有限的，但指数$k$是可以无穷的，所以必然存在循环。这个循环周期为$2^w-1$。

生成元：在域$GF\left(2^w\right)$中的生成元为2，当$k$大于等于$2^w-1$时，$g^k=g^{\left[k\%\left(2^w-1\right)\right]}$。

$g$为域$GF\left(2^w\right)$上的生成元，集合$\left\{g^0,g^1,...,g^{\left(2^w-2\right)}\right\}$包含了域$g$为域$GF\left(2^w\right)$上所有非零元素。

对于$g^k=a$，有正过程和逆过程。指数$k$求$a$是正过程，元素$a$求$k$是逆过程。

---

例：$GF\left(2^3\right)$的一个生成多项式为$x^3+x+1$

生成元对应的非零元素集合$\left\{2^0,2^1,2^2,2^3,2^4,2^5,2^6\right\}$

$2^0\to 1$

$2^1\to x$

$2^2\to x^2$

$2^3\to x+1$

$2^4\to x^2+x$

$2^5\to x^2+x+1$

$2^6\to x^2+1$

$nextloop\to$

$2^7\to 1$

$\cdots \cdots$

生成元$2^k$中$k$值	元素${}_{gf}$	对应二进制	十进制
-	0	000	0
0	1	001	1
1	$x$	010	2
2	$x^2$	100	4
3	$x+1$	011	3
4	$x^2+x$	110	6
5	$x^2+x+1$	111	7
6	$x^2+1$	101	5
7	$1$	001	1
…	…	…	…

clear all;
close all;
clc;

index = 4;
len = pow2(index)-1;

gfnums1 = gf(randi([0 len],[1,3]),index);
gfnums2 = gf(randi([0 len],[1,3]),index).';

nums1 = double(gfnums1.x);
nums2 = double(gfnums2.x);

mid = gf(zeros(size(nums1)),index);
fprintf('\t伽罗华域对位乘:\n');
for i=1:size(nums1,2)
   y = gfnums1(i)*gfnums2(i);
   mid(i) = y;
   disp(double(y.x));
end

fprintf('\t伽罗华域将对位乘相加得到伽罗华域结果:\n');
y = gf(0,index);
for i=1:size(nums1,2)
   y = y+mid(i);
   disp(double(y.x));
end

temp = gfnums1*gfnums2;
fprintf('\t伽罗华域的乘积:\n');
disp(double(temp.x));
fprintf('\t自然域的乘积:\n');
disp(nums1*nums2);
fprintf('\t自然域的乘积模:\n');
disp(mod(nums1*nums2,pow2(index)));

	伽罗华域对位乘:
     3

     7

     1

	伽罗华域将对位乘相加得到伽罗华域结果:
     3

     4

     5

	伽罗华域的乘积:
     5

	自然域的乘积:
   194

	自然域的乘积模:
     2