| For All Solvers
NF杯2024

NF杯2024(P)

 特に断らない限り合同式は 55 を法とする.f(1)=1f(1)=155 の倍数ではない.n2n\geq 2 において,

f(n)=k=1nnkgcd(n,k)=ndn1kngcd(n,k)=dkd=ndn1mn/dgcd(n/d,m)=1m=ndn1mdgcd(d,m)=1m\begin{aligned} f(n) & = \sum_{k=1}^n\frac{nk}{\text{gcd}(n,k)}\\ & = n\sum_{d|n} \sum_{\substack{1\leq k\leq n\\ \text{gcd}(n,k)=d}}\frac{k}{d}\\ &=n\sum_{d|n}\sum_{\substack{1\leq m\leq n/d\\ \text{gcd}(n/d,m)=1}}m\\ &=n\sum_{d|n}\sum_{\substack{1\leq m\leq d\\ \text{gcd}(d,m)=1}}m \end{aligned}

ここで,dd 以下の dd と互いに素な正の整数の個数を φ(d)\varphi(d) とすると,

1mdgcd(d,m)=1m={1(d=1)d2φ(d)(d>1)\sum_{\substack{1\leq m\leq d\\ \text{gcd}(d,m)=1}}m =\begin{cases} 1&(d=1)\\ \dfrac{d}{2}\varphi(d)&(d\gt1) \end{cases}

が成り立つ.実際,d=1d=1 のときは明らかで,d>1d\gt1 のとき,mmdd 以下の dd と互いに素な正の整数ならば dmd-m もそうであるので,dd 以下の dd と互いに素な整数の値の平均は d2\dfrac{d}{2} である.よってそのような整数の総和は d2φ(d)\dfrac{d}{2}\varphi(d) となる.

 恒等写像およびオイラー関数 φ\varphi は乗法的であるので,関数 nφ(n)n\varphi(n) も乗法的である.したがって素数 p1,,pp_1,\ldots,p_{\ell} と正の整数 e1,,ee_1,\ldots,e_{\ell} により n=p1e1pen=p_1^{e_1}\cdots p_\ell^{e_{\ell}} と素因数分解すれば,

f(n)=n(1+d>1,dnd2φ(d))=n(12+dnd2φ(d))=n2(1+dndφ(d))=n2(1+i=1(k=1eipikφ(pik)))=n2(1+i=1(1+pi(pi1)++pi2ei1(pi1)))=n2(1+i=1(1pi+pi2pi2ei1+pi2ei))=n2(1+i=1pi2ei+1+1pi+1)\begin{aligned} f(n)&=n\Bigg(1+\sum_{d\gt1,d|n}\frac{d}{2}\varphi(d)\Bigg)\\ &=n\Bigg(\frac{1}{2}+\sum_{d|n}\frac{d}{2}\varphi(d)\Bigg)\\ &=\frac{n}{2}\Bigg(1+\sum_{d|n}d\varphi(d)\Bigg)\\ &=\frac{n}{2}\Bigg(1+\prod_{i=1}^{\ell}\left(\sum_{k=1}^{e_i}p_i^k\varphi(p_i^k)\right)\Bigg)\\ &=\frac{n}{2}\Bigg(1+\prod_{i=1}^{\ell}(1+p_i(p_i-1)+\cdots+p_i^{2e_i-1}(p_i-1))\Bigg)\\ &=\frac{n}{2}\Bigg(1+\prod_{i=1}^{\ell}(1-p_i+p_i^2-\cdots-p_i^{2e_i-1}+p_i^{2e_i})\Bigg)\\ &=\frac{n}{2}\Bigg(1+\prod_{i=1}^{\ell}\frac{p_i^{2e_i+1}+1}{p_i+1}\Bigg) \end{aligned}

よって,f(n)f(n)55 の倍数となるのは,nn55 の倍数となるときか,

i=1pi2ei+1+1pi+11(i)\prod_{i=1}^{\ell}\frac{p_i^{2e_i+1}+1}{p_i+1}\equiv -1\tag{i}

となるときのいずれかである.N2024N^{2024} の正の約数のうち 55 の倍数であるものの個数は,55 の指数が 11 以上であるものの個数であるので,100100 以下の素数の個数が 2525 であることから,

20252420242025^{24}\cdot 2024

と計算される.以下,nn55 の倍数でないとする.

 nn の素因数 pip_i について, pi2ei+1+1pi+1\dfrac{p_i^{2e_i+1}+1}{p_i+1}55 の倍数となるとすると,pi2ei+11p_i^{2e_i+1}\equiv-1 より pi2(2ei+1)1p_i^{2(2e_i+1)}\equiv1 となる.一方フェルマーの小定理より pi41p_i^4\equiv1 であるので,gcd(2(2ei+1),4)=2\gcd(2(2e_i+1),4)=2 より pi21p_i^2\equiv1 すなわち pi±1p_i\equiv \pm1 でなくてはならない.pi2ei+11p_i^{2e_i+1}\equiv-1 と合わせると,pi2ei+1+1pi+1\dfrac{p_i^{2e_i+1}+1}{p_i+1}55 の倍数となるのは pi1p_i\equiv-1 のときに限られることが分かる.一方 pi1p_i\equiv-1 のとき,

pi2ei+1+1pi+1=1pi+pi2+pi2ei2ei+1\dfrac{p_i^{2e_i+1}+1}{p_i+1}=1-p_i+p_i^2-\cdots+p_i^{2e_i}\equiv2e_i+1 となる.

 よって,100100 以下の素数のうち 55 で割った余りが 1,2,31,2,3 となるものを q1,,qxq_1,\ldots,q_x,余りが 44 となるものを r1,,ryr_1,\ldots,r_y とおき,それらに対応する指数をそれぞれ f1,,fx,g1,,gyf_1,\ldots,f_x,g_1,\ldots,g_y とすると,

i=1xqi2fi+1+1qi+1≢0,i=1yri2gi+1+1ri+1i=1y(2gi+1)\prod_{i=1}^x\frac{q_i^{2f_i+1}+1}{q_i+1}\not\equiv0, \quad\prod_{i=1}^y\frac{r_i^{2g_i+1}+1}{r_i+1}\equiv\prod_{i=1}^y(2g_i+1)

が成り立つ.特に,55 の倍数でない N2024N^{2024} の正の約数 nn のうち,

1i,pirypi2ei+1+1pi+1=(i=1xqi2fi+1+1qi+1)(j=1y1rj2gj+1+1rj+1)≢0\prod_{1\leq i\leq \ell,p_i\neq r_y}\frac{p_i^{2e_i+1}+1}{p_i+1}=\Bigg(\prod_{i=1}^x\frac{q_i^{2f_i+1}+1}{q_i+1}\Bigg)\Bigg(\prod_{j=1}^{y-1}\frac{r_j^{2g_j+1}+1}{r_j+1}\Bigg)\not\equiv0

となるものの個数は,j=1,,y1j=1,\ldots,y-1 に対して 2gj+1≢02g_j+1\not\equiv0 すなわち gj≢2g_j\not\equiv 2 となるような 00 以上 20242024 以下の整数 gjg_j の選び方を考えることで,

2025x1620y12025^x\cdot 1620^{y-1}

と計算できる.この条件が成り立っているとき,(i) が成り立つような gyg_y の選び方は 55 を法としてちょうど 11 つ存在するので,55 の倍数でない N2024N^{2024} の正の約数 nn のうち (i) が成り立つものの個数は,

2025x1620y14052025^x\cdot 1620^{y-1}\cdot405

 100100 以下の素数のうち 55 で割って 44 余るのは,下一桁が 99 である 19,29,59,79,8919,29,59,79,8955 個であるので,x=19,y=5x=19,y=5 となる.従って,N2024N^{2024} の正の約数のうち f(n)f(n)55 の倍数となるものの割合は,

ab=2025242024+20251916204405202525=4055(552024+44)20256=552024+44202555=63252566328125\begin{aligned} \frac{a}{b}&=\frac{2025^{24}\cdot2024+2025^{19}\cdot 1620^4\cdot 405}{2025^{25}}\\ &=\frac{405^5\cdot(5^5\cdot2024+4^4)}{2025^6}\\ &=\frac{5^5\cdot2024+4^4}{2025\cdot 5^5}\\ &=\frac{6325256}{6328125} \end{aligned}

特に解答すべき値は 12653381\bold{12653381} である.

解説YouTube

解説YouTubeが存在しません.