Encontre a enésima Fibohexaprime


23

O desafio neste momento é encontrar o n º Fibohexaprime . A definição de um Fibohexaprime é a seguinte:

Primeiro observamos uma lista com números de Fibonacci:

N  | Fibonacci number
1  | 1 
2  | 1 
3  | 2 
4  | 3 
5  | 5 
6  | 8 
7  | 13 
8  | 21 
9  | 34 
10 | 55 
11 | 89 
12 | 144 
13 | 233 
14 | 377 
15 | 610
16 | 987 
17 | 1597

Depois disso, convertemos os números em hexadecimal:

N  | Fib  | Hex 
1  | 1    | 1
2  | 1    | 1
3  | 2    | 2
4  | 3    | 3
5  | 5    | 5
6  | 8    | 8
7  | 13   | D
8  | 21   | 15
9  | 34   | 22
10 | 55   | 37
11 | 89   | 59
12 | 144  | 90
13 | 233  | E9
14 | 377  | 179
15 | 610  | 262
16 | 987  | 3DB
17 | 1597 | 63D

A partir dos números hexadecimais, filtramos as letras. Tudo o que nos resta são números. Precisamos verificar se esses números são primos:

hex |  filtered |  is prime? |  N =
1   >  1        >  false
1   >  1        >  false
2   >  2        >  true         1
3   >  3        >  true         2
5   >  5        >  true         3
8   >  8        >  false
D   >  0        >  false
15  >  15       >  false
22  >  22       >  false
37  >  37       >  true         4
59  >  59       >  true         5
90  >  90       >  false
E9  >  9        >  false
179 >  179      >  true         6
262 >  262      >  false
3DB >  3        >  true         7
63D >  63       >  false

Se o número filtrado é primo, chamamos isso de Fibohexaprime . Você pode ver que N = 7, o número de fibonacci relacionado é 987.

A tarefa é simples: quando receber uma entrada usando STDIN ou uma alternativa aceitável, escreva um programa ou uma função que produza a enésima Fibohexaprime usando STDOUT ou uma alternativa aceitável.

Casos de teste

Input - Output
1     - 2
2     - 3
3     - 5
4     - 55
5     - 89
6     - 377
7     - 987
8     - 28657
9     - 75025
10    - 121393
11    - 317811
12    - 5702887
13    - 9227465
14    - 39088169
15    - 102334155
16    - 32951280099
17    - 4052739537881
18    - 806515533049393
19    - 7540113804746346429

As regras:

  • Dado um número inteiro entre 1e 19(os valores acima 20excedem o valor máximo para um número inteiro assinado de 64 bits), imprima o valor correspondente.
  • Você pode escrever uma função ou um programa.
  • Isso é , então a submissão com a menor quantidade de bytes ganha!

Do jeito que está redigido, parece que as funções também devem ler de STDIN e gravar em STDOUT. Isso está correto? Normalmente, permitimos que as funções aceitem argumentos e retornem valores como conveniente.
Alex A.

2
@AlexA. Ambos são alternativas aceitáveis. A leitura do STDIN e o uso do STDOUT não são obrigatórios.
Adnan

Respostas:


4

Pitão, 27 bytes

Leu,eGsGbU2ye.fq1lPs-.HyZGQ

Demonstração

ycalcula o enésimo número de Fibonacci. Um .floop encontra o fibohexaprime de acordo com a entrada.


12

MATL , 28 bytes

Isso usa o MATL versão 1.0.0 , publicada em Esolangs em 12 de dezembro, anterior a esse desafio.

1Hi:"`tb+t16YAt58<)YtZp~]]1$

Exemplo

>> matl 1Hi:"`tb+t16YAt58<)YtZp~]]1$
> 10
121393

Explicação

O código é semelhante ao da resposta de Martin Büttner .

1           % number literal
H           % paste from clipboard H. Initial contents: 2
i:          % vector of equally spaced values from 1 to input value           
"           % for                      
  `         % do...while         
    t       % duplicate                           
    b       % bubble up element in stack          
    +       % addition 
    t       % duplicate                   
    16YA    % convert integer to string representation in base 16
    t       % duplicate             
    58      % number literal: first ASCII code after '9'           
    <       % is less than? (element-wise)    
    )       % reference () indexing with logical index from previous comparison
    Yt      % convert string to number 
    Zp      % true for prime numbers                                
    ~       % logical 'not'
  ]         % end                                                   
]           % end                                                   
1$          % input specification for final implicit display function

4
A primeira resposta MATL do mundo! Bom trabalho, Luis!
copo

1
Viva o MATL! Bem-vindo ao mundo do código de golfe!
rayryeng - Restabelece Monica

8

CJam, 28 bytes

TXri{{_@+_Gb{A<},Abmp!}g}*p;

Teste aqui.

Explicação

TX        e# Push 0 and 1 to initialise Fibonacci computation.
ri        e# Read input and convert to integer N.
{         e# Run this block N times...
  {       e#   While the condition on top of the stack is truthy...
    _@+   e#     Compute next Fibonacci number (dropping the second-to-last one).
    _Gb   e#     Duplicate and convert to base 16.
    {A<}, e#     Keep only digits less than 10.
    Ab    e#     Convert from base 10.
    mp!   e#     Check that it's not a prime.
  }g
}*
p;        e# Print the last number we found and discard the one before.

7

Perl 6 , 62 bytes

Meu primeiro passo de fazê-lo funcionar foi:

{(grep *[1].is-prime,map {$_,+[~] .base(16)~~m:g/\d/},(1,1,*+*...*))[$_-1;0]} # 77

Ao combinar o grepe map, posso remover 10 bytes

{(map {$_ if is-prime [~] .base(16)~~m:g/\d/},(1,1,*+*...*))[$_-1]} # 67

Se eu usar em grepvez de map, economizo mais 5 bytes:

{(grep {is-prime [~] .base(16)~~m:g/\d/},(1,1,*+*...*))[$_-1]} # 62

uso:

# give it a name
my &code = {...}

say code $_ for 1..^20;

2
3
5
55
89
377
987
28657
75025
121393
317811
5702887
9227465
39088169
102334155
32951280099
4052739537881
806515533049393
7540113804746346429

3

Mathematica 111 bytes

Ainda pode haver espaço para golfe adicional.

t=Table[Fibonacci@k,{k,1600}];f@n_:=PrimeQ@FromDigits[Select[n~IntegerDigits~16,#<10&]];
g@k_:=Select[t,f][[k]]

g[7]

987


g[19]

7540113804746346429


3

Julia, 123 bytes

n->(a=[];i=1;while endof(a)<n b=([1 1;1 0]^i)[1];(s=filter(isdigit,hex(b)))>""&&isprime(parse(s))&&push!(a,b);i+=1end;a[n])

Esta é uma função anônima que aceita um número inteiro e retorna um número inteiro. Para chamá-lo, dê um nome, por exemplo f=n->....

Ungolfed:

function f(n::Integer)
    # Initialize an array and an index
    a = []
    i = 1

    # Loop while we've generated fewer than n fibohexaprimes
    while endof(a) < n
        # Get the ith Fibonacci number
        b = ([1 1; 1 0]^i)[1]

        # Filter the hexadecimal representation to digits only
        s = filter(isdigit, hex(b))

        # If there are digits to parse, parse them into an
        # integer, check primality, and push the Fibonacci
        # number if prime
        s > "" && isprime(parse(s)) && push!(a, b)

        # Next
        i += 1
    end

    # Return the last generated
    return a[n]
end

3

GAP , 204 bytes

Essa resposta é bastante comum, exceto que o GAP é legal o suficiente para encontrar os próximos dois fibohexaprimes (e mais frio ainda, encontra-os em milissegundos com o código fornecido).

gap>f(20);                                                                    
31940434634990099905
gap> f(21);
12776523572924732586037033894655031898659556447352249
gap> f(22);
971183874599339129547649988289594072811608739584170445
gap> f(23);
1324695516964754142521850507284930515811378128425638237225
gap> f(24);
187341518601536966291015050946540312701895836604078191803255601777

Observe que f (24) está entre 2 ^ 216 e 2 ^ 217.

Aqui está o código:

f:=function(n)local c,i,x;c:=1;i:=0;while c<=n do x:=HexStringInt(Fibonacci(i));RemoveCharacters(x,"ABCDEFGHIJKLMNOPQRSTUVWXYZ");x:=Int(x);if IsPrime(x) then c:=c+1;fi;i:=i+1;od;Print(Fibonacci(i-1));end;

Provavelmente ainda há algum golfe que poderia ser feito. Eu acho que a implementação é bem direta.

Ungolfed:

f:=function(n)
    local counter,i,x;
    counter:=1;i:=0;
    while counter<=n do
        x:=HexStringInt(Fibonacci(i));
        RemoveCharacters(x,"ABCDEFGHIJKLMNOPQRSTUVWXYZ");
        x:=Int(x);
        if IsPrime(x) then
            counter:=counter+1;
        fi;
        i:=i+1;
    od;
    Print(Fibonacci(i-1));
end;

3

C, 186 183 bytes

#include<stddef.h>
size_t a,b,c,d,m,x;size_t F(n){a=0,b=1;while(n){x=b;b+=a;a=x;c=0,m=1;while(x)d=x%16,m*=d<10?c+=m*d,10:1,x/=16;d=c>1;x=2;while(x<c)if(c%x++==0)d=0;d&&--n;}return a;}

O teste de primalidade é muito ineficiente, de modo que a computação se esforçará um pouco n > 16e se tornará dolorosamente longa paran = 19 . No entanto, funciona e fornece os resultados esperados.

O código pressupõe que size_té um tipo de 64 bits, que é verdadeiro para Linux e Windows de 64 bits.


Bônus: infelizmente somos obrigados a usar tipos de 64 bits, o que leva a uma sobrecarga de 33 bytes. A versão seguinte funciona para n <= 15usar inte é de 150 bytes de comprimento:

a,b,c,d,m,x;F(n){a=0,b=1;while(n){x=b;b+=a;a=x;c=0,m=1;while(x)d=x%16,m*=d<10?c+=m*d,10:1,x/=16;d=c>1;x=2;while(x<c)if(c%x++==0)d=0;d&&--n;}return a;}

Teste principal:

#include <stdio.h>

int main() {
  printf("Input - Output\n");
  for (int i = 1; i < 20; ++i) {
    printf("%2d    - %ld\n", i, F(i));
  }
}

Você poderia economizar um pouco usando size_te soltando a inclusão? É específico da implementação, mas parece ser de 64 bits no Linux e no Windows gcc de 64 bits (e desde quando nos preocupamos com a portabilidade no codegolf?). (nota: %ldnão é de 64 bits no Windows de 64 bits; necessidades %lld)
Bob

@ Bob eu pensei sobre isso, mas size_tnão é um builtin, é definido em stddef.h(que por sua vez é direta ou indiretamente incluído por praticamente qualquer outro cabeçalho). De um jeito ou de outro, eu preciso de um #include. Eu ainda salvar 2 bytes por pode usar size_tem vez de uint64_t, embora :)
Stefano Sanfilippo

Também graças à lldpouco, eu não ter a oportunidade de testá-lo no Windows (mas portabilidade não importa, certo?)
Stefano Sanfilippo

Hum, deve ter vindo stdio.henquanto eu estava testando. De qualquer forma - você ainda pode salvar um casal incluindo, em math.hvez de stddef.h.
Bob

math.hnão faz o truque para mim (GCC 4.9 com GNU libc)
Stefano Sanfilippo

2

Python 2, 127 bytes

N=input();a,b=0,1
while N:a,b=b,a+b;t=int(''.join(c for c in hex(b)if ord(c)<65));N-=(t>1)*all(t%x for x in range(2,t))
print b

O algoritmo pode ser muito mais eficiente. Em particular, a verificação de primalidade (t>1)*all(t%x for x in range(2,t))verifica os fatores em potencial até o t-1momento, quando, na verdade, seria necessário apenas verificar o nível da raiz quadrada . Como rangearmazena uma lista inteira na memória no Python 2, isso leva a um MemoryErrorat N=17(na minha máquina usando as configurações padrão).


2

Ruby, 160 bytes

->i{t,o,k=[],[],0;f=->n{t[n]||=n<3?1:f[n-2]+f[n-1]};(r=('%x'%f[k]).scan(/\d/).join.to_i;(r>1&&(r==2||(2...r).none?{|j|r%j==0}))&&o<<r;k+=1)while !o[i-1];t[k-1]}

Ungolfed:

-> i {
  t, o, k = [], [], 0
  f = -> n {
    t[n] ||= n < 3 ? 1 : f[n-2] + f[n-1]
  }
  while !o[i-1] do
    r=('%x'%f[k]).scan(/\d/).join.to_i
    o << r if (r > 1 && (r == 2 || (2...r).none?{|j| r%j == 0 }))
    k+=1
  end
  t[k-1]
}

Uso:

# Assign the anonymous function to a variable
m = ->i{t,o,k=[],[],0;f=->n{t[n]||=n<3?1:f[n-2]+f[n-1]};(r=('%x'%f[k]).scan(/\d/).join.to_i;(r>1&&(r==2||(2...r).none?{|j|r%j==0}))&&o<<r;k+=1)while !o[i-1];t[k-1]}

m[2]
=> 3
m[19]
=> 7540113804746346429

2

R, 164 bytes

g=function(n){f=function(m)ifelse(m<3,1,f(m-1)+f(m-2));p=0;while(n){p=p+1;x=gsub("\\D","",sprintf("%x",f(p)));x[x==""]=1;y=1:x;if(sum(!tail(y,1)%%y)==2)n=n-1};f(p)}

Recuado, com novas linhas:

g=function(n){
    f = function(m)ifelse(m<3,1,f(m-1)+f(m-2)) #Fibonacci function
    p = 0
    while(n){
        p = p+1
        x = gsub("\\D","",sprintf("%x",f(p))) #To Hex, and get rid of non-digits
        x[x==""] = 1 #If x is empty string
        y = 1:x #Converts to integer(!) and save the 1-to-x sequence to a variable
        if(sum(!tail(y,1)%%y)==2) n = n-1 #If prime, decrements counter
        }
    f(p)
    }

Exemplos:

> g(1)
[1] 2
> g(5)
[1] 89
> g(10)
[1] 121393
> g(12)
[1] 5702887
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.