Question

我正在研究JavaScript以使三种不同的算法在同一代码中工作，我为每种算法设置了一个函数。我正在尝试获得三种Primality测试方法：试验部门，Fermat Primality测试和Miller-Rabin Primality测试。前两个工作正常，但Miller-Rabin不是。我对JavaScript和编程一般都很陌生，所以如果你能找到我出错的地方或想办法让它发挥作用，请告诉我！谢谢！

// 1531 6389 68819 688889 6388819
// 68883889 688838831 1000000009
// 561 is a Carmichael number; a Fermat pseudoprime with the property a^n-1 = 1 mod n, for any "a" coprime to 561.
input = 5491763;
numTrials = 2000;

document.getElementById("input").innerHTML = input;

function TrialDiv(n) {

  if (n === 1) {
    return false;
  } else if (n === 2) {
    return true;
  } else {
    for (var x = 2; x < n; x++) {
      if (n % x === 0) {
        return false;
      }
    }
    return true;
  }
}

if ((TrialDiv(input)) === true) {
  a = "Prime"
} else if ((TrialDiv(input)) === false) {
  a = "Composite"
}
//---------------------------------------------------------------------------
function gcd(x, y) {
  while (y !== 0) {
    var z = x % y;
    x = y;
    y = z;
  }
  return x;
}

function getRndInteger(max) {
  return Math.floor(Math.random() * (max - 2)) + 2;
}
//--------------------------------------------------------------------------
function Fermat(n) {

  for (var t = 0; t = numTrials; t++) {
    m = getRndInteger(input);
    if (gcd(m, n) !== 1) {
      return false;
    }
  }
  return (Math.pow(m, n - 1) % n !== 1);
}

if ((Fermat(input)) === true) {
  b = "Prime";
} else if ((Fermat(input)) === false) {
  b = "Composite";
}
//---------------------------------------------------------------------------
function genD(n) { // Generates "d" such that "n-1 = 2^s * d"
  var p = n - 1;
  var d = p / 2;
  while (d % 2 === 0) {
    d = d / 2;
  }
  return d;
}

function genS() { // Generates "s" such that "n-1 = 2^s * d"
  var s = Math.log2(p / d);
  return s;
}
//---------------------------------------------------------------------------
function MillerRabin(n) {

  for (var t = 0; t < numTrials; t++) {
    m = getRndInteger(input);
    if (gcd(m, n) !== 1) {
      return false;
    } else {
      for (var r = 0; r < genS(); r++) {
        power = (Math.pow(2, r) * genD(input));
        if (Math.pow(m, genD(input)) % n === 1 || Math.pow(m, power) % n === -1) {
          return true;
        } else {
          return false;
        }
      }
      return true;
    }
    return true;
  }
}

if ((MillerRabin(input)) === true) {
  c = "Prime";
} else if ((MillerRabin(input)) === false) {
  c = "Composite";
}

<body>
  <button type="button" onclick='document.getElementById("TrialDivision").innerHTML = a;  document.getElementById("FermatTest").innerHTML = b; document.getElementById("MillerRabinTest").innerHTML = c; '>Show</button>

  <hr>
  <b style="color:rgb(0,0,120)">PRIMALITY TESTS</b>
  <p></p>
  Input:
  <l id="input"></l>

  <hr>
  <h5 style="color:rgb(160,0,0)">TRIAL DIVISION</h5>
  <p></p>
  Output:
  <i id="TrialDivision"></i>

  <hr>
  <h5 style="color:rgb(160,0,0)">FERMAT PRIMALITY TEST</h5>
  <p></p>
  Output:
  <i id="FermatTest"></i>

  <hr>
  <h5 style="color:rgb(160,0,0)">MILLER-RABIN PRIMALITY TEST</h5>
  <p></p>
  Output:
  <i id="MillerRabinTest"></i>
</body>
<script>

这就是我写它的方式，这纯粹是由我从每个测试的原始数学算法创建的。当输入数为素数时，Miller-Rabin输出没有显示任何内容;算法无法识别它。但它确实能够正确识别复合材料。

请告诉我你想到的任何改进！

Answer 1

循环未运行时遇到一些问题。我不确定该算法应该如何工作，所以我不知道它是否正常工作，但是直接问题是由于泄漏变量和未在函数中声明的变量（这使得它们全局化）。

解决这个直接问题的方法是在函数中声明局部变量，这样你就不依赖于通过其他函数设置局部变量了。

function genD(n) { // Generates "d" such that "n-1 = 2^s * d"
    var p = n - 1;
    var d = p / 2;
    while (d % 2 === 0) {
        d = d / 2;
    }
    return d;
}

function genS(n) { // Generates "s" such that "n-1 = 2^s * d"
    var p = n - 1
    var d = p / 2
    var s = Math.log2(p / d);
    return s;
}

一旦你开始修复崩溃的错误，还有一些其他大问题。例如，您的fermat()可能会永远运行：

for (var t = 0; t = numTrials; t++) {

应该是：

for (var t = 0; t == numTrials; t++) {

每个循环

=将t设置为numTrails

我认为你应该从头开始，将所有功能放在一起，以及所有其他逻辑代码一起调用这些功能，这样你就可以看到发生了什么。

JavaScript上的Miller Rabin Primality测试错误

1 个答案: