Количество положительных рациональных

14

Можно показать, что положительные рациональные числа нумеруются с помощью следующего процесса:

  1. Ноль имеет порядковый номер 0
  2. Расположите другие числа в сетке так, чтобы строка a, столбец b содержала a / b
  3. Постройте диагональный зигзаг сверху справа внизу слева
  4. Ведите счет уникальных чисел, встречающихся вдоль зигзага

Вот картина зигзага:

Начните с 1/1, начальное движение вправо

Итак, числа встречаются в порядке

1/1, 2/1, 1/2, 1/3, 2/2, 3/1, 4/1, 3/2, 2/3, 1/4, 1/5, 2/4, 3/3, 4/2, 5/1, 6/1, 5/2, 4/3, 3/4, 2/5, 1/6, 1/7, 2/6, 3/5, 4/4, 5/3 ...

И упрощенные, уникальные числа встречаются

1, 2, 1/2, 1/3, 3, 4, 3/2, 2/3, 1/4, 1/5, 5, 6, 5/2, 4/3, 3/4, 2/5, 1/6, 1/7, 3/5, 5/3, ...

Вызов:

  • Если даны два целых числа больше нуля p и q , выведите порядковое число p / q
  • р и q не обязательно взаимно просты
  • Самый короткий код выигрывает
  • Стандартные лазейки запрещены

Тестовые случаи:

Вот первые 24 встреченных рациональных числа и желаемый результат для каждого:

1/1: 1
2/1: 2
1/2: 3
1/3: 4
2/2: 1
3/1: 5
4/1: 6
3/2: 7
2/3: 8
1/4: 9
1/5: 10
2/4: 3
3/3: 1
4/2: 2
5/1: 11
6/1: 12
5/2: 13
4/3: 14
3/4: 15
2/5: 16
1/6: 17
1/7: 18
2/6: 4
3/5: 19

И, для дальнейших тестов, вот 200 первых положительных рациональных чисел в порядке:

1, 2, 1/2, 1/3, 3, 4, 3/2, 2/3, 1/4, 1/5, 
5, 6, 5/2, 4/3, 3/4, 2/5, 1/6, 1/7, 3/5, 5/3, 
7, 8, 7/2, 5/4, 4/5, 2/7, 1/8, 1/9, 3/7, 7/3, 
9, 10, 9/2, 8/3, 7/4, 6/5, 5/6, 4/7, 3/8, 2/9, 
1/10, 1/11, 5/7, 7/5, 11, 12, 11/2, 10/3, 9/4, 8/5, 
7/6, 6/7, 5/8, 4/9, 3/10, 2/11, 1/12, 1/13, 3/11, 5/9, 
9/5, 11/3, 13, 14, 13/2, 11/4, 8/7, 7/8, 4/11, 2/13, 
1/14, 1/15, 3/13, 5/11, 7/9, 9/7, 11/5, 13/3, 15, 16, 
15/2, 14/3, 13/4, 12/5, 11/6, 10/7, 9/8, 8/9, 7/10, 6/11, 
5/12, 4/13, 3/14, 2/15, 1/16, 1/17, 5/13, 7/11, 11/7, 13/5, 
17, 18, 17/2, 16/3, 15/4, 14/5, 13/6, 12/7, 11/8, 10/9, 
9/10, 8/11, 7/12, 6/13, 5/14, 4/15, 3/16, 2/17, 1/18, 1/19, 
3/17, 7/13, 9/11, 11/9, 13/7, 17/3, 19, 20, 19/2, 17/4, 
16/5, 13/8, 11/10, 10/11, 8/13, 5/16, 4/17, 2/19, 1/20, 1/21, 
3/19, 5/17, 7/15, 9/13, 13/9, 15/7, 17/5, 19/3, 21, 22, 
21/2, 20/3, 19/4, 18/5, 17/6, 16/7, 15/8, 14/9, 13/10, 12/11, 
11/12, 10/13, 9/14, 8/15, 7/16, 6/17, 5/18, 4/19, 3/20, 2/21, 
1/22, 1/23, 5/19, 7/17, 11/13, 13/11, 17/7, 19/5, 23, 24, 
23/2, 22/3, 21/4, 19/6, 18/7, 17/8, 16/9, 14/11, 13/12, 12/13, 
11/14, 9/16, 8/17, 7/18, 6/19, 4/21, 3/22, 2/23, 1/24, 1/25 

Прикоснитесь к обратному вопросу , где первый ход опущен, поэтому вы не можете использовать ответы для создания дополнительных тестовых случаев.

гема
источник
Интересно, существуют ли альтернативные схемы нумерации, которые делают код более коротким?
qwr
1
Числители фракций: oeis.org/A157807 Знаменатели: oeis.org/A157813 Нет совпадений для порядковой последовательности: oeis.org/…
qwr
Понимаю. Вы должны уменьшить дробь, а затем считать. это не только зигзаг
Дон Яркий

Ответы:

4

Желе ,  21  20 байт

Вероятно, побиваемый количеством байтов, используя некоторую умную математику ...

:g/
ǵSRRUĖ€UÐeẎÇ€Qi

Монадическая ссылка, принимающая список, [p,q]который возвращает назначенное натуральное число p/q.

Попробуйте онлайн! Или посмотрите набор тестов .

Как?

Прежде всего обратите внимание, что встречающаяся N- я диагональ содержит все рациональные числа сетки, для которых сумма числителя и знаменателя равна N + 1 , поэтому, учитывая функцию, которая сводит [p,q]пару к простейшей форме ( [p/gcd(p,q),q/gcd(p,q)]), мы можем построить диагонали до нужно быть *, уменьшить все записи, дедуплицировать и найти индекс упрощенного ввода.

* на самом деле еще один здесь, чтобы сохранить байт

:g/ - Link 1, simplify a pair: list of integers, [a, b]
  / - reduce using:
 g  - Greatest Common Divisor -> gcd(a, b)
:   - integer division (vectorises) -> [a/gcd(a,b), b/gcd(a,b)]

ǵSRRUĖ€UÐeẎÇ€Qi - Main Link: list of integers, [p, q]
Ç                - call last Link as a monad (simplify)
 µ               - start a new monadic chain (call that V)
  S              - sum -> the diagonal V will be in plus one
   R             - range -> [1,2,3,...,diag(V)+1]
    R            - range (vectorises) -> [[1],[1,2],[1,2,3],...,[1,2,3,...,diag(V)+1]]
     U           - reverse each       -> [[1],[2,1],[3,2,1],[diag(V)+1,...,3,2,1]]
      Ė€         - enumerate €ach     -> [[[1,1]],[[1,2],[2,1]],[[1,3],[2,2],[3,1]],[[1,diag(V)+1],...,[diag(V)-1,3],[diag(V),2],[diag(V)+1,1]]]
         Ðe      - apply only to the even indexed items:
        U        -   reverse each     -> [[[1,1]],[[2,1],[1,2]],[[1,3],[2,2],[3,1]],[[4,1],[3,2],[2,3],[1,4]],...]
           Ẏ     - tighten            -> [[1,1],[2,1],[1,2],[1,3],[2,2],[3,1],[4,1],[3,2],[2,3],[1,4],...]
            Ç€   - for €ach: call last Link as a monad (simplify each)
                 -                    -> [[1,1],[2,1],[1,2],[1,3],[1,1],[3,1],[4,1],[3,2],[2,3],[1,4],...]
              Q  - de-duplicate       -> [[1,1],[2,1],[1,2],[1,3],[3,1],[4,1],[3,2],[2,3],[1,4],...]
               i - index of V in that list
Джонатан Аллан
источник
3

Perl 6 ,  94  90 байт

->\p,\q{(({|(1…($+=2)…1)}…*)Z/(1,{|(1…(($||=1)+=2)…1)}…*)).unique.first(p/q,:k)+1}

Проверь это

{(({|(1…($+=2)…1)}…*)Z/(1,{|(1…(($||=1)+=2)…1)}…*)).unique.first($^p/$^q):k+1}

Проверь это

Это в основном генерирует всю последовательность значений и останавливается, когда находит совпадение.

Expanded:

{  # bare block lambda with placeholder parameters $p,$q

  (
      ( # sequence of numerators

        {
          |( # slip into outer sequence (flatten)

            1      # start at one
            
            (
              $    # state variable
              += 2 # increment it by two each time this block is called
            )
            
            1      # finish at one
          )

        }
         * # never stop generating values
      )


    Z/   # zip using &infix:« /  » (generates Rats)


      ( # sequence of denominators

        1,  # start with an extra one

        {
          |( # slip into outer sequence (flatten)

            1
            
            (
              ( $ ||= 1 ) # state variable that starts with 1 (rather than 0)
              += 2        # increment it by two each time this is called
            )
            
            1
          )
        }
         * # never stop generating values
      )


  ).unique            # get only the unique values
  .first( $^p / $^q ) # find the first one that matches the input
  :k                  # get the index instead (0 based)
  + 1                 # add one               (1 based)
}

({1…($+=2)…1}…*)генерирует бесконечную последовательность числителей ( |(…)используется выше для выравнивания)

(1 2 1)
(1 2 3 4 3 2 1)
(1 2 3 4 5 6 5 4 3 2 1)
(1 2 3 4 5 6 7 8 7 6 5 4 3 2 1)
(1 2 3 4 5 6 7 8 9 10 9 8 7 6 5 4 3 2 1)

(1,{1…(($||=1)+=2)…1}…*) генерирует бесконечную последовательность знаменателей

1
(1 2 3 2 1)
(1 2 3 4 5 4 3 2 1)
(1 2 3 4 5 6 7 6 5 4 3 2 1)
(1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1)
(1 2 3 4 5 6 7 8 9 10 11 10 9 8 7 6 5 4 3 2 1)
Брэд Гилберт b2gills
источник
3

Python 2 , 157 144 137 134 126 125 байт

def f(p,q):a=[((j-i)/(i+1.))**(j%-2|1)for j in range(p-~q)for i in range(j)];return-~sorted(set(a),key=a.index).index(p*1./q)

Попробуйте онлайн!

4 байта сохранены благодаря г-ну Xcoder ; 1 байт от Джонатона Фреха .

Как отметил г-н Xcoder, в Python 3 мы можем сделать немного лучше, поскольку, среди прочего, целочисленное деление по умолчанию выводит результаты, и мы можем легче распаковать lists:

Python 3 , 117 байт

def f(p,q):a=[((j-i)/-~i)**(j%-2|1)for j in range(p-~q)for i in range(j)];return-~sorted({*a},key=a.index).index(p/q)

Попробуйте онлайн!

Час Браун
источник
128 байтов (-6) путем замены дроби и использования **(j%-2|1)и p-~q.
г-н Xcoder
@Г-н. Xcoder: Вы все о негативном модуле сегодня! :) Я думаю, что мне все еще нужно +1в конце, так как 1,1"должен" дать 1, а не 0.
Час Браун
125 байт .
Джонатан Фрех
124 байта :) Да, отрицательный модуль оказывается действительно полезным!
г-н Xcoder
На самом деле 117 байтов в Python 3
г-н Xcoder
3

Python 3 ,157, 146, 140133 байта

def f(p,q):a=[(i+i-abs(j-i-i))/(abs(j-i-i+.5)+.5)for i in range(p+q)for j in range(4*i)];return sorted(set(a),key=a.index).index(p/q)

Попробуйте онлайн!

выиграл 11 байтов благодаря Джонатану Фреху

выиграл еще 6 байтов, а затем 7 благодаря Chas Brown

bobrobbob
источник
1
146 байт .
Джонатан Фрех
140 байт .
Час Браун
@bobrobbob: Добро пожаловать в PPCG! Я не совсем уверен, как работает ваш алгоритм (хотя он явно работает); но экспериментально кажется, что вы можете сохранить еще несколько байтов, заменив их range(max(p,q)+1)на range(p+q).
Час Браун
1
Вы можете сохранить еще несколько байтов, используя {*a}вместо set(a).
г-н Xcoder
2

J, 41 , 35 , 30 байтов

-11 байт благодаря FrownyFrog

%i.~0~.@,@,[:]`|./.[:%/~1+i.@*

Попробуйте онлайн!

оригинальный 41 байт пост с объяснением

%>:@i.~[:([:~.@;[:<@|.`</.%"1 0)~(1+i.@*)

ungolfed

% >:@i.~ [: ([: ~.@; [: <@|.`</. %"1 0)~ 1 + i.@*

объяснение

                  +
                  | Everything to the right of this
                  | calculates the list
p (left arg)      |                                      create the
divided by q      |                                      diagonals.  yes,
      (right)     |                                     +this is a         +create this list
   |              |        ++ finally rmv ^alternately  |primitive ADVERB  |1..(p*q), and pass
   |   + the index          | the boxes,  |box, or      |in J              |it as both the left
   |   | of that  |         | and remove  |reverse and  |                  |and right arg to the
   |   | in the   |         | any dups    |box, each    |                  |middle verb, where each
   |   | list on  |         |             |diagonal     |                  |element will divide the
   |   | the right|         |             |             |       +----------+entire list, producing
   |   | plus 1   |         |             |             |       |          |the undiagonalized grid
   |   |          |         |             |             |       |          |
   |   |          |         |             |             |       |          |
   |   |          +         |             |             |       |          |
  ┌+┬──|──────────┬─────────|─────────────|─────────────|───────|──────────|─────────────┐
  │%│┌─+───────┬─┐│┌──┬─────|─────────────|─────────────|───────|────────┬─|────────────┐│
  │ ││┌──┬─┬──┐│~│││[:│┌────|─────────────|─────────────|───────|─────┬─┐│┌+┬─┬────────┐││
  │ │││>:│@│i.││ │││  ││┌──┬|───────┬─────|─────────────|───────|────┐│~│││1│+│┌──┬─┬─┐│││
  │ ││└──┴─┴──┘│ │││  │││[:│+──┬─┬─┐│┌──┬─|─────────────|─┬─────|───┐││ │││ │ ││i.│@│*││││
  │ │└─────────┴─┘││  │││  ││~.│@│;│││[:│┌|───────────┬─+┐│┌─┬─┬+──┐│││ │││ │ │└──┴─┴─┘│││
  │ │             ││  │││  │└──┴─┴─┘││  ││+────────┬─┐│/.│││%│"│1 0││││ ││└─┴─┴────────┘││
  │ │             ││  │││  │        ││  │││┌─┬─┬──┐│<││  ││└─┴─┴───┘│││ ││              ││
  │ │             ││  │││  │        ││  ││││<│@│|.││ ││  ││         │││ ││              ││
  │ │             ││  │││  │        ││  │││└─┴─┴──┘│ ││  ││         │││ ││              ││
  │ │             ││  │││  │        ││  ││└────────┴─┘│  ││         │││ ││              ││
  │ │             ││  │││  │        ││  │└────────────┴──┘│         │││ ││              ││
  │ │             ││  │││  │        │└──┴─────────────────┴─────────┘││ ││              ││
  │ │             ││  ││└──┴────────┴────────────────────────────────┘│ ││              ││
  │ │             ││  │└──────────────────────────────────────────────┴─┘│              ││
  │ │             │└──┴──────────────────────────────────────────────────┴──────────────┘│
  └─┴─────────────┴──────────────────────────────────────────────────────────────────────┘

Попробуйте онлайн!

Ион
источник
35:1+~.@;@([:<@|.`</.%~/~)@(1+i.@*)i.%
FrownyFrog
Спасибо очень мило! Я полностью обновлю его позже, так как это потребует повторения объяснения ...
Иона
И 30:%i.~0~.@,@,[:]`|./.[:%/~1+i.@*
FrownyFrog
это умно, используйте 0 и ~. чтобы избежать бокса и прироста, мне это нравится
Иона
2

Python 3, 121 байт

import math
def o(x,y):
 p=q=n=1
 while x*q!=p*y:a=p+q&1;p+=1+a*~-~-(p<2);q+=1-~-a*~-~-(q<2);n+=math.gcd(p,q)<2
 return n
orlp
источник
2

Ржавчина, 244 байта

Я создал простую формулу, чтобы найти «простой» порядковый номер «простого» зигзага без ограничений головоломки, используя формулы треугольных чисел: https://www.mathsisfun.com/algebra/triangular-numbers.html . Это было изменено по модулю 2, чтобы учесть, что зигзаги переворачивают свое направление в каждом диагональном ряду головоломки. Это функция h ()

Тогда для основного трюка этой головоломки: как не считать определенные повторяющиеся значения, такие как 3/3 против 1/1, 4/2 против 2/1, в зигзагообразной тропе. Я запустил 1-200 примеров и заметил разницу между простым зигзагообразным треугольным счетчиком и рисунком, который желает загадка. Шаблон «пропущенных» чисел - 5, 12, 13, 14, 23 и т. Д., Что привело к попаданию в OEIS. Это один из описанных Робертом Пнем, в https://oeis.org/A076537 , чтобы «дедуплицировать» числа типа 3/3, 4/2 и 1/1, вы можете проверить, является ли GCD> 1 для х, у всех «предыдущих» ординалов в зигзаге. Это циклы for, а g () - это gcd.

Я предполагаю, что с некоторыми встроенными gcd это было бы короче, но я вроде не мог найти его очень легко (я немного новичок в Rust, и Integer меня смутил), и мне нравится тот факт, что этот использует прямую целочисленную арифметику, и нет встроенных или библиотек любого рода.

fn f(x:i64,y:i64)->i64 {
        fn h(x:i64,y:i64)->i64 {let s=x+y;(s*s-3*s+4)/2-1+(s+1)%2*x+s%2*y}
        fn g(x:i64,y:i64)->i64 {if x==0 {y} else {g(y%x,x)}}
        let mut a=h(x,y);
        for i in 1..x+y {for j in 1..y+x {if h(i,j)<h(x,y) && g(i,j)>1 {a-=1;}}}
        a
}
не яркий
источник
1

JavaScript (ES6), 86 байт

Принимает ввод в синтаксисе карри (p)(q).

p=>q=>(g=x=>x*y?x*q-y*p?g(x+d,g[x/=y]=g[x]||++n,y-=d):n:g(x+!~d,y+=!~(d=-d)))(d=n=y=1)

Попробуйте онлайн!

Arnauld
источник
0

Javascript, 79 байт

a=(p,q)=>p*q==1?1:1+((p+q)%2?q==1?a(p-1,q):a(p+1,q-1):p==1?a(p,q-1):a(p-1,q+1))

(Я новичок в коде игры в гольф, так что это, вероятно, можно легко улучшить)

объяснение

let a = (p, q) => p * q == 1 // If they are both 1
    ? 1
    // Do a recursive call and increment the result by 1
    : 1 + (
        (p + q) % 2 // If on an even-numbered diagonal
        ? q == 1 // If at the beginning of the diagonal
            ? a(p - 1, q) // Go to previous diagonal
            : a(p + 1, q - 1) // Go one back
        : p == 1 // rougly the same
            ? a(p, q - 1)
            : a(p - 1, q + 1)
    )
Bary12
источник
4
(3,5)должно привести к 19(не 24) , так как (1,1)==(2,2)==(3,3), (2,4)==(1,2), (4,2)==(2,1)и (2,6)==(1,3). (т. е. (2,2)должно привести к 1не 5, и т. д ...)
Джонатан Аллан