%I#9 2017年3月9日09:57:59

%S 0,1,2,3,4,5,6,1,8,1,10,12,1,14,11,16,18,1,20,12,2,1,24,14,19，

%电话28,1,30,1,16,1,34,29,36,1,20,27,40,1,42,1,22,1,46,1,48,49,26,51,52,1，

%U 54,51,28,1,58,1,60,1,32,57,64,65,66,1,34,1,70,1,72,1,38,51,76,1,78,1,40,1,82,1,84,1,44,59,88,1,90,1,46

%N a（N）=N-A280704（N）。

%C问题：所有术语都是非负的吗？它们发生在哪里？

%H Antti Karttunen，n的表，n=1..13098的a（n）</a>

%F a（n）=n-A280704（n）=n-（A250469（n）/gcd（A003961（n），A250468（n）））。

%t f[n_]：=f[n]=其中[n==1，1，PrimeQ@n，NextPrime@n，True，Times@@Replace[FactorInteger[n]，{p_，e_}:>f[p]^e，1]]；g[n_]：=如果[n==1，0，PrimePi@FactorInteger[n][[1，1]]]；函数[s，MapIndexed[Function[t，First@#2-t/GCD[t，f@First@@#2]][Lookup[s，g[First@#2]+1][[#1]]-Boole[First@#2==1]]&，#]&@Map[Position[Lookup[s，g@#]，#][1，1]&，Range@120]]@PositionIndex@Array[g，10^4]（*Michael De Vlieger_，2017年3月8日，版本10*）

%o（方案）（定义（A280701 n）（-n（A28070 4 n）））

%Y参考A003961、A250469、A280702和A280704。

%K nonn公司

%氧1,3

%2017年3月8日，安蒂·卡图内