외로운 프라임 찾기


21

외로운 소수 (내가 부르는 것처럼)는 소수이며 너비가있는 숫자 그리드가 주어지면 w ≥ 3직교 또는 대각선으로 인접한 다른 소수가없는 소수입니다.

예를 들어,이 그리드를 여기서 w = 12(프라임은 굵은 체로 강조 표시됨)

1   2   3   4   5   6   7   8   9   10  11  12
13  14  15  16  17  18  19  20  21  22  23...
 ...86  87  88  89  90  91  92  93  94  95  96
97  98  99  100 101 102 103 104 105 106 107 108
109 110 111 112 113 114 115 116 117 118 119 120

두 소수 ( 103 , 107) 에만 직교 또는 대각선으로 어떤 소수도 소수가 없음을 알 수 있습니다 . 외로운 소수가 없기 때문에 섹션을 건너 뛰었습니다. (실제로 37을 제외하고)

당신의 임무는, 두 개의 입력이 주어 w ≥ 3지고 i ≥ 1, 폭이있는 숫자 그리드에서 첫 번째 외로운 소수를 결정하는 것입니다 w. 여기서 외로운 소수는보다 크거나 같아야합니다 i. 입력은 문자열로 취하는 것을 포함하여 합리적인 형식으로 가져올 수 있습니다. 너비에 대한 외로운 소수가 보장됩니다 w.

그리드가 줄 바꿈되지 않습니다.

예 :

w  i   output
11 5   11
12 104 107
12 157 157
9  1   151
12 12  37

이것이 이므로 가장 짧은 코드가 승리합니다!


이유입니다 w=12하지 37외로운 주요? 주변의 숫자 중 어느 것도 {25, 26, 38, 49, 50}소수입니다.
Jonathan Frech

@JonathanFrech 예, 테스트 사례에 포함됩니다.
Okx

답변:


8

C (GCC) , 159 (158) 149 바이트

P(n,d,b){for(d=b=1<n;n>++d;)b*=n%d>0;n=b;}F(w,i){w=P(i)&!(P(i-w)|P(i+w)|i%w>1&(P(~-i)|P(i+~w)|P(i+~-w))|i%w>0&(P(-~i)|P(-~i-w)|P(i-~w)))?i:F(w,++i);}

온라인으로 사용해보십시오!


줄 바꿈을 건너 뛰고 1 바이트를 저장할 수 있습니다. 온라인으로 사용해보십시오!
xanoetux

@ceilingcat 훌륭한 제안, 감사합니다.
Jonathan Frech

5

자바 스크립트 (ES6) 116 104 바이트

카레 구문으로 입력을 (w)(i)받습니다.

w=>g=i=>!(C=(k,n=d=i+k)=>n>0?n%--d?C(k,n):d>1:1)(0)&[i,x=1,i-1].every(j=>C(x-w)&C(w+x--)|j%w<1)?i:g(i+1)

테스트 사례

댓글

w =>                    // main function, taking w
  g = i =>              // g = recursive function, taking i
    !(                  //
      C = (             // define C:
        k,              //   a function taking an offset k
        n = d = i + k   //   and using n and d, initialized to i + k
      ) =>              //
        n > 0 ?         //   if n is strictly positive:
          n % --d ?     //     decrement d; if d does not divide n:
            C(k, n)     //       do a recursive call
          :             //     else:
            d > 1       //       return true if d > 1 (i.e. n is composite)
        :               //   else:
          1             //     return true (n is beyond the top of the grid)
    )(0) &              // !C(0) tests whether i is prime (or equal to 1, but this is safe)
    [                   // we now need to test the adjacent cells:
      i,                //   right side: i MOD w must not be equal to 0
      x = 1,            //   middle    : always tested (1 MOD w is never equal to 0)
      i - 1             //   left side : (i - 1) MOD w must not be equal to 0
    ]                   // for each value j defined above,
    .every(j =>         // and for x = 1, 0 and -1 respectively:
      C(x - w) &        //   test whether i - w + x is composite
      C(w + x--) |      //            and i + w + x is composite
      j % w < 1         //   or j MOD w equals 0, so that the above result is ignored
    ) ?                 // if all tests pass:
      i                 //   return i
    :                   // else:
      g(i + 1)          //   try again with i + 1

2

파이썬 2 , 144 바이트

f=lambda w,i,p=lambda n:all(n%j for j in range(2,n))*(n>1):i*(any(map(p,~-i%w*(i+~w,i-1,i+w-1)+(i-w,i+w)+i%w*(i-w+1,i+1,i-~w)))<p(i))or f(w,i+1)

온라인으로 사용해보십시오!

순서대로 인수 : w, i.

여기에는 외부 모듈이 사용되지 않습니다.

Python 2 + Sympy, 127 바이트

import sympy
f=lambda w,i,p=sympy.isprime:i*(any(map(p,~-i%w*(i+~w,i-1,i+w-1)+(i-w,i+w)+i%w*(i-w+1,i+1,i-~w)))<p(i))or f(w,i+1)

온라인으로 사용해보십시오!

여기서 유일한 차이점 sympy.isprime은 수동으로 구현 된 프라임 체크 기능 대신 사용한다는 점 입니다.


2

MATL , 38 바이트

xx`@1G*:5MeZpt3Y6Z+>3LZ)ft2G<~)X<a~}2M

온라인으로 사용해보십시오! 또는 모든 테스트 사례를 확인하십시오 .

설명

코드는 본질적으로 각 반복에서 한 행씩 챌린지에 설명 된대로 그리드를 확대하는 루프로 구성됩니다.

각 반복에서 그리드를 생성 한 후 마지막 행이 제거되고 (프라임이 외로운 지 여부를 알 수 없음) 나머지 숫자는 하나 이상의 외로운 소수가 있는지 테스트합니다. 이것은 2D 컨볼 루션을 통해 이루어집니다.

외로운 소수가 있으면 루프를 종료하고 첫 번째 소수를 출력합니다. 그렇지 않으면 다음 반복을 진행하여 더 큰 그리드를 시도합니다.

(코드는 실제로 행이 아닌 열로 확대 된 바뀐 버전의 격자를 사용합니다.)

xx        % Take two inputs (implicit): w, i. Delete them. They get copied
          % into clipboard G
`         % Do...while
  @       %   Push iteration index (1-based)
  1G      %   Push w
  *       %   Multiply
  :       %   Range from 1 to that
  5M      %   Push w again (from automatic clipboard M)
  e       %   Reshape into a matrix with w rows in column-major order
  Zp      %   Is prime? Element-wise
  t       %   Duplicate
  3Y6     %   Push neighbour mask: [1 1 1; 1 0 1; 1 1 1]
  Z+      %   2D convolution, maintaining size
  >       %   Greater than? Element-wise. Gives true for lonely primes
  3LZ)    %   Remove the last column
  f       %   Find linear indices of nonzeros
  t       %   Duplicate
  2G      %   Push i
  <~      %   Not less than?
  )       %   Use as logical index: this removes lonle primes less than i
  X<      %   Minimum. This gives either empty or a nonzero value
  a~      %   True if empty, false if nonzero. This is the loop condition.
          %   Thus the loop proceeds if no lonely prime was found
}         % Finally (execute on loop exit)
  2M      %   Push the first found lonely prime again
          % End (implicit). Display (implicit)

1

줄리아 0.6, 135 바이트

using Primes
f(w,i,p=isprime)=findfirst(j->(a=max(j-1,0);b=min(j+1,w);c=a:b;!any(p,v for v=[c;c+w;c-w]if v>0&&v!=j)&&p(j)&&j>=i),1:w*w)

TIO에는 Primes패키지 가 없습니다 . 내가 모든 외로운 소수를 (반환 할 수있어 경우 5 바이트 짧다 findfirst된다 find). Julia의 기능을 벗어나려는 Julia의 시도 Base는 골프를 아프게한다 (Julia의 목표가 아님) Primes는 0.4에 포함되었습니다.

ungolfed (주로)

function g(w,i)
    for j=i:w*w
        a,b=max(j-1,0),min(j+1,w)
        c=a:b
        !any(isprime,v for v=[c;c+w;c-w]if v>0&&v!=j)&&isprime(j)&&return j
    end
end

1

젤리 , 20 바이트

+‘ÆRœ^ḷ,ḷ’dạ/Ṁ€ṂḊð1#

온라인으로 사용해보십시오!

작동 원리

+‘ÆRœ^ḷ,ḷ’dạ/Ṁ€ṂḊð1#  Main link. Left argument: i. Right argument: w.

                 ð    Combine the links to the left into a chain and begin a new,
                      dyadic chain with arguments i and w.
                  1#  Call the chain to the left with left argument n = i, i+1, ...
                      and right argument w until 1 of them returns a truthy value.
                      Return the match.
+                       Yield n+w.
 ‘                      Increment, yielding n+w+1.
  ÆR                    Yield all primes in [1, ..., n+w+1].
      ḷ                 Left; yield n.
    œ^                  Multiset OR; if n belongs to the prime range, remove it; if
                        it does not, append it.
       ,ḷ               Wrap the resulting array and n into a pair.
         ’              Decrement all involved integers.
          d             Divmod; map each integer k to [k/w, k%w].
           ạ/           Reduce by absolute difference, subtracting [n/w, n%w] from
                        each [k/w, k%w] and taking absolute values.
             Ṁ€         Take the maximum of each resulting pair.
                        A maximum of 0 means that n is not prime.
                        A maximum of 1 means that n has a prime neighbor.
               Ṃ        Take the minimum of the maxima.
                Ḋ       Dequeue; map the minimum m to [2, ..., m].
                        This array is non-empty/truthy iff m > 1.


0

클린 , 181 ... 145 바이트

import StdEnv
@w i=hd[x+y\\y<-[0,w..],x<-[1..w]|x+y>=i&&[x+y]==[a+b\\a<-[y-w,y,y+w]|a>=0,b<-[x-1..x+1]|0<b&&b<w&&all((<)0o(rem)(a+b))[2..a+b-1]]]

온라인으로 사용해보십시오!

언 골프 드 :

@ w i
    = hd [
        x+y
        \\ y <- [0, w..]
        ,  x <- [1..w]
        | x+y >= i && [x+y] == [
            a+b
            \\ a <- [y-w, y, y+w]
            | a >= 0
            ,  b <- [x-1..x+1]
            | 0 < b && b < w && all ((<) 0 o (rem) (a+b)) [2..a+b-1]
            ]
        ]

0

젤리 ,  30  29 바이트

내 생각에 이것은 아마도 공정한 마진으로 이길 수 있습니다.

ÆPŒR+€×¥+©⁸’:⁹Ġ®ṁLÞṪFÆPS’¬ð1#

i왼쪽과 w오른쪽 에있는 2 차원 링크가 외로운 소수를 반환합니다.

온라인으로 사용해보십시오!

방법?

ÆPŒR+€×¥+©⁸’:⁹Ġ®ṁLÞṪFÆPS’¬ð1# - Link: i, w                     e.g. 37, 12
                           1# - find the 1st match starting at i and counting up of...
                          ð   - ...everything to the left as a dyadic link
                              - (n = i+0; i+1; ... on the left and w on the right):
ÆP                            -   is i prime: 1 if so, 0 if not     1
  ŒR                          -   absolute range: [-1,0,1] or [0]   [-1,0,1]
       ¥                      -   last two links as a dyad (w on the right):
      ×                       -     multiply (vectorises)           [-12,0,12]
    +€                        -     add for €ach       [[-13,-1,11],[-12,0,12],[-11,1,13]]
                              -     - i.e. the offsets if including wrapping
          ⁸                   -   chain's left argument, i
        +                     -   add                  [[24,36,48],[25,37,49],[26,38,50]]
                              -     - i.e. the adjacents if including wrapping
         ©                    -   copy to the register
           ’                  -   decrement            [[23,35,47],[24,36,48],[25,37,49]]
             ⁹                -   chain's right argument, w
            :                 -   integer division               [[1,2,3],[2,3,4],[2,3,4]]
              Ġ               -   group indices by value         [[1],[2,3]]
                              -     - for a prime at the right this would  be [[1,2],[3]]
                              -     - for a prime not at an edge it would be [[1,2,3]]
               ®              -   recall from register [[24,36,48],[25,37,49],[26,38,50]]
                ṁ             -   mould like           [[24,36,48],[[25,37,49],[26,38,50]]]
                  Þ           -   sort by:
                 L            -     length             [[24,36,48],[[25,37,49],[26,38,50]]]
                   Ṫ          -   tail                             [[25,37,49],[26,38,50]]
                              -     - i.e the adjacents now excluding wrapping
                    F         -   flatten                          [25,37,49,26,38,50]
                     ÆP       -   is prime? (vectorises)           [0,1,0,0,0,0]
                       S      -   sum                              1
                        ’     -   decrement                        0
                         ¬    -   not                              1            

내 생각 엔 이것이 아마도 아마도 공정한 마진으로 이길 수 있을까요? 골프 언어는 쉽지 않습니다.
Outgolfer Erik

아니,하지만 내의 추측 이 (도) 젤리에서 (그렇지 않으면 껍질!) 내 방법, 또는 더 나은 방법을 절약 할 수있는 몇 가지 방법이있을 수있다.
Jonathan Allan

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.