%I#29 2019年12月7日12:18:22
%S 1,1,2,1,4,1,4,4,3,2,8,3,7,7,9,2,8,5,10,10,8,6,19,6,12,9,8,22,9,12,
%电话:12,15,10,31,9,11,14,24,13,23,9,24,17,16,10,35,15,23,25,20,12,40,17,
%U 34,21,18,14,37,17,24,25,41,20,39,14,31,34,33,18,42,16,32,37,41,18,44,25
%N抗-phi(N)。
%C anti-phi(n)=不可被n的任何反除数整除的整数数。
%C这个序列的旧定义是:anti-phi(n)=与n的反除数互素的整数个数<=n,但这与条目不匹配。
%C关于反除数的定义,请参见A066272。
%H Nathaniel Johnston,n的表,n的a(n)表示n=2..10000
%H Jon Perry,<a href=“http://www.users.globalnet.co.uk/~perry/methemats/antidivisorother2.htm“>反钓鱼功能
%H Jon Perry,反歧视软件[Cached copy]
%H Jon Perry,反驱动程序:更多反驱动程序
%e 10有反除数3、4、7。不可被3、4、7和小于10中的任何一个整除的数字是1、2、5。因此,抗φ(10)=3。
%p#需要A066272中实现的antidivisors()
%p A066452:=proc(n)local ad,isad,j,k,totad:ad:=antidivisors(n):totad:=0:对于j从1到n-1,do isad:=1:对于k从1到nops(ad)do如果(j mod ad[k]=0),则isad:=0:break:fi:od:totad:=totad+isad:od:return totad:end:
%p序列(A066452(n),n=2..50);#_纳撒尼尔·约翰斯顿,2011年4月20日
%o(Python)
%o定义A066452(n):
%o。。。。return len([x表示范围(1,n)中的x,如果全部([x%d表示范围(2,n)中的d,如果(n%d)和(2*n)%d表示[d-1,0,1]])])#_Chai Wah Wu_,2014年8月7日
%o(PARI)antidiv(n)={my(v=[]);对于(k=2,n-1,如果(abs((n%k)-k/2)<1,v=concat(v,k)););v;}
%o a(n)={my(vad=antidiv(n));my(nbad=0);对于(j=1,n-1,isad=1;对于(k=1,#vad,if((j%vad[k])==0,isad=0;break););nbad+=isad;)
%Y参考A058838、A066241。
%K nonn,简单
%氧2,3
%乔恩·佩里,2001年12月29日
%E更好的定义和更多术语摘自Diana L.Mecum_,2007年7月1日
|