数的本原根

我已尝试实现中描述的算法 here 找到素数的本原根。

它适用于小素数,但当我尝试大数字时,它不再返回正确答案。

然后我注意到 a^(p-1)/pi 往往是一个很大的数字,它会返回 inf 在MATLAB中,所以我认为 (p-1) 可能会有帮助,但我不知道怎么做。

我用MATLAB写了一小段代码,就在这里。

clear all
clc

%works with prime =23,31,37,etc.

prime=761; %doesn't work for this value

F=factor(prime-1); % the factors of prime-1
 for i = 2: prime-1
        a=i;
        tag =1;
        for j= 1 :prime-1
            if (isprime(j))
                 p_i = j;
                 if(mod(a^((prime-1)/p_i),prime)== 1)
                     tag=0;
                     break 
                else
                    tag = tag +1;
                 end
            end
        end
        if (tag > 1 )
            a %it should only print the primitive root
            break
        end
end

欢迎任何意见。 谢谢