%I#12 2014年10月19日11:48:10

%第14295827741196011260118457200112052124919259992859129947页，

%电话：3321133493603738149412274264935794530746099495952057，

%电话：523875431959107591975962967891709517347776175037811579204193607114889

%N素数p，这样p^4-p^3+1和p^4-p ^3-1也是素数。

%H K.D.Bajpai，n的表格，n=1..2918的a（n）</a>

%e 1429在序列中，因为1429，（1429^4-1429^3+1）和（1429*4-1429*3-1）都是素数。

%p KD:=proc（）局部a，b，d；a： =ithprime（n）；b： =a^4-a^3+1；d： =a^4-a^3-1；如果是isprime（b）和isprim（d），则返回（a）；fi；结束：seq（KD（），n=1..20000）；

%t选择[Prime[Range[3000]]，PrimeQ[#^4-#^3+1]&&PrimeQ[#^4-#^3-1]&]

%t c=0；a=2；Do[k=素数[n]；如果[PrimeQ[k^4-k^3+1]和&PrimeQ[k^4-k ^3-1]，c=c+1；打印[c，“”，k]]，{n，12000000}]；

%t pQ[n_]：=模[{c=n^4-n^3}，AllTrue[c+{1，-1}，PrimeQ]]；选择[Prime[Range[11000]]，pQ]（*程序使用Mathematica版本10中的AllTrue函数*）（*Harvey P.Dale_，2014年10月19日*）

%o（PARI）s=[]；对于素数（p=2120000，if（isprime（p^4-p^3+1）&isprime，s=concat（s，p））；2014年2月18日，科林·巴克

%Y参见A000040、A237639、A237641和A237642。

%K nonn公司

%O 1,1号机组

%A _K.D.Bajpai，2014年2月18日