%I#14 2019年8月5日03:03:35

%S 0,0,1,1,2,0,3,2,3,2,1,3,3,3,1,3,4,4,1,4,12,2,5,8,5,1,5,3,6,6,7,6,6，

%温度4,2,4,3,6,11,6,4,3,17,6,8,3,2,10,9,6,12,8,9,5,2,5,2,3,13,5,14,8，

%U 12,7,8,9,6,13,9,4,10,3,13,12,4,8,4

%N素数p<N，使得p*（N-p）是一个本原根模素数（N）。

%C猜想：对于所有n>6，a（n）>0。

%C我们已经验证了所有n=7。。。，2*10^5.

%孙志伟，n的表，n的a（n）=1..10000</a>

%孙志伟，<a href=“http://arxiv.org/abs/1405.0290“>关于模素数本原根的新观察</a>，arXiv:1405.0290[math.NT]，2014。

%e a（18）=1，因为17是素数，17*（18-17）=17 a本原根模素数（18）=61。

%e a（20）=1，因为11是素数，11*（20-11）=99是一个本原根模素数（20）=71。

%e a（26）=1，因为2是素数，而2*（26-2）=48是原根模素数（26）=101。

%e a（27）=1，因为17是素数，17*（27-17）=170 a本原根模素数（27）=103。

%t dv[n_]：=除数[n]

%t Do[m=0；Do[Do[If[Mod[（素数[k]*（n-素数[k]））^（部分[dv[Prime[n]-1]，i]），素数[n]==1，转到[aa]]，{i，1，长度[dv[素数[n]-1]-1}]；m=m+1；标签[aa]；继续，{k，1，PrimePi[n-1]}]；

%t打印[n，“”，m]；继续，{n，1,80}]

%Y参考A000040、A000720、A243164。

%K nonn公司

%O 1,5型

%A _孙志伟，2014年6月4日