\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 
A Multidimensional Ces\'{a}ro Type Identity and Applications \\
}
\vskip 1cm
\large
Olivier Bordell\`es\\
2, all\'{e}e de la Combe\\
43000 Aiguilhe\\
France\\
\href{mailto:borde43@wanadoo.fr}{\tt borde43@wanadoo.fr}
\end{center}

\vskip .2 in

\newcommand{\N}{\mathbb {N}}
\newcommand{\Z}{\mathbb {Z}}
\newcommand{\Q}{\mathbb {Q}}
\newcommand{\R}{\mathbb {R}}
\newcommand{\C}{\mathbb {C}}
\newcommand{\K}{\mathbb {K}}


\begin{abstract}
We generalize a well-known identity due to Ces\'{a}ro and get some asymptotics on some multivariate gcd-sum functions. Several examples when the functions are completely multiplicative are given.
\end{abstract}


\section{Introduction}
\label{s1}

In 1885, Ces\'{a}ro \cite{ces} proved the following identity.
\begin{lemma}
\label{le1}
Let $f$ be an arithmetic function and $\varphi$ denote Euler's totient function. For every positive integer $n$, we have
$$\sum_{i=1}^{n} f(\gcd(i,n)) = (f \star \varphi)(n).$$
\end{lemma}


This follows from
$$\sum_{i=1}^{n} f(\gcd(i,n)) = \sum_{d \mid n} f(d) \sum_{\substack{k \leqslant n/d \\ \gcd(k,n/d)=1}} 1 = \sum_{d \mid n} f(d) \varphi \left ( \frac{n}{d} \right ) = (f \star \varphi)(n)$$
where $(f \star g)(n) = \sum_{d \mid n} f(d) g(n/d)$ is the usual Dirichlet convolution product. It should be mentioned that such an identity also occurs with some other convolution products where the summation is over some subset of the set of the divisors of $n$. For instance, T\'{o}th \cite{tot} showed that
$$\sum_{i \in \mathrm{Reg}(n)} f(\gcd(i,n)) = (f \odot \varphi)(n)$$
where $f \odot g$ is the unitary convolution product and where the notation $i \in \mathrm{Reg}(n)$ means that $1 \leqslant i \leqslant n$ and there exists an integer $x$ such that $i^2x \equiv i \pmod n$.


In this work, we intend to generalize Lemma~\ref{le1} to higher dimensions and present some applications to various common number-theoretic functions.

The following result is easily deduced in a similar fashion to Lemma~\ref{le1}.

\begin{proposition}
\label{pro2}
Let $r \in \Z_{\geqslant 1}$ and $f : \left( \Z_{\geqslant 1} \right)^r \longrightarrow \C$ be an arithmetic function of $r$ variables. For any $n \in \Z_{\geqslant 1}$
$$\sum_{i_1=1}^n \dotsb \sum_{i_r=1}^n f \left( \gcd(i_1,n),\dotsc,\gcd(i_r,n) \right) = \sum_{d_1 \mid n} \dotsb \sum_{d_r \mid n} f \left( d_1,\dotsc,d_r \right) \varphi(n/d_1) \dotsb \varphi(n/d_r).$$
\end{proposition}

\begin{proof}
Similarly as in Lemma~\ref{le1}, we have
\begin{eqnarray*}
   & & \sum_{i_1=1}^n \dotsb \sum_{i_r=1}^n f \left( \gcd(i_1,n),\dotsc,\gcd(i_r,n) \right) \\
   &=& \sum_{d_1 \mid n} \dotsb \sum_{d_r \mid n} f \left( d_1,\dotsc,d_r \right) \sum_{\substack{k_1 \leq n/d_1 \\ \gcd(k_1,n/d_1)=1}} \dotsb \sum_{\substack{k_r \leq n/d_r \\ \gcd(k_r,n/d_r)=1}} 1\\
    &=& \sum_{d_1 \mid n} \dotsb \sum_{d_r \mid n} f \left( d_1,\dotsc,d_r \right) \varphi(n/d_1) \dotsb \varphi(n/d_r)
\end{eqnarray*}
as announced.
\end{proof}

\section{Applications}
\label{s2}

\subsection{Main result}
\label{ss1}

\begin{theorem}
\label{t3}
Let $r \in \Z_{\geqslant 1}$ and $f : \Z_{\geqslant 1} \longrightarrow \C$ be a completely multiplicative function such that, for any prime number $p$, $|f(p)| \leqslant 1$. Set $\theta := \underset{p \, \mathrm{prime}}{\max} \left| 1 - f(p) \right|$ and assume $\theta r \geqslant 1$. Then, for any $x$ sufficiently large
$$\sum_{n \leqslant x} \sum_{i_1=1}^n \dotsb \sum_{i_r=1}^n f \left( \gcd(i_1,n) \right) \dotsb f \left( \gcd(i_r,n) \right) = C_r(f) \, x^{r+1} + O \left( \left( x \log^{\theta} x \right)^r \right)$$
where 
\begin{equation}
C_r(f) := \frac{1}{r+1} \prod_p \left( 1 - \frac{1}{p} \right) \left( 1 + \sum_{\alpha = 1}^\infty \frac{1}{p^\alpha} \left( 1 - \left ( 1 - \left( \frac{f(p)}{p} \right)^\alpha \right)\frac{1-f(p)}{p-f(p)} \right)^r \right). \label{e1}
\end{equation}
\end{theorem} 

This result immediately yields the following examples.

\begin{corollary}
\label{cor4}
Let $r \in \Z_{\geqslant 1}$. 
\begin{enumerate}
   \item[$1.$] For any $x$ sufficiently large
   $$\sum_{n \leqslant x} \sum_{i_1=1}^n \dotsb \sum_{i_r=1}^n \frac{1}{\gcd(i_1,n)\dotsb\gcd(i_r,n)} = C_r  x^{r+1} + O \left( (x \log x)^r \right),$$
where
$$C_r = \frac{1}{r+1} \prod_p \left( 1 - \frac{1}{p} \right) \left( 1 + \sum_{\alpha = 1}^\infty \frac{1}{p^\alpha} \left( 1 - \frac{1-p^{-2\alpha}}{p+1} \right)^r \right).$$
   \item[$2.$] Let $\lambda$ be the Liouville function. Then for any $x$ sufficiently large
   $$\sum_{n \leqslant x} \sum_{i_1=1}^n \dotsb \sum_{i_r=1}^n \lambda \left( \gcd(i_1,n) \right) \dotsb \lambda \left( \gcd(i_r,n) \right) = C_r  x^{r+1} + O \left( x^r (\log x)^{2r} \right),$$
where
$$C_r = \frac{1}{r+1} \prod_p \left( 1 - \frac{1}{p} \right) \left( 1 + \sum_{\alpha = 1}^\infty \frac{1}{p^\alpha} \left( 1 - \frac{2 \left ( 1 - (-1)^\alpha p^{-\alpha} \right)}{p+1} \right)^r \right).$$   
   \item[$3.$] Let $\ell$ be an odd prime number, $(\cdot / \ell)$ be the Legendre symbol modulo $\ell$ and define {\rm \cite{bor}}
   $$\Omega_\ell (n) := \sum_{\substack{p \mid n \\ (p/\ell) = -1}} 1 \quad \text{and} \quad \lambda_\ell(n) := (-1)^{\Omega_\ell (n)}.$$
   Then for any $x$ sufficiently large
   $$\sum_{n \leqslant x} \sum_{i_1=1}^n \dotsb \sum_{i_r=1}^n \lambda_\ell \left( \gcd(i_1,n) \right) \dotsb \lambda_\ell \left( \gcd(i_r,n) \right) = C_r  x^{r+1} + O \left( x^r (\log x)^{2r} \right),$$
where
$$C_r = \frac{1}{r+1} \prod_{(p/\ell) = -1} \left( 1 - \frac{1}{p} \right) \left( 1 + \sum_{\alpha = 1}^\infty \frac{1}{p^\alpha} \left( 1 - \frac{2 \left ( 1 - (-1)^\alpha p^{-\alpha} \right)}{p+1} \right)^r \right).$$   
   \item[$4.$] Let $q > 1$ be an integer and $\chi$ be a real non-principal Dirichlet character modulo $q$. Then for any $x$ sufficiently large
$$\sum_{n \leqslant x} \sum_{i_1=1}^n \dotsb \sum_{i_r=1}^n \chi \left( \gcd(i_1,n) \right) \dotsb \chi \left( \gcd(i_r,n) \right) = C_r x^{r+1} + O \left( x^{r} (\log x)^{2r} \right),$$
where 
\begin{eqnarray*}
   C_r &=& \frac{1}{r+1} \prod_{\chi(p)=-1} \left( 1 - \frac{1}{p} \right) \left( 1 + \sum_{\alpha = 1}^\infty \frac{1}{p^\alpha} \left( 1 - \frac{2 \left ( 1 - (-1)^\alpha p^{-\alpha} \right)}{p+1} \right)^r \right) \\
   & & \times \prod_{\chi(p)=0} \left( 1 - \frac{1}{p} \right) \left( 1 + \frac{1}{p} \left( 1 - \frac{1}{p} \right)^{r-1} \right).
\end{eqnarray*}
\end{enumerate}
\end{corollary}


\subsection{Average order of certain multiplicative functions}
\label{ss2}

The proof of Theorem~\ref{t3} uses the following tool, which is of independent interest.

\begin{proposition}
\label{pro5}
Let $F$ be a multiplicative function for which there exists $\kappa \geqslant 1$ such that, for all prime powers $p^\alpha$
$$|(F \star \mu) (p^\alpha)| \leqslant \kappa p^{-\alpha}.$$
Then
$$\sum_{n \leqslant x} F(n) = x \prod_p \left( 1 - \frac{1}{p} \right) \left( 1 + \sum_{\alpha = 1}^\infty \frac{F(p^\alpha)}{p^\alpha} \right) + O \left( (\log x)^{\kappa} \right).$$
\end{proposition}

\begin{proof}
Let $G := F \star \mu$. The function $G$ is multiplicative and satisfies 
$$|G(n)| \leqslant \frac{\kappa^{\omega(n)}}{n} \quad \left(n \in \Z_{\geqslant 1} \right).$$
For all $x \geqslant 1$, we then have
$$\sum_{d \leqslant x} \frac{|G(d)|}{d} \leqslant  \sum_{d \leqslant x} \frac{\kappa^{\omega(d)}}{d^2} \ll 1,$$
so that the series $\sum_{d \geqslant 1} G(d)/d$ is absolutely convergent and one may check that
$$1 + \sum_{\alpha=1}^\infty \frac{G(p^\alpha)}{p^\alpha} = \left( 1 - \frac{1}{p} \right) \left( 1 + \sum_{\alpha = 1}^\infty \frac{F(p^\alpha)}{p^\alpha} \right) \quad \left( p \ \textrm{prime} \right).$$
Hence
\begin{eqnarray*}
   \sum_{n \leqslant x} F(n) &=& \sum_{n \leqslant x} (G \star \mathbf{1})(n) = \sum_{d \leqslant x} G(d) \left \lfloor \frac{x}{d} \right \rfloor \\
   &=& x \sum_{d \leqslant x} \frac{G(d)}{d} + O \left( \sum_{d \leqslant x} \left | G(d) \right | \right) \\
   &=& x \sum_{d=1}^{\infty} \frac{G(d)}{d} + O \left( \sum_{d \leqslant x} \left | G(d) \right | + x \sum_{d > x} \frac{\left | G(d) \right |}{d} \right) \\
   &=& x \prod_p \left( 1 - \frac{1}{p} \right) \left( 1 + \sum_{\alpha = 1}^\infty \frac{f(p^\alpha)}{p^\alpha} \right) + O \left( R(x)\right)
\end{eqnarray*}
say, and by partial summation, we finally get
\begin{eqnarray*}
   R(x) & \ll & \sum_{d \leqslant x} |G(d)| + x \int_x^\infty \left( \sum_{d \leqslant t} |G(d)| \right) \frac{\textrm{d}t}{t^2} \\
   & \ll & \sum_{d \leqslant x} \frac{\kappa^{\omega(d)}}{d} + x \int_x^\infty \left( \sum_{d \leqslant t} \frac{\kappa^{\omega(d)}}{d} \right) \frac{\textrm{d}t}{t^2} \\
   & \ll & (\log x)^{\kappa} + x \int_x^\infty \left( \log t \right)^{\kappa} \, \frac{\textrm{d}t}{t^2} \ll (\log x)^{\kappa} 
\end{eqnarray*}
as required.
\end{proof}

\subsection{Proof of Theorem~\ref{t3}}
\label{ss3}

Let $f$ satisfy the hypotheses of Theorem~\ref{t3} and define
$$F_r (n) := \left( \frac{(\varphi \star f)(n)}{n} \right)^r \quad \left( n,r \in \Z_{\geqslant 1} \right).$$
We start with the following lemma.

\begin{lemma}
\label{le6}
Let $r \in \Z_{\geqslant 1}$ and $f : \Z_{\geqslant 1} \longrightarrow \C$ be a completely multiplicative function such that, for any prime $p$, $|f(p)| \leqslant 1$. Set $\theta := \underset{p \, \mathrm{prime}}{\max} \left| 1 - f(p) \right|$. For any prime powers $p^\alpha$ and any $r \in \Z_{\geqslant 1}$
$$\left | \left( F_r \star \mu \right) \left( p^\alpha \right) \right | \leqslant \frac{\theta r}{p^\alpha}.$$
\end{lemma}

\begin{proof}
We first note that
\begin{equation}
   \left | 1 - \left( 1 - \left( \frac{f(p)}{p} \right)^\alpha \, \right)\frac{1-f(p)}{p-f(p)} \right | \leqslant 1 \quad \left( p \ \mathrm{prime}, \ \alpha \in \Z_{\geqslant 0} \right). \label{e2}
\end{equation}
Now for any prime powers $p^\alpha$, straightforward computations show that
$$F_r \left( p^\alpha \right) = \left( 1 - \left( 1 - \left( \frac{f(p)}{p} \right)^\alpha \right)\frac{1-f(p)}{p-f(p)} \right)^r$$
so that
\begin{eqnarray*}
   & & \left | \left( F_r \star \mu \right) \left( p^\alpha \right) \right | \\
   &=& \left| \left( 1 - \left( 1 - \left( \frac{f(p)}{p} \right)^\alpha \, \right)\frac{1-f(p)}{p-f(p)} \right)^r - \left( 1 - \left( 1 - \left( \frac{f(p)}{p} \right)^{\alpha-1} \right)\frac{1-f(p)}{p-f(p)} \right)^r \right |.
\end{eqnarray*}
Note that
$$\left | (1-a)^r-(1-b)^r \right | = \left | (b-a) \sum_{k=0}^{r-1}(1-a)^k(1-b)^{r-k-1} \right |  \leqslant r \rho^{r-1} |b-a|$$
whenever $|1-a|, \; |1-b| \leqslant \rho$ ($\rho > 0$) so that, using \eqref{e2}
$$\left | \left( F_r \star \mu \right) \left( p^\alpha \right) \right | \leqslant r \left( \frac{|f(p)|}{p} \right)^\alpha \left | 1 - \frac{1}{f(p)} \right | \leqslant  \frac{\theta r}{p^\alpha} \left| f(p) \right|^{\alpha-1} \leqslant \frac{\theta r}{p^\alpha}$$
as asserted. 
\end{proof}

We now are in a position to prove Theorem~\ref{t3}.

\begin{proof}[Proof of Theorem~\ref{t3}]
From Lemma~\ref{le6}, we may apply Proposition~\ref{pro5} to the function $F_r$ and get
$$\sum_{n \leqslant x} F_r(n) = (r+1) \, C_r(f) \, x + O \left( (\log x)^{\theta r} \right)$$
where $C_r(f)$ is defined in \eqref{e1}. Now from Proposition~\ref{pro2}, we get
$$\sum_{n \leqslant x} \sum_{i_1=1}^n \dotsb \sum_{i_r=1}^n f\left( \gcd(i_1,n) \right) \dotsb f \left( \gcd(i_r,n) \right) = \sum_{n \leqslant x} \left( \varphi \star f \right)^r (n) = \sum_{n \leqslant x} n^r F_r(n)$$
and we conclude the proof using partial summation.
\end{proof}

\section{Acknowledgments}
 I express my gratitude to the anonymous referee for his careful reading of the manuscript and the many valuable suggestions and corrections he made.
 

\begin{thebibliography}{9}
   \bibitem{bor} P. Borwein, S. K. K. Choi, and M. Coons, Completely multiplicative functions taking values in $\{-1,1\}$, {\it Trans. Amer. Math. Soc.} {\bf 362} (2010), 6279--6291. 
   \bibitem{ces} E. Ces\'{a}ro, \'{E}tude moyenne du plus grand commun diviseur de deux nombres, {\it Ann. Mat. Pura Appl.} {\bf 13} (1885), 235--250. 
   \bibitem{tot} L. T\'{o}th, A gcd-sum function over regular integers modulo $n$, \textit{J. Integer Sequences} {\bf 12} (2009), 
   \href{https://cs.uwaterloo.ca/journals/JIS/VOL12/Toth/toth3.html}{Article 09.2.5}. 
\end{thebibliography}

\bigskip
\hrule
\bigskip

\noindent 2010 {\it Mathematics Subject Classification}:
Primary 11A25; Secondary 15A15, 15A18, 11C20.

\noindent \emph{Keyword: } Multivariate gcd-sum function, Ces\'{a}ro's identity.

\bigskip
\hrule
\bigskip

\noindent
(Concerned with sequence 
\seqnum{A018804}.)

\bigskip
\hrule
\bigskip

\vspace*{+.1in}
\noindent
Received November 30 2014; revised versions received January 29 2015; 
February 11 2015.
Published in {\it Journal of Integer Sequences}, February 15 2015.

\bigskip
\hrule
\bigskip

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


\end{document}