Подсчет от 1 до n без последовательных чисел

19

Цель

Вам дано целое число n( n > 1). Вы должны вывода , сколько перестановок чисел 1на nесть , которые начинаются с 1, заканчиваются на n, и не имеют два последовательных целых чисел , которые отличаются на 1.

В качестве альтернативы, если вы возьмете полный граф K_nи удалите ребра пути, 1-2-3-...-nвы должны посчитать гамильтоновы пути от 1до nв оставшемся графе.

В примерах будет использоваться f(n)функция, которая принимает nи выводит количество допустимых перестановок, но ваша отправка может быть функцией или программой.


Примеры

Для n = 6возможного решения1-3-5-2-4-6

Тем 1-3-5-2-6-4не менее, это не является действительным решением, поскольку оно не заканчивается 6.

На самом деле, для n = 6, есть только 2 решения ( 1-4-2-5-3-6есть другое).

Отсюда f(6) = 2.


Для n = 4единственных перестановок, которые начинаются 1и заканчиваются в 4являются 1-2-3-4и 1-3-2-4. В обоих из них 2смежно с 3, давая последовательные целые числа, которые отличаются на 1. Следовательно f(4) = 0.


Контрольные примеры

f(6) = 2
f(4) = 0
f(8) = 68
f(13) = 4462848

Критерий победы

Это код-гольф, самый короткий ответ выигрывает.

Philippe
источник
7
Видите ли, дети, вы не можете просто проверить, сколько перестановок [2..n-1]не содержат дельт 1или -1, вы также должны проверить, что ни одна из них не начинается 2или не заканчивается n-1...
ETHproductions
1
Список должен начинаться с 1 и заканчиваться номером?
Okx
3
Может быть, ОП означает «смежный», а не «последовательный»?
Стилз
6
Причудливая последовательность здесь: algo.inria.fr/libraries/autocomb/graphs99.ps, где на странице 6 написано, Q_ser:=z + 2 z^6 + 10 z^7 + 68 z^8 + 500 z^9 + 4174 z^10 + 38774 z^11 + 397584z^12 + 4462848 z^13 + 54455754 z^14я сейчас провожу некоторое время, пытаясь использовать формулы, но я не могу составить тот, который генерирует последовательность. Удивительно, что показатель степени z является входом формулы, а результат - коэффициентом умножения. Тот, кто может вывести формулу из этого, может быть один с самым коротким ответом в байтах
Кристиан
1
@ChristiaanWesterbeek, которая называется производящей функцией для последовательности. Существует много последовательностей с генерирующей функцией, которая имеет более замкнутую форму, чем сама последовательность, это классно!
Кармейстер

Ответы:

6

MATL , 16 байт

qtq:Y@0&Yc!d|qAs

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

Для входов, превышающих 12его, не хватает памяти.

объяснение

q      % Implicitly input n. Push n-1
tq     % Duplicate and subtract 1: pushes n-2
:      % Range [1 2 ... n-2]
Y@     % Matrix with all permutations, each in a row
0      % Push 0
&Yc    % Append n-1 and predend 0 to each row
!      % Tranpose
d      % Consecutive differences along each column
|      % Absolute value
q      % Subtract 1
A      % All: true if all values in each column are non-zero
s      % Sum. Implicitly display
Луис Мендо
источник
1
Работает хорошо, молодец :)
Филипп
1
Хотя в этой проблеме произошел действительно хороший прогресс, ваше решение по-прежнему самое короткое. Это также быстрее, чем желе. Congratz!
Филипп
19

Mathematica, 58 байт, полиномиальное ( n ) время

Abs[Sum[(k-1)Hypergeometric2F1[k,k-#,2,2](#-k)!,{k,#}]-1]&

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

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

Пусть S - множество всех перестановок из [1,…, n] с σ 1 = 1, σ n = n , и пусть S i - множество перестановок σ ∈ S, таких что | σ i - σ i + 1 | = 1. Тогда искомый счет

| S | - | S 1 ∪ ⋯ ∪ S n - 1 | = ∑ 2 ≤ kn + 1; 1 ≤ i 2 <⋯ < i k - 1 < n (−1) k - 2 | S i 2 ∩ ⋯ ∩ S i k - 1 |.

Теперь | S i 2 ∩ ⋯ ∩ S i k - 1 | зависит только от k и от числа j серий последовательных индексов в [ i 1 , i 2 ,…, i k - 1 , i k ], где для удобства мы фиксируем i 1 = 0 и i k = n . В частности,

| S i 2 ∩ ⋯ ∩ S i k - 1 | = 2 Дж - 2 ( н - К ) !, для 2 ≤ jkn ,
| S i 2 ∩ ⋯ ∩ S i k - 1 | = 1, для j = 1, k = n + 1.

Число таких наборов индексов [ i 1 , i 2 ,…, i k - 1 , i k ] с j пробегами равно

( k - 1 C j - 1 ) ( n - k C j - 2 ), для 2 ≤ jkn ,
1, для j = 1, k = п + 1.

Результат тогда

(−1) n - 1 + ∑ 2 ≤ kn2 ≤ jk (−1) k - 2 ( k - 1 C j - 1 ) ( n - k C j - 2 ) 2 j - 2 ( н - к )!

Внутренняя сумма по J можно записать с помощью гипергеометрический 2 F 1 функцию :

(−1) n - 1 + ∑ 2 ≤ kn (−1) k ( k - 1) 2 F 1 (2 - k , k - n ; 2; 2) ( n - k )!

к которому мы применяем преобразование Пфаффа, которое позволяет нам убрать степени -1, используя абсолютное значение:

(−1) n - 1 + ∑ 2 ≤ kn (−1) n ( k - 1) 2 F 1 ( k , k - n ; 2; 2) ( n - k )!
= | −1 + ∑ 1 ≤ kn ( k - 1) 2 F 1 ( k , k - n ; 2; 2) ( n - k )! |.

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

In[1]:= Table[Abs[Sum[(k-1)Hypergeometric2F1[k,k-#,2,2](#-k)!,{k,#}]-1]&[n],{n,50}]

Out[1]= {1, 0, 0, 0, 0, 2, 10, 68, 500, 4174, 38774, 397584, 4462848, 

>    54455754, 717909202, 10171232060, 154142811052, 2488421201446, 

>    42636471916622, 772807552752712, 14774586965277816, 297138592463202402, 

>    6271277634164008170, 138596853553771517492, 3200958202120445923684, 

>    77114612783976599209598, 1934583996316791634828454, 

>    50460687385591722097602304, 1366482059862153751146376304, 

>    38366771565392871446940748410, 1115482364570332601576605376898, 

>    33544252621178275692411892779180, 1042188051349139920383738392594332, 

>    33419576037745472521641814354312790, 

>    1105004411146009553865786545464526206, 

>    37639281863619947475378460886135133496, 

>    1319658179153254337635342434408766065896, 

>    47585390139805782930448514259179162696722, 

>    1763380871412273296449902785237054760438426, 

>    67106516021125545469475040472412706780911268, 

>    2620784212531087457316728120883870079549134420, 

>    104969402113244439880057492782663678669089779118, 

>    4309132147486627708154774750891684285077633835734, 

>    181199144276064794296827392186304334716629346180848, 

>    7800407552443042507640613928796820288452902805286368, 

>    343589595090843265591418718266306051705639884996218154, 

>    15477521503994968035062094274002250590013877419466108978, 

>    712669883315580566495978374316773450341097231239406211100, 

>    33527174671849317156037438120623503416356879769273672584588, 

>    1610762789255012501855846297689494046193178343355755998487686}
Андерс Касеорг
источник
3
Мой ум взорван, хорошая работа
Филипп
6

Желе , 17 16 байт

ṖḊŒ!ð1;;⁹IỊṀðÐḟL

Монадическая ссылка.

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

Как?

ṖḊŒ!ð1;;⁹IỊṀðÐḟL - Link: number n
Ṗ                - pop (implicit range build) -> [1,n-1]
 Ḋ               - dequeue -> [2,n-1]
  Œ!             - all permutations of [2,n-1]
    ð       ðÐḟ  - filter discard those entries for which this is truthy:
     1;          -   1 concatenated with the entry
       ;⁹        -   ...concatenated with right (n)
         I       -   incremental differences
          Ị      -   is insignificant (absolute value <=1)
           Ṁ     -   maximum
               L - length (the number of valid arrangements)
Джонатан Аллан
источник
Извините, но это не соответствует контрольным случаям
Филипп
1
Да, вы сделали ту же ошибку, Оккс, и я сделал сначала. Вы должны учитывать тот факт, что второе число не может быть 2, а число от второго до последнего не может быть n-1
ETHproductions
@Philippe исправил это.
Джонатан Аллан
Я не думаю, что использование IỊṀявляется действительным. В частности, что, если, например, -2есть одна из дельт? Вы можете исправить с помощью IAỊṀ+1.
Эрик Outgolfer
1
@JonathanAllan Ох, я думал, что это вернулось x <= 1.
Эрик Outgolfer
5

Japt , 19 18 байт

o2 á è_pU äÉ m²e>1

Проверьте это онлайн!Я бы не рекомендовал тестировать что-либо большее, чем10 .

объяснение

o2 á è_  pU äÉ  m²  e>1
o2 á èZ{ZpU ä-1 mp2 e>1}
                          : Implicit: U = input integer
o2                        : Create the range [2..U-1].
   á                      : Generate all permutations of this range.
     èZ{               }  : Check how many permutations Z return a truthy value:
        ZpU               :   Push U to the end of Z.
            ä-1           :   Push 1 to the beginning of Z, then take the difference
                          :   of each pair of items.
                m         :   Map each item X to
                 p2       :     X ** 2. This gives a number greater than 1 unless the
                          :     item is 1 or -1.
                    e>1   :   Return whether every item in this list is greater than 1.
                          :   This returns `true` iff the permutation contains no
                          :   consecutive pairs of numbers.
                          : Implicit: output result of last expression
ETHproductions
источник
Молодец! Забавно, что мой код грубой силы не может быть больше n = 13, ахах
Филипп
@Philippe Я бы не рекомендовал принимать так быстро, я уверен, что это будет короче в 05AB1E или Jelly ;-)
ETHproductions
Сбои в тесте 1.
Okx
2
@Okx OP указал, что мы можем предположить n > 1.
ETHproductions
5

05AB1E , 17 байт

L¦¨œʒ¹1Š)˜¥Ä1å_}g

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

Okx
источник
Это не дает правильных результатов, извините
Филипп
@Philippe На каком тестовом примере?
Okx
@Philippe Исправлено.
Okx
¹1Š)˜сохраняет байт.
Emigna
5

Haskell, 76 65 байт

Сохранено 11 байт благодаря @xnor.

Используя результат Q_recпоиска на странице 7 из @ ChristiaanWesterbeek, мы получим

f 1=1
f n|n<6=0
f n=sum$zipWith((*).f)[n-5..][n-4,1,10-2*n,4,n-2]

Я не понимаю, как их следующий результат ha связан с этим, но после ускорения (сначала по запоминанию, см. Более ранние версии, затем, как показано ниже), я получаю их числа.

Хотя вышесказанное вполне подходит n=20, это, в сущности, пример того, как не делать рекурсию. Вот более быстрая версия (только для n>=6), которая также будет нуждаться только в постоянной памяти - если бы только цифры не продолжали увеличиваться ...

f n=last$foldl(#)[1,0,0,0,0][6..n]
l#n=tail l++[sum$zipWith(*)l[n-4,1,10-2*n,4,n-2]]

Что дает

Prelude> f 50
1610762789255012501855846297689494046193178343355755998487686
Prelude> f 500
659178618863924802757920269977240274180092211041657762693634630044383805576666007245903670780603497370173231423527767109899936008034229541700392144282505597945561328426013937966521561345817045884498867592832897938083071843810602104434376305964577943025310184523643816782047883794585616331928324460394146825636085453532404319881264974005968087265587062691285454120911586459406436421191277596121471930913837355151842093002557978076653884610826296845041929616496533544124347765641367732716560025553179112645454078955409181466212732427071306363820080109636358537270466838558068527692374178581063316309789026101221004745226182671038004326069705775312654329754698423385241664984156235692539255677944294995403233446243315371404887473868003155621849544566385172835597260848972758443874423271017007843907015007416644383573987606586308556317833384896267539628278571497402655322562624217658332870157802254043614726316296058329670971054977099155788604175817828380564156329839201579006169173002756295957371639199917376529472990059986681882194726437566769717959443857298155265292535858523609764515938314672724480762724541633037484152303637096

Это тоже не проблема, f 5000но я не хочу вставлять результат ...


Кстати, можно не использовать причудливую математику и все же не использовать (ультра) грубую силу. Во-первых, вместо просмотра всех перестановок, посмотрите на частичные перестановки и расширяйте их только тогда, когда они еще не являются недействительными. Бесполезно смотреть на все перестановки, начиная с 1 6 5. Во-вторых, некоторые частичные перестановки похожи 1 3 5 7и 1 5 3 7имеют одинаковые действительные продолжения, поэтому обрабатывайте их вместе. Используя эти идеи, я мог вычислить значения до n=16 0,3 с.

Кристиан Сиверс
источник
Вы можете написать рекурсивное выражение короче как точка, извлекая из коэффициентов: f n=sum$zipWith((*).f)[n-5..][n-4,1,10-2*n,4,n-2].
xnor
@xnor Хорошо, спасибо!
Кристиан Сиверс
Это хорошая работа, я удивлен результатами, которые придумали это сообщество! Жаль, что это гольф ^^
Филипп
4

Python, 125 байт

from itertools import*
lambda n:sum(p[-1]-p[0]==n-1and all(~-abs(x-y)for x,y in zip(p,p[1:]))for p in permutations(range(n)))
shooqie
источник
Выглядит довольно быстро, хорошая работа!
Филипп
2
117 байт
овс
3

Mathematica, 66 байт

Count[Permutations@Range@#,x:{1,__,#}/;FreeQ[Differences@x,1|-1]]&

объяснение

Functionс первым аргументом #.

Count[                                                             (* Count the number of *)
      Permutations@                                                (* permutations of *)
                   Range@#,                                        (* the list {1, ..., #} *)
                           x:{1,__,#}                              (* of the form {1, __, #} *)
                                     /;                            (* such that *)
                                             Differences@x,        (* the list of differences of consecutive elements *)
                                       FreeQ[                      (* is free of elements of the form *)
                                                           1|-1    (* 1 or -1 *)
                                                               ]]&
ngenisis
источник
3

Javascript (ES6), 100 74 72 60 байт

f=n=>n--<6?!n|0:f(n)*--n+4*f(n--)-2*f(n--)*--n+f(n)*++n+f(n)

Ниже приведена версия до гольф-мастерства @PeterTaylor

f=n=>n<6?n==1|0:(n-4)*f(n-5)+f(n-4)-2*(n-5)*f(n-3)+4*f(n-2)+(n-2)*f(n-1)

Благодаря ответу @ChristianSievers, которому удалось написать решение Haskell из бумаги которую я нашел после поиска в Google '0, 2, 10, 68, 500, 4174, 38774, 397584', вот версия Javascript, которая также не переставляет.

использование

for (i=1; i<=20; i++) {
  console.log(i, f(i))
}

1 1 
2 0 
3 0 
4 0 
5 0 
6 2 
7 10 
8 68 
9 500 
10 4174 
11 38774 
12 397584 
13 4462848 
14 54455754 
15 717909202 
16 10171232060 
17 154142811052 
18 2488421201446 
19 42636471916622 
20 772807552752712
Кристиан Вестербик
источник
1
В описании задачи спрашивается только f(n)когда n>1, поэтому не имеет значения, за что вы возвращаетесь n=1. Тоже считаю f(1)=1правильным.
Кристиан Сиверс
Вы можете комбинировать особые случаи как n<6?n==1|0:для дальнейшего сохранения двух символов.
Питер Тейлор
Отлично. Я с поправкой на эти 2 комментария.
Кристиан
1
И переупорядочив условия и полагаясь на порядок оценки, можно снизить до 60:f=n=>n--<6?!n|0:f(n)*--n+4*f(n--)-2*f(n--)*--n+f(n)*++n+f(n)
Питер Тейлор
1

Брахилог , 26 байт

{⟦₁pLh1&~tLs₂ᶠ{-ȧ>1}ᵐ}ᶜ|∧0

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

объяснение

{                    }ᶜ       Output = count the number of outputs of:
 ⟦₁pL                           L is a permutation of [1, …, Input]
    Lh1                         The head of L is 1
       &~tL                     The tail of L is the Input
          Ls₂ᶠ                  Find all sublists of length 2 of L
              {    }ᵐ           Map on each sublist:
               -ȧ>1               The elements are separated by strictly more than 1
                       |      Else (no outputs to the count)
                        ∧0    Output = 0
Fatalize
источник
1

Python 3 , 109 107 102 байта

q=lambda s,x,n:sum(q(s-{v},v,n)for v in s if(v-x)**2>1)if s else x<n;f=lambda n:q({*range(2,n)},1,n-1)

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

Удалил четыре байта, не пытаясь поставить в одну строку функцию (как предложено @shooqie) и другой байт, заменив absквадрат. (Требуется Python 3.5+)

RICi
источник
103 байта
shooqie
0

Mathematica, 134 байта

(s=Permutations@Range[2,#-1];g=Table[Join[Prepend[s[[i]],1],{#}],{i,Length@s}];Length@Select[Union@*Abs@*Differences/@g,FreeQ[#,1]&])&


контрольные примеры от 2 до 12

{0, 0, 0, 0, 2, 10, 68, 500, 4174, 38774, 397584}

J42161217
источник
0

Python 2 , 105 байт

lambda n:reduce(lambda a,i:a+[i*a[-5]+a[-4]+2*(1-i)*a[-3]+4*a[-2]+(i+2)*a[-1]],range(2,n),[0,1]+4*[0])[n]

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

Это основано на статье Филиппа Флажоле, обнаруженной @Christiaan Westerbeek ; это намного быстрее и на два байта короче моего решение Python 3, которое перечисляет возможные перестановки. (В Python 3 reduceбыл досадно перенесен вfunctools .)

Существует гораздо более короткая версия, использующая точечный продукт numpy, но она переполняется довольно быстро и требует импорта numpy. Но для чего это стоит:

lambda n:reduce(lambda a,i:a+[dot([i,1,2-2*i,4,i+2],a[-5:])],range(2,n),[0,1]+4*[0])[n]
RICi
источник