cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 21-22 of 22 results.

A290560 Generalized Lucas-Carmichael numbers for D=9697.

Original entry on oeis.org

1, 35, 143, 323, 385, 455, 595, 665, 899, 935, 1045, 1295, 1547, 1729, 2639, 2737, 2821, 2915, 3289, 3689, 4355, 4465, 5005, 5183, 5291, 6479, 6721, 8855, 8911, 9215, 9361, 10153, 10439, 10465, 11305, 11663, 11951, 15841, 17119, 18095, 19981, 20909, 22607
Offset: 1

Views

Author

Keywords

Comments

On the set Lc(Z/NZ,D) = {(x,y) in (Z/NZ)^2 : x^2 - Dy^2 = 1 (mod N)}, define an operation as follows: (x,y)x(z,w) = (xz+Dyw, xw+zy) (mod N). The set Lc(Z/NZ, D) endowed with this operation is a group. Moreover, the set of Lucas numbers endowed with this operation is a subgroup of Lc(Z/NZ, D).
The following results appear in Babinkostova, et al.: If q is a prime, then #Lc(Z/(q^e)Z, D) = (q-(D|q))q^(e-1).
The group Lc(Z/(q^e)Z, D) is cyclic for e > 0. This result was proven in Hinkel, 2007 for the case when e = 1. We showed that the statement is true for e > 1 (Babinkostova, et al.).
The following notions are introduced in Babinkostova, et al.: A composite integer N is a generalized Lucas pseudoprime (or Lucas pseudoprime in Babinkostova, et al.) to base P in Lc(Z/NZ, D) and integer D if (N-(D|N))P = O, where O is the identity of the group.
We define a composite integer N to be a generalized Lucas-Carmichael number if for all P in Lc(Z/NZ, D) it is true that (N-(D|N))P = O.
The following Korselt-like criterion holds for a generalized Lucas-Carmichael number: A composite number N is a generalized Lucas-Carmichael number if and only if N is squarefree and for every prime factor q of N, (q-(D|q)) divides (N-(D|N)).
This sequence is a list of generalized Lucas-Carmichael numbers for D=9697.
For prime values of D less than 10000 and odd nonprime values of N less than 1000000, this is the longest sequence of generalized Lucas-Carmichael numbers.
The resulting sequence of generalized Lucas-Carmichael numbers is based on work done by L. Babinkostova, B. Bentz, M. I. Hassan, and H. Kim.

Examples

			We will illustrate an example using the Korselt criterion for generalized Lucas pseudoprimes. Let us observe the second term, 35. Note that 35 = 5*7, so that it is squarefree. Now note that (5-(9697|5)) = 6 and (7-(9697|7)) = 6, both of which divide (35-(9697|35)) = 36. Therefore, by the Korselt criterion for generalized Lucas pseudoprimes, we have that 35 is a generalized Lucas Carmichael number for D = 9697.
		

Crossrefs

Programs

  • SageMath
    # A program in SageMath is given in the links section.

A371805 Composite numbers k that divide A001644(k) - 1.

Original entry on oeis.org

182, 25201, 54289, 63618, 194390, 750890, 804055, 1889041, 2487941, 3542533, 3761251, 6829689, 12032021, 12649337, 18002881
Offset: 1

Views

Author

Robert FERREOL, Apr 06 2024

Keywords

Comments

If k is prime, k divides A001644(k) - 1; and since A001644(k) satisfies a tribonacci recurrence, these numbers could be called tribonacci pseudoprimes.

Examples

			(A001644(182)-1)/182 = 8056145960961609628091266244940745410646318417.
		

Crossrefs

Cf. A001644.
Cf. A005845 (composite k that divide Lucas(k) - 1).
Cf. A013998 (composite k that divide Perrin(k) - 1).

Programs

  • Maple
    A001644:=proc(n) option remember: if n=0 then 3 elif n=1 then 1 elif n=2 then 3 else A001644(n-1)+A001644(n-2)+A001644(n-3) fi end:
    test:=n->A001644(n) mod n = 1:select(test and not isprime, [seq(n, n=1..100000)]);
  • Mathematica
    seq[kmax_] := Module[{x = 1, y = 3, z = 7, s = {}, t}, Do[t = x + y + z; If[Mod[t, k] == 1 && CompositeQ[k], AppendTo[s, k]]; x = y; y = z; z = t, {k, 4, kmax}]; s]; seq[200000] (* Amiram Eldar, Apr 06 2024 *)
  • Python
    from sympy import isprime
    from itertools import count, islice
    def agen(): # generator of terms
        t0, t1, t2 = 3, 1, 3
        for k in count(1):
            t0, t1, t2 = t1, t2, t0+t1+t2
            if k > 1 and not isprime(k) and (t0-1)%k == 0:
                yield k
    print(list(islice(agen(), 5))) # Michael S. Branicky, Apr 07 2024

Extensions

a(13)-a(15) from Amiram Eldar, Apr 07 2024
Previous Showing 21-22 of 22 results.