Наименьшее положительное целое число, взаимно простое с двумя последними предшественниками и еще не появившееся; a (1) = 1, a (2) = 2

10

Определение

  • Два целых числа взаимно просты, если они не имеют общих положительных делителей, кроме 1.
  • a(1) = 1
  • a(2) = 2
  • a(n)наименьшее целое положительное число , которое копростое к a(n-1)и a(n-2)и еще не появилось, для целого числа n >= 3.

задача

  • Учитывая положительное целое число n, вывод / печать a(n).

пример

  • a(11) = 6потому что 6взаимно прост с последними двумя предшественниками (а именно, 11и 13) и 6не появился раньше.

Ноты

  • Обратите внимание, что последовательность не является возрастающей, что означает, что элемент может быть меньше, чем его предшественник.

Спекуляции

  • Вы должны использовать 1-индексированный.

Testcases

n      a(n)
1      1
2      2
3      3
4      5
5      4
6      7
7      9
8      8
9      11
10     13
11     6
12     17
13     19
14     10
15     21
16     23
17     16
18     15
19     29
20     14
100    139
1000   1355
10000  13387
100000 133361

счет

  • Поскольку взаимно простое означает, что эти два числа имеют только один делитель ( 1) и 1является небольшим числом, ваш код должен быть как можно меньше с точки зрения количества байтов.

Ссылки

Дрянная Монахиня
источник
4
Эти "причины" для короткого кода ...
Луис Мендо
1
Интересно, почему это было понижено. Конечно, не из-за ужасного обоснования?
Конор О'Брайен
@ Конор Не я. На самом деле я проголосовал. Я надеюсь, что люди увидят и обоснование, и мой комментарий как шутки
Луис Мендо
3
Проблема с этими «забавными» обоснованиями кодового гольфа состоит в том, что мне нужно прочитать плохую шутку, охватывающую четыре строки, чтобы выяснить, что это стандартный кодовый гольф. Это просто затеняет правила вызова без уважительной причины.
Мартин Эндер
1
@ ConorO'Brien Не все браузеры всегда показывают заголовок (а также мобильное приложение), и мы обычно описываем выигрыш в посте в дополнение к использованию тега, потому что сам тег ничего не значит для новичков на сайт. Несмотря на то, что я знаком с нашими тегами типа задачи, я никогда не читаю их, чтобы выяснить, как оценивается задача, но пытаюсь найти ее в теле задачи. Тэг предназначен для категоризации, поиска и конкретной информации о типе задачи в вики-теге.
Мартин Эндер

Ответы:

5

Python 3.5, 160 141 126 124 121 109 байт

Это простая реализация определения последовательности. Предложения по игре в гольф приветствуются.

Изменить: -17 байт благодаря Leaky Nun. -9 байт благодаря Питеру Тейлору. -6 байт благодаря Sp3000 и переходу на Python 3.5.

import math;f=lambda n,r=[2,1],c=3:n<2and r[1]or(c in r)+math.gcd(c,r[0]*r[1])<2and f(n-1,[c]+r)or f(n,r,c+1)

Ungolfing:

import math
def f(n, r=[2,1], c=3):
    if n<2:
        return r[1]
    elif (c in r) + math.gcd(c,r[0]*r[1]) < 2:
        return f(n-1, [c]+r)
    else:
        return f(n, r, c+1)
Sherlock9
источник
Для Python 3.5+ import mathон g=math.gcdдолжен быть короче, чем определение вашего собственного g. Для до 3.5 вы можете сделать from fractions import*для gcd.
Sp3000
Если вы инициализируете c=3внутри цикла, вам нужно сделать это только один раз. По моим подсчетам вы экономите 3 символа.
Питер Тейлор
Существует также экономия в 2 символа от построения массива наоборот: вы должны использовать r=[c]+rвместо +=, но три отрицательных индекса становятся положительными. И еще есть 2-символьная экономия от перезаписи как лямбды, хотя это довольно радикальное изменение: from fractions import*;F=lambda n,r=[2,1],c=3:n<2and r[1]or(c in r)+gcd(r[0]*r[1],c)<2and F(n-1,[c]+r)or F(n,r,c+1)и не нужно, printпотому что это больше не полная программа.
Питер Тейлор
2

MATL , 28 27 байт

2:i:"`@ym1MTF_)Zdqa+}@h]]G)

Код медленный, но дает правильный результат.

Попробуйте онлайн! Или проверьте первые десять случаев .

Небольшая модификация кода создает график последовательности:

2:i:"`@ym1MTF_)Zdqa+}@h]]G:)XG

Рассматривайте это как искусство ASCII или с графическим выводом в автономном компиляторе:

введите описание изображения здесь

введите описание изображения здесь

объяснение

2:         % Push [1 2] to initiallize the sequence
i:         % Input n. Push [1 2 ... n]
"          % For loop: repeat n times
  `        %   Do while loop
    @      %     Push iteration index, starting at 1. This is the candidate number
           %     to extend the sequence
    y      %     Duplicate vector containing the sequence so far
    m      %     Is member? Gives true if the candidate is in the sequence
    1M     %     Push candidate and vector again
    TF_)   %     Get last two elements of the vector
    Zd     %     GCD between the candidate and those two elements. Produces a
           %     two-element vector
    qa     %     True if any of the two results exceeds 1, meaning
           %     the candidate is not coprime with the latest two sequence values
    +      %     Add. This corresponds to logical "or" of the two conditions, namely
           %     whether the candidate is member of the sequence so far, and
           %     whether it is not coprime with the latest two. In either case
           %     the do...while must continue with a next iteration, to try a new
           %     candidate. Else the loop is exited, and the current candidate
           %     is the new value of the sequence
  }        %   Finally (execute when the loop is exited)
    @h     %     Push current candidate and concatenate to the sequence vector
  ]        %   End do...while
]          % End for
G)         % Get n-th value of the sequence. Implicitly display
Луис Мендо
источник
1

C 185 байтов

G(a,b){return a%b?G(b,a%b):b;}
i,j,k;f(n){int a[n+2];for(i=0;i++<n;){a[i]=i<3?i:0;for(j=2;!a[i];++j){for(k=i;--k;){if(a[k]==j)++j,k=i;}a[G(a[i-1],j)*G(a[i-2],j)<2?i:0]=j;}}return a[n];}
orlp
источник
1

На самом деле , 38 37 35 33 31 30 байт

Это простая реализация определения функции. Предложения по игре в гольф приветствуются. Попробуйте онлайн!

Изменить: -3 байта благодаря Leaky Nun.

2R#╗,;`1";2±╜tπg@╜í+Y"£╓╖`nD╜E

Ungolfing:

2R#╗    Push [1,2] and store it in register 0
,;      Take input and duplicate
`1      Start function, push 1
  "       Start string
  ;       Duplicate i
  2±╜t    Push (list in register 0)[-2:]
  πg      gcd(i, product of list[-2:])
  @╜í     Rotate the gcd and bring up i, check for i in list (0-based, -1 if not found)
  +Y      Add the gcd and the index, negate (1 if coprime and not found in list, else 0)
  "£      End string, turn into a function
╓       Push first (1) values where f(x) is truthy, starting with f(0)
╖`      Append result to the list in register 0, end function
n       Run function (input) times
D╜E     Return (final list)[n-1]
Sherlock9
источник
1
Много манипуляций со стеком
Leaky Nun
0

Haskell, 81 73 байта

c l@(m:n:_)=m:c([x|x<-[1..],gcd(m*n)x<2,all(/=x)l]!!0:l)
((0:1:c[2,1])!!)

Пример использования: ((0:1:c[2,1])!!) 12-> 17.

Создайте список всего a(n), начиная с 0исправления индекса на основе 1 1и затем c[2,1]. cзанимает заголовок списка аргументов, lза которым следует рекурсивный вызов с последующим добавлением следующего номера (совмещенный, не замеченный ранее) l. Выберите nй элемент этого списка.

Ними
источник
0

R, 141 байт

 f=Vectorize(function(n)ifelse(n>3,{c=3;a=f(n-1);b=f(n-2);d=f(4:n-3);while(!c%%which(!a%%1:a)[-1]||!c%%which(!b%%1:b)[-1]||c%in%d)c=c+1;c},n))

ungolfed

f=Vectorize( function(n)     #build a recursive function. Vectorize allows
    if(n>3) {                #the function to be called on vectors.
        c=3                  #Tests size. Builds some frequent variables.
        a=f(n-1)
        b=f(n-2)
        d=f(4:n-3)           #Should really golf this out, but its horribly slow.
        while(!c%%which(!a%%1:a)[-1]||!c%%which(!b%%1:b)[-1]||c%in%d)
              c=c+1          #If we are coprime and not already seen. add.
        c
     } else n)               #First three are 1,2,3.
user5957401
источник
0

Mathematica, 97 90 байт

a@1=1;a@2=2;a@n_:=SelectFirst[Range[2n],GCD[a[n-1]a[n-2],#]<2&&!MemberQ[a/@Range[n-1],#]&]

Исходя из моей догадки, что a(n) < 2n для всех n.

Чтобы ускорить запуск, добавьте a@n=после оригинала, :=чтобы функции не нужно было пересчитывать предыдущие значения .

Сохранено 7 байт благодаря Sherlock9 (если gcd(a,b)=1потом gcd(ab,m) = gcd(a,m)*gcd(b,m))


источник
Это не гипотеза, так как на странице OEIS написано, что " ABS(a(n)-n) < n"
Дрянная Монахиня
@ LeakyNun Спасибо. Страница OEIS была недоступна несколько минут назад, и меня беспокоил возможный контрпример для большого числа n.
0

Pyth, 23 байта

eu+Gf&-TGq1iT*F>2G1tQ]1

Тестирование

Довольно простая реализация, но с некоторыми приятными уловками игры в гольф.

eu+Gf&-TGq1iT*F>2G1tQ]1
 u                 tQ]1    Apply the following function input - 1 times,
                           where G is the current state (List of values so far)
  +G                       Add to G
    f             1        The first number, counting up from 1
      -TG                  That has not been seen so far
     &                     And where
               >2G         The most recent two numbers
             *F            Multiplied together
           iT              Gcd with the current number being checked
         q1                Equals 1
e                          Output the final element of the list.
isaacg
источник