A103876 - OEIS

A103876

A test for divisibility by the n-th prime.

2, 1, 9, 5, 17, 16, 26, 3, 11, 4, 30, 14, 37, 53, 6, 20, 7, 51, 71, 58, 80, 29, 10, 72, 32, 98, 79, 38, 13, 41, 125, 134, 15, 47, 114, 50, 121, 161, 18, 19, 135, 59, 179, 21, 156, 68, 206, 163, 215, 24, 25, 77, 184, 242, 27, 83, 28, 198, 205, 92, 31, 219, 95, 33, 101, 104

(list; graph; refs; listen; history; text; internal format)

OFFSET

4,1

COMMENTS

Given a number M, remove its last digit d, then subtract d*a(n). If the result is divisible by prime(n), then M is also divisible by prime(n). This process may be repeated.

Values of a(n) can be quickly calculated by finding the smallest multiple of prime(n) which ends in a 1, and removing this last digit. E.g., 7 -> 21 -> 2, 11 -> 11 -> 1, 13 -> 91 -> 9, 17 -> 51 -> 5, 19 -> 171 -> 17.

REFERENCES

Alfred S. Posamentier, Math Charmers, Tantalizing Tidbits for the Mind, Prometheus Books, NY, 2003, page 76-81.

LINKS

Table of n, a(n) for n=4..69.

FORMULA

a(n) = p - (10 mod p)^(-1) where p = prime(n). - Joerg Arndt, Jan 23 2023

MATHEMATICA

a[n_] := Block[{p = Prime[n], k = 1}, While[ Mod[10k + 1, p] != 0, k++ ]; k]; Table[ a[n], {n, 4, 69}]

PROG

(Python)

import sympy

[pow(-10, -1, p) for p in sympy.primerange(7, 300)]