====== 합성곱 (Convolution) ======

  * [[https://ko.wikipedia.org/wiki/%ED%95%A9%EC%84%B1%EA%B3%B1|위키피디아]]에 있는 합성곱에 대한 요약설명은 "하나의 함수와 또 다른 함수를 반전 이동한 값을 곱한 다음, 구간에 대해 적분하여 새로운 함수를 구하는 수학 연산자" 이다.
  * 이것을 식으로 정의하면, 함수 f와 g의 합성곱 f*g는 $ (f  * g )(t) = \int_{-\infty}^\infty f(\tau) g(t - \tau)\, d\tau $ 으로 쓸수 있다.
  * 하지만 이렇게 연속함수를 대상으로 하는 것은 수학이나 신호처리에서 주로 쓰이는 형태이고, PS에서 필요한 것은 이산함수 또는 수열간의 합성곱이다. 다음처럼 표기할수 있다
    * $ (f  * g)(m) = \sum_n {f(n) g(m - n)} $
    * 수열의 경우는 f(i)를 수열의 i번째 항을 나타내는 이산함수가 된다고 생각하면 된다.
  * 이러한 이산합성곱을 구하는 것은 이산 푸리에 변환을 이용해서 효율적으로 구할수 있다. [[ps:FFT]]를 사용하면 길이가 n인 두개의 수열에 대해서 O(nlogn)에 합성곱을 구할수 있다.
    * 이것의 응용은 [[ps:FFT]]를 참고. 대표적으로는 두개의 다항식을 곱한 결과식의 계수, 큰수의 곱셈, 두 수열의 원소들의 모든 pairwise sum 등등이 있다.
  * 이제 원래의 합성곱에서 중간에 들어가는 연산자를 일반화해서 확장시켜보자.
    * 원래 식을 바꿔쓰면, $(f  * g)(m) = \sum_n{f(n) g(m-n)} = \sum_{i+j=m} {f(i) g(j)} $ 이 되는데, i+j=m 부분에서  +를 다른 연산자로 치환해서 $(f  * g)(m) = \sum_{i★j=m} {f(i) g(j)}$ 으로 확장하는 것이다. 
    * 보통은 사용되는 연산자 이름을 넣어서 XXX convollution 이라고 부른다. 예를 들어 ★ 에 xor 연산자가 들어가는 경우는 xor convolution, GCD 가 들어가는 경우는 GCD convolution 등으로 부른다. (PS분야 밖에서 이렇게 부르는 것을 잘 보지는 못했다.. PS에서만 쓰이는 용어일수도 있다)
  * 비트 연산자를 사용하는 컨볼루션 (xor/or/and convolution) 의 경우는 Walsh-Hadamard transform을 이용해서 효율적으로 구할수 있다.
    * [[ps:FWHT]]를 사용하면 길이가 n인 두개의 수열에 대해서 O(nlogn)에 convolution을 구할수 있다.
    * FFT와 같은 방식으로 이것을 활용하면 두 수열의 원소들에 대한 모든 pairwise xor/or/and 의 합을 구할 수 있다
  * 연산 결과가, 어떤 Poset에서의 least common ancestor나 highest common descendent 로 정의될수 있는 연산자에 대해서는,  Zeta transform / Mobius transform 을 이용해서 효율적으로 구할수 있다.
    * or, and, gcd, lcm 등의 연산자가 여기에 해당된다
    * SOS DP 를 이용해서 zeta/mobius transform을 계산하면, 길이가 n인 두개의 수열에 대해서 O(nlogn)에 convolution을 구할수 있다.