%I#17 2022年12月12日16:03:13
%S 1,1,2,2,3,2,3,3,4,2,4,3,4,4,3,4,5,5,5,5,4,5,5,5,5,5,5,5,5,6,6,5,
%T 6,5,6,5,1,6,6,6,6,6,1,7,6,6-6,7,6,7,6,16,7,16,8,6,7~6,6~8,8,7,9,
%U 7,9,7,10,7,10,7,10,10,7,10,7,11
%N a(N)=小于等于N且与N互素的正整数个数!!(n!!=A006882(n))。
%H Robert Israel,n的表,a(n)表示n=1..10000</a>
%F From _Robert Israel_,2022年12月8日:(开始)
%F如果n是奇数,a(n)=A070939(n)。
%F如果n>2是偶数,则a(n)=1+A056171(n)。(结束)
%p f:=进程(n)局部t;
%p如果n::奇数,则ilog2(n)+1
%p else 1+数字理论:-pi(n)-numtheory:-pi(n/2)
%功率因数
%p结束过程:
%p f(2):=1:
%p映射(f,[1..100]美元);#_罗伯特·伊斯雷尔,2022年12月8日
%t a[n_]:=模[{},co=0;对于[i=1,i<n+1,i++,如果[GCD[n!!,i]==1,co++]];co];表[a[n],{n,1,80}](*Stefan Steinerberger_,2007年6月5日*)
%t表[Total[Boole[CoprimQ[n!!,Range[n]]],{n,80}](*哈维·P·戴尔,2022年12月12日*)
%Y参考A006882、A056171、A070939。
%K nonn公司
%氧1,3
%A _罗伊查询,2007年6月3日
%E更多条款,来自斯特凡·斯坦纳伯格,2007年6月5日
|