20210903下午

20210903下午

难度上来了，好耶

	T1	T2	T3	T4	T5	T6	T7	T8
预测	100	100	100	100	100	100	0	100
一测	100	100	0	100	0	100	0	100

T1：
先找因数和，再找因数和的因数和，判断是否相等与越界。

T2：
模拟，注意别在这轮把下轮出局的数删了。

T3：
和前几天一样构造二叉树，边造边统计即可，数组开小了寄了，但数据范围在哪呢？怎么回事呢？

T4：
如题目下提示所示，bfs即可。

T5：
矩阵转移即可，设$f_i$为数列第$i$项，$S_i$为前$i$项和。
$$\begin{gathered} \left[\begin{array}{c}{S}_i\\ f_i\\ f_{i-1}\end{array}\right]=\left[\begin{array}{c}{S}_{i-1}+f_i\\ f_i\\ f_{i-1}\end{array}\right]=\left[\begin{array}{c}1\times S_{i-1}+1\times f_{i-1}+1\times f_{i-2}\\ 0\times S_{i-1}+1\times f_{i-1}+1\times f_{i-2}\\ 0\times S_{i-1}+1\times f_{i-1}+0\times f_{i-2}\end{array}\right]=\left[\begin{array}{c}{S}_{i-1}\\ f_{i-1}\\ f_{i-2}\end{array}\right]\times \left[\begin{array}{ccc}1& 1& 1\\ 0& 1& 1\\ 0& 1& 1\end{array}\right]=\cdots =\left[\begin{array}{c}{S}_1\\ f_1\\ f_0\end{array}\right]\times {\left[\begin{array}{ccc}1& 1& 1\\ 0& 1& 1\\ 0& 1& 1\end{array}\right]}^{n-1}= \\ \left[\begin{array}{c}1\\ 1\\ 0\end{array}\right]\times {\left[\begin{array}{ccc}1& 1& 1\\ 0& 1& 1\\ 0& 1& 1\end{array}\right]}^{n-1} \end{gathered}$$
矩阵快速幂即可。考试自信不开long long，他真的。。。怎么这么猪比啊，我真的哭死。
T6：
很有难度的几何。
首先确定思路用总的方案数减共线的方案数。可以考虑出总的减横着竖着的情况为$C_{nm}^3-C_n^3\times m-C_m^3\times n$，接下来就是考虑斜线的情况。

考虑枚举$(1,1)$与$(n+1,m+1)$之间的每一个点，它与原点$(0,0)$便确定了两个点，那在这条线段上的每一点都与这两个点构成三点共线（确定在线段上就可以确定不重不漏），那么这条线段上除两端点有多少个端点呢，答案是$(i,j)$到$(0,0)$之间有$gcd(i,j)-1$个点，也很好证明。易得这条线段从下往上除$(0,0)$外第一个点为$(i/gcd(i,j),j/gcd(i,j))$（即坐标互质的那个点），且从0开始横向每$i/gcd(i,j)$就有一个点，所以点数即为$\frac{i}{i/gcd(i,j)}=gcd(i,j)$，去掉$(i,j)$就是$gcd(i,j)-1$个点。（上图可验证，$gcd(2,4)=2$，$(0,0)-(2,4)$上也正好有$2-1=1$个点）所以每条线段的贡献即为$gcd(i,j)-1$。

然后发现线段可以在图内平移，因为横纵坐标为$(n+1,m+1)$，所以可以纵向平移$n+1-i$位，横向平移$m+1-j$位，总共$(n+1-i)(m+1-j)$条线段，总贡献就是$(n+1-i)(m+1-j)(gcd(i,j)-1)$，注意这样是考虑左下角的情况，右下角的情况是相同的，所以总的情况数为
${\sum }_{i=1}^n{\sum }_{j=1}^m(n+1-i)(m+1-j)(gcd(i,j)-1)$。
总复杂度为$O(n^2)$

优化
前置知识：
$[\cdots ]$表示当$[]$内语句判定为真时为1，否则为0。
欧拉函数$\varphi (n)$：1-n中与n互质的函数，即
$\varphi (n)={\sum }_{i=1}^n[gcd(i,n)==1]$
$\varphi (n)$可用欧拉筛在$O(n)$时间内求出$1-n$的函数值，这个可以去自行了解，在此不多加证明。
$\varphi (n)$具有性质$n={\sum }_{d|n}\varphi (d)$，可用莫比乌斯反演证明，而另一种证明如下：
定义$f(x)$代表$1-n$中与$n$最大公约数为x的数的个数。即$f(x)={\sum }_{i=1}^n[gcd(n,i)==x]$，则显然有$n={\sum }_{d|n}f(d)$，即枚举n的每个约数，$1-n$中每个数与n的最大公约数一定是n的一个约数。
因为若$gcd(i,j)=x$，则$gcd(\frac{i}{x},\frac{j}{x})=1$，所以
$f(x)={\sum }_{i=1}^n[gcd(n,i)==x]={\sum }_{i=1}^{\frac{n}{x}}[gcd(\frac{n}{x},i)==1]=\varphi (\frac{n}{x})$
所以得$n={\sum }_{d|n}f(d)={\sum }_{d|n}\varphi (\frac{n}{d})$，当枚举每个$d|n$时$\frac{n}{d}$等价于$d$，所以$n={\sum }_{d|n}\varphi (d)$得证。

好吧，其实你看不懂上面的证明也是没关系的，你只需要记住结论$n={\sum }_{d|n}\varphi (d)$就行。
大家可能都知道当看到$[gcd(i,j)==1]$的时候我们会用莫比乌斯函数$\mu (n)$来反演，事实上当遇到求$gcd(i,j)$时也有一种常见转换思路即利用上述公式将$n$代换为$gcd(i,j)$，即$gcd(i,j)={\sum }_{d|gcd(i,j)}\varphi (d)$。
我们再来看看原式变为
${\sum }_{i=1}^n{\sum }_{j=1}^m(n+1-i)(m+1-j)({\sum }_{d|gcd(i,j)}\varphi (d)-1)$。
因为$\varphi (1)=1$，所以我们可以把这一项提出来消去后面的$1$，变为
${\sum }_{i=1}^n{\sum }_{j=1}^m(n+1-i)(m+1-j){\sum }_{d|gcd(i,j)}^{d\ne 1}\varphi (d)$。
再来看这个枚举的d，我们是在枚举$i,j$最大公约数的约数，那就是枚举$i,j$的所有公约数，那就可以换一个思路，一个数$d$作为公约数被计算到，只需要将其分别扩大$i$倍和$j$倍，即枚举所有的$d$，再将其扩大$i$倍与$j$倍，将d作为其公约数进行计算，所以对式子进行变形
${\sum }_{d=2}^{\min (n,m)}{\sum }_{i=1}^{\lfloor \frac{n}{d}\rfloor }{\sum }_{j=1}^{\lfloor \frac{m}{d}\rfloor }(n+1-id)(m+1-jd)\varphi (d)$。
将互不干扰的部分拆开，得到
${\sum }_{d=2}^{\min (n,m)}\varphi (d){\sum }_{i=1}^{\lfloor \frac{n}{d}\rfloor }(n+1-id){\sum }_{j=1}^{\lfloor \frac{m}{d}\rfloor }(m+1-jd)$。
现在前一项可以$O(n)$预处理得到，那么就对后两项进行处理，发现后两项显然是等差数列形式，以$n$项为例，对其求和。
首项：$n-d+1$
末项：$n+1-\lfloor \frac{n}{d}\rfloor \times d=n\mathrm{m}\mathrm{o}\mathrm{d}d+1$
项数：$\lfloor \frac{n}{d}\rfloor$
故${\sum }_{i=1}^{\lfloor \frac{n}{d}\rfloor }(n+1-id)=\frac{(n-d+1+n\mathrm{m}\mathrm{o}\mathrm{d}d+1)\times \lfloor \frac{n}{d}\rfloor }{2}$
对m同理，式子最终化为
$\frac{{\sum }_{d|n}^{\min (n,m)}\varphi (d)(n-d+n\mathrm{m}\mathrm{o}\mathrm{d}d+2)(m-d+m\mathrm{m}\mathrm{o}\mathrm{d}d+2)}{4}$
成功优化至$O(n)。$

T7：
贪心，考试的时候没写（去推上面那个题了），考完发现我三年前的代码竟然在没spj的aoj上过了，震惊。
对连续两项$(a_1,b_1)，(a_2,b_2)$分析，前项在前的时间为$a_1+b_2+\max (a_2,b_1)$，后项在前时间为$a_2+b_1+\max (a_1,b_2)$，若前项在前更优，则有$a_1+b_2+\max (a_2,b_1)<a_2+b_2+\max (a_1,b_2)$，移项得$a_1+b_2-\max (a_1,b_2)<a_2+b_1-\max (a_2,b_1)$，发现最大项被减掉了，所以得到$\min (a_1,b_2)<\min (a_2,b_1)$。我们以这个来排序就做完了。。。。。。
但据dl所说，这个式子不满足不等比性的传递性，所以不能直接用这个式子排序。（详见浅谈邻项交换排序的应用以及需要注意的问题）
所以我们对这个式子进一步分析，发现这个式子可以根据$a,b$大小进行分类：
1.$a_1<b_1$且$a_2<b_2$时，根据$a_1<a_2$排序。
2.$a_1=b_1$且$a_2=b_2$时，可随意排序。
3.$a_1>b_1$且$a_2>b_2$时，根据$b_1>b_2$排序。
因为优先找$a$比$b$小的，所以设计判定函数$d_i$：
$d_i=\{\begin{array}{ll}-1& a_i<b_i\\ 0& a_i=b_i\\ 1& a_i>b_i\end{array}$
以$d_i$为第一关键字，相等再按上述$a,b$条件排序，模拟计算出时间即可。
T8：
数学题，据说有用三分的，tqlorz。
设人距离左边墙长度为$d$，影长为$l$。
分类，首先是影子全在地上的情况，即$d\times \frac{H}{H-h}<=D$时，根据相似得$l=d\times \frac{H}{H-h}-d=d\times \frac{h}{H-h}$，随$d$单调上升，所以当影子抵到右墙底时，即$d=D\times \frac{H-h}{H}$时，影长最长$l_{\max }=D\times \frac{h}{H}$。

然后是影子挂墙上的情况，即$d\times \frac{H}{H-h}>D$且$d<=D$时，延长出去，根据相似可得$d^{\prime }=d\times \frac{H}{H-h}$，再根据相似得$\frac{L}{h}=\frac{d^{\prime }-D}{d^{\prime }-d}$，化简得$L=H-\frac{D(H-h)}{d}$，加上地上的影子得$l=D+H-\frac{D(H-h)}{d}-d$，欲求$l_{\max }$，则使后两项最小，显然后两项为对钩函数，取最低点得当$d=\sqrt{D(H-h)}$时，$l_{\max }=D+H-2\sqrt{D(H-h)}$，这个值显然是比第一种情况的值大的，代入计算就做出这道题啦。

然后你去算就会发现3个样例两个过不了。

然而你的d是有范围的。超范围就要分类讨论。
当$\sqrt{D(H-h)}>D\Rightarrow H-h>D$时，也就是人站墙里面了，那这时$d$应该取$D$，此时$l_{\max }=h$。
当$d^{\prime }<D\Rightarrow \sqrt{D(H-h)}>H$时，也就是影子到墙根里去了，那这时应该取情况1的答案，即$l_{\max }=D\times \frac{h}{H}$。
答案完备
$l_{\max }=\{\begin{array}{ll}h& H-h>D\\ D\times \frac{h}{H}& \sqrt{D(H-h)}>H\\ D+H-2\sqrt{D(H-h)}& \text{else}\end{array}$

总结：难度上升了，摩多摩多。