A371035 a(n) = A086330(prime(n)).

0, 2, 7, 18, 43, 73, 113, 159, 203, 334, 496, 706, 863, 874, 1097, 1124, 1560, 2033, 2073, 2409, 2462, 3336, 3345, 3634, 3958, 4657, 5198, 5284, 5186, 6096, 7801, 8594, 9270, 9167, 10659, 10578, 12375, 12227, 13221, 13769, 15958, 16458, 18820, 17919, 18722

Offset

1,2

Comments

The sequence sometimes decreases, as for example at a(29) = 5186 < 5284 = a(28).

Links

Table of n, a(n) for n=1..45.

Example

For n = 3, a(n) = A086330(prime(3)) = A086330(5) = (2! mod 5) + (3! mod 5) + (4! mod 5) = 2 + 1 + 4 = 7.

Prog

(Python)
from sympy import isprime
l = []
for i in range(2, 185):
    if isprime(i):
        sum = 0
        reminder = 1
        for j in range(2, i):
            reminder = (reminder * j) % i
            sum += reminder
        l.append(sum)
print(l)
(Python)
from sympy import prime
def A371035(n):
    a, c, p = 0, 1, prime(n)
    for m in range(2, p):
        c = c*m%p
        a += c
    return a # Chai Wah Wu, Apr 16 2024
(PARI) a(n) = my(p=prime(n)); sum(m=2, p, m! % p); \\ Michel Marcus, Apr 11 2024

Crossrefs

Cf. A086330, A000040.

Sequence in context: A147611 A007991 A037294 * A076857 A243717 A337482

Adjacent sequences: A371032 A371033 A371034 * A371036 A371037 A371038

Keyword

easy,nonn

Author

Alexandre Herrera, Apr 10 2024

Status

approved