由于matlab中double类型的数据小数点后面老会自动加一些乱七八糟的。
所以在模逆过程中把double类型的数据先换成int类型,最后再换回来。
代码如下:
function [key,with] = modinv(F,m)
format long g
i=1;
a=int64(mod(det(F),m));
while(1)
if (mod((a*i),m)==1)
break;
else
i=i+1;
end
end
with=int64(mod(det(F)*inv(F),m));
key=mod(double(i*with+m*100),m);
end
参数F是输入矩阵,m是模数。模逆后的矩阵是key。
如果你能保证m是质数,就可以直接使用该方法。
否则,在while循环前,加入一个判断,判断a是否和m互质,代码网上都有,可以自己找了进行添加。