bzoj4174，大力推导

传送门
原式=$\sum_{n=1}^N \sum_{m=1}^M \sum_{k=1}^m \frac{n*k-n+x}{m}$
=$\sum_{n=1}^N \sum_{m=1}^M \sum_{k=1}^m (n*(k-1)+x-(n*(k-1)+x)\pmod m )/m$
然后对于和式的前半部分，大力推一波代数公式，就行了
后一半部分，需要点数论知识。
利用完全剩余系的理论，设$d=gcd(n,m)$，我们知道$(n*(k-1)+x)\pmod m )$刚好分成$\frac{m}{d}$段，每段由0,d,2*d,一直到m-d，然后$(n*(k-1)+x)\pmod m$就等于$\frac{m^2-m*d}{2}+x\pmod d*m$
把$m＾２$项利用代数方法取出计算后，剩下的就是莫比乌斯反演模板了，复杂度应该是线性的，比网上那些O(nlogn)的高到不知道哪里去了，在bzoj上是rank13。
然而好久没打莫比乌斯反演，一开始莫比乌斯函数给筛错了，我太菜了。

#include<iostream>
#include<cstdio>
using namespace std;
typedef long long ll;
const ll mo=998244353;
const int N=500010;
ll n,m,x,i,j,ans,p[N],mu[N],xb,u,v,inv2=499122177,f[N],s[N];
bool b[N];
double xx;
inline int min(int a,int b){
    return a>b?b:a;
}
inline ll got(int n,int m){
    ll ans=0;
    for(ll i=1,j;i<=n && i<=m;i=j+1){
        j=min(n/(n/i),m/(m/i));
        ans=(ans+(s[j]-s[i-1])*(n/i)%mo*(m/i)%mo+mo)%mo;
    }
    return ans;
}
int main(){
    cin>>n>>m>>xx;
    mu[1]=1;
    for(i=2;i<=n;++i){
        if(!b[i])mu[p[++xb]=i]=-1;
        for(j=1;j<=xb;++j){
            x=i*p[j];
            if(x>n)break;
            b[x]=1;
            if(i%p[j]==0)break;
            mu[x]=-mu[i];
        }
    }
    for(i=1;i<=n;++i)s[i]=s[i-1]+mu[i];
    x=xx;
    ans=n*m%mo*x%mo;
    for(i=1;i<=m;++i)
        ans=(ans+(i-1)*(n*(n+1)%mo*inv2%mo)%mo*inv2%mo)%mo;
    ans=(ans+mo-m*(m+1)%mo*inv2%mo*inv2%mo*n%mo)%mo;
    for(i=1;i<=n;++i){
        u=n/i;
        v=m/i;
        ans=(ans+(i-(x%i<<1)+mo)*got(u,v)%mo*inv2%mo)%mo;
    }
    printf("%lld\n",ans);
}