%I#15 2020年10月12日11:14:27

%S 1,1,1,2,1,1,1,3,4,2,1,6,5,2,9,11,12,5,7,9,8,8,17,12,11,16,12,23,20，

%电话：16,17,17,23,17,26,18,19,25,26,32,38,21,21,18,27,25,24,27,52,30,44,33，

%U 19，44，54，45，57，14，29，27，39，58，28，41，39，62，26，25，69，48，51，61，44，47，37，63，77，55，55，41

%N本原根g mod prime（N）的个数，其中g^x==x（mod p）与2<=x<=prime（N）-2没有解。

%H Robert Israel，n的表，n=1..1000的a（n）</a>

%H M.Levin，C.Pomerance，K.Soundararajan，<a href=“https://doi.org/10.1007/978-3-642-14518-6_5“>离散对数的不动点</a>。摘自：Hanrot G.，Morain F.，ThoméE.（eds）算法数论。ANTS 2010。计算机科学讲义，第6197卷。施普林格，柏林，海德堡（2010）。

%H数学溢出，<a href=“https://mathoverflow.net/q/269368“>g^x的不动点（模a素数）</a>

%对于n>=2，F a（n）=A008330（n）-A174407（n）。

%p f：=进程（n）局部p，r，S，r，x；

%p p：=i素数（n）；

%p r：=数量理论：-primroot（p）；

%pS:=选择（t->igcd（t，p-1）=1，{$1..p-1}）；

%p R：=地图（s->R&^s mod p，s）；

%x的p从2到p-2 do

%p R:=删除（t->（t&^x-x mod p=0），R）；

%p od；

%p nops（R）；

%p端程序；

%p映射（f，[1..100]美元）；

%t连接[{1}，表[p=素数[n]；EulerPhi[EulerPhi[p]]-长度[Select[PrimitiveRootList[p]，MemberQ[PowerMod[#，Range[p-1]，p]-Range[p1]，0]&]]，{n，2，100}]]（*Jean-François Alcover_，2020年10月11日，A174407中的_T.D.Noe_之后*）

%Y参考A008330，A174407。

%K nonn公司

%O 1,5型

%以色列罗贝尔特，2017年5月10日