project_euler.problem_686.sol1

Project Euler Problem 686: https://projecteuler.net/problem=686

2^7 = 128 is the first power of two whose leading digits are “12”. The next power of two whose leading digits are “12” is 2^80.

Define p(L,n) to be the nth-smallest value of j such that the base 10 representation of 2^j begins with the digits of L.

So p(12, 1) = 7 and p(12, 2) = 80.

You are given that p(123, 45) = 12710.

Find p(123, 678910).

Functions

log_difference(→ float)

This function returns the decimal value of a number multiplied with log(2)

solution(→ int)

This function calculates the power of two which is nth (n = number)

Module Contents

project_euler.problem_686.sol1.log_difference(number: int) float

This function returns the decimal value of a number multiplied with log(2) Since the problem is on powers of two, finding the powers of two with large exponents is time consuming. Hence we use log to reduce compute time.

We can find out that the first power of 2 with starting digits 123 is 90. Computing 2^90 is time consuming. Hence we find log(2^90) = 90*log(2) = 27.092699609758302 But we require only the decimal part to determine whether the power starts with 123. So we just return the decimal part of the log product. Therefore we return 0.092699609758302

>>> log_difference(90)
0.092699609758302
>>> log_difference(379)
0.090368356648852
project_euler.problem_686.sol1.solution(number: int = 678910) int

This function calculates the power of two which is nth (n = number) smallest value of power of 2 such that the starting digits of the 2^power is 123.

For example the powers of 2 for which starting digits is 123 are: 90, 379, 575, 864, 1060, 1545, 1741, 2030, 2226, 2515 and so on. 90 is the first power of 2 whose starting digits are 123, 379 is second power of 2 whose starting digits are 123, and so on.

So if number = 10, then solution returns 2515 as we observe from above series.

We will define a lowerbound and upperbound. lowerbound = log(1.23), upperbound = log(1.24) because we need to find the powers that yield 123 as starting digits.

log(1.23) = 0.08990511143939792, log(1,24) = 0.09342168516223506. We use 1.23 and not 12.3 or 123, because log(1.23) yields only decimal value which is less than 1. log(12.3) will be same decimal value but 1 added to it which is log(12.3) = 1.093421685162235. We observe that decimal value remains same no matter 1.23 or 12.3 Since we use the function log_difference(), which returns the value that is only decimal part, using 1.23 is logical.

If we see, 90*log(2) = 27.092699609758302, decimal part = 0.092699609758302, which is inside the range of lowerbound and upperbound.

If we compute the difference between all the powers which lead to 123 starting digits is as follows:

379 - 90 = 289 575 - 379 = 196 864 - 575 = 289 1060 - 864 = 196

We see a pattern here. The difference is either 196 or 289 = 196 + 93.

Hence to optimize the algorithm we will increment by 196 or 93 depending upon the log_difference() value.

Let’s take for example 90. Since 90 is the first power leading to staring digits as 123, we will increment iterator by 196. Because the difference between any two powers leading to 123 as staring digits is greater than or equal to 196. After incrementing by 196 we get 286.

log_difference(286) = 0.09457875989861 which is greater than upperbound. The next power is 379, and we need to add 93 to get there. The iterator will now become 379, which is the next power leading to 123 as starting digits.

Let’s take 1060. We increment by 196, we get 1256. log_difference(1256) = 0.09367455396034, Which is greater than upperbound hence we increment by 93. Now iterator is 1349. log_difference(1349) = 0.08946415071057 which is less than lowerbound. The next power is 1545 and we need to add 196 to get 1545.

Conditions are as follows:

1) If we find a power whose log_difference() is in the range of lower and upperbound, we will increment by 196. which implies that the power is a number which will lead to 123 as starting digits. 2) If we find a power, whose log_difference() is greater than or equal upperbound, we will increment by 93. 3) if log_difference() < lowerbound, we increment by 196.

Reference to the above logic: https://math.stackexchange.com/questions/4093970/powers-of-2-starting-with-123-does-a-pattern-exist

>>> solution(1000)
284168
>>> solution(56000)
15924915
>>> solution(678910)
193060223