JZOJ4872. 太阳神

题目大意

给定一个$n$，求有多少对$(a,b)$满足$a\leq n,b\leq n,lcm(a,b)>n$

Data Constraint
$n\leq10^10$

题解

直接计算不太好算。考虑用所有$n^2$的方案减去不合法的方案。
现在要求$\sum_{i=1}^{n}\sum_{j=1}^{n}[\frac{i\times j}{(i,j)}\leq n]$
令$d=(i,j)$
易得$\sum_{d=1}^n\sum_{i=1}^{n}\sum_{j=1}^n[(i,j)=1][i\times j\times d\leq n]$
莫比乌斯反演一波：$\sum_{d=1}^{n}\sum_{k=1}^{n}\mu(k)\sum_{i=1}^n\sum_{j=1}^n[i\times j\times d\leq \frac{n}{k^2}]$
转换主体：$\sum_{k=1}^{\sqrt{n}}\mu(k)\sum_{d=1}^{n}\sum_{i=1}^n\sum_{j=1}^n[i\times j\times d\leq \frac{n}{k^2}]$
前面是根号的，现在要求后面那一部分，也就是求$\sum_{a=1}\sum_{b=1}\sum_{c=1}abc\leq n$
分类讨论即可。假设$a\leq b\leq c$，那么$a\in[1,n^{\frac{1}{3}}]$，$b\in[1,\sqrt{\frac{n}{a}}]$，暴力枚举$a,b$统计即可。通过积分可以算出，这个部分的复杂度是$O(n^{\frac{1}{6}})$

时间复杂度：$O(n^{\frac{2}{3}})$

SRC

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std ;

#define N 100000 + 10
typedef long long ll ;
const int MO = 1e9 + 7 ;

bool flag[N] ;
int Pri[N] , miu[N] ;
ll n , ans ;

void Pre() {
    miu[1] = 1 ;
    for (int i = 2 ; i < N ; i ++ ) {
        if ( !flag[i] ) {
            Pri[++Pri[0]] = i ;
            miu[i] = -1 ;
        }
        for (int j = 1 ; j <= Pri[0] ; j ++ ) {
            if ( i * Pri[j] >= N ) break ;
            flag[i*Pri[j]] = 1 ;
            if ( i % Pri[j] == 0 ) { miu[i*Pri[j]] = 0 ; break ; }
            else miu[i*Pri[j]] = -miu[i] ;
        }
    }
}

ll Solve( ll c ) {
    ll ret = 0 ;
    for (int a = 1 ; (ll)a * a * a <= c ; a ++ ) {
        ret = (ret + 3ll * (c / a / a - a) % MO + 1ll) % MO ;
        int UP = sqrt( c / a ) ;
        for (int b = a + 1 ; b <= UP ; b ++ ) {
            ret = (ret + 6ll * (c / a / b - b) % MO + 3ll) % MO ;
        }
    }
    return ret ;
}

int main() {
    freopen( "ra.in" , "r" , stdin ) ;
    freopen( "ra.out" , "w" , stdout ) ;
    Pre() ;
    scanf( "%lld" , &n ) ;
    ans = 0 ;
    int UP = sqrt(n) ;
    for (int i = 1 ; i <= UP ; i ++ ) {
        ans = (ans + (ll)miu[i] * Solve( n / i / i ) % MO + MO) % MO ;
    }
    n %= MO ;
    ans = (n * n % MO - ans + MO) % MO ;
    printf( "%lld\n" , ans ) ;
    return 0 ;
}

以上.