Визначення послідовності

Побудуйте послідовність натуральних чисел a(n)наступним чином:

1. a(0) = 4
2. Кожен доданок a(n), окрім першого, є найменшим числом, яке задовольняє наступному:
   а) a(n)складене число,
   б) a(n) > a(n-1),
   в) a(n) + a(k) + 1- складене число для кожного 0 <= k < n.

Отже, ми починаємо з a(0) = 4. Наступний запис a(1)повинен бути 9. Це не може бути 5або 7так як ті НЕ композит, і воно не може бути 6або 8тому , що 6+4+1=11це не є складовим і 8+4+1=13не складова. Нарешті, 9+4+1=14що є складовим, так a(1) = 9.

Наступний запис, a(2)має бути 10, так як це найменше число більше , ніж 9з 10+9+1=20і 10+4+1=15як композит.

Для наступного запису, 11і 13обидва, тому що вони не композит. 12є поза тим, 12+4+1=17що не є складовим. 14є поза тим, 14+4+1=19що не є складовим. Таким чином, 15це наступний член послідовності , тому що 15це композитний і 15+4+1=20, 15+9+1=25і 15+10+1=26все кожен композит, тому a(3) = 15.

Ось перші 30 термінів у цій послідовності:

4, 9, 10, 15, 16, 22, 28, 34, 35, 39, 40, 46, 52, 58, 64, 70, 75, 76, 82, 88, 94, 100, 106, 112, 118, 119, 124, 125, 130, 136

Це OEIS A133764 .

Виклик

З огляду на ціле вхідне число n, виведіть цей nтермін у цій послідовності.

Правила

• Ви можете вибрати індексацію на основі 0 або 1. Будь ласка, вкажіть, хто у ваших поданнях.
• Можна вважати, що вхід і вихід відповідають цілому цілому типу вашої мови.
• Введення та вихід можуть бути надані будь-яким зручним способом .
• Прийнятна або повна програма, або функція. Якщо функція, ви можете повернути вихід, а не надрукувати його.
• Стандартні лазівки заборонені.
• Це код-гольф, тому застосовуються всі звичайні правила гольфу, і найкоротший код (у байтах) виграє.

Лушпиння , 11 байт

!üȯṗ→+fotpN

1-індексований. Спробуйте в Інтернеті!

Пояснення

!üȯṗ→+fotpN  Implicit input, a number n.
          N  The list of positive integers [1,2,3,4,..
      f      Keep those
         p   whose list of prime factors
       ot    has a nonempty tail: [4,6,8,9,10,12,..
 ü           De-duplicate wrt this equality predicate:
     +       sum
    →        plus 1
  ȯṗ         is a prime number.
             Result is [4,9,10,15,16,..
!            Get n'th element.

Perl 6 , 70 байт

{(4,->+_{first {none($^a X+0,|(_ X+1)).is-prime},_.tail^..*}...*)[$_]}

Спробуйте його 0-індексувати

Розширено:

{  # bare block lambda with implicit parameter $_

  (  # generate the sequence

    4, # seed the sequence

    -> +_ { # pointy block that has a slurpy list parameter _ (all previous values)

      first

      {  # bare block with placeholder parameter $a

        none(                 # none junction
            $^a               # placeholder parameter for this inner block
          X+                
            0,                # make sure $a isn't prime
            |( _ X+ 1 )       # check all a(k)+1
        ).is-prime            # make sure none are prime
      },

      _.tail ^.. *            # start looking after the previous value
    }

    ...                       # keep generating values until

    *                         # never stop

  )[$_]                       # index into the sequence
}

Python 2 , 112 107 байт

дякую містеру Xcoder за байт.

n=-1,4;v=5
exec"while any(all((v-~k)%i for i in range(2,v))for k in n):v+=1\nn+=v,;v+=1\n"*input()
print~-v

Спробуйте в Інтернеті!

Python 2 , 115 109 байт

n=-1,4;v=4;x=input()
while x:v+=1;k=1^any(all((v-~k)%i for i in range(2,v))for k in n);n+=(v,)*k;x-=k
print v

Спробуйте в Інтернеті!

JavaScript (ES6), 83 байти

1-індексований

f=(n,a=[-1,p=4])=>a[n]||f(n,a.some(x=>(P=n=>n%--x?P(n):x<2)(x-=~p),p++)?a:[...a,p])

Демо

f=(n,a=[-1,p=4])=>a[n]||f(n,a.some(x=>(P=n=>n%--x?P(n):x<2)(x-=~p),p++)?a:[...a,p])

for(n = 1; n <= 30; n++) {
  console.log('a(' + n + ') = ' + f(n))
}

Прокоментував

Допоміжна функція P () , повертаючи істину, якщо n є простим, або помилковим в іншому випадку:

P = n => n % --x ? P(n) : x < 2

Примітка: його потрібно викликати з x = n .

Основна функція f () :

f = (               // given:
  n,                //   n = target index
  a = [-1, p = 4]   //   a = computed sequence with an extra -1 at the beginning
) =>                //   p = last appended value
  a[n] ||           // if a[n] exists, stop recursion and return it
  f(                // otherwise, do a recursive call to f() with:
    n,              //   n unchanged
    a.some(x =>     //   for each value x in a[]:
      P(x -= ~p),   //     rule c: check whether x + p + 1 is prime
                    //     rule a: because a[0] = -1, this will first compute P(p)
      p++           //     rule b: increment p before the some() loop starts
    ) ?             //   end of some(); if truthy:
      a             //     p is invalid: use a[] unchanged
    :               //   else:
      [...a, p]     //     p is valid: append it to a[]
  )                 // end of recursive call

05AB1E , 21 байт

0-індексований

®4Iµ)˜D¤N‹sN+>p_P*iN¼

Спробуйте в Інтернеті!

Мова Вольфрама (Mathematica) , 65 байт

±0=-1;±1=4;±n_:=±(n-1)+1//.a_/;PrimeQ/@Array[a+1+±#&,n,0,Or]:>a+1

Використовує кодування CP-1252 (Windows за замовчуванням). 1-індексований.

Спробуйте в Інтернеті!

Java 8, 186 173 байт

n->{int a[]=new int[n+1],r=a[n]=4;a:for(;n>0;)if(c(++r)<2){for(int x:a)if(x>0&c(r-~x)>1)continue a;a[--n]=r;}return r;}int c(int n){for(int i=2;i<n;n=n%i++<1?0:n);return n;}

0-індексований.
На жаль, прості чеки (або анти-прості / складені чеки в цьому випадку) не такі дешеві на Яві.

Пояснення:

Спробуйте в Інтернеті.

n->{                     // Method with integer as both parameter and return-type
  int a[]=new int[n+1],  //  Integer-array of size `n+1`
      r=a[n]=4;          //  Start the result and last item at 4
  a:for(;n>0;)           //  Loop as long as `n` is larger than 0
    if(c(++r)<2){        //   Raise `r` by 1, and if it's a composite:
      for(int x:a)       //    Inner loop over the array
        if(x>0           //     If the item in the array is filled in (non-zero),
           &c(r-~x)>1)   //     and if `r+x+1` is a prime (not a composite number):
          continue a;}   //      Continue the outer loop
      a[--n]=r;}         //    Decrease `n` by 1, and put `r` in the array
  return r;}             //  Return the result

// Separated method to check if a given number is a composite number
// (It's a composite number if 0 or 1 is returned, otherwise it's a prime.)
int c(int n){for(int i=2;i<n;n=n%i++<1?0:n);return n;}

Ruby + -rprime, 85 75 байт

->n{*a=x=4
n.times{x+=1;!x.prime?&&a.none?{|k|(x+k+1).prime?}?a<<x:redo}
x}

Спробуйте в Інтернеті!

Лямбда, що повертає n-індексований n-й елемент.

-10 байт: Використовуйте redoі потрійний оператор замість loop... breakі умовний ланцюг

Безголівки:

->n{
  *a=x=4                         # x is the most recent value: 4
                                 # a is the list of values so far: [4]
  n.times{                       # Repeat n times:
    x += 1                       # Increment x
    !x.prime? &&                 # If x is composite, and
      a.none?{|k|(x+k+1).prime?} #   for all k, a(n)+x+1 is composite,
      ? a<<x                     # Add x to a
      : redo                     # Else, restart the block (go to x+=1)
  }
  x                              # Return the most recent value
}

C (gcc) , 170 байт

P(n,d,b){for(b=d=n>1;++d<n;)b=b&&n%d;n=b;}h(n,N,b,k){if(!n)return 4;for(b=N=h(n-1);b;)for(b=k=!N++;k<n;b|=P(h(k++)-~N));n=N;}f(n,j){for(j=0;n--;)if(P(h(++j)))j++;n=h(j);}

Спробуйте в Інтернеті!

C (gcc) ,  140  138 байт

Дякуємо @Jonathan Frech за збереження двох байтів!

c(n,i){for(i=1;++i<n;)i=n%i?i:n;i=i>n;}f(n){int s[n],k,j,i=0;for(*s=k=4;i++-n;i[s]=k)for(j=!++k;j-i;)2-c(k)-c(k-~s[j++])?j=!++k:f;n=n[s];}

0-індексований

Спробуйте в Інтернеті!