%I#21 2015年7月27日05:46:44

%序号1,2,7,9,10,13,17,19,25,26,34,37,41,43,49,50,62,67,73,74,82,87,89,94，

%电话：97,99105106109113121122129130133137139145146157，

%U 162167173178181185187193194206214217218223237241249电话

%N a（1）=1，a（N-1）是平方模a（N），a（N。

%C a（n）是大于a（n-1_R.J.Mathar，2015年7月27日

%H T.D.Noe，n表，n=1..1000的a（n）</a>

%e对于n=3，我们有一个（2）=2。2不是二次剩余mod 3，因为二次剩余mod 3是{0,1}，参见A011655。2不是二次剩余mod 4，因为二次剩余mod 4是{0,1}，参见A000035。2不是二次剩余mod 5，因为二次剩余mod 5是{0,1,4}，参见A070430。2不是二次剩余mod 6，因为二次剩余mod 6是{0,1,3,4}，参见A070431。2是二次剩余mod 7，因为二次剩余mod 7是{0,1,2,4}，参见A053879。因此a（3）=7_R.J.Mathar，2015年7月27日

%p A034796：=程序（n）

%p选项记忆；

%p如果n=1，则

%第1页；

%p其他

%p表示来自procname（n-1）+1 do的a

%p如果numtheory[quadres]（进程名（n-1），a）=1，则

%p返回a；

%p end if；

%p端do：

%p end if；

%p结束程序：#_R.J.Mathar_，2015年7月27日

%t剩余Q[n_，k_]：=长度[Select[Range[Floor[k/2]]^2，Mod[#，k]==n&，1]]==1；a[1]=1；a[n_]：=a[n]=对于[k=a[n-1]+1，真，k++，如果[resideQ[a[n-1'，k]，返回[k]]]；表[a[n]，{n，1，60}]（*Jean-François Alcover_，2013年8月13日*）

%K诺恩，不错

%O 1,2号机组

%A·热心的W·威尔逊_

%E澄清定义，Joerg Arndt_，2013年8月14日