%I#28 2020年5月26日06:27:25

%S 1,1,2,1,1,2,3,1,2,1,1,1,2,1,2,1,3,4,1,1,2,1,1,1,3,1,1,2,1,3,1,4,5,1,2,1，

%温度3,2,1,1,2,1,1,1,3,1,1,4,1,2,3,1,2,1,2,1,2,2,1,1,4,1,5,6,1,1,2,1,1,1,1，

%U 2,3,1,2,1,1,1,1,4,1,3,2,1,1,1,2,1,1,3,1,1,1,1,1,1,1,4,3,1,1,5,1,2,3,2,1,2,1,2,1,1,4,4

%N a（N）是2^k-1除以N的最大整数k。

%C A161788（n）=2 ^a（n）-1。a（A161790（n））=1。

%C猜想：对于n>0和m>0，gcd（n，m）=a（2^n+2^m-2）_Velin Yanev，2017年8月24日

%H Robert Israel，n的表，a（n）表示n=1..10000</a>

%p A161789：=从ilog2（n+1）到0的k的proc（n）乘以-1 do，如果n mod（2^k-1）=0，则返回（k）；fi；od：结束：序列（A161789（n），n=1..120）；#_R.J.Mathar，2009年6月27日

%p#备选方案：

%p N:=200:#对于a（1）。。a（否）

%p V：=矢量（N，1）：

%p表示k从2到ilog2（N）do

%p t：=2^k-1；

%p V[[序列（i，i=t..N，t）]]：=k

%日期：

%p转换（V，列表）；#_罗伯特·伊斯雷尔，2020年5月12日

%tkn[n_]：=模[{k=Floor[Log[2，n]]+1}，而[！可除[n，2^k-1]，k--]；k] ；阵列[kn，110]（*哈维·P·戴尔，2012年3月26日*）

%o（PARI）a（n）=forstep（k=logint（n+1,2），1，-1，if（n%（2^k-1）==0，return（k）））\\_Charles R Greathouse IV_，2017年8月25日

%Y参考A000225、A161788、A161790。

%K nonn，简单

%氧1,3

%A _罗伊查询，2009年6月19日

%E由R.J.Mathar_于2009年6月27日延期