%I#9 2014年7月12日22:43:37
%S 2,13,10,17,6,37,12,13,16,27,24,71,16,31,64,43,18,43,26,23,32,29,24,
%电话:79,32,53,34,61,92,47,40,33,34,57,36,47,40.53,40,79,44,43,68,91,68,57,
%U 66,61,60,53,58,83,60,91,94,61,82,61,70101,82,71,68145,82,67,76,69100
%N使k^8+N^8为素数的最小数k>N。
%如果n在A153504中,则C a(n)=n+1。
%H Jens Kruse Andersen,n的表,n的a(n)=1..10000</a>
%e 13^8+14^8=2291519777不是素数,13^8+1 5^8=3378621346不是素数。13^8+16^8=5110698017是质数。因此a(13)=16。
%o(Python)
%o导入交响乐
%o来自sympy import isprime
%o定义a(n):
%o。。对于范围(n+1,10**4)中的k:
%o。。。。如果isprime(k**8+n**8):
%o。。。。。。返回k
%o n=1
%o当n<100时:
%o。。打印(a(n),结束=',')
%o。。n+=1
%o(PARI)a(n)=对于(k=n+1,10^4,如果(i素数(k^8+n^8),返回(k))
%o n=1;而(n<100,打印1(a(n),“,”);n++)
%Y参见A158979、A089489、A242553。
%K nonn公司
%O 1,1号机组
%2014年7月8日,A _Derek Orr_
|