【每日一题】ABC243G-Sqrt | 数据范围暴力 | 中等

题目内容

原题链接

$T$ 组数据，
每组数据给定一个整数 $x$ ，构造一个下标从 $1$ 开始的，长度为 $10^{100}$ 的数组 $a$ 。
满足 $a_1=x$ ，$a_i\times a_i\le a_{i-1}(2\le i\le 10^{100})$

问可以构造输出多少种不同的数组 $a$ 。

数据范围

$1\le x\le 9\times 10^{18}$

题解

一个关键的点在于，后一个数 $a_i$ 至多是前一个数 $a_{i-1}$ 的平方根，即 $a_i\le \sqrt{a_{i-1}}$

以下用 $\sqrt{x}$ 表示数 $x$ 开平方根后下取整

所以当每个 $a_i=\sqrt{a_{i-1}}$
考虑极端情况，
$a_1=9\times 10^{18}$
$a_2=3\times 10^9$
$a_3=54772$
$a_4=233$
$a_5=15$
$a_6=3$
$a_7=1$

故至多 $6$ 次开根，$a_i$ 就必须为 $1$ 了。

定义 $f_{i,j}$ 表示数组长度为 $i$ ，第 $i$ 个数为 $j$ 的 方案数

可以知道的是，$f_{2,j}=1$ ，其中 $j^2\le x$
那么对于， $f_{3,k}$ 来说，$f_{3,k}=\sum _{j=k^2}^{\sqrt{x}}{f}_{2,j}=\sqrt{x}-k^2+1$

对于 $f_4,f_5,f_6,f_7$ 来说，取值范围太小，可以直接暴力。

代码

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;

const int N = 55000;
ll dp[N];
ll ndp[N];

void solve() {
    ll a;
    cin >> a;

/*
    dp[0][a] = 1;

    b = sqrtl(a)
    dp[1][b] = 1;
    dp[1][b-1] = 1
    dp[1][b-2] = 1
    dp[1][b-3] = 1
    ...
    dp[1][1] = 1

    c = sqrtl(b)
    dp[2][c] = dp[1][c^2]+dp[1][c^2+1]+...+dp[1][b]
    dp[2][c-1] = dp[1][(c-1)^2]+dp[1][(c-1)^2+1]+...dp[1][c]
*/
    ll b = sqrtl(a);
    ll c = sqrtl(b);


    // dp[2]
    dp[c + 1] = 0;
    for (ll i = c; i >= 1; --i) {
        dp[i] = b - i * i + 1;
        dp[i] += dp[i + 1];
    }

    // dp[3...]
    while (c > 1) {
        ll d = sqrtl(c);

        ndp[d + 1] = 0;
        for (ll i = d; i >= 1; --i) {
            ndp[i] = dp[i * i];
            ndp[i] += ndp[i + 1];
        }

        for (int i = 1; i <= d; ++i) dp[i] = ndp[i];

        c = d;
    }

    cout << dp[1] << "\n";
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);

    int T = 1;
    cin >> T;
    while (T--) {
        solve();
    }

    return 0;
}