[Tutorial] GCD Convolution

Правка en1, от szaranczuk, 2023-02-04 20:59:01

On today's POI training camp round I have learnt a nice technique that could possibly be useful in some number theory problems. I couldn't find any CF article on it, so I think it's fair enough to share it on my own.

Remark on used notation

In some sums I will use an Iverson notation

Problem: Squarefree Function

Let's define a Squarefree Function $$$b(x)$$$ for any positive integer $$$x$$$ as $$$x$$$ divided by a greatest perfect square, that divides $$$x$$$.

For example: $$$b(1) = 1$$$, $$$b(2) = 2$$$, $$$b(4) = 1$$$, $$$b(6) = 6$$$, $$$b(27) = 3$$$, $$$b(54) = 6$$$, $$$b(800) = 2$$$

Given an array $$$a$$$ of $$$n \leq 10^5$$$ positive integers, where each $$$a_i \leq 10^5$$$ compute sum

\begin{gather} \sum\limits_{1 \leq i < j \leq n}b(a_i\cdot a_j) \end{gather}

Technique: GCD Convolution

You might probably heard about a sum convolution. For two arrays $$$b$$$, and $$$c$$$, it is defined as an array $$$d$$$ such that \begin{gather} d_k = \sum\limits_{i,j}a_i\cdot a_j[i + j = k] \end{gather} If not, it's basically the same thing as a polynomial multiplication. If $$$B(x) = b_0 + b_1x + b_2x^2 + ... + b_nx^n$$$, and $$$C(x) = c_0 + c_1x + c_2x^2 + ... + c_mx^m$$$, then $$$(B \cdot C)(x) = d_0 + d_1x + d_2x^2 + ... + d_{n+m}x^{n+m}$$$

Теги number theory, gcd, convolution, tutorial

История

 
 
 
 
Правки
 
 
  Rev. Язык Кто Когда Δ Комментарий
en16 Английский szaranczuk 2023-02-05 16:37:54 33
en15 Английский szaranczuk 2023-02-05 01:19:42 0 (published)
en14 Английский szaranczuk 2023-02-05 01:17:32 108
en13 Английский szaranczuk 2023-02-05 01:01:55 653
en12 Английский szaranczuk 2023-02-05 00:10:36 15
en11 Английский szaranczuk 2023-02-04 23:51:22 151
en10 Английский szaranczuk 2023-02-04 23:37:19 54
en9 Английский szaranczuk 2023-02-04 23:33:55 7
en8 Английский szaranczuk 2023-02-04 23:31:29 118
en7 Английский szaranczuk 2023-02-04 22:50:30 340
en6 Английский szaranczuk 2023-02-04 22:44:31 682
en5 Английский szaranczuk 2023-02-04 22:15:46 1037
en4 Английский szaranczuk 2023-02-04 21:39:05 627
en3 Английский szaranczuk 2023-02-04 21:18:33 233
en2 Английский szaranczuk 2023-02-04 21:06:12 325
en1 Английский szaranczuk 2023-02-04 20:59:01 1344 Initial revision (saved to drafts)