%I#17 2021年5月7日09:10:37
%S 0,1,2,1,4,1,6,1,2,1,4,1,12,1,2,1,41,18,1,2,4,1,4,1,2,3,1,1,
%温度2,1,4,1,6,1,2,1,4,1,12,1,2,4,1,18,1,2,1,4,1,24,1,2,2,4,1,60,1,2,1,
%U 4,1,6,1,2,1,4,1,12,1,2,1,41,18,1,2,1,4,1,24,1,2,4,1,90,1,2,2,1,4,6,12,4,12,1,1,2,4,18,1
%N a(N)=N模A002110(A257993(N)),a(0)=0。
%H Antti Karttunen,n表,n=0..2310的a(n)</a>
%F a(0)=0,对于n>=1,a(n)=n模A002110(A257993(n))。
%F或a(n)=A276088(n)*A002110(A276084(n))。
%F其他身份。对于所有n>=0:
%F a(n)=n-A276093(n)。
%t{0}~Join~表[k=1;While[!互质Q[Prime@k,n],k++];Mod[n,Product[Prime@i,{i,k}]],{n,79}](*_Michael De Vlieger_,2017年6月22日*)
%o(方案,两个版本)
%o(定义(A276094 n)(如果(零?n)n(让回路((n n)(i 1)(pr 1)))(让*((p(A000040 i))(d(模n p)))
%o(定义(A276094 n)(如果(零?n)n(模n(A002110(A257993 n))))
%o(Python)
%o从sympy导入nextprime,primepi,primarial
%o定义a053669(n):
%o p=2
%o为True时:
%o如果n%p:返回p
%o其他:p=下一素数(p)
%o定义a257993(n):返回primepi(a053669(n))
%o def a002110(n):如果n<1,则返回1
%o定义a(n):如果n==0,则返回0,否则返回n%a002110(a257993(n))
%o打印([a(n)代表范围(101)内的n)]#_Indranil Ghosh,2017年6月22日
%Y参见A000040、A002110、A257993、A276084、A2760088、A276093。
%K nonn,基础
%0、3
%A _Antti Karttunen,2016年8月22日
|