找质数(素数)

shiyiming

最近天气比较冷,人气也不旺,现在弄个小题目热热身:
寻找质数(素数)集合,一直到最小的两质数之差为62.

shiyiming

用Ｃ，Ｍａｔｌａｂ等都很好编，用ＳＡＳ不会，现在只会用ＳＡＳ读读数据整理整理再ｍｅａｎｓ，ｇｌｍ这种分析一下。。。期待高手解答～～

Qiong

不明白呀，跟集合有什么关系？
是由小到大的找质数，找到间隔大于等于62？

shiyiming

Does this work? <!-- s:-) --><img src="{SMILIES_PATH}/icon_smile.gif" alt=":-)" title="Smile" /><!-- s:-) -->

[code:2eemt1d1]data primes(keep = prime);
  retain difference testOn previousPrime;

  prime = 2;
  output;

  difference = 0;
  testOn = 3;

  do while (difference &lt;= 62);
    isPrime = 1;

    do divisor = 3 to int(sqrt(testOn)) by 2;
      if mod(testOn, divisor) = 0 then do;
        isPrime = 0;
        leave;
      end;
    end;
   
    if isPrime = 1 then do;
      prime = testOn;
      difference = prime - previousPrime;
      previousPrime = Prime;
      output;
    end;

    testOn = testOn + 2;
  end;
run;[/code:2eemt1d1]