P6810 「MCOI-02」Convex Hull 凸包（数论）

题意：

给定n,m,p，要求计算：
$ans={\sum }_{i=1}^n{\sum }_{j=1}^{m}d(i)d(j)d(gcd(i,j))$
答案对p取模。
其中d(x)=x的约数个数。

数据范围：n,m<=2e6，p<=1e9

解法：

$$\begin{gathered} ans=\sum _{i=1}^n\sum _{j=1}^{m}d(i)d(j)d(gcd(i,j)) \\ 因为d(gcd(i,j))=\sum _{k|i,k|j}1,所以: \\ ans=\sum _{i=1}^n\sum _{j=1}^{m}d(i)d(j)\sum _{k|i,k|j}1 \\ 考虑枚举k: \\ ans=\sum _{k=1}^{min(n,m)}\sum _{k|i,i<=n}\sum _{k|j,j<=m}d(i)d(j) \\ ans=\sum _{k=1}^{min(n,m)}\sum _{k|i,i<=n}d(i)\sum _{k|j,j<=m}d(j) \\ O(n\ast logn)预处理一下: \\ f[k]=\sum _{k|i,i<=n}d(i) \\ g[k]=\sum _{k|j,j<=m}d(j) \end{gathered}$$

code：

#include <bits/stdc++.h>
using namespace std;
const int maxm=2e6+5;
int d[maxm];
int f[maxm];
int g[maxm];
int n,m,mod;
signed main(){
    scanf("%d%d%d",&n,&m,&mod);
    for(int i=1;i<maxm;i++){
        for(int j=i;j<maxm;j+=i){
            d[j]=(d[j]+1)%mod;
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=i;j<=n;j+=i){
            f[i]=(f[i]+d[j])%mod;
        }
    }
    for(int i=1;i<=m;i++){
        for(int j=i;j<=m;j+=i){
            g[i]=(g[i]+d[j])%mod;
        }
    }
    int ans=0;
    for(int k=1;k<=min(n,m);k++){
        ans=(ans+1ll*f[k]*g[k]%mod)%mod;
    }
    printf("%d\n",ans);
    return 0;
}

---