So, I was watching an interview with Terence Tao which mentioned his famous Problem Six from the 1988 International Math Olympiad, which he only managed to get one point out of seven for. In the competition, the final problem he had to solve was:
Let a and b be positive integers such that ab + 1 divides a2 + b2. Show that a2 + b2Ā / ab + 1 is the square of an integer.
At first I thought this just meant you had to find a single solution at which it is true that:
(a2 + b2) mod (ab + 1) ā 0
AND
sqrt((a2 + b2) / (ab + 1)) = [some integer]
Which is actually a really easy problem. If you do the obvious thing when testing number-by-number, you get it immediately. If a = 1 and b = 1, then (12Ā + 12)/(1*1 + 1) = 2/2 = 1, and the root of 1 is of course 1.
But, like, this looked wayĀ tooĀ easy, so I wanted to know if there were other instances of this forming a perfect square. Then I found that a = 2 and b = 8 does it. Noticing a potential pattern, I checked a = 3 and b = 27. Sure enough, another solution.
At this point, I didnāt want to keep checking these by hand, so I tried the first 25 iterations of b = a3Ā using a Python script.
a = 1
while(a Ā <= 25):
Ā b = a**3
Ā if((a**2 + b**2)%(a*b + 1) == 0):
Ā Ā x = (a**2 + b**2)/(a*b + 1)
Ā Ā print("A =", a)
Ā Ā print("B =", b)
Ā Ā print("The result,", x, "is a square of", x**0.5)
Ā a +=1
Upon running it, I found that for all instances of b = a3, the result of the equation is a perfect square. Specifically, the result is a perfect square of a. ThatĀ was the cool part. (I already suspected this after trying a = 2 and a = 3, but double checking with a computer was trivial enough.)
So then I tried to figure out whyĀ that was the case. Firstly, if a3Ā = b, and a and b are both squared in the numerator, then the numerator is effectively a2Ā + a6. Next, if the denominator is a*(a3) + 1, then itās a4Ā + 1. Finally, the result of the equation is a2. Therefore:
Thus, it is proven that for any positive integers b and a such that b = a3, ab + 1 divides a2 + b2Ā and a2 + b2Ā / ab + 1 is the square of an integer.
...Except, that wasnāt what I needed to prove. I needed to prove that for any Z+ a and bĀ such that ab + 1 divides a2 + b2,Ā a2 + b2 / ab + 1 is a perfect square.
Aaaaand I donāt yet know how to do that. But Iām working on it! So far, I know that not all instances ofĀ āab + 1 divides a2 + b2ā³ are instances of b = a3, because a = 8 and b = 30 also does it, with a result of 4, which is the square of 2. Meanwhile, a isnāt even a factor of b here: 8 = 2*2*2, while 30 = 2*3*5.
Iām not sure what Iāll need to do to prove the general case, but Iāll work on it some more when I get back from dinner this evening. In the meantime, if anyone has hints/suggestions for how I can tackle this or things I should test, please let me know! Just please donāt write a proof in the notes unless you write it in words and rot13 it.
Anya is live and ready to show you everything. Watch her strip, dance, and perform exclusive shows just for you. Interact in real-time and make your fantasies come true.
ā Live Streamingā Interactive Chatā Private Showsā HD Quality
Anya is LIVE right now
FREE
Free to watch ⢠No registration required ⢠HD streaming
OK, so, yesterday I said I was investigating Problem 6 of the 1988 International Math Olympiad. The problem as stated is:
Let a and b be positive integers such that ab + 1 divides a2 + b2. Show that a2 + b2 / ab + 1 is the square of an integer.
I interpreted this to mean:
For any positive integersĀ a and b such that ab + 1 divides a2 + b2, show thatĀ a2 + b2 / ab + 1 must be a perfect square.
In the process of investigating it, I would up proving that for any positive integers b and a such that b = a3, ab + 1 divides a2 + b2and a2 + b2 / ab + 1 is the square of an integer. This is related, but not quite what Iām looking for. My specific reasoning was:
Firstly, if a3 = b, and a and b are both squared in the numerator, then the numerator is effectively a2 + a6. Next, if the denominator is a*(a3) + 1, then itās a4 + 1. Finally, the result of the equation is a2. Therefore:
However, b = a3Ā isnāt the only time ab + 1 divides a2Ā + b2, nor the only time the result is a perfect square. For example, I noticed yesterday that a = 8 and b = 30 does it.
Anyway, thatās what happened yesterday, so now youāre all caught up. Now, the thing that nerdsniped me for the past two hours: What is the pattern of theseĀ āirregularā solutions, such as a=8 & b=30?
I fired up my compiler and wrote a script to search for pairs of numbers like this, and I discovered a pattern. The equation is solved when a=8 & b=30, orĀ when a=30 & b=112, orĀ when a=112 & b=418...
That is, thereās a cycle. The b of one solution is the a of the next, with some larger number being b next time. This goes on for as far up as I checked*, which means thereās a long (I expect infinite) sequence of integers where you can plug in any adjacent pair and find a solution.
However, this wasnāt the only sequence like this. In fact, after I stopped thinking of these asĀ āthe irregular solutionsā, I realised there must be a bunchĀ of sequences like this. Let me explain:
The sequence I just discussed ({8, 30, 112, 418...}) doesnāt start at 8. It starts at 2. Because, as you might recall, a=2 & b=8 is a solution. So we have a sequence {2, 8, 30, 112, 418...}.
(The pairs you can form have interesting properties - the first has a ratio of 4 (ie: 8/2 =4), while every subsequent pair has a slightly smaller ratio [3.75, 3.73333 (recurring), 3.7321 (truncated), 3.7320 (truncated)]. Because the difference between adjacent ratios becomes smaller as well, I assume this converges on some limit. However, at present, Iām not sure how Iād find it.)
However, if youāre like me, when you saw that 2 & 8 was the start of a sequence, you thoughtĀ āHold on, does that mean 3 & 27Ā is also the start of a sequence?ā If so, you would be right! Thereās a wholeĀ ānotherĀ āirregularā sequence {3, 27, 240, 2133, 18957...}. This also consists of pair-wise solutions, and the ratios also start at a perfect square (32) and then decline toward some (supposed) lower bound.
Unsurprisingly, thereās a sequence like this for everyĀ āregularā pair I checked (ie: pairs where b = a3), and I assume thereās one for all integer values of a. (Unproven, but conjectured strongly enough that Iād be pretty surprised if I were wrong.)
The result of the equation when given inputs from any of the sequences ends up being a square, of course. Specifically, a square of the first number in the sequence. So, for the first sequence I talked about, the equation always outputs 4. For this reason, I consider the first element of the sequence to be the defining element, so Iāve been calling each sequence Sfirst element. So, the first sequence discussed would be S2. (I expect there isĀ a sequence S1, but that it simply contains {1, 1, 1, 1, 1, 1, 1...}).
The prime factorisations of the numbers in each sequence havenāt told me much useful stuff yet. So far, I canāt see a pattern that would let me predict** the next number in the sequence using its factorisation. Only two result really jump out at me:
The first thing is that every number in a given sequence has the first number of that sequence as a factor. This was intuitive to me so it wasnāt very interesting. What wasĀ interestingĀ is that the third, sixth, ninth and so on elements of every single sequence tested are all divisible by 10. I have no idea whyĀ this keeps happening, but it does, and thatās cool enough that it might merit more research.
However, if anyone reading this discovers anything more interesting in the factorisations, please let me know.
Anyway, this has almost certainly all been a massive sidetrack from actually solving the problem. But, hey, I started working on the problem because it was interesting, and this sidetrack was interesting, so I donāt consider this lost in any way.
*I checked for the existence of all the sequences S2...S10, and several members of each set, and the factorisations of some members, and the ratios of some pairs. All those fiddly numbers and my obscure hypotheses are under the cut for anyone who thought this post wasnāt already long enough, or wants to use my results to investigate further.
**Edit: Shit, I only realised after writing this that I should just look these sequences up in the G-d damn OEIS. Sure enough, every sequence Iāve checked is already there, and has an equation for it. Hereās S2. Iām really kicking myself now. Itās a bit late to go diving into them, but Iāll let you guys know what I uncover.
Hi, extremely strange people who want to read even more math. Here is some more stuff from my exploration:
Firstly, here are my guesses about what might be going on here. Each conjecture has varying amounts of confidence, based on how intuitively correct the pattern seems to me and how confident I am that I understand the relevant math.
Conjectures with >95% confidence:
There are no natural numbers a such that b = a3does not yield a solution.
:. The set So is countably infinite.
:. There is a countably infinite number of sets S1...Sn.
All members of the set SnĀ will have n as a factor.
Conjectures with >70% confidence:
Each set S1...Sn is countably infinite
The set of all the sets S1..n (hence forth Sab) is countably infinite (for the same reason the rationals are countable)
Every third element (the third, the sixth, the ninth, etc) of any sequence S2...Sn will be divisible by 10
The limit for the ratios for a given sequence will be greater than n2-1, for any sequence Sn.
Conjectures with >50% confidence:
There are no solutions to the equation which are not in the super set Sab (bearing in mind that So is a subset of Sab., since the listing the first two elements of each set S1...Sn gives you the first 2n terms of So.)
For each set Sn, every prime is a factor of some element of that set. (Strictly conditional on all sets S1...SnĀ being infinite.)
If you can prove, disprove, or otherwise contribute to me thinking about any of these conjectures, you will win the prize of me expressing my gratitude by sending emoji to your inbox. Thatās basically as valuable as a Fieldās Medal, right? š”šÆš”
Next up, my method of investigation thus far:
The first thing I did was iteratively test loooooots of possible {a, b} pairs. Like, so many. Literally a million. Hereās the code that did that:
This worked down from very high numbers to very low ones, testing lots of possibilities. The list generated included the original sequence of āregularā pairsĀ (So, which contains {{1, 1}, {2, 8}, {3, 27}, {4, 64}...}. Unlike the SnĀ sequences, only the pairs shown as subsets are valid solutions.), Ā where b = a3,Ā up to {17, 4913}.
Additionally, it generated the first few terms of S2Ā and S3, as well as the first term of each of S4, and S5. This is the point at which it was 100% clear to me that there was a sequence of sequences S2..n. (I then assumed a sequence S1Ā which is as described before the cut.)
In testing the numbers in each sequence, I found that the ratios of the pairs declined as I described before the cut. The ratios for the first couple elements of S2Ā are:
As you can see, it appears to converge. Iām not sure how Iād figure out what the limit is, but Iām pretty confident itās >3.7. Furthermore, you can see that each ratio is lower than the one before it, which means I could discount the possibility of a pair such that their ratio would be larger than the ratio of the pair before it. This allowed me to greatly restrict the search space when checking for eachĀ āfamilyā of numbers (ie: numbers in the same sequence).
My method for searching within eachĀ āfamilyā was to take its initial member āaā (eg: In S2Ā that would beĀ ā2ā²) and then test whether a*a2Ā was a solution. Given that it was, I would enter it into my array, assign b to be that number, and then take the ratio of b/a to be the largest possible ratio of a pair of solutions.
Then I would use the number Iād just discovered as my new āaā, multiply it by the ratio Iād found last time, and then check if that number was a solution. If not, Iād decrement it by one over and over again until I found a solution, and then Iād save that and repeat the cycle, with the ratio becoming smaller each time.
This algorithm, converted to Python, was:
stopper = 1
ceiling = 25
family = [5]
for a in family:
Ā b = int(a*ceiling)
Ā while(b > a):
Ā Ā if((a**2 + b**2)%(a*b + 1) == 0):
Ā Ā Ā x = (a**2 + b**2)/(a*b + 1)
Ā Ā Ā ceiling = b/a
Ā Ā Ā family.append(b)
Ā Ā Ā stopper +=1
Ā Ā Ā print("The ratio is", ceiling)
Ā Ā Ā break
Ā Ā b -=1
Ā if(stopper > 9):
Ā Ā break
print(family)
#This will print the first 10 elements of S5. If you want a different sequence Sn, then replace theĀ ā5ā² in familyās array with n and the 25 for ceiling with n**2
Once the difference between the ratios of pairs became small enough, b ended up being pretty much the same as a*ceiling rounded down to the nearest integer.
OK, if youāre still reading this, now is when I paste all the numbers I found in case you want the values. First, Iāll list the first elements of the the sets I found, with the ratio of the last two elements in each list. Iāll give the first ten for S2Ā to S4, and the first 5 for the others up to S10. Underneath that, Iāll give you some prime factorisations, in case youāre into those.
Colson Corp. had $700,000 net income in 2015. On January 1, 2015 there were 200,000 shares of common stock outstanding. On April 1, 20,000 shares were issued and on September 1, Colson bought 30,000 shares of treasury stock. There are 30,000 options to buy common stock at $40 a share outstanding. The market price of the common stock averaged $50ā¦