数论杂题笔记

CF338D GCD Table

先考虑如果会出现，那么会出现在哪一行，因为行数是所有 $a$ 的倍数，显然直接 $lcm(a)$ 最优。

a[i] 的列数设为 $y+i$ ，那么有 $a[i]|y+i$ 那么就有 $n $个同余方程，解出来最小的 $y$ 显然最优。最后验算一下就做完了。

51nod 1042

$$\begin{array}{rl}ans& =\sum _{i=a}^{b}ib/\gcd (i,b)\\ & =\sum _{i=a}^{b}ib\sum _{t}1/t[\gcd (i,b)=t]\\ & =\sum _{i=a}^{b}ib\sum _{t|i,t|b}f(t)\\ & =b\sum _{t|b}f(t)\sum _{t|i}i\end{array}$$

${\sum }_{t|i}i$ 就很好求，就是 $a$ 到 $b$ 中 $t$ 的倍数的和。

思考我们临时设的函数 $f$ ，他要满足 ${\sum }_{t|s}f(t)=1/s$

这个其实就挺简单了，当然可以直接大力算。

最后枚举 $t$ 加答案即可。

51nod 1227

$$\begin{array}{rl}f(n)& =\frac{1}{n}\sum _{i=1}^{n}lcm(i,n)\\ & =\sum _{d|n}\sum _{i}i/d[\gcd (i,n)=d]\\ & =\sum _{d|n}\sum _{i=1}^{\lfloor n/d\rfloor }i[\gcd (id,n)=d]\\ & =\sum _{d|n}\sum _{i=1}^{\lfloor n/d\rfloor }i[\gcd (i,n/d)=1]\\ & 因为有\sum _{i=1}^{n}i[\gcd (i,n)=1]=\frac{n\phi (n)+[n=1]}{2}\\ f(n)& =\frac{1}{2}(1+\sum _{d|n}d\phi (d))\end{array}$$

求前缀和的话就是 $ans={\sum }_{d=1}^{n}d\phi (d)\lfloor n/d\rfloor$。对求 $d\phi (d)$ 前缀和就是杜教筛模板了。

project euler 530

$$\begin{array}{rl}ans& =\sum _{i=1}^n\sum _{d|i}\gcd (d,i/d)\\ & =\sum _{i=1}^n\sum _{d|i}\sum _{h|d,i/d}\phi (h)\\ & =\sum _h\phi (h)\sum _{a,b\ge 1}[ab\le n/h^2]\end{array}$$

暴力枚举 $h$ ，对于每个 $h$ 右边根号算。然后就做完了，时间复杂度不会算，反正它跑了519366769次，快得一p。

51nod 1847

（对这道题真的很无语

（反正是我蒻就对了2333

给出一个结论：若 $i$ , $j$ 不互质 ，则 $sgcd(i,j)=\gcd (i,j)/minp(gcd(i,j))$ 。

那么答案就是

${\sum }_{i=2}^n{(\frac{i}{minp(i)})}^k(2{\sum }_{i=1}^{n/i}\phi (i)-1)$

后面杜教筛，前面min25筛，结束。

cometoj 1098

”即便它再像一个能筛的函数，也要先推它的定义式“。——未来哲学家yxy

若 $x=\prod p_i^{e_i}$
$$\begin{array}{rl}f(x)& =\prod p_i^{\lfloor e^i/2\rfloor }\\ & =\frac{x}{\prod p_i^{\lceil e^i/2\rceil }}\\ & =\sum _{j=1}^x[\prod _{i=1}^k{p}_i^{\lceil e^i/2\rceil }|j]\\ & =\sum _{j=1}^x[x|j^2]\end{array}$$
第二部到第三步用了的除法的定义 $\frac{a}{b}={\sum }_{i=1}^a[b|i]$ ，意思是 $a$ 里面有多少个 $b$ 的倍数。

最后一步意思是当 $a|b$ ，当且仅当 $a^2|b^2$ （就**离谱。

考虑求前缀和，我们知道 $f(x)$ 是积性函数，且 $f(p)=1,p\in Prime$ 。

powerful number筛，结束。