[bzoj4173]数学

题目大意

懒得写

推导

$[n\ mod\ k+m\ mod\ k\geq k]=[n-\lfloor\frac{n}{k}\rfloor k+m-\lfloor\frac{m}{k}\rfloor k>=k]$
于是等价于
$[\lfloor\frac{n+m}{k}\rfloor-\lfloor\frac{n}{k}\rfloor-\lfloor\frac{m}{k}\rfloor>=1]$
注意这个式子的值也只可能是0或1，因此可以写成
$\lfloor\frac{n+m}{k}\rfloor-\lfloor\frac{n}{k}\rfloor-\lfloor\frac{m}{k}\rfloor$
代入原式
$\sum_{k=1}^{n+m}(\lfloor\frac{n+m}{k}\rfloor-\lfloor\frac{n}{k}\rfloor-\lfloor\frac{m}{k}\rfloor)\phi(k)$
展开来
$\sum_{k=1}^{n+m}\lfloor\frac{n+m}{k}\rfloor \phi(k)-\sum_{k=1}^n\lfloor\frac{n}{k}\rfloor \phi(k)-\sum_{k=1}^m\lfloor\frac{m}{k}\rfloor phi(k)$
注意到$\sum_{i=1}^n \lfloor\frac{n}{i}\rfloor \phi(i)=\sum_{i=1}^n\sum_{d|i}\phi(d)=\sum_{i=1}^ni$
于是你就奥妙重重的发现刚刚那条式子等于nm。
至于前面那个phi(n)和phi(m)随便求一下。

#include<cstdio>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int mo=998244353;
ll i,j;
ll n,m,l,k,t,ans;
int main(){
    scanf("%lld%lld",&n,&m);
    ans=(n%mo)*(m%mo)%mo;
    k=l=n;
    t=floor(sqrt(n));
    fo(i,2,t){
        if (k%i==0){
            l=l/i*(i-1);
            while (k%i==0) k/=i;
        }
    }
    if (k>1) l=l/k*(k-1);
    ans=ans*(l%mo)%mo;
    k=l=m;
    t=floor(sqrt(m));
    fo(i,2,t){
        if (k%i==0){
            l=l/i*(i-1);
            while (k%i==0) k/=i;
        }
    }
    if (k>1) l=l/k*(k-1);
    ans=ans*(l%mo)%mo;
    printf("%lld\n",ans);
}