%I#13 2017年3月9日09:58:30

%S 1,1,1,1,1,1,1,7,1,9,1,11,13,1,5,17,1,19,1,21,1,23,1,25,13,9,1,29，

%T 1,31,17,33,1,7,1,37,19,13,1,41,1,43,23,45,1,47,1,1,25,1,53,1,5,29，

%U 57,1,59,1,61,31,7,1,1,67,35,69,1,71,1,73,37,25,1,77,1,79,41,81,1,83,1,85,43,29,1,89,1,91,47,93,1,19,1,97,49,33,1

%N a（N）=A250469（N）/A280702（N）=A250468（N）/gcd（A003961（N），A250466（N））。

%C注：尽管A280703（352）=3，但a（352）=1为A003961（352）=3159=3^5*13，而A250469（352）=1053=3^4*13。（因此，A266645（352）=176=352/2。）问题：对于某些整数k>1，A003961（n）=k*A250469（n）是否有更多的n？参见A280703中的评论。

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

%F a（n）=A250469（n）/A280702（n）=A250468（n）/gcd（A003961（n），A250466（n））。

%F A280701（n）=n-a（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，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（方案）（定义（A280704 n）（/（A250469 n）（A280702 n）））

%Y参考A003961、A250469、A266645、A280701、A280702、A280703。

%K nonn公司

%O 1.8个

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