题目描述

小易总是感觉饥饿,所以作为章鱼的小易经常出去寻找贝壳吃。最开始小易在一个初始位置x_0。对于小易所处的当前位置x,他只能通过神秘的力量移动到 4 * x + 3或者8 * x + 7。因为使用神秘力量要耗费太多体力,所以它只能使用神秘力量最多100,000次。贝壳总生长在能被1,000,000,007整除的位置(比如:位置0,位置1,000,000,007,位置2,000,000,014等)。小易需要你帮忙计算最少需要使用多少次神秘力量就能吃到贝壳。

输入描述:

输入一个初始位置x_0,范围在1到1,000,000,006

输出描述:

输出小易最少需要使用神秘力量的次数,如果使用次数使用完还没找到贝壳,则输出-1

示例1

输入

125000000

输出

1

分析:

  这道题我们只能把每步都分为两种情况,使用神秘力量1(4 * x + 3)和使用神秘力量2(8 * x + 7)。从出发点开始枚举,使用广度优先遍历算法(BFS)。由于贝壳出现在能被1,000,000,007整除的位置,所以我们只需要考虑%1000000007后的结果。我们要记录初次到达某个位置时使用了几次神秘力量。

第一种方法:

from collections import deque
mod = 1e9+7
n = int(raw_input().strip())
currentPos = n%mod
power = {}
power[currentPos] = 0
d = deque()
d.append(currentPos)
flag = False
while len(d):
currentPos = d.popleft()
if power[currentPos] > 100000:
break
if currentPos == 0:
flag = True
break
nextPos = (4*currentPos+3)%mod
if nextPos not in power:
power[nextPos] = power[currentPos]+1
d.append(nextPos)
nextPos = (8*currentPos+7)%mod
if nextPos not in power:
power[nextPos] = power[currentPos]+1
d.append(nextPos)
if flag:
print(power[currentPos])
else:
print(-1)

第二种方法:

观察变换形式,并做变形:

4x+3=4(x+1)-1

8x+7=8(x+1)-1

如果多层嵌套呢?

y=4x+3

8y+7=8((4(x+1)-1)+1)-1=8(4(x+1))-1=32(x+1)-1

如果你多枚举一些,就会发现,能变换出的数的形式都是:

a(x+1)-1,其中a是2的>=2的幂次数(4、8、16、32、64、……)

我们可以利用这个特点

考虑直接枚举那个a,从2^2一直到……等等,最大是2的多少次?

答:直接考虑最大情况,每次变换都选择8x+7那种,也就是,每次a乘上8,也就是说,最坏是(2^3)^100000=2^300000次

所以,枚举a,从2^2次,一直到2^300000次

然后,对每个a检查一下,乘起来结果%1e9+7是不是0,如果是0,说明100000次之内有解

——问:那最小要执行几次变换?

答:我们直接贪心,尽量让a乘8(乘2次8和乘3次4一样大,当然是乘8越多,变换次数越少)

——问:如果我发现a==2^5或a==2^4的时候满足要求,但是5和4才不能表示成3的倍数,怎么办?

答:别忘了你手上还有4x+3的变换(就是a乘4的变换)

对5这种情况,除以3余2,那刚好,用一次乘4的变换就行了

对4这种情况,除以3余1,我们考虑,消去一个乘8的变换,用2个乘4的变换代替并补足。

n = int(raw_input().strip())
mod = int(1e9+7)
ans = -1
time = 4
for i in range(1,300001):
x = (n*time+time-1)%mod
if x == 0:
ans = (i+1)/3
if (i+1)%3:
ans += 1
break
time = (time*2)%mod
print(ans)

第二种方法要比第一种方法高效一点

参考博客:

http://blog.csdn.net/fcxxzux/article/details/52138964#t0

05-04 03:22