تابع فی اویلر¶
تابع فی اویلر، که با نام تابع فی $\phi (n)$ نیز شناخته میشود، تعداد اعداد صحیح در بازه ۱ تا $n$ (شامل خود ۱ و $n$) را میشمارد که نسبت به $n$ اول هستند. دو عدد را نسبت به هم اول (متباین) میگوییم هرگاه بزرگترین مقسومعلیه مشترک آنها برابر با ۱ باشد ($۱$ نسبت به هر عددی اول در نظر گرفته میشود).
در ادامه مقادیر $\phi(n)$ برای چند عدد صحیح مثبت اول آمده است:
ویژگیها¶
ویژگیهای زیر برای محاسبه تابع فی اویلر برای هر عددی کافی هستند:
- اگر $p$ یک عدد اول باشد، آنگاه برای تمام $q$ های در بازه $1 \le q < p$ داریم $\gcd(p, q) = 1$. بنابراین:
- اگر $p$ یک عدد اول و $k \ge 1$ باشد، آنگاه دقیقاً $p^k / p$ عدد بین ۱ و $p^k$ وجود دارد که بر $p$ بخشپذیر هستند. این به ما نتیجه زیر را میدهد:
-
اگر $a$ و $b$ نسبت به هم اول باشند، آنگاه:
$$\phi(a b) = \phi(a) \cdot \phi(b).$$اثبات این رابطه بدیهی نیست و از قضیه باقیمانده چینی نتیجه میشود. قضیه باقیمانده چینی تضمین میکند که برای هر $0 \le x < a$ و هر $0 \le y < b$، یک $z$ یکتای $0 \le z < a b$ وجود دارد که $z \equiv x \pmod{a}$ و $z \equiv y \pmod{b}$ است. به راحتی میتوان نشان داد که $z$ نسبت به $ab$ اول است اگر و تنها اگر $x$ نسبت به $a$ و $y$ نسبت به $b$ اول باشند. بنابراین، تعداد اعداد صحیح اول نسبت به $ab$ برابر با حاصلضرب تعداد اعداد اول نسبت به $a$ و $b$ است.
-
به طور کلی، برای $a$ و $b$ که نسبت به هم اول نیستند، معادله زیر برقرار است:
$$\phi(ab) = \phi(a) \cdot \phi(b) \cdot \dfrac{d}{\phi(d)}$$که در آن $d = \gcd(a, b)$ است.
بنابراین، با استفاده از سه ویژگی اول، میتوانیم $\phi(n)$ را از طریق تجزیه $n$ به عوامل اول (تجزیه $n$ به حاصلضرب عوامل اولش) محاسبه کنیم. اگر $n = {p_1}^{a_1} \cdot {p_2}^{a_2} \cdots {p_k}^{a_k}$ که در آن $p_i$ ها عوامل اول $n$ هستند،
پیادهسازی¶
در اینجا پیادهسازی با استفاده از تجزیه به عوامل اول با پیچیدگی $O(\sqrt{n})$ آمده است:
int phi(int n) {
int result = n;
for (int i = 2; i * i <= n; i++) {
if (n % i == 0) {
while (n % i == 0)
n /= i;
result -= result / i;
}
}
if (n > 1)
result -= result / n;
return result;
}
محاسبه تابع فی اویلر از ۱ تا $n$ در $O(n \log\log{n})$¶
اگر به مقدار فی برای تمام اعداد بین ۱ تا $n$ نیاز داشته باشیم، تجزیه تمام $n$ عدد کارآمد نیست. میتوانیم از ایدهای مشابه غربال اراتوستن استفاده کنیم. این روش هنوز بر اساس ویژگی نشان داده شده در بالا است، اما به جای بهروزرسانی نتیجه موقت برای هر عامل اول برای هر عدد، تمام اعداد اول را پیدا میکنیم و برای هر کدام، نتایج موقت تمام اعدادی را که بر آن عدد اول بخشپذیر هستند، بهروزرسانی میکنیم.
از آنجایی که این رویکرد اساساً با غربال اراتوستن یکسان است، پیچیدگی زمانی آن نیز یکسان خواهد بود: $O(n \log \log n)$
void phi_1_to_n(int n) {
vector<int> phi(n + 1);
for (int i = 0; i <= n; i++)
phi[i] = i;
for (int i = 2; i <= n; i++) {
if (phi[i] == i) {
for (int j = i; j <= n; j += i)
phi[j] -= phi[j] / i;
}
}
}
ویژگی مجموع مقسومعلیهها¶
این ویژگی جالب توسط گاوس اثبات شده است:
در اینجا مجموع روی تمام مقسومعلیههای مثبت $d$ از $n$ است.
به عنوان مثال، مقسومعلیههای 10 عبارتند از ۱، ۲، ۵ و ۱۰. از این رو $\phi{(1)} + \phi{(2)} + \phi{(5)} + \phi{(10)} = 1 + 1 + 4 + 4 = 10$.
یافتن فی از ۱ تا $n$ با استفاده از ویژگی مجموع مقسومعلیهها¶
ویژگی مجموع مقسومعلیهها همچنین به ما امکان محاسبه فی تمام اعداد بین ۱ تا $n$ را میدهد. این پیادهسازی کمی سادهتر از پیادهسازی قبلی مبتنی بر غربال اراتوستن است، اما پیچیدگی زمانی آن کمی بدتر است: $O(n \log n)$
void phi_1_to_n(int n) {
vector<int> phi(n + 1);
phi[0] = 0;
phi[1] = 1;
for (int i = 2; i <= n; i++)
phi[i] = i - 1;
for (int i = 2; i <= n; i++)
for (int j = 2 * i; j <= n; j += i)
phi[j] -= phi[i];
}
کاربرد در قضیه اویلر¶
مشهورترین و مهمترین ویژگی تابع فی اویلر در قضیه اویلر بیان شده است:
در حالت خاصی که $m$ اول باشد، قضیه اویلر به قضیه کوچک فرما تبدیل میشود:
قضیه اویلر و تابع فی اویلر اغلب در کاربردهای عملی ظاهر میشوند، به عنوان مثال هر دو برای محاسبه وارون ضربی به پیمانه استفاده میشوند.
به عنوان یک نتیجه فوری، به همارزی زیر نیز میرسیم:
این رابطه امکان محاسبه $x^n \bmod m$ را برای $n$ های بسیار بزرگ فراهم میکند، به خصوص اگر $n$ نتیجه یک محاسبه دیگر باشد، زیرا به ما امکان میدهد $n$ را تحت یک پیمانه محاسبه کنیم.
نظریه گروهها¶
$\phi(n)$ برابر است با مرتبه گروه ضربی به پیمانه n $(\mathbb Z / n\mathbb Z)^\times$، یعنی گروه یکهها (عناصری که وارون ضربی دارند). عناصری که وارون ضربی دارند دقیقاً همانهایی هستند که نسبت به $n$ اول هستند.
مرتبه ضربی یک عنصر $a$ به پیمانه $n$ که با $\operatorname{ord}_n(a)$ نمایش داده میشود، کوچکترین $k>0$ است به طوری که $a^k \equiv 1 \pmod m$. $\operatorname{ord}_n(a)$ اندازه زیرگروه تولید شده توسط $a$ است، بنابراین طبق قضیه لاگرانژ، مرتبه ضربی هر عنصر $a$ باید $\phi(n)$ را بشمارد. اگر مرتبه ضربی $a$ برابر با $\phi(n)$ باشد، که بزرگترین مقدار ممکن است، آنگاه $a$ یک ریشه اولیه است و گروه طبق تعریف، دوری است.
تعمیم¶
نسخه کمتر شناخته شدهای از همارزی آخر وجود دارد که امکان محاسبه کارآمد $x^n \bmod m$ را برای $x$ و $m$ که نسبت به هم اول نیستند، فراهم میکند. برای $x$ و $m$ دلخواه و $n \geq \log_2 m$:
اثبات:
فرض کنید $p_1, \dots, p_t$ مقسومعلیههای اول مشترک $x$ و $m$ باشند و $k_i$ توان آنها در $m$ باشد. با اینها $a = p_1^{k_1} \dots p_t^{k_t}$ را تعریف میکنیم، که باعث میشود $\frac{m}{a}$ نسبت به $x$ اول باشد. و فرض کنید $k$ کوچکترین عددی باشد که $a$ عدد $x^k$ را بشمارد. با فرض $n \ge k$ میتوانیم بنویسیم:
همارزی بین خط سوم و چهارم از این واقعیت ناشی میشود که $ab \bmod ac = a(b \bmod c)$. در واقع اگر $b = cd + r$ با $r < c$ باشد، آنگاه $ab = acd + ar$ با $ar < ac$.
از آنجایی که $x$ و $\frac{m}{a}$ نسبت به هم اول هستند، میتوانیم قضیه اویلر را به کار ببریم و فرمول کارآمد زیر را به دست آوریم (کارآمد است چون $k$ بسیار کوچک است؛ در واقع $k \le \log_2 m$):
به کار بردن این فرمول دشوار است، اما میتوانیم از آن برای تحلیل رفتار $x^n \bmod m$ استفاده کنیم. میتوان دید که دنباله توانها $(x^1 \bmod m, x^2 \bmod m, x^3 \bmod m, \dots)$ پس از $k$ عنصر اول (یا کمتر) وارد یک دور به طول $\phi\left(\frac{m}{a}\right)$ میشود. $\phi\left(\frac{m}{a}\right)$ عدد $\phi(m)$ را میشمارد (زیرا $a$ و $\frac{m}{a}$ نسبت به هم اول هستند و داریم $\phi(a) \cdot \phi\left(\frac{m}{a}\right) = \phi(m)$)، بنابراین میتوانیم بگوییم که دوره تناوب طولی برابر با $\phi(m)$ دارد. و از آنجایی که $\phi(m) \ge \log_2 m \ge k$ است، میتوانیم فرمول مطلوب و بسیار سادهتر زیر را نتیجه بگیریم:
مسائل تمرینی¶
- SPOJ #4141 "Euler Totient Function" [سختی: بسیار آسان]
- UVA #10179 "Irreducible Basic Fractions" [سختی: آسان]
- UVA #10299 "Relatives" [سختی: آسان]
- UVA #11327 "Enumerating Rational Numbers" [سختی: متوسط]
- TIMUS #1673 "Admission to Exam" [سختی: زیاد]
- UVA 10990 - Another New Function
- Codechef - Golu and Sweetness
- SPOJ - LCM Sum
- GYM - Simple Calculations (F)
- UVA 13132 - Laser Mirrors
- SPOJ - GCDEX
- UVA 12995 - Farey Sequence
- SPOJ - Totient in Permutation (easy)
- LOJ - Mathematically Hard
- SPOJ - Totient Extreme
- SPOJ - Playing with GCD
- SPOJ - G Force
- SPOJ - Smallest Inverse Euler Totient Function
- Codeforces - Power Tower
- Kattis - Exponial
- LeetCode - 372. Super Pow
- Codeforces - The Holmes Children
- Codeforces - Small GCD