JZOJ4919. 神炎皇

原创于 2016-12-13 15:10:23 发布 · 832 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

题解同时被 3 个专栏收录

161 篇文章

订阅专栏

----------数学----------

20 篇文章

订阅专栏

欧拉函数

2 篇文章

订阅专栏

题目大意

给定 $n$ 求有多少个数对 $(a,b)$ 满足 $a+b\leq n$ 且 $(a+b)$ 是 $ab$ 的因子。

Data Constraint
$n\leq10^{14}$

题解

显然要求 $\frac{ab}{a+b}$ 是一个整数。
先提出一个公因数 $d$ ，得

d a ' b ' a ' + b '

$\frac{da'b'}{a'+b'}$

∵ (a', a' + b') = 1, (b', a' + b') = 1

$\because (a',a'+b')=1,(b',a'+b')=1$

∴ (a' b', a' + b') = 1

$\therefore (a'b',a'+b')=1$
所以要求

d $d$ 是

a′+b′ $a'+b'$ 的倍数。
令

k=a′+b′ $k=a'+b'$
则

d $d$ 可写为

d=k×c $d=k\times c$ 的形式。
又

k×d≤n $k\times d\leq n$

∴c≤nk2 $\therefore c\leq\frac{n}{k^2}$
所以

d $d$ 的取值有

nk2 $\frac{n}{k^2}$ 个

a′b′ $a'b'$ 的取值有

ϕ(k) $\phi(k)$ 个。
枚举

k $k$ 然后统计即可。

时间复杂度： $O(\sqrt{n})$

SRC

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

#define N 10000000 + 10
typedef long long ll ;

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

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

int main() {
    freopen( "uria.in" , "r" , stdin ) ;
    freopen( "uria.out" , "w" , stdout ) ;
    Pre() ;
    scanf( "%lld" , &n ) ;
    ll UP = sqrt(n) ;
    for (ll k = 2 ; k <= UP ; k ++ ) ans += n / (k * k) * Phi[k] ;
    printf( "%lld\n" , ans ) ;
    return 0 ;
}