CodeForces 1K-1600R-1423K

1423K-Lonely Numbers

Description

Link

Analysis

题意简述

定义：

孤独数：在一个给定的数集 $A$ 中的元素 $a$，不存在 $b\in A$，使得 $gcd(a,b),\frac{a}{gcd(a,b)},\frac{b}{gcd(a,b)}$ 能构成三角形的三边长

$t$ 组询问（$1\le t\le 10^6$），每组询问给出一个整数 $n$ （$1\le n\le 10^6$），求 $[1,n]$ 内孤独数的个数

分析

从简化条件的形式入手，若 $a,b,c$ 三边能构成三角形，则 $ka,kb,kc(k\ne 0)$ 三边也能构成三角形

由此，判断 $gcd(a,b),\frac{a}{gcd(a,b)},\frac{b}{gcd(a,b)}$ 能否构成三角形，可以等价于 $gc{d}^2(a,b),a,b$ 能否构成三角形

-

显然，对于互质的数 $a,b$，$gcd(a,b)=1$，则 $1,a,b$ 必然不可能构成三角形

（实数的三岐性：$a,b$ 必然存在序关系 max ⁡ { a , b } ≥ min ⁡ { a , b } + 1 \max\{a,b\}\ge\min\{a,b\}+1 max{a,b}≥min{a,b}+1）

由此，$1$ 一定是孤独数

-

设 $gcd(a,b)=g,a=mg,b=ng$，进一步化简条件，可得 $g$ 必须满足 max ⁡ { m , n } − min ⁡ { m , n } < g < m + n \max\{m,n\}-\min\{m,n\}< g< m+n max{m,n}−min{m,n}<g<m+n

-

对于质数 $x$：

其配对的元素一定是其倍数（否则互质，$g=1$，由上述结论，不可能构成三角形）

讨论 $x$ 的配对情况：

从 $a=1\cdot x$ 出发（ $m=1$），观察其最小的倍数 $b=n\cdot x$（即需满足 $g=x$ ）

若满足 $n>g-m$，则 $n_{min}=g-m+1=x-1+1=x$

则质数 $x$ 不是孤独数，当且仅当 $x^2$ 处于当前范围中（对于 $\forall kx,1\le k<x,k+1\le x$）

-

对于合数 $y$：

若 $y$ 为质数 $\sqrt{y}$ 的平方，根据上述结论，$y$ 与 $\sqrt{y}$ 配对

否则，根据唯一分解定理，$y$ 可以表示为 $p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}..p_k^{{\alpha }_k}$ 的形式，其中 $p_1<p_2<...<p_k$，$a_i\ge 1$

取 $a=y=p_1^{{\alpha }_1}{p}_2^{{\alpha }_2}..p_k^{{\alpha }_k}$，$b=p_1(p_1^{{\alpha }_1-1}{p}_2^{{\alpha }_2}..p_k^{{\alpha }_k}-1)$

则 $g=p_1,m=p_1^{{\alpha }_1-1}{p}_2^{{\alpha }_2}..p_k^{{\alpha }_k},n=p_1^{{\alpha }_1-1}{p}_2^{{\alpha }_2}..p_k^{{\alpha }_k}-1$

$m-n=1<g$，$m+n>g$ ，则 $y$ 必然存在小于 $y$ 的数与其配对

即：合数一定不是孤独数

-

由此，$[1,n]$ 中孤独数的个数，即 $(\sqrt{n},n]$ 中质数的个数 $+1$（包含 $1$）

只需要筛法预处理 $[1,n]$ 中的质数个数，利用前缀和即可 $O(1)$ 求解每个询问

时间复杂度：$O(n)$ （线性筛法）

Code

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

constexpr const int N = 1e6 + 15;

int t, n, p[N], f[N], cnt;
bool st[N];

void primes(int x)
{
    for(int i = 2; i <= x; i++)
    {
        if(!st[i])
        {
            p[++cnt] = i;
            f[i] = f[i - 1] + 1;
        }
        else f[i] = f[i - 1];
        for(int j = 1; p[j] <= x / i; j++)
        {
            st[i * p[j]] = 1;
            if(i % p[j] == 0) break;
        }
    }
}

int main()
{
    cin.tie(0);
    cout.tie(0);
    ios::sync_with_stdio(0);

    primes(1e6 + 10);
    cin >> t;
    while(t--)
    {
        cin >> n;
        cout << f[n] - f[(int)sqrt(n)] + 1 << '\n';
    }

    return 0;
}

Tag

筛法、前缀和、$gcd$、唯一分解定理