A078606 Constant c(p) used in determining divisibility by the n-th prime, p=A000040(n), for n>=4.
-2, -1, 4, -5, 2, 7, 3, -3, -11, -4, 13, -14, 16, 6, -6, -20, -7, 22, 8, 25, 9, -29, -10, 31, -32, 11, 34, -38, -13, -41, 14, 15, -15, -47, 49, -50, 52, 18, -18, -19, 58, -59, 20, -21, 67, -68, 23, 70, 24, -24, -25, -77, 79, 27, -27, -83, -28, 85, 88, -92, -31, 94, -95, -33, -101, -104, 35, 106, 36, -110, 112, 38
Offset: 4
Examples
The first few terms are c(7)=-2, c(11)=-1, c(13)=4. To find out if a number is divisible by 7, take the last digit, double it and subtract the result from the rest of the number. If you get an answer divisible by 7 (including 0), then the original number is divisible by 7. If you do not know the new number's divisibility, you can apply the rule again. Example: If you had 203, you would subtract 2*3=6 from 20 to get 14; since 14 is divisible by 7, so is 203.
Links
- Ethan Magness, Steven Sinnott and Robert L. Ward, Divisibility Rules
Programs
-
Mathematica
c[p_] := If[(v=PowerMod[10, -1, p])>p/2, v-p, v]; c/@Prime/@Range[4, 100]
-
PARI
a(n) = centerlift(Mod(1,prime(n))/10); \\ Kevin Ryde, Feb 18 2023
-
Python
import sympy [(pow(10, -1, p))-p*(p%10%6==1) for p in sympy.primerange(7, 300)] # Nicholas Stefan Georgescu, Jan 18 2023
Extensions
Edited by Dean Hickerson, Dec 23 2002
Comments