order.m

来自「Generates Alist format files for SFT LDP」· M 代码 · 共 23 行

function ord = order(x,m)
%Note: not sure how well this program works
over = 0;
pow = 1;

if x < m
    for i=1:m
        %pow = mod(x^i,m);
        pow = mod(pow*x,m);
        if pow == 1
            break;
        end
    end
    if pow == 1
        ord = i;
    else
        ord = -1;
    end
else
    ord = -1;
    
end