ddidderr
04bb89375c
The binary now derives command-line parsing from clap instead of maintaining a custom parser. SEARCH_LIMIT remains an optional positional argument and keeps its default of 1_000_000_000, while clap now owns usage errors, --help, and --version output. The parser stores the limit as NonZeroU64 so zero is rejected before the search starts. The existing CLI parsing tests now exercise clap directly, and the README documents the generated help/version flags plus the top-level program structure. Test Plan: - cargo clippy - cargo clippy --benches - cargo clippy --tests - cargo test - cargo run -- --help Trailer: Refs: local request to replace custom argument parsing with clap Dependencies: none
2.0 KiB
hcn
Find highly composite numbers up to a search limit.
cargo run --release -- [SEARCH_LIMIT]
If no limit is given, the default is 1_000_000_000.
Use --help to see the generated command-line help and --version to print
the package version.
Structure
src/main.rs- Command-line interface: parses the optional
SEARCH_LIMITargument withclapand rejects non-positive limits before the search starts. - Search setup: defines the default limit and the fixed prime list used for candidate generation.
- Candidate generation: recursively builds only exponent sequences that can produce highly composite record candidates.
- Record filtering: sorts candidates and prints each new divisor-count record.
- Command-line interface: parses the optional
Algorithm
A number with prime factorization
n = 2^a * 3^b * 5^c * ...
has
(a + 1) * (b + 1) * (c + 1) * ...
divisors.
The exact primes do not affect the divisor count; only the exponents do. For a fixed exponent list, the smallest possible number is made by putting the largest exponent on the smallest prime:
a >= b >= c >= ...
If a larger exponent appears on a larger prime, swapping those two exponents keeps the divisor count unchanged but makes the number smaller.
So a highly composite number cannot have that shape. If it did, the smaller swapped number would already have the same number of divisors, so the original number would not be the first record.
That means record candidates are only numbers like:
2^a * 3^b * 5^c * ... where a >= b >= c >= ...
All other numbers are redundant for finding new records.
The program:
- Recursively generates numbers from the first primes.
- Only tries exponent sequences where each exponent is no larger than the previous one.
- Stops a branch as soon as the number exceeds the search limit.
- Computes the divisor count directly from the exponents.
- Sorts candidates by number and prints each new divisor-count record.
This avoids factoring every number in the range.