观察s(n)=60的前几个例子，除以2的序列都是 ***32**, ***128***, ***256***这种结构，到底是哪种，取决于n的大小，2^x总是小于n/2.可以把小于n的公比为2的等比数列都列举出来，去凑能用+（n-1)连上的结构，如果最后总长度=60就成了

因为满足s(n)=x的n上限是2^x,所以x很小的可以求出全部，x很大就无能为力了，有的x对应的n还很多 >>> def gen(x): ... n=x ... c=0 ... while n>1: ... if n%2==0:n//=2 ;c+=1 #;print(n) ... else: n+=(x-1) ... return c ... >>> gen(18) 8 >>> [i for i in range(2,300,2) if gen(i)==8] [18, 52, 86, 256] >>> [i for i in range(2,3000,2) if gen(i)==60] [62, 144, 176, 184, 226, 288, 306, 326, 370, 386, 404, 428, 430, 496, 526, 534, 550, 672, 716, 756, 794, 862, 916, 976, 1002, 1086, 1156, 1210, 1282, 1288, 1322, 1396, 1436, 1526, 1576, 1600, 1656, 1846, 1892, 1926, 1964, 2014, 2016, 2136, 2146, 2266, 2276, 2380, 2476, 2502, 2584, 2666, 2746, 2822, 2926] >>> len([i for i in range(2,3000,2) if gen(i)==60]) 55

s(n)==n-2的特点。每次除以2的商都不同，如果相同那个就死循环了 >>> gens(20) 10 5 12 6 3 11 15 17 18 9 14 7 13 16 8 4 2 1

有可能奇数的下一步仍然是奇数 >>> gens(18) 9 13 15 16 8 4 2 1

实际上就是求不定方程组的整数解 2^x*奇数=n 奇数+n-1=2^y 其中第1、2个方程可以合并成一个，消去奇数 n/(2** x)+n-1==2**y x+y=洗牌次数