Bzoj 4173 数学

原题链接：
http://www.lydsy.com/JudgeOnline/problem.php?id=4173
貌似题目是付费才可以看的。这里阐述一下题意吧：

定义整数集合$S(n,m)$

它表示所有满足$n\ mod\ k+m\ mod\ k\geq k$的所有的k组成的集合。

要求你计算 $$\varphi(n)\varphi(m)\sum_{i\in S(n,m)}\varphi(i)\ \ mod \ \ 998244353$$

<script type="math/tex; mode=display" id="MathJax-Element-56"></script>

题意有： $$\frac{answer}{\varphi(n)\varphi(m)}=\sum_{i=1}^{n+m}\varphi(i)[n\ mod\ i+m\ mod\ i\geq i]$$

$n\ mod\ i+m\ mod\ i\geq i$ 的对立条件： $$n\ mod\ i+m\ mod\ i< i$$

也就是说: $$n\ mod\ i+m\ mod \ i=(n+m)mod\ i$$

则： $$n-\Big\lfloor\frac{n}{i}\Big\rfloor i+m-\Big\lfloor\frac{m}{i}\Big\rfloor i=(n+m)-\Big\lfloor\frac{n+m}{i}\Big\rfloor i$$

那么也就是说：

$$\Big\lfloor\frac{n}{i}\Big\rfloor +\Big\lfloor\frac{m}{i}\Big\rfloor=\Big\lfloor\frac{n+m}{i}\Big\rfloor$$

对于： $$\Big\lfloor\frac{n+m}{i}\Big\rfloor=\Big\lfloor\frac{n}{i}\Big\rfloor +\Big\lfloor\frac{m}{i}\Big\rfloor+\Big\lfloor\Big\{\frac{n}{i}\Big\}+\Big\{\frac{m}{i}\Big\}\Big\rfloor$$

$$0\leq\Big\{\frac{n}{i}\Big\}+\Big\{\frac{m}{i}\Big\}<2$$

所以：

$$\Big\lfloor\frac{n+m}{i}\Big\rfloor-\Big\lfloor\frac{n}{i}\Big\rfloor -\Big\lfloor\frac{m}{i}\Big\rfloor=0\ \ or\ \ 1$$

则$n\ mod\ i+m\ mod\ i\geq i$ 等价于：

$$\Big\lfloor\frac{n+m}{i}\Big\rfloor-\Big\lfloor\frac{n}{i}\Big\rfloor -\Big\lfloor\frac{m}{i}\Big\rfloor=1$$

所以： $$\sum_{i=1}^{n+m}\varphi(i)[n\ mod\ i+m\ mod\ i\geq i]\\=\sum_{i=1}^{n+m}\varphi(i)\Big(\Big\lfloor\frac{n+m}{i}\Big\rfloor-\Big\lfloor\frac{n}{i}\Big\rfloor -\Big\lfloor\frac{m}{i}\Big\rfloor\Big)\\=\sum_{i=1}^{n+m}\varphi(i)\Big\lfloor\frac{n+m}{i}\Big\rfloor-\sum_{i=1}^{n}\varphi(i)\Big\lfloor\frac{n}{i}\Big\rfloor -\sum_{i=1}^{m}\varphi(i)\Big\lfloor\frac{m}{i}\Big\rfloor$$

对于: $$\sum_{i=1}^{n}\varphi(i)\Big\lfloor\frac{n}{i}\Big\rfloor=\sum_{i=1}^{n}\varphi(i)\sum_{i|d,d\leq n}1\\=\sum_{d=1}^n \sum_{i|d}\varphi(i)=\sum_{d=1}^{n}d=\frac{n(n+1)}{2}$$

所以： $$answer=\varphi(n)\varphi(m)\Big(\frac{(n+m)(n+m+1)-n(n+1)-m(m+1)}{2}\Big)$$

代码：

#include <algorithm>
#include <vector>
#include <stdio.h>
#include <cmath>

using namespace std;
typedef long long LL;

const LL P=998244353;

LL phi(LL n)
{
    LL tmp=n;
    for(LL k=2;k*k<=n;k++)
    {
        if(n%k)continue;
        tmp=tmp-tmp/k;
        while(n%k==0)n/=k;
    }
    if(n>1)tmp=tmp-tmp/n;
    return tmp;
}
LL Pow(LL a,LL b)
{
    LL tmp=1;
    while(b)
    {
        if(b&1)
            tmp=tmp*a%P;
        a=a*a%P;
        b>>=1;
    }
    return tmp;
}
int main ()
{
    LL n,m;
    scanf("%lld %lld",&n,&m);
    LL ans=(phi(n)%P)*(phi(m)%P)%P;
    n%=P;
    m%=P;
    LL Iv2=Pow(2,P-2);
    ans=ans*(((n+m)*(n+m+1)%P*Iv2%P-n*(n+1)%P*Iv2%P-m*(m+1)%P*Iv2%P+P*3)%P)%P;
    printf("%lld\n",ans);
    return 0;
}