%I#24 2021年2月24日03:51:34
%S 0,0,0,1,1,2,2,2,2,3,3,33,3,1,3,4,4,4,5,5,5,5,5,1,5,5,
%T 5,5,5,1,5,5,5,5和5,6,6,6,7,7,7,7,7,1,7,8,8,8,8,9,9,9,
%U 9,9,9,10,10,10,10,10,11,11,11,11,1,11,1,11,11,11,11,11,11,11,11,12
%N前N个无平方数中的计数因子7。
%H Reinhard Zumkeller,n的表,n=1..10000的a(n)</a>
%如果A005117(n)在A008589中,则F a(n)-a(n-1)=1,否则为0_罗伯特·伊斯雷尔(Robert Israel),2015年8月23日
%传真(n)~n/8.-_Amiram Eldar,2021年2月24日
%e前10个无平方数是:1,2,3,5,6=2*3,7,10=2*5,11,13和14=2*7:7和14可以被7整除,因此a(10)=2。
%p N:=1000:#使用无平方数<=N
%p M:=map(proc(t)if numtheory:-issqrfree(t)then if t mod 7=0 then 1 else 0 fi end proc,[$1..N]):
%p列表工具:-部分和(M);#_罗伯特·伊斯雷尔(Robert Israel),2015年8月23日
%t使用[{sf=Select[Range[200],SquareFreeQ]},Accumulate[If[Divisible[#,7],1,0]&/@sf]](*_哈维·P·戴尔,2013年3月21日*)
%o(哈斯克尔)
%o a072750 n=a072750_列表!!(n-1)
%o a072750_list=scanl1(+)$map((0^)。(`mod`7))005117_list
%o--_Reinhard Zumkeller,2013年3月25日
%o(PARI)
%o n=94;k=0;bag=列表();a=矢量(n);
%o直到(n==0,k++;如果(issquarefree(k),listput(bag,k));n——);
%o表示(i=2,#bag,a[i]=a[i-1]+(bag[i]%7==0);
%o打印(a);\\_Gheorghe Coserea,2015年8月23日
%Y参见A005117、A008589、A072747、A07274、A0727.49、A072 751。
%K nonn公司
%O 1,10号
%A _Reinhard Zumkeller,2002年7月8日
%E姓名由_Gheorghe Coserea_澄清,2015年8月23日
|