Project Euler Lab - Problem 752

#752 - Powers of $1+\sqrt 7$

● AdvancedOfficial difficulty: 45%Naive enumeration is infeasible; requires a mathematical reductionTier C - reduced scale in browser; full scale in notebookNot viewed
↖ Euler Lab

When \((1+\sqrt 7)\) is raised to an integral power, \(n\), we always get a number of the form \((a+b\sqrt 7)\).
We write \((1+\sqrt 7)^n = \alpha(n) + \beta(n)\sqrt 7\).

For a given number \(x\) we define \(g(x)\) to be the smallest positive integer \(n\) such that: \[\begin{align} \alpha(n) &\equiv 1 \pmod x\qquad \text{and }\\ \beta(n) &\equiv 0 \pmod x\end{align} \] and \(g(x) = 0\) if there is no such value of \(n\). For example, \(g(3) = 0\), \(g(5) = 12\).

Further define \[ G(N) = \sum_{x=2}^N g(x)\] You are given \(G(10^2) = 28891\) and \(G(10^3) = 13131583\).

Find \(G(10^6)\).

This problem is taken from Project Euler, Problem 752.
Problem text © Project Euler, licensed under CC BY-NC-SA 4.0. Original: projecteuler.net/problem=752. Published Sunday, 21st March 2021, 01:00 am. Solved by 800 members at time of mirroring.

Why this is useful

Optimization. The transferable skill is replacing infeasible enumeration with a mathematical reduction - the core move in calibration and large-scale computation (Phases 10, 13).

We classify relevance honestly - not every Euler problem is a trading application.

Prerequisites

Lessons that prepare you:
19.14 Computational Complexity, Feasibility Estimation, and Proving Algorithms Correct · 19.10 Search: Backtracking, Branch-and-Bound, Binary Search, Meet-in-the-Middle

Recommended stepping-stone problems: #209 · #637 · #749

Concepts: brute-force-reduction

Likely techniques: binary-search

Learning mode

Pick how much scaffolding you want. Your choice is remembered per problem.

Scratchpad

Mathematical notes, formulas, pseudocode, hypotheses, complexity notes. Saved automatically with your progress.

Python workbench

Tier C - reduced scale in browser; full scale in notebook
Browser runs a reduced, clearly-labelled educational scale; the original scale is provided in a local notebook.

Real Python (Pyodide) in a sandboxed Web Worker - no network, no filesystem, no DOM access. Ctrl/Cmd+Enter runs. Escape leaves the editor. Stop terminates the worker.

Python runtime not loaded (it boots on first run - a one-time local load).

Check your answer

Answers are checked against a salted hash held in a separate file - not printed in this page. This prevents accidental spoilers; it is not cryptographic protection (see the build notes).

Progressive hints

Confidence

Low confidence schedules this problem for spaced review, even if you solved it.

Reference solution

Spoiler
The complete original explanation (interpretation, naive approach, insight, proof, complexity, Python implementation, tests, common mistakes, alternatives) is hidden and lazy-loaded. Reveal it only after a meaningful attempt - the struggle is where the learning happens.