Home
About
Resume
Projects
Links
Blog
Download notebook
{ "cells": [ { "cell_type": "markdown", "id": "a469358e-2f57-4bd6-9c8d-c59f3d77920b", "metadata": {}, "source": [ "### Q44\n", "Find the pair of pentagonal numbers, \\(P_j\\) and \\(P_k\\), for which their sum and difference are pentagonal and \\( D = |P_k − P_j| \\) is minimised; what is the value of D?" ] }, { "cell_type": "code", "execution_count": 1, "id": "820c6f59-d5d1-4f67-a379-b41d8bb337ac", "metadata": {}, "outputs": [], "source": [ "def pentagonal_num(n):\n", " return n*(3*n-1)//2\n", "\n", "def min_D():\n", " pentagonal_num_dict = {}\n", " found_pentagonal_num = set()\n", " i = 1\n", " while True: \n", " if i not in pentagonal_num_dict:\n", " pentagonal_num_dict[i] = pentagonal_num(i)\n", " found_pentagonal_num.add(pentagonal_num_dict[i])\n", " for j in range(1,i)[::-1]:\n", " if 2*pentagonal_num_dict[j] >= pentagonal_num_dict[i]:\n", " continue\n", " if (pentagonal_num_dict[i] - pentagonal_num_dict[j] in found_pentagonal_num and\n", " pentagonal_num_dict[i] - 2*pentagonal_num_dict[j] in found_pentagonal_num):\n", " return pentagonal_num_dict[i] - 2*pentagonal_num_dict[j]\n", " i += 1" ] }, { "cell_type": "code", "execution_count": 2, "id": "1ef770af-4b63-4d07-abef-0cbe9483e888", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "CPU times: user 250 ms, sys: 0 ns, total: 250 ms\n", "Wall time: 246 ms\n" ] }, { "data": { "text/plain": [ "5482660" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "%%time\n", "min_D()" ] }, { "cell_type": "code", "execution_count": null, "id": "5ba24417-6761-490e-8ea3-524a29a4ebef", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.10" } }, "nbformat": 4, "nbformat_minor": 5 }
Previous Notebook:
Project Euler Q43
Next Notebook:
Project Euler Q45
Loading