「SDOI2015」约数个数和

「SDOI2015」约数个数和

设$d(x)$为$x$的约数个数，求

$$\sum_{i=1}^n \sum_{j=1}^m d(ij)$$

$n, m <= 50000$
数据组数 $<= 50000$

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

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

设$ij$唯一分解后为

$$\prod_k q_k^{qi_k + qj_k}$$
则 $ij$约数个数为
$$\prod_k (qi_k + qj_k + 1)$$
所以对于一个质数 $p$，它对$d(ij)$的贡献就是
$$qi + qj + 1 = \sum_{x=0}^{qi} \sum_{y=0}^{qj} [gcd(p^x,p^y) = 1]$$
$$\therefore d(ij) = \prod_k \bigg( \sum_{x=0}^{qi_k} \sum_{y=0}^{qj_k}[gcd(p_k^x, p_k^y) = 1] \bigg)$$
令 $x$为$i$， $y$为$j$所有可能的质因数排列方式，可以枚举地的求出 $d(ij)$的值
$$d(ij) = \sum_{x|i} \sum_{y|j} [gcd(x, y) = 1]$$

所以我们得到

$$ans = \sum_{i=1}^n \sum_{j=1}^m \sum_{x|i} \sum_{y|j} [gcd(x, y) = 1]$$
可以看出 $[1, n]$中，一个数 $x$的倍数有$\lfloor\dfrac{n}{x}\rfloor$个
$$\therefore ans = \sum_{x=1}^n \sum_{y=1}^m \lfloor \dfrac{n}{x} \rfloor \lfloor \dfrac{m}{y} \rfloor [gcd(x, y) = 1]$$

$$\because \sum_{d|n} \mu(d) = [ n =1 ]$$

$$\therefore ans = \sum_{x=1}^n \sum_{y=1}^m \sum_{k|gcd(x, y)} \mu(k)\lfloor \dfrac{n}{x} \rfloor \lfloor \dfrac{m}{y} \rfloor$$

$$\therefore ans = \sum_{x=1}^n \sum_{y=1}^m \sum_{k|x 且k|y} \mu(k)\lfloor \dfrac{n}{x} \rfloor \lfloor \dfrac{m}{y} \rfloor$$

我们发现现在在枚举同时为$x, y$的因数$k$，
我们可以转换一下，枚举$k$的倍数$x, y$

$$ans = \sum_{k=1}^{lim} \sum_{x=1}^{\big\lfloor \dfrac{n}{k} \big\rfloor} \sum_{y=1}^{\big\lfloor \dfrac{m}{k} \big\rfloor} \mu(k)\lfloor \dfrac{n}{kx} \rfloor \lfloor \dfrac{m}{ky} \rfloor$$

$lim = min(n, m)$
此时$x, y$的意义是$k$的多少倍，这样原来的$x, y$就对应现在的$kx, ky$
我们设$n = n/k, m = m/k$，交换一下可得

$$ans = \sum_{k=1}^{lim} \mu(k) \sum_{x=1}^n \sum_{y=1}^m \lfloor \dfrac{n}{x} \rfloor \lfloor \dfrac{m}{y} \rfloor$$

我们发现

$$\sum_{x=1}^n \sum_{y=1}^m \lfloor \dfrac{n}{x} \rfloor \lfloor \dfrac{m}{y} \rfloor = \sum_{x=1}^n \lfloor \dfrac{n}{x} \rfloor * \sum_{y=1}^m \lfloor \dfrac{m}{y} \rfloor$$

这样的话，我们明显可以分块$O(n\sqrt{n})$预处理出

$$f(n) = \sum_{i=1}^n \lfloor \dfrac{n}{i} \rfloor$$

那么

$$ans = \sum_{k=1}^{lim} \mu(k) f(\lfloor \dfrac{n}{k} \rfloor)f(\lfloor \dfrac{m}{k} \rfloor)$$

我们知道$\lfloor \dfrac{n}{k} \rfloor$只有$\sqrt{n}$种结果，所以我们可以在$O(\sqrt{n})$的时间内求出单次的$ans$

所以我们在$O(n\sqrt{n})$内做出这道题。

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

code：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 70000;

bool is_prime[N];
int prime[N], tot = 0, mobius[N]; 
int n, m; ll f[N];

void init_mobius( int size )
{
    memset( is_prime, true, sizeof( is_prime ) );
    is_prime[1] = true; mobius[1] = 1;

    for( int i = 2; i <= size; i ++ )
    {
        if( is_prime[i] )
            prime[++tot] = i, mobius[i] = -1;
        for( int j = 1; j <= tot; j ++ )
        {
            int k = i * prime[j];
            if( k > size ) break;

            is_prime[k] = false;
            if( i % prime[j] == 0 )
                { mobius[k] = 0; break; }
            mobius[k] = -mobius[i];  
        }
    }

//  for( int i = 1; i <= 20; i ++ )
//      printf( "%d ", mobius[i] ); printf( "mobius\n" ); 

    for( int i = 2; i <= size; i ++ )
        mobius[i] += mobius[i-1];
} 

inline ll calc_f( int x )
{
    ll ret = 0; int next = 0;
    for( int i = 1; i <= x; i = next+1 )
    {
//      printf( "x:%d, n:%d\n", x, n );
        next = x/(x/i); if( next > x ) next = x;
        ret += (ll)(next - i + 1) * ( x/i );
    }
//  printf( "f(%d) = %d\n", x, ret );
    return ret;
} 

inline void solve()
{
    scanf( "%d%d", &n, &m );

    ll ans = 0; int next = 0, lim = min( n, m );
    for( int i = 1; i <= lim; i = next + 1 )
    {
        next = min( n/(n/i), m/(m/i) ); if( next > lim ) next = lim;
        ans += 1ll * f[n/i] * f[m/i] * (mobius[next] - mobius[i-1]);  

//      printf( "i:%d, %lld\n", i, ans );
    }

    printf( "%lld\n", ans );
}

int main()
{
//  freopen( "output.txt", "w", stdout );

    init_mobius( 60000 );
    for( int i = 1; i <= 60000; i ++ )
        f[i] = calc_f( i );

//  for( int i = 1; i <= 20; i ++ )
//      printf( "%d ", mobius[i] );

    int ttt; scanf( "%d", &ttt );
    while( ttt -- ) 
        solve();

    return 0;
}

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

莫比乌斯函数的性质：

$$\sum_{d|n} \mu(d) = [n = 1]$$