%I#13 2015年8月10日16:35:19

%S 1,0,21,1,45,79120305484697404186135460361227579453051，

%电话：35524989501186265717816403919280351927218127136972621，

%电话：101393476108788608226092102824550190311128523426277933427911543

%N a（N）=（2^p+1）^（p-1）模p^2，其中p是素数（N）。

%C A260507给出了a（n）==1的素数。

%H Felix Fröhlich，n的表，n=1..10000时的a（n）</a>

%F a（n）=A098640（n）^（A000040（n）-1）模数A000040。

%tf[n_]：=块[{p=素数@n}，PowerMod[2^p+1，p-1，p^2]；数组[f，46]（*_Robert G.Wilson v_，2015年7月29日*）

%o（PARI）a（n）=升力（Mod（2^prime（n）+1，prime（n^2）^（prime（n-1）））

%Y参考A000040、A098640、A260507。

%K nonn公司

%氧1,3

%2015年7月28日，A_Felix Fröhlich