Nが現れる素数をJuliaでやってみた

この記事は@kjunichiの2017年パーソナルアドベントカレンダーの1日目の記事です。

integers.hatenablog.com

これはすごい、どうやってるんだと、パッと見は思った。

swdrsker.hatenablog.com

を読んで、なんとなくJulia言語で試してたくなり、やってみた。

5が現れる素数を昼休みに恥ずかしいバグに悩まされながらも計算で来た #julialng / “this five is prime number · GitHub” https://t.co/gJKRpEViM2
— kjunichi (@kjunichi) 2017年11月30日

Pkg.add("Primes")
using Primes
cnt=0
isPrime=false
for a in 1:9
    #println("a =", a)
    for b in 0:9      
        for c in 0:9
            for d in 0:9    
                for ee in 0:9
                    for f in 0:9                                   
                        for g in 0:9
                            for h in 0:9                                               
                                for i in [1,3,5,7,9]
                                    str5 =
                                    "ABC000000000000000"*
                                    "005555555555555500"*
                                    "005555555555555500"*
                                    "005555000000000000"*
                                    "005555000000000000"*
                                    "000555555555500000"*
                                    "000555555555555000"*
                                    "000000000000555500"*
                                    "000000000000555500"*
                                    "005555555555555000"*
                                    "005555555555000000"*
                                    "DEF000000000000GHI"
                                    str5=replace(str5,"A",a)
                                    str5 = replace(str5, "B",b)
                                    str5 = replace(str5, "C",c)
                                    str5 = replace(str5, "D",d)
                                    str5 = replace(str5, "E",ee)
                                    str5 = replace(str5, "F",f)
                                    str5 = replace(str5, "G",g)
                                    str5 = replace(str5, "H",h)
                                    str5 = replace(str5, "A",a)
                                    str5 = replace(str5, "I",i)
                                    if isprime(parse(BigInt,str5))
                                        #println(a,b,c,d,ee,f,g,h,i)
                                        println(str5)
                                        isPrime=true
                                        break
                                    end
                                    cnt=cnt+1
                                end
                                if isPrime 
                                    break
                                end
                            end
                            if isPrime 
                                break
                            end
                        end
                        if isPrime 
                            break
                        end
                    end
                    if isPrime 
                         break
                    end
                end
                if isPrime 
                    break
                end
            end
            if isPrime 
                break
            end
        end
        if isPrime 
            break
        end
    end
    if isPrime 
        break
    end
end
#println(cnt)

100000000000000000
005555555555555500
005555555555555500
005555000000000000
005555000000000000
000555555555500000
000555555555555000
000000000000555500
000000000000555500
005555555555555000
005555555555000000
000000000000000437