<html><head><meta name="color-scheme" content="light dark"></head><body><pre style="word-wrap: break-word; white-space: pre-wrap;">\documentclass[12pt,reqno]{article}

\usepackage[usenames]{color}
\usepackage{amssymb}
\usepackage{graphicx}
\usepackage{amscd}

\usepackage[colorlinks=true,
linkcolor=webgreen,
filecolor=webbrown,
citecolor=webgreen]{hyperref}

\definecolor{webgreen}{rgb}{0,.5,0}
\definecolor{webbrown}{rgb}{.6,0,0}

\usepackage{color}
\usepackage{fullpage}
\usepackage{float}

\usepackage{psfig}
\usepackage{graphics,amsmath,amssymb}
\usepackage{amsthm}
\usepackage{amsfonts}
\usepackage{latexsym}
\usepackage{epsf}

\setlength{\textwidth}{6.5in}
\setlength{\oddsidemargin}{.1in}
\setlength{\evensidemargin}{.1in}
\setlength{\topmargin}{-.1in}
\setlength{\textheight}{8.4in}

\newcommand{\seqnum}[1]{\href{http://oeis.org/#1}{\underline{#1}}}

\begin{document}

\begin{center}
\epsfxsize=4in
\leavevmode\epsffile{logo129.eps}
\end{center}

\theoremstyle{plain}
\newtheorem{theorem}{Theorem}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}

\theoremstyle{definition}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{example}[theorem]{Example}
\newtheorem{conjecture}[theorem]{Conjecture}

\theoremstyle{remark}
\newtheorem{remark}[theorem]{Remark}

\begin{center}
\vskip 1cm{\LARGE\bf Generalized Multiple Counting Jacobsthal \\
\vskip .12in Sequences of Fermat Pseudoprimes}
\vskip 1cm
\large
M. H\"{u}srev C{\i}lasun\\
Emerging Circuits and Computation Group\\
Electrical and Electronics Faculty\\
Istanbul Technical University\\
Maslak, Istanbul\\
Turkey\\
\href{mailto:cilasun@itu.edu.tr}{\tt cilasun@itu.edu.tr}
\end{center}

\vskip 0.2 in

\begin{abstract}
 This study involves definitions of regular and representational multiple-counting Jacobsthal sequences of Carmichael numbers. We introduce recurrence relations for multiple-counting Jacobsthal sequences and show their association with Fermat's little theorem. We also provide matrix representations and generalized Binet formulas for defined sequences. This leads to a better understanding of how certain composite numbers are distributed among consecutive powers.
\end{abstract}

\section{Introduction}
The Jacobsthal sequence \cite{ref1} is defined as follows:
\begin{equation}
j_{n+2}=j_{n+1}+2j_{n}, j_{1}=1, j_{0}=0, n\geq 0  \label{eq:1}
\end{equation}
Horadam and Hoggatt have written many papers \cite{ref1,ref4,ref9} about the Jacobsthal sequence, which has many remarkable properties, such as counting microcontroller skip instructions \cite{ref11} and counting the number of ways to tile a $3\times(n-1)$ rectangle with $2\times 2$ and $1\times 1$ tiles \cite{ref10}. The Jacobsthal sequence can be expressed in floor function notation \cite{ref2}. It is given in the description of sequence \seqnum{A001045} that the elements of the Jacobsthal sequence count the multiples of $3$ in between successive powers of $2$, as illustrated in Table~\ref{tablo1}. 

\begin{table}[h]
 \centering \footnotesize {
\begin{tabular}{|c|c c c c c c c c c c c|}
\hline
Interval boundaries &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;\\
of powers of 2 &amp; $2^0$ &amp;  &amp; $2^1$ &amp;  &amp; $2^2$ &amp;  &amp; $2^3$ &amp;  &amp; $2^4$ &amp;  &amp; $%
2^5$ \\ \hline
Number of multiples of 3 &amp;  &amp; 0 &amp;  &amp; 1 &amp;  &amp; 1 &amp;  &amp; 3 &amp;  &amp; 5&amp; \\
 in each interval. &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;  &amp;\\
\hline
\end{tabular}
}
\caption{Number of multiples of $3$ in between powers of $2$. The sequence is identical to the original Jacobsthal sequence.}
\label{tablo1}
\end{table}

As the initial motivation of this paper, it is notable that there are other Jacobsthal-like recurrence relations that hold between some other exponents and bases. For instance, as it can be understood from a superficial glance at the sequence \seqnum{A007910}, when double of any term is added to three consecutive terms, the very next term is obtained. This study investigates various aspects of those recurrence sequences. 

Another concept that will be discussed in this paper is Fermat's little theorem  \cite{ref3}, which can be summarized as $x^{\rho -1}\equiv 1$ (mod $\rho$). The theorem states that the given condition is satisfied when $\rho$ and $x$ are coprime and $\rho$ is prime. However, vice versa is not always true. The condition is still valid for some composite $\rho$s called \emph{Fermat pseudoprimes} \cite{ferm}, \emph{Carmichael numbers} \cite{ref7} or \emph{$K_1$ Kn\"{o}del numbers} \cite{knodel}. In particular for the $x=2$ case, the $\rho$s are called \emph{Poulet numbers} \cite{pouletsarrus} or \emph{Fermatians} \cite{ref7} if they satisfy Fermat's given condition \cite{ref3}.

\section{Generalized multiple-counting Jacobsthal sequences}

\begin{definition}
	Let $x$ and $\rho$ be elements of $\{2,3,4,5,\dots \}$ and $n$ be an element of $\{0,1,2,3,\dots \}$. The \textit{multiple-counting Jacobsthal sequence} $J(x,\rho)$ is defined by setting $J_{n}$ equal to the number of multiples of $\rho $ which are greater than $x^{n}$ and less than $x^{n+1}$.
\end{definition}
	In floor function notation, the number of multiples of $\rho$ less than $x^{n+1}$ is $\left \lfloor \frac{x^{n+1}}{\rho} \right \rfloor$. Therefore, the number of multiples of $\rho$ greater than $x^{n}$ and less than $x^{n+1}$ is found as follows:
\begin{equation} \label{floornot}
J_n = \left \lfloor \frac{x^{n+1}}{\rho} \right \rfloor - \left \lfloor \frac{x^{n}}{\rho} \right \rfloor
\end{equation}


\begin{theorem} Let $x$ and $\rho $ be elements of $\{1,2,3,\dots
\}$. Let $x $ and $\rho $ be relatively prime and satisfy
Fermat's $x^{\rho -1}\equiv  1$ (mod $\rho$) condition. For any number $n$ from the set $\{0,1,2,\dots \}$, the sequence with starting terms $J_{0},J_{1},\dots J_{\rho -1}$ which counts the multiples of $\rho $ between $x^{n}$ and $x^{n+1}$ also satisfies
the following recurrence relation for $\rho \geq 3$:
\begin{equation}
J_{n+\rho -1}=(x-1)\sum_{i=1}^{\rho -2}{J_{n-i+\rho -1}}+xJ_{n}
\label{eq:2}
\end{equation}
\end{theorem}%


\begin{proof}
	
	The expression in \eqref{eq:2} can be transformed to
	\begin{equation}\label{eq:3}
	J_{n+\rho-1}=(x-1)\sum^{\rho-1}_{i=1}{J_{n-i+\rho-1}}+J_n
	\end{equation}
Replace $J_n$ with its floor function notation \eqref{floornot}. The following expression is obtained when intermediate elements of the sigma notation are eliminated.
\begin{equation}\label{eq:e}
 \left \lfloor \frac{x^{n+\rho}}{\rho}\right\rfloor-\left\lfloor\frac{x^{n+\rho-1}}{\rho}\right\rfloor= (x-1)\left ( \left\lfloor\frac{x^{n+\rho-1}}{\rho}\right\rfloor -\left\lfloor\frac{x^n}{\rho}\right\rfloor\right ) +\left\lfloor \frac{x^{n+1}}{\rho}\right\rfloor-\left\lfloor\frac{x^n}{\rho}\right\rfloor
\end{equation}
Let $a$ be an arbitrary constant. Assume $x^{\rho-1}$ is equal to $a\rho+1$ as a result of the $x^{\rho - 1} \equiv 1$ (mod $\rho$) condition. When $a\rho+1$ is substituted into \eqref{eq:e}, the floor function takes the $a$-terms out due to its integer exclusion property. Both sides of the equation turn into $a J_n$. The right and left hand sides then simplify, and the proof is complete. 
\end{proof}

\begin{example}
Table~\ref{tablo4} contains several examples of multiple-counting sequences for different $\rho$s and $x$s

\begin{table}[h]
\centering
\begin{tabular}{|c|c|c|}
\hline 
$\rho$ &amp; $x$ &amp;\textbf{Notes on sequences} \\ 
\hline 
3 &amp; 2 &amp; {\footnotesize Identical to Jacobsthal sequence: $J_n=J_{n-1}+2J_{n-2},n\geq 2, J_0=0, J_1=1$} \\
\hline 
5 &amp; 3 &amp; {\footnotesize $J_n=2J_{n-1}+2J_{n-2}+2J_{n-3}+3J_{n-4},n\geq 4, J_0=0, J_1=1, J_2=4, J_3=11$}\\
\hline
3 &amp; 10 &amp; {\footnotesize $J_n=9J_{n-1}+10J_{n-2},n\geq 2, J_0=3, J_1=30$}\\
&amp;&amp;  \footnotesize $J_n$ counts $(n+1)$-digit numbers which are divisible by 3 \\
\hline
\end{tabular}
\caption{Some multiple-counting sequences} \label{tablo4}
\end{table}

In their paper, Cook and Bacon \cite{ref8} defines various sequences $J_n$ for various $\rho$s, called \textit{higher order Jacobsthal sequences}, when $x=2$. 
\end{example}

\begin{theorem} \label{binet}
For any element of $J_n$, the Binet expression
\begin{equation}  \label{eq:7}
J_n=\sum^n_{i=1}{\lambda_i^{n-i+2}v_{1,i}(\sum_{j=1}^n{v_{i,j}^{-1}u_j^{%
\rho-2}})}
\end{equation}
is always satisfied when the key matrix

\[
K_{(\rho-1) \times (\rho-1)}= \left(
\begin{array}{cccccc}
x-1 &amp; x-1 &amp; \dots &amp; x-1 &amp; x-1 &amp; x \\
1 &amp; 0 &amp; \dots &amp; 0 &amp; 0 &amp; 0 \\
0 &amp; 1 &amp; \ddots &amp; 0 &amp; 0 &amp; 0 \\
\vdots &amp; \ddots &amp; \ddots &amp; \ddots &amp; \vdots &amp; \vdots \\
0 &amp; 0 &amp; \ddots &amp; 1 &amp; 0 &amp; 0 \\
0 &amp; 0 &amp; \dots &amp; 0 &amp; 1 &amp; 0 \\
&amp;  &amp;  &amp;  &amp;  &amp;
\end{array}
\right)
\]

satisfies $Ku^n=u^{n+1}$ for the vector

\[
u^n= \left(
\begin{array}{c}
u_1 \\
u_2 \\
\vdots \\
u_{p-1} \\
\end{array}
\right) = \left(
\begin{array}{c}
J_n \\
J_{n-1} \\
\vdots \\
J_{n-p+2} \\
\end{array}
\right)
\]
if $K$ is diagonalizable and $K$'s eigenvector matrix $V$ is non-singular.
Note that elements $u_{i}$,$v_{ij}$ belong to $u$,$V$.
\end{theorem}


\begin{proof}
The Binet form for $J_n$ can be obtained using Kalman's method \cite{ref5}. After decomposing $K=V\Lambda V^{-1}$, the equation $Ku^n=u^{n+1}$  is transformed to $V\Lambda V^{-1}u^n=u^{n+1}$ where $\Lambda$ is the diagonal eigenvalue matrix of $K$. We only need the top elements of vectors, thus multiply the expression by
$$
\left(
\begin{array}{cccc}
1&amp;0&amp;\dots&amp;0\\
\end{array}
\right)
$$
a $(\rho-1)\times 1$ matrix from the left. Any particular term of the sequence
$J_n$ satisfies the equation
\begin{equation} \label{regmat}
J_n= \left(
\begin{array}{cccc}
1&amp;0&amp;\dots&amp;0\\
\end{array}
\right) (V\Lambda^{n-\rho+2}(V^{-1}u^{\rho-2}))
\end{equation}
Since
$$
(V^{-1}u^{\rho-2}) = \left(
\begin{array}{cccc}
v^{-1}_{1,1}&amp;v^{-1}_{1,2}&amp;\dots&amp;v^{-1}_{1,\rho-1}\\
v^{-1}_{2,1}&amp;v^{-1}_{2,2}&amp;\dots&amp;v^{-1}_{2,\rho-1}\\
\dots&amp;\dots&amp;\ddots&amp;\dots\\
v^{-1}_{{\rho-1 },1}&amp;v^{-1}_{{\rho-1},2}&amp;\dots&amp;v^{-1}_{\rho-1,\rho-1}\\
\end{array}
\right) \left(
\begin{array}{c}
J_{\rho-2}\\
J_{\rho-3}\\
\vdots\\
J_{0}\\
\end{array}
\right)
$$
 and
$$
 V =
\left(
\begin{array}{cccc}
v_{1,1}&amp;v_{1,2}&amp;\dots&amp;v_{1,\rho-1}\\
v_{2,1}&amp;v_{2,2}&amp;\dots&amp;v_{2,\rho-1}\\
\dots&amp;\dots&amp;\ddots&amp;\dots\\
v_{{\rho-1 },1}&amp;v_{{\rho-1},2}&amp;\dots&amp;v_{\rho-1,\rho-1}\\
\end{array}\\
\right)
$$
and
$$
\Lambda^{n-\rho+2} =
\left(
\begin{array}{cccc}
\lambda_{1}^{n-\rho+2}&amp;0&amp;\dots&amp;0\\
0&amp;\lambda_{2}^{n-\rho+2}&amp;\ddots&amp;\vdots\\
\vdots&amp;\ddots&amp;\ddots&amp;0\\
0&amp;\dots&amp;0&amp;\lambda_{\rho-1}^{n-\rho+2}\\
\end{array}
\right)\\
$$
are known, we can compute \eqref{regmat} after the required substitution.
Recall that $\lambda_1,\lambda_2, \ldots, \lambda_n$ 
denote the eigenvalues of $K$. After the outcome of the computation is simplified, we obtain the initial expression \eqref{eq:7}. 
\end{proof}

\begin{example}
	Table~\ref{tablo5} contains several Binet formulas of multiple-counting sequences for different $\rho$s and $x$s
	
	\begin{table}[h]
		\centering
		\begin{tabular}{|c|c|c|c|}
			\hline 
			$\rho$ &amp; $x$ &amp;\textbf{Sequence}&amp; \textbf{Binet form} \\ 
			\hline 
			3 &amp; 2 &amp;  $J_n=J_{n-1}+2J_{n-2}$&amp; $\frac{2^n}{3} - \frac{(-1)^n}{3}$\\ 
			 &amp;  &amp; {\footnotesize $n\geq 2, J_0=0, J_1=1$} &amp;\\
			\hline 
			3 &amp; 4 &amp; {\footnotesize $J_n=3J_{n-1}+4J_{n-2}$}&amp;{\footnotesize$4^n$}\\
			&amp; &amp; {\footnotesize $n\geq 2, J_0=1, J_1=4$}&amp;\\
			\hline
			3 &amp; 10 &amp; {\footnotesize $J_n=9J_{n-1}+10J_{n-2}$}&amp; {\footnotesize$	3 \frac{10^n}{10}$}\\
			 &amp; &amp; {\footnotesize $n\geq 2, J_0=3, J_1=30$}&amp;\\	
			\hline
		\end{tabular}
		\caption{Some Binet forms of multiple-counting sequences when $\rho$ is 3.} \label{tablo5}
	\end{table}
\end{example}

\vskip 0.7cm
\section{Generalized multiple-counting representational Jacobsthal sequences}

\begin{definition}
Let $x$ and $\rho$ be elements of  $\{2,3,4,5, \ldots\}$ and $n$ be an element of $\{0, 1, 2, 3, \dots\}$. Then $S(x,\rho)$ will be called a \textit{multiple-counting representational Jacobsthal sequence} when $S_n$ is equal to the number of integer multiples of $\rho$, that are greater than $0$ and less than $x^{n+1}$. Representational sequences can be compactly expressed as
\begin{equation}
S_n =\left \lfloor \frac{x^{n+1}}{\rho} \right \rfloor
\end{equation}
The strong interrelation between regular multiple-counting sequences and multiple-counting representational Jacobsthal sequences is expressed by

\begin{equation}\label{eq:11}
S_n=\sum^{n}_{i=0}{J_i}
\end{equation}

\end{definition}

The summation idea of a multiple-counting Jacobsthal sequence is very similar to the \textit{Jacobsthal representation sequence} which Horadam \cite{ref4} explores in detail. The representation sequence satisfies a similar recurrence relation. However, the constant term is different.

\begin{theorem}
Let $x$ and $\rho$ be elements of $\{1,2,3,\dots\}$ and suppose they satisfy Fermat's $x^{\rho-1}\equiv 1$ (mod $\rho$) condition. For any element $n$ of $\{0,1,2,\dots\}$, the
sequence with starting terms $S_0,S_1, \dots, S_{p-1}$, which counts multiples of $\rho$ between $0$ and
$x^{n+1}$, satisfies the following recurrence relation:
\begin{equation}  \label{eq:8}
S_{n+\rho-1}=(x-1)\sum^{\rho-2}_{i=1}{S_{n-i+\rho-1}}+xS_n+\phi
\end{equation}
or
\begin{equation}  \label{eq:9}
S_{n+\rho-1}=(x-1)\sum^{\rho-1}_{i=1}{S_{n-i+\rho-1}}+S_n+\phi,
\end{equation}
where $\phi$ is a unique constant derived from the particular combination of
$\rho$, $x$ and the  initial conditions $J_0,J_1,\dots, J_{p-2}$. Note that the formula for $\phi$ is
\begin{equation}  \label{eq:10}
\phi=\sum^{\rho-1}_{i=0}{J_i}-(x-1)\sum^{\rho-1}_{i=1}{\sum^{i-2}_{j=0}{J_j}
}
\end{equation}
\end{theorem}

\begin{proof}
Rewrite \eqref{eq:9} using property \eqref{eq:11}:
\begin{equation}\label{eq:12}
\sum^{n+\rho-1}_{i=0}{J_i}=(x-1)\sum^{\rho-1}_{i=1}{\sum^{n-i+\rho-1}_{j=0}{J_j}}+\sum^{n}_{i=0}{J_i}+\phi
\end{equation}
Separate \eqref{eq:12} into $n$-dependent and $n$-independent parts:

\begin{equation}\label{eq:13}
\sum^{\rho-2}_{i=0}{J_i}+\sum^{n+\rho-1}_{i=\rho-1}{J_i}=(x-1)\sum^{\rho-1}_{i=1}{(\sum^{i-2}_{j=0}{J_j}+\sum^{n-i+\rho-1}_{j=i-1}{J_j})}+\sum^{n}_{i=0}{J_i}+\phi
\end{equation}

\begin{equation}\label{eq:14}
\sum^{\rho-2}_{i=0}{J_i}+\underbrace{\sum^{n+\rho-1}_{i=\rho-1}{J_i}}=(x-1)\sum^{\rho-1}_{i=1}{\sum^{i-2}_{j=0}{J_j}}+\underbrace{(x-1)\sum^{\rho-1}_{i=1}{\sum^{n-i+\rho-1}_{j=i-1}{J_j}}}+\underbrace{\sum^{n}_{i=0}{J_i}}+\phi
\end{equation}

The $n$-dependent underbraced parts of Eq.~\eqref{eq:14} are
eliminated by induction on \eqref{eq:3}. Since the $n$-independent parts
of the expression above are constant for the whole sequence, $\phi$ also
depends on a constant value and the required proof is provided.
\end{proof}

\begin{theorem}
To obtain a matrix-based Binet representation for the sequence $S$, we
use a similar approach as for the sequence $J$. Since the only
difference between the definitions of the two sequences $J$ and $S$ is a
constant $\phi$, we can easily modify the key matrix to insert the $\phi$
dependency with notational changes to get a simple formula.
\begin{equation}  \label{eq:15}
S_n=\sum^{n}_{i=1}{\mu^{n-i+2}_{i}w_{1,i}(\sum^{n}_{j=1}{w^{-1}_{i,j}%
\tau_{j}^{\rho-1}})}
\end{equation}
Assuming $L$ is the key matrix
\[
L_{\rho \times \rho}= \left(
\begin{array}{ccccccc}
x-1 &amp; x-1 &amp; \dots &amp; x-1 &amp; x-1 &amp; x &amp; 1 \\
1 &amp; 0 &amp; \dots &amp; 0 &amp; 0 &amp; 0 &amp; 0 \\
0 &amp; 1 &amp; \ddots &amp; 0 &amp; 0 &amp; 0 &amp; 0 \\
\vdots &amp; \ddots &amp; \ddots &amp; \ddots &amp; \vdots &amp; \vdots &amp;  \\
0 &amp; 0 &amp; \ddots &amp; 1 &amp; 0 &amp; 0 &amp; 0 \\
0 &amp; 0 &amp; \dots &amp; 0 &amp; 1 &amp; 0 &amp; 0 \\
0 &amp; 0 &amp; \dots &amp; 0 &amp; 0 &amp; 0 &amp; 1 \\
&amp;  &amp;  &amp;  &amp;  &amp;  &amp;
\end{array}
\right)
\]
of $L\tau^n=\tau^{n+1}$ where $\tau$ satisfies
\[
\tau^n= \left(
\begin{array}{c}
\tau_1 \\
\tau_2 \\
\vdots \\
\tau_{p-1} \\
\phi \\
\end{array}
\right) = \left(
\begin{array}{c}
S_n \\
S_{n-1} \\
\vdots \\
S_{n-p+2} \\
\phi \\
\end{array}
\right)
\]
if $L$ is diagonalizable and $L$'s eigenvector matrix $W$ is non-singular.
Note that the elements $w_{ij}$ belong to $W$ and $\mu_1,\mu_2,
\dots, \mu_n$ denote the eigenvalues of $L$. 
\end{theorem}
\begin{proof}
Proof is similar to the proof of Theorem~\ref{binet} and therefore omitted.
\end{proof}

\section{Concluding remarks}
\hfill
   
Several sequences in OEIS \cite{ref2} are identical to some multiple-counting Jacobsthal sequences. \seqnum{A007910} is a multiple-counting Jacobsthal sequence when ($x=2$, $\rho=5$) and \seqnum{A077947} when ($x=2$, $\rho=7$), which are also called fourth and sixth order Jacobsthal sequences \cite{ref8}. Additionally, \seqnum{A000302} and \seqnum{A093138} are multiple-counting Jacobsthal sequences for ($x=4$, $\rho=3$) and ($x=10$, $\rho=3$).

Unlike the prime-counting function $\pi(x)$ (which is \seqnum{A006880} of OEIS \cite{ref2}), the sequences $J_n$ and $S_n$ count the composites of a given prime $\rho$. There is only one exception for both sequences: their initial terms. If $\rho$ is prime, $J_n$ and $S_n$ count all composites of the prime $\rho$.

The outcomes of this research can be used in various fields, such as electronics and computer science. The strong connection between the described sequences and Fermat's little theorem might lead to a better understanding of number fields. The distribution of prime and composite numbers is a significant field in present day number theory. Considering this, the future of multiple-counting Jacobsthal sequences is promising.
 
\section{Acknowledgement}
I would like to thank Mustafa Sefa Y\"{o}r\"{u}ker, Mustafa Zengin, Adem \c{S}ahin, and Ayhan Yenilmez for helpful advice.

\begin{thebibliography}{99}

\bibitem{ref2} 
N. Sloane, The On-Line Encyclopedia of Integer Sequences, 2015,
\url{http://oeis.org}.

\bibitem{ref1} 
A. F. Horadam, 
Jacobsthal and Pell curves, 
\textit{Fibonacci Quart.} \textbf{26} (1988), 
77--83.

\bibitem{ref3} 
G. H. Hardy and E. M. Wright,
\textit{An Introduction to the Theory of Numbers},
Oxford Univ. Press, 
1979.

\bibitem{ref4}
A. F. Horadam,
Jacobsthal representation numbers,
\textit{Fibonacci Quart.} \textbf{34} (1996), 
40--54.

\bibitem{ref5} 
D. Kalman, 
Generalized Fibonacci numbers by matrix methods, 
\textit{Fibonacci Quart.} \textbf{20} (1982), 
73--76.

\bibitem{ref7} 
D. Shanks, 
\textit{Solved and Unsolved Problems in Number Theory},  
AMS Chelsea, 
2001.

\bibitem{ref8} 
C. K. Cook and M. R. Bacon, 
Some identities for Jacobsthal and Jacobsthal-Lucas numbers satisfying higher order recurrence relations, 
\textit{Ann. Math. Inform.} \textbf{41} (2013), 
27--39.

\bibitem{ref9}
G. E. Bergum, L. Bennett, A. F. Horadam, and S. D. Moore, 
Jacobsthal polynomials and a conjecture concerning Fibonacci-like matrices, 
\textit{Fibonacci Quart.} \textbf{23} (1985), 
240--248.

\bibitem{ref10}
S. Heubach, 
Tiling an $m$-by-$n$ area with squares of size up to $k$-by-$k$ $(m\leq 5)$,
\textit{Congr. Numer.} \textbf{140} (1999), 
43--64.

\bibitem{ref11}
A. Da\c{s}demir, On the Jacobsthal numbers by matrix method, 
\textit{Fen Derg.} \textbf{7} (2012), 
69--76.


\bibitem{ferm}
R. Crandall and C. B. Pomerance, 
\textit{Prime Numbers: a Computational Perspective}, 
Springer Science Business Media, 
2006.

\bibitem{knodel}
A. Makowski, 
Generalization of Morrow's $D$ numbers, 
\textit{Bull. Belg. Math. Soc. Simon Stevin} \textbf{36} (1962), 
71.

\bibitem{pouletsarrus}
D. H. Lehmer, 
Errata for Poulet's table, 
\textit{Math. Comp.} \textbf{25} (1971), 
944--945.
 
\end{thebibliography}

\bigskip
\hrule
\bigskip

\noindent 2010 {\it Mathematics Subject Classification}:
Primary 11Bxx; Secondary 11Y55, 11A15.

\noindent \emph{Keywords: } 
Carmichael number, Fermat's little theorem, Binet formula, floor
function, multiple-counting sequence, Fermat pseudoprime, Jacobsthal
sequence.

\bigskip
\hrule
\bigskip

\noindent (Concerned with sequences
\seqnum{A000302},
\seqnum{A001045},
\seqnum{A006880},
\seqnum{A007910},
\seqnum{A077947}, and
\seqnum{A093138}.)

\bigskip
\hrule
\bigskip

\vspace*{+.1in}
\noindent
Received  September 26 2015;
revised versions received  October 23 2015; December 1 2015; December
19 2015.
Published in {\it Journal of Integer Sequences}, January 10 2016.

\bigskip
\hrule
\bigskip

\noindent
Return to
\htmladdnormallink{Journal of Integer Sequences home page}{http://www.cs.uwaterloo.ca/journals/JIS/}.
\vskip .1in


\end{document}

                                                                                

                                                                                


</pre></body></html>