找出第n个六六元


23

这次的挑战是找到第nFibohexaprimeFibohexaprime的定义如下:

我们首先观察一个带有斐波那契数的列表:

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

之后,我们将数字转换为十六进制:

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

从十六进制数字中,我们过滤出字母。我们只剩下数字。我们需要检查以下数字是否为质数:

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

如果过滤后的数字是素数,我们称其为Fibohexaprime。您可以看到N = 7的相关斐波那契数是987。

任务很简单,当使用STDIN或可接受的替代方法进行输入时,编写程序或函数,使用STDOUT或可接受的替代方法输出第n个Fibohexaprime。

测试用例

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

规则:

  • 给定一个介于1和之间的整数19(上述值20超过64位有符号整数的最大值),输出相应的值。
  • 您可以编写函数或程序。
  • 这是,因此以最少的字节提交为准!

听起来,函数听起来也必须同时从STDIN读取并写入STDOUT。那是对的吗?通常,我们允许函数方便地接受参数和返回值。
Alex A.

2
@AlexA。两者都是可接受的替代方案。从STDIN读取并使用STDOUT不是强制性的。
阿德南(Adnan)2015年

Answers:


4

Pyth,27个字节

Leu,eGsGbU2ye.fq1lPs-.HyZGQ

示范

y计算第n个斐波那契数。甲.f环根据输入发现fibohexaprime。


12

MATL,28字节

它使用的MATL版本1.0.0早于此挑战,于12月12日在Esolangs中发布。

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

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

说明

该代码类似于MartinBü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
世界上第一个MATL答案!干得好,路易斯!
烧杯

1
祝您万事如意!欢迎来到代码高尔夫世界!
rayryeng-恢复莫妮卡2015年

8

CJam,28个字节

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

在这里测试。

说明

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个字节

我刚开始工作的第一步是:

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

通过结合grepmap,我可以删除10个字节

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

如果我使用grep而不是map,则会再节省5个字节:

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

用法:

# 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字节

仍有可能有其他打高尔夫球的空间。

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

朱莉娅123字节

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])

这是一个接受整数并返回整数的匿名函数。要给它起个名字,例如f=n->...

取消高尔夫:

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字节

这个答案几乎没有什么特别的,除了GAP非常酷,可以找到下几个fibohexaprimes(更酷的是,它用给定的代码在几毫秒内找到了它们)。

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

请注意,f(24)在2 ^ 216和2 ^ 217之间。

这是代码:

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;

可能还有一些打高尔夫球的事情。我认为实施非常简单。

取消高尔夫:

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字节

#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;}

素数测试效率很低,因此计算会费时费力,n > 16并且会很痛苦n = 19。但是,它可以正常工作并给出预期的结果。

该代码假定它size_t是64位类型,对于64位Linux和Windows均适用。


奖励:不幸的是,我们必须使用64位类型,这导致33个字节的开销。以下版本适用于n <= 15使用int,长度为150字节:

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;}

测试主体:

#include <stdio.h>

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

您能否通过使用size_t和删除包含项来节省一点时间?它是特定于实现的,但是在64位Linux和Windows gcc中似乎都是64位的(而且从什么时候开始我们关心代码高尔夫的可移植性?)。(旁注:%ld不是在64位Windows中为64位;需要%lld
Bob

我考虑过@Bob,但size_t它不是内置的,它是在stddef.h(它实际上直接或间接地包含在任何其他标头中)定义的。一种方式,我需要一个#include。我仍然可以使用size_t而不是来保存2个字节uint64_t:)
Stefano Sanfilippo

同样感谢您lld,我没有机会在Windows上进行测试(但可移植性并不重要,对吗?)
Stefano Sanfilippo

嗯,一定是stdio.h我测试时产生的。无论如何,您仍然可以通过包含math.h而不是来节省一对夫妇stddef.h
鲍勃

math.h对我没有帮助(使用GNU libc的GCC 4.9)
Stefano Sanfilippo

2

Python 2,127字节

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

该算法可能更加有效。尤其是,素数检查会一直(t>1)*all(t%x for x in range(2,t))检查所有潜在因素t-1,而实际上只需要检查平方根底数即可。由于range将整个列表存储在Python 2中的内存中,因此会导致MemoryErrorat N=17(在我的机器上使用默认设置)。


2

Ruby,160个字节

->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]}

取消高尔夫:

-> 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]
}

用法:

# 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字节

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)}

缩进,换行:

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)
    }

例子:

> g(1)
[1] 2
> g(5)
[1] 89
> g(10)
[1] 121393
> g(12)
[1] 5702887
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.