Trova i numeri mancanti nella sequenza K di Fibonacci


20

Ispirato da questa domanda Math.SE .

sfondo

La sequenza di Fibonacci (chiamata F) è la sequenza, che inizia in modo 0, 1tale che ogni numero ( F(n)) (dopo i primi due) sia la somma dei due prima di esso ( F(n) = F(n-1) + F(n-2)).

Una sequenza K di Fibonacci mod (chiamata M) è la sequenza dei numeri di Fibonacci mod K ( M(n) = F(n) % K).

Si può dimostrare che la sequenza K di Fibonacci mod è ciclica per tutte le K, poiché ogni valore è determinato dalla coppia precedente, e ci sono solo K 2 possibili coppie di numeri interi non negativi entrambi inferiori a K. Perché la sequenza di Fibonacci mod K è ciclico dopo la prima coppia ripetuta di termini, un numero che non appare nella sequenza K di Fibonacci mod prima che la prima coppia ripetuta di termini non compaia mai.

Per K = 4

0 1 1 2 3 1 0 1 ...

Per K = 8

0 1 1 2 3 5 0 5 5 2 7 1 0 1 ...

Si noti che per K = 8, 4 e 6 non compaiono prima del ripetuto 0 1, quindi 4 e 6 non appariranno mai nella sequenza di Fibonacci mod 8.

Sfida

Dato un intero K strettamente maggiore di 0, genera tutti gli interi non negativi inferiori a K che non compaiono nella sequenza di Fibonacci mod K.

Regole

  • Sono vietate le scappatoie predefinite .

  • I / O predefinito .

  • Programmi o funzioni sono accettabili .

  • Puoi presumere che K si adatterà al tuo tipo intero nativo ( entro limiti ragionevoli ).

  • Se ci sono numeri non negativi inferiori a K che non compaiono nella sequenza K di Fibonacci mod K, il tuo programma / funzione dovrebbe generare tutti questi numeri in modo ragionevole.

  • Se non ci sono numeri interi non negativi inferiori a K che non compaiono nella sequenza K di Fibonacci mod, il tuo programma / funzione può indicarlo restituendo un elenco vuoto, stampando nulla, producendo un errore, ecc.

  • L'ordine non ha importanza.

  • Questo è , quindi vince la risposta più breve in ogni lingua.

Casi test

Genera casi di test online!

Casi di prova non vuoti

  8 [4, 6]
 11 [4, 6, 7, 9]
 12 [6]
 13 [4, 6, 7, 9]
 16 [4, 6, 10, 12, 14]
 17 [6, 7, 10, 11]
 18 [4, 6, 7, 9, 11, 12, 14]
 19 [4, 6, 7, 9, 10, 12, 14]
 21 [4, 6, 7, 9, 10, 11, 12, 14, 15, 16, 17, 19]
 22 [4, 6, 7, 9, 15, 17, 18, 20]
 23 [4, 7, 16, 19]
 24 [4, 6, 9, 11, 12, 14, 15, 18, 19, 20, 22]
 26 [4, 6, 7, 9, 17, 19, 20, 22]
 28 [10, 12, 14, 16, 18, 19, 23]
 29 [4, 6, 7, 9, 10, 11, 12, 14, 15, 16, 17, 18, 19, 20, 22, 23, 24, 25, 27]
 31 [4, 6, 9, 12, 14, 15, 17, 18, 19, 22, 25, 29]
 32 [4, 6, 10, 12, 14, 18, 20, 22, 26, 28, 30]
 33 [4, 6, 7, 9, 15, 17, 18, 20, 24, 26, 27, 28, 29, 31]
 34 [4, 6, 7, 9, 10, 11, 12, 14, 15, 16, 17, 18, 19, 20, 22, 23, 24, 25, 27, 28, 30]
 36 [4, 6, 7, 9, 10, 11, 12, 14, 16, 18, 20, 22, 23, 24, 25, 26, 27, 29, 30, 31, 32]
 37 [9, 10, 14, 17, 20, 23, 27, 28]
 38 [4, 6, 7, 9, 10, 11, 12, 14, 15, 16, 18, 19, 20, 22, 23, 24, 25, 26, 27, 28, 29, 31, 32, 33, 36]
 39 [4, 6, 7, 9, 15, 17, 19, 20, 22, 24, 30, 32, 33, 35]
...
200 [4, 6, 12, 14, 20, 22, 28, 30, 36, 38, 44, 46, 52, 54, 60, 62, 68, 70, 76, 78, 84, 86, 92, 94, 100, 102, 108, 110, 116, 118, 124, 126, 132, 134, 140, 142, 148, 150, 156, 158, 164, 166, 172, 174, 180, 182, 188, 190, 196, 198]
...
300 [6, 18, 30, 42, 54, 66, 78, 90, 102, 114, 126, 138, 150, 162, 174, 186, 198, 210, 222, 234, 246, 258, 270, 282, 294]
...
400 [4, 6, 10, 12, 14, 20, 22, 26, 28, 30, 36, 38, 42, 44, 46, 52, 54, 58, 60, 62, 68, 70, 74, 76, 78, 84, 86, 90, 92, 94, 100, 102, 106, 108, 110, 116, 118, 122, 124, 126, 132, 134, 138, 140, 142, 148, 150, 154, 156, 158, 164, 166, 170, 172, 174, 180, 182, 186, 188, 190, 196, 198, 202, 204, 206, 212, 214, 218, 220, 222, 228, 230, 234, 236, 238, 244, 246, 250, 252, 254, 260, 262, 266, 268, 270, 276, 278, 282, 284, 286, 292, 294, 298, 300, 302, 308, 310, 314, 316, 318, 324, 326, 330, 332, 334, 340, 342, 346, 348, 350, 356, 358, 362, 364, 366, 372, 374, 378, 380, 382, 388, 390, 394, 396, 398]
...

Casi di test vuoti (nessun output, errore, elenco vuoto, ecc. È un output accettabile)

1, 2, 3, 4, 5, 6, 7, 9, 10, 14, 15, 20, 25, 27, 30, 35 ... 100 ...

Relazionato:

Conteggio delle orbite di Fibonacci

Trova il periodo pisano


Sandbox (cancellato).
pizzapants184

Risposte:


6

Gelatina , 9 8 byte

²RÆḞ%ḟ@Ḷ

Provalo online!

Basato sul periodo pisano p(n) <= 6nda A001175 . Inoltre, p(n) <= 6n <= n^2per n >= 6e p(n) <= n^2per n < 6. Ho salvato questo byte grazie a Dennis.


²dovrebbe funzionare invece di ×6.
Dennis,

6

Haskell , 70 byte

Una certa quantità di byte salvati grazie a Esolanging Fruit

8 byte salvati grazie a Laikoni

a=1:scanl(+)1a
f x=[u|u<-[2..x-1],and[mod b x/=u|(_,b)<-zip[1..x^2]a]]

Provalo online!


@EsolangingFruit Ah grazie! Stavo arrivando a una conclusione simile anch'io.
Wheat Wizard

read$showfunziona invece che fromIntegerin questo caso e salva due byte.
Laikoni,

L'uso zip[1..x^2]per il troncamento consente di risparmiare qualche byte in più: provalo online!
Laikoni,

@Laikoni Ci è voluto un po 'ma ho apportato la modifica. Grazie, è una buona idea.
Wheat Wizard

5

Perl 6 ,  43 42 39  32 byte

{^$_ (-)(1,1,(*+*)%$_...->\a,\b{!a&&b==1})}

Provalo

{^$_∖(1,1,(*+*)%$_...->\a,\b{!a&&b==1})}

Provalo

{^$_∖(1,1,(*+*)%$_...{!$^a&&$^b==1})}

Provalo

{^$_∖(1,1,(*+*)%$_...!*&*==1)}

Provalo

Allargato:

{  # bare block lambda with implicit parameter 「$_」

  ^$_               # Range upto and excluding the input

                   # set minus (U+2216)

  (                 # generate the Fibonacci sequence mod k

    1, 1,           # seed the sequece (can't be 0,1)

    ( * + * ) % $_  # add two values and modulus the input (lambda)

    ...             # keep doing that until

                    # it matches 0,1
    !*              #   negate the first param (1 when 0)
    &               #   and Junction
    *               #   second param
    == 1            #   both match 1

  )
}

3

> <> , 48 byte

01\
?!\:&+{:}%:1$0p&$:
v0\~:1=?
>?!;1-::0g?!nao:

Provalo online!

Riceve input tramite il flag -v.

Stampa molte nuove righe in eccesso, ma esegue il lavoro. Questo utilizza fondamentalmente la prima riga per memorizzare l'insieme di numeri che sono apparsi finora nella sequenza.

Come funziona:

01\    Input is already on the stack
...... Initialises the sequence with 1 and 0
...... Goes to the second line
......

......
..\:&+{:}% Gets the next number in the modded Fibonacci sequence while preserving the previous number
......
......

......
..........:1$0p&$: Puts a 1 at that cell number on the first line
.......
.......

......             If the number is a 0 go to the third line
?!\..............: Check if the next number is a 1, meaning we've reached the end of the sequence
v0\~:1=?           Go to the fourth line if so
>.....             Re-add the 0 and go back to the second line if not

......           While input:
......             Get the cell from the first line
......             If not 0: print the number
>?!;1-::0g?!nao:   Finally, print a newline and decrement the input


3

MATL , 19 18 byte

0lbU:"yy+]vG\G:qX~

Provalo online!

-1 byte grazie a Guiseppe.

  bU:"   ]         % Do K^2 (>6K) times.
0l    yy+          %  Fibbonaci
                X~ % Set exclusive difference between
          vG\      %  the fibonacci numbers mod K
             G:q   %  and 0...K-1

18 byte ; riorganizzare recupera il tuo uso di X~!
Giuseppe,

@Giuseppe Grazie! Ancora molto a lungo però ...
Sanchises il


2

Buccia , 13 12 10 byte

Grazie @Zgarb per -2 byte!

-U2m%⁰İfŀ⁰

Stampa un elenco vuoto nel caso in cui vengano visualizzati tutti i numeri interi, provalo online!

Spiegazione

-U2m%⁰İfŀ⁰  -- named argument ⁰, example with: 8
-           -- difference of
        ŀ⁰  -- | lowered range: [0,1,2,3,4,5,6,7]
            -- and
      İf    -- | Fibonacci sequence: [1,1,2,3,5,8,13,21,34,55,89,144,233,377…
   m%⁰      -- | map (modulo ⁰): [1,1,2,3,5,0,5,5,2,7,1,0,1,1…
 U2         -- | keep longest prefix until 2 adjacent elements repeats: [1,1,2,3,5,0,5,5,2,7,1,0,1]
            -- : [4,6]

È possibile utilizzare U2per ottenere il prefisso più lungo in cui nessuna coppia adiacente si ripete.
Zgarb,


2

R, 92 86 byte

Grazie a @Giuseppe per aver salvato 6 byte!

function(k,n=!!0:2){while(any((z=tail(n,2))-n[1:2]))n=c(n,sum(z)%%k);setdiff(1:k-1,n)}

Provalo online!

Implementazione piuttosto semplice ( versione precedente , ma stesso concetto):

function(k,
         K=1:k-1,      #Uses default arguments to preset variables for legibility 
         n=c(0,1,1)){  #(wouldn't change byte-count to put them in the body of the function)
    while(any((z=tail(n,2))!=n[1:2])) #Do as long as first 2 elements are not identical to last 2 elements
        n=c(n,sum(z)%%k) #Built the fibonacci mod k sequence
    K[!K%in%n] #Outputs integers < k if not in sequence.
}


@Giuseppe ah setdiff, buona idea!
plannapus,

70 byte porting l' 1:k^2approccio che tutti usano
Giuseppe

2

Python 3, 173 152 143 131 byte

f=lambda n,m,a=0,b=1:a%m if n<=0else f(n-1,m,b,a+b)
p=lambda n,i=2,y={0}:y^{*range(n)}if f(i,n)==1>f(i-1,n)else p(n,i+1,y|{f(i,n)})

Un ringraziamento speciale a @ovs.

Provalo online

Come funziona?

La prima funzione accetta due parametri m e n, e restituisce l'ennesimo numero di Fibonacci mod m. La seconda funzione scorre attraverso i numeri di Fibonacci mod k e verifica se si ripetono 0 e 1. Memorizza i numeri in un elenco e li confronta con un elenco contenente i numeri 1-n. I numeri duplicati vengono rimossi e vengono restituiti i numeri rimanenti.


Fa parte dell'intestazione e non è obbligatorio includerlo nel codice.
Manish Kundu,

Va bene. @ovs Grazie per avermelo detto, non ne ero a conoscenza.
Manish Kundu,

1
131 byte creando set con parentesi graffe anziché set()e confronti concatenati.
Ovs


2

Rubino , 47 byte

->n{a=b=1;[*1...n]-(1..n*n).map{a,b=b,a+b;a%n}}

Provalo online!

Mentre utilizza una parte della stessa logica, questo non è basato sulla risposta di GB .

Spiegazione:

->n{
  a=b=1;   # start sequence with 1,1
  [*1...n] # all the numbers from 1 to n-1 as an array
           # 0 is excluded as it should never be in the final answer 
  -  # set operation; get all items in the first set and not in the second
  (1..n*n).map{ # n squared times
    a,b=b,a+b;  # assign next fibonacci numbers 
    a%n         # return a fibonacci number mod n
  }    # Map to an array
}

2

Lisp comune, 106 byte

(lambda(k)(do((a 1 b)c(b 1(mod(+ a b)k)))((=(1- b)0 a)(dotimes(i k)(or(member i c)(print i))))(push a c)))

Provalo online!



1

Elisir , 148 144 byte

 fn x->Enum.to_list(1..x-1)--List.flatten Enum.take_while Stream.chunk(Stream.unfold({1,1},fn{p,n}->{rem(p,x),{n,p+n}}end),2),&Enum.sum(&1)!=1end

Provalo online!

Non è una risposta particolarmente competitiva, ma è stato davvero divertente giocare a golf! L'elisir è una lingua abbastanza leggibile, ma segue una spiegazione per il caos dei personaggi nel mezzo.


Questa spiegazione è in due sezioni, il mod-fibonacci e il funzionamento su di esso

Mod-fib:

Stream.unfold({1,1},fn{p,n}->{rem(p,x),{n,p+n}}end)

Questo restituisce un flusso infinito di mod fibonacci x. Inizia con un accumulatore {1,1}e applica all'infinito la seguente operazione: dato accumulatore {p,n}, uscita p mod xallo stream. Quindi, impostare l'accumulatore su {n,p+n}.

Il riposo:

fn x->                              Define a fxn f(x) that returns
  Enum.to_list(1..x-1)--            The numbers from 1..x-1 that are not in
  List.flatten                      The flattened list constructed by
    Enum.take_while                 Taking from mod-fib until
      Stream.chunk(                 A 2-size chunk
        Stream.unfold({1,1},fn{p,n}->{rem(p,x),{n,p+n}}end) (of mod fib)
        ,2)
      ,&Enum.sum(&1)!=1             sums to 1, representing [0,1] or [1,0]
end



1

JavaScript (ES6), 84 byte

f=(n,a=0,b=1,q=[...Array(n).keys()])=>a*b+a-1?f(n,b,(a+b)%n,q,q[b]=0):q.filter(x=>x)

1

Python 3, 76 byte

def t(n,r=[1]):
 while n*n>len(r):r+=[sum(r[-2:])%n]
 return{*range(n)}-{*r}

Questo semplicemente esamina il ciclo più lungo possibile di numeri di Fibonnaci (n ^ 2) e crea un elenco di tutti i numeri che si verificano in quel momento. Per semplificare la logica i numeri sono memorizzati modulo n.

Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.