[논문 리뷰] Local Euler-Maclaurin formula for polytopes

Nicole Berline, Michèle Vergne|arXiv (Cornell University)|2005. 07. 13.

Pharmacological Effects of Medicinal Plants참고 문헌 14인용 수 8

한 줄 요약

이 논문은 유리 볼록 다면체에 대한 국소적이고 계산 가능한 오일러-매클로린 공식을 제시하며, 정수 점 위의 이산 합을 면 위의 다항함수에 작용하는 미분 연산자의 합으로 표현한다. 주요 기여는 각 면에서의 횡단 원뿔에만 의존하는 유리수 계수를 가진 미분 연산자 $ D(\frak{p}, \frak{f}) $ 를 구성한 것으로, 고정된 차원에서 다항식 시간 알고리즘을 통해 에르하르트 쿼asi다항식과 이산 합을 효율적으로 계산할 수 있게 한다.

ABSTRACT

We give a local Euler-Maclaurin formula for rational convex polytopes in a rational euclidean space . For every affine rational polyhedral cone C in a rational euclidean space W, we construct a differential operator of infinite order D(C) on W with constant rational coefficients, which is unchanged when C is translated by an integral vector. Then for every convex rational polytope P in a rational euclidean space V and every polynomial function f (x) on V, the sum of the values of f(x) at the integral points of P is equal to the sum, for all faces F of P, of the integral over F of the function D(N(F)).f, where we denote by N(F) the normal cone to P along F.

연구 동기 및 목표

• 유리 볼록 다면체에 대한 국소적이고 계산 가능한 오일러-매클로린 공식을 구축하여, 이러한 연산자가 존재한다는 오랫동안 남아있던 추측을 해결한다.
• 다항함수 $ h $ 에 대해 $ \text{Card}(\frak{p} \bigcap \bbZ^d) $ 와 $ \text{sum}_{x \bigcap \frak{p} \bigcap \bbZ^d} h(x) $ 를 계산하는 방법을 제공하며, 차원과 차수를 고정할 경우 입력 크기와 다항식 시간 내에 실행 가능하도록 한다.
• 바빈ок의 다각형 내 정수 점 수를 세는 다항식 시간 알고리즘을, 국소적이고 유리수 계수를 가진 미분 연산자를 사용하여 격자 점 위의 가중치 합으로 확장한다.
• 등변 Riemann-Roch 정리에 기반하여 토릭 다양체의 토드 클래스에 대한 표준적이고 국소적인 공식을 수립한다.
• 국소 공식의 구조적이고 알고리즘적인 실현을 제공하여, 작은 차원에서 에르하르트 쿼اسي다항식과 이산 합을 실질적으로 계산할 수 있도록 한다.

제안 방법

• 각 유리 애파인 원뿔 $ \frak{a} $ 와 관련된 생성함수 $ \nu(\frak{a})(\theta) $ 를 정의하며, 이는 $ \frak{a} $ 내의 격자 점 위의 지수 합에 기반하고, 격자 이동에 의해 정규화된다.
• 미분 연산자 $ D(\frak{p}, \frak{f}) $ 의 기호를 함수 $ \nu(\frak{t}(\frak{p}, \frak{f})) $ 를 통해 구성하며, 여기서 $ \frak{t}(\frak{p}, \frak{f}) $ 는 면 $ \frak{f} $ 沿해 $ \frak{p} $ 의 횡단 원뿔이며, $ \bbR^d $ 에서의 스칼라 곱을 사용한다.
• 함수 $ \nu(\frak{a}) $ 가 $ \nu(\frak{a}_1 \bigcup \frak{a}_2) = \nu(\frak{a}_1) + \nu(\frak{a}_2) - \nu(\frak{a}_1 \bigcap \frak{a}_2) $ 를 만족하는 치환성 성질을 이용하여 원뿔 분해에서의 일致성을 확보한다.
• 지수 생성함수와 하위 차원 면에 대한 적분을 통한 공식을 이용해 $ D(\frak{p}, \frak{f}) $ 를 재귀적으로 정의함으로써 국소성과 유리수 계수를 보장한다.
• 바빈ок의 원뿔을 단순원뿔로 부호 있는 분해하는 방법을 사용하여 알고리즘을 구현함으로써, 차원과 차수가 고정된 경우 $ \text{sum}_{x \bigcap \frak{p} \bigcap \bbZ^d} h(x) $ 를 다항식 시간 내에 계산할 수 있다.
• 각 면의 기여를 합산하여 에르하르트 쿼اسي다항식을 계산하기 위해 방법을 적용하며, 계수는 모듈로 함수 $ \text{mod}(at, m) $ 로 표현된다.

실험 결과

연구 질문

• RQ1 유리 다면체에 대해, 각 면에서의 횡단 원뿔에만 의존하는 국소적이고 계산 가능한 오일러-매클로린 공식을 구성할 수 있는가?
• RQ2 국소성과 계산 가능성 조건을 만족하는 표준적이고 유리수 계수를 가진 미분 연산자 $ D(\frak{p}, \frak{f}) $ 가 존재하는가? 즉, $ \text{sum}_{x \bigcap \frak{p} \bigcap \bbZ^d} h(x) = \text{sum}_{\frak{f} \bigcap \frak{F}(\frak{p})} \text{int}_{\frak{f}} D(\frak{p}, \frak{f}) \bullet h $ 가 성립하는가?
• RQ3 국소적 미분 연산자를 기반으로 한 면별 분해를 통해 유리 다면체의 에르하르트 쿼اسي다항식을 효율적으로 계산할 수 있는가?
• RQ4 이러한 연산자를 사용하여 국소적이고 알고리즘적인 방법으로 토릭 다양체의 등변 Riemann-Roch 정리를 실현할 수 있는가?
• RQ5 다항함수 $ h $ 에 대해 $ \text{sum}_{x \bigcap \frak{p} \bigcap \bbZ^d} h(x) $ 를 계산하는 데 필요한 계산 복잡도는 무엇이며, 차원과 차수가 고정된 경우 다항식 시간 내에 수행할 수 있는가?

주요 결과

• 논문은 유리수 계수를 가진 미분 연산자 $ D(\frak{p}, \frak{f}) $ 를 구성하여 오일러-매클로린 공식 (1) 을 만족시키며, 국소성과 계산 가능성까지 확보함으로써 [3] 에서 제기된 추측을 해결한다.
• 연산자 $ D(\frak{p}, \frak{f}) $ 는 격자 이동에 대해 불변하는 횡단 원뿔 $ \frak{t}(\frak{p}, \frak{f}) $ 에만 의존하며, 이는 국소성과 격자 보존 등급변환에 대한 불변성을 보장한다.
• 상수 함수 $ h(x) = 1 $ 인 경우 공식은 $ \text{Card}(\frak{p} \bigcap \bbZ^d) = \text{sum}_{\frak{f}} \nu_0(\frak{p}, \frak{f}) \text{vol}(\frak{f}) $ 로 단순화되며, 유리수 계수 $ \nu_0(\frak{p}, \frak{f}) $ 를 가지며, Danilov, Morelli, McMullen 의 결과를 확인한다.
• 이 방법을 통해 유리 다면체의 에르하르트 쿼اسي다항식을 각 면의 기여를 합산하여 계산할 수 있으며, 계수는 모듈로 함수 $ \text{mod}(at, m) $ 로 표현된다. 이는 예제 36 에서 확인된다.
• 차원 $ d $ 와 함수 $ h $ 의 차수가 고정된 경우, 합 $ \text{sum}_{x \bigcap \frak{p} \bigcap \bbZ^d} h(x) $ 는 다항식 시간 내에 계산 가능하며, 바빈ок의 결과를 확장한다.
• 이 방법은 치수 11^5 의 확대된 삼각형에서 $ x_1^{48}x_2^{48} $ 의 합을 계산하는 데 성공적으로 적용되었으며, 약 한 시간 내에 300자리 숫자를 도출하여 작은 차원에서 기존 소프트웨어를 능가했다.