Рассчитать MU-номера


19

Первые два MU-числа - это 2 и 3. Любое другое MU-число - это наименьшее число, которое еще не появилось, и которое может быть выражено как произведение двух ранее различных MU-номеров точно одним способом.

Вот первые 10

2, 3, 6, 12, 18, 24, 48, 54, 96, 162

задача

При заданном положительном числе рассчитать и вывести n- ное MU-число.

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

OEIS A007335


1
0-индексация или 1-индексация?
HyperNeutrino

1
@HyperNeutrino Либо в порядке.
Пшеничный волшебник

2
Есть идеи, почему они называются MU-номерами? (

Ответы:


5

Pyth, 22 21 байт

@u+Gfq2l@GcLTGheGQhB2

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

0 индексированные.

объяснение

@u+Gfq2l@GcLTGheGQhB2Q    Implicitly append Q and read+eval input to it.
                  hB2     Take the list [2, 2 + 1].
 u               Q        Put the list in G and apply this Q times:
               eG           Get last number in G.
              h             Add one.
    f                       Starting from that, find the first T such that:
          cLTG                Divide T by each of the numbers in G.
        @G                    Find the quotients that are also in G.
       l                      Get the number of such quotients.
     q2                       Check that it equals 2.
  +G                        Append that T to G.
@                    Q    Get the Q'th number in G.

@Знак на последней строке смещена. Я не могу сделать предлагаемое редактирование, так как это двухзначное изменение.
user2357112 поддерживает Монику

@ user2357112 Исправлено.
PurkkaKoodari

4

Haskell, 80 77 байтов

l#(a:b)|[x]<-[a|i<-l,j<-l,i<j,i*j==a]=a:(a:l)#b|1<2=l#b
((2:3:[2,3]#[4..])!!)

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

Как это устроено

2:3:             -- start the list with 2 and 3 and append a call to # with
    [2,3]        -- the list so far and
         #[4..]  -- list of candidate elements

l # (a:b)        -- l -> list so far, a -> next candidate element, b -> rest c.el.
  | [x]<-[...]   -- if the list [...] is a singleton list
    =a:(a:l#b) -- the result is a followed by a recursive call with l extended
                    by a and b
  | 1<2=l#b      -- if it's not a singleton list, drop a and retry with b

                 -- the [...] list is
 [ i<-l,j<-l,    -- loop i through l and j through l and whenever   
       i<j,      -- i<j and
       i*j==a]   -- i*j==a
  a|             -- add a to the list              

3

Желе , 22 байта

ŒcP€ḟ⁸ṢŒgLÞḢḢṭ
2,3Ç¡ị@

Монадическая ссылка, 1-индексированная.

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

Как?

ŒcP€ḟ⁸ṢŒgLÞḢḢṭ - Link 1, add the next number: list, a  e.g. [2,3,6,12,18,24]
Œc             - unordered pairs                            [[2,3],[2,6],[2,12],[2,18],[2,24],[3,6],[3,12],[3,18],[3,24],[6,12],[6,18],[6,24],[12,18],[12,24],[18,24]]
  P€           - product of €ach                            [6,12,24,36,48,18,36,54,72,72,108,144,216,288,432]
     ⁸         - chain's left argument, a                   [2,3,6,12,18,24]
    ḟ          - filter discard                             [36,48,36,54,72,72,108,144,216,288,432]
      Ṣ        - sort                                       [36,36,48,54,72,72,108,144,216,288,432]
       Œg      - group runs of equal elements               [[36,36],[48],[54],[72,72],[108],[144],[216],[288],[432]]
          Þ    - sort by:
         L     -   length                                   [[48],[54],[108],[144],[216],[288],[432],[36,36],[72,72]]
           Ḣ   - head                                       [48]
            Ḣ  - head                                       48
             ṭ - tack to a                                  [2,3,6,12,18,24,48]

2,3Ç¡ị@ - Link: number, i                              e.g. 7
2,3     - literal [2,3]                                     [2,3]
    ¡   - repeat i times:
   Ç    -   call last link (1) as a monad                   [2,3,6,12,18,24,48,54,96]
     ị@ - index into with swapped @rguments (with i)        48

3

R , 127 118 111 108 105 100 98 90 байтов

8 байтов благодаря Джузеппе.

r=3:2;for(i in 1:scan())r=c(min((g=(r%o%r)[i:-1<i])[colSums(g%o%g==g*g)+g%in%r<3]),r);r[3]

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


Мне потребовалось целую вечность, чтобы осознать, что это <имеет более низкий приоритет, чем +я, поэтому я не мог понять, что, черт возьми, +g%in%r<3делал, и пока я делал это, вы играли в гольф две части, которые я собирался предложить ... +1
Джузеппе

@ Giuseppe Я только начал изучать R сегодня ... приятно познакомиться с приличным игроком в R.
Утренняя монахиня

Я собирался сказать тебе то же самое .............
Джузеппе

Ах, еще одна вещь, которую вы можете использовать n=scan()вместо определения функции для чтения из stdin; это
Джузеппе

Ошибка ввода:0
Rift

2

CJam (32 байта)

4,{_2m*{~>},::*1$-$e`$0=|}qi*-2=

Онлайн демо с 0-индексацией.

Я не уверен, что можно сделать что-то еще, кроме тривиального перевода спецификации, с одним исключением: начиная со списка [0 1 2 3](вместо [2, 3]), я сразу же сохраняю один байт при инициализации и еще два, когда могу сделать 0=|(добавляя только новый элемент, потому что его частота есть 1и уже есть в списке), но не вводите никаких ложных элементов, потому что для каждого xв списке 0*xи 1*xуже есть в списке.



1

Mathematica, 154 байта

простая модификация кода, найденная по ссылке oeis

(s={2,3};Do[n=Select[Split@Sort@Flatten@Table[s[[j]]s[[k]],{j,Length@s},{k,j+1,Length@s}],#[[1]]>s[[-1]]&&Length@#==1&][[1,1]];AppendTo[s,n],{#}];s[[#]])&

1

PHP , 130 байт

0 индексированные

for($r=[2,3];!$r[$argn];$r[]=$l=min($m)/2){$m=[];foreach($r as$x)foreach($r as$y)($p=$x*$y)<=$l|$y==$x?:$m[$p]+=$p;}echo$r[$argn];

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

расширенный

for($r=[2,3];!$r[$argn]; #set the first to items and loop till search item exists
$r[]=$l=min($m)/2){ # add the half of the minimum of found values to the result array
  $m=[]; # start with empty array
  foreach($r as$x) # loop through result array
    foreach($r as$y) # loop through result array
      ($p=$x*$y)<=$l|$y==$x? # if product is greater as last value and we do multiple two distinct values
        :$m[$p]+=$p; # add 2 times or more the product to array so we drop 36 cause it will be 144  
}
echo$r[$argn]; # Output 

PHP , 159 байт

0 индексированные

for($r=[2,3];!$r[$argn];$r[]=$l=min(array_diff_key($m,$d))){$d=$m=[];foreach($r as$x)foreach($r as$y)$x<$y?${dm[$m[$p=$x*$y]<1&$p>$l]}[$p]=$p:0;}echo$r[$argn];

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

PHP , 161 байт

0 индексированные

for($r=[2,3];!$r[$argn];$r[]=$l=min(array_diff($m,$d))){$d=$m=[];foreach($r as$x)foreach($r as$y)$x<$y?${dm[!in_array($p=$x*$y,$m)&$p>$l]}[]=$p:0;}echo$r[$argn];

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


1

Mathematica, 140 байт

(t=1;s={2,3};While[t<#,s=AppendTo[s,Sort[Select[First/@Select[Tally[Times@@@Permutations[s,{2}]],#[[2]]==2&],#>Last@s&]][[1]]];t++];s[[#]])&

1

MATL , 25 байт

3:i:"t&*9B#u2=)yX-X<h]2_)

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

объяснение

3:     % Push [1 2 3]. Initial array of MU numbers, to be extended with more numbers
i:     % Input n. Push [1 2 ... n]
"      % Do this n times
  t    %   Duplicate array of MU numbers so far
  &*   %   Matrix of pair-wise products
  9B   %   Push 9 in binary, that is, [1 0 0 1]
  #    %   Specify that next function will produce its first and fourth ouputs
  u    %   Unique: pushes unique entries (first output) and their counts (fourth)
  2=   %   True for counts that equal 2
  )    %   Keep only unique entries with count 2
  y    %   Duplicate (from below) array of MU numbers so far
  X-   %   Set difference
  X<   %   Minimum. This is the new MU number
  h    %   Concatenate vertically horizontally to extend the array
]      % End
2_     % Push 2 negated, that is, -2
)      % Get entry at position -2, that is, third-last. Implicitly display

1

Perl 6 , 96 байт

{(2,3,{first *∉@_,@_.combinations(2).classify({[*]
$_}).grep(*.value==1)».key.sort}...*)[$_]}

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

  • 2, 3, { ... } ... *является бесконечной последовательностью, где каждый элемент, начинающийся с третьего, вычисляется с помощью блока кода с разделителями-скобками. Поскольку блок кода принимает свои аргументы через @_массив slurpy , он получает всю текущую последовательность в этом массиве.
  • @_.combinations(2) является последовательностью всех 2-элементных комбинаций @_ .
  • .classify({ [*] $_ }) классифицирует каждый 2-кортеж по его продукту, создавая хеш, где продукты - это ключи, а значения - список 2-кортежей, имеющих этот продукт.
  • .grep(*.value == 1) выбирает эти пары ключ-значение из хэша, где значение (т. е. список пар, имеющих этот ключ в качестве продукта) имеет размер 1.
  • ».keyвыбирает только ключи каждой пары. Это список продуктов, которые возникают только из одной комбинации факторов текущей последовательности.
  • .sort сортирует продукты численно.
  • first * ∉ @_, ... находит первый из тех товаров, которые еще не появились в последовательности.

1

JavaScript (ES6), 119 118 117 байт

Рекурсивная функция, которая принимает индекс на основе 0.

f=(n,a=[2,m=3])=>a[n]||a.map(c=>a.map(d=>c<d&(d*=c)>m?b[d]=b[d]/0||d:0),b=[])|f(n,a.push(m=b.sort((a,b)=>a-b)[0])&&a)

Как?

На каждой итерации функции f () мы используем последний член m последовательности и первоначально пустой массив b для идентификации следующего члена. Для каждого произведения d> m из двух ранее различных MU-чисел мы делаем:

b[d] = b[d] / 0 || d

а затем сохранить минимальное значение б .

Вышеупомянутое выражение оценивается следующим образом:

b[d]               | b[d] / 0  | b[d] / 0 || d
-------------------+-----------+--------------
undefined          | NaN       | d
already equal to d | +Infinity | +Infinity
+Infinity          | +Infinity | +Infinity

Это гарантирует, что продукты, которые можно выразить более чем одним способом, никогда не будут выбраны.

Отформатировано и прокомментировано

f = (n, a = [2, m = 3]) =>           // given: n = input, a[] = MU array, m = last term
  a[n] ||                            // if a[n] is defined, return it
  a.map(c =>                         // else for each value c in a[]:
    a.map(d =>                       //   and for each value d in a[]:
      c < d &                        //     if c is less than d and
      (d *= c) > m ?                 //     d = d * c is greater than m:
        b[d] = b[d] / 0 || d         //       b[d] = either d or +Infinity (see 'How?')
      :                              //     else:
        0                            //       do nothing
    ),                               //   end of inner map()
    b = []                           //   initialization of b[]
  ) |                                // end of outer map()
  f(                                 // do a recursive call:
    n,                               //   - with n
    a.push(                          //   - push in a[]:
      m = b.sort((a, b) => a - b)[0] //     m = minimum value of b[]
    ) && a                           //     and use a[] as the 2nd parameter
  )                                  // end of recursive call

демонстрация



0

Python 3 2 , 167 139 136 133 123 121 120 118 байт

a=[2,3];exec'p=[x*y for x in a for y in a if x-y];a+=min(q for q in p if p.count(q)+(q in a)<3),;'*input();print a[-2]

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


Спасибо @ Mr.Xcoder и @LeakyNun за улучшения!


159 байт , просто удалив ненужные пробелы и скобки.
г-н Xcoder

@ Mr.Xcoder Спасибо за улучшения. Я не уверен, что изменение p.count(q)==1на p.count(q)>0допустимо, потому что это код, который обеспечивает «точно одно условие» задачи.
Погоня Фогели

p.count(q)-~(q in a)<=3эквивалентноp.count(q)+(q in a)<3
Leaky Nun

@ LeakyNun спасибо!
Погоня Фогели
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.