Aladdin and the Flying Carpet（数论）

描述

问有多少种矩形，边 $x,y$，满足 $x\ast y=a$ 且 $b\le x<y$。

链接

思路

这道题是在刷 $kuangbin$ 基础题题单时遇到的。至于为什么要写这篇题解，是因为网上很多题解都是不对的，所以我想写一篇我认为合适的做法来给别人作为参考。

他们的做法是：预处理 $10^6$ 以内所有素数，用算术基本定理计算 $a$ 所有约数和，通过素数定理可知这个过程复杂度是 $O(\frac{\sqrt{a}}{log\sqrt{a}})$，然后又暴力枚举 $1$ ~ $min(b,\sqrt{a})$ 中 $a$ 的约数，减去这些约数的个数。当 $b$ 和 $\sqrt{(}a)$ 差不多大时，这个复杂度是 $O(\sqrt{a})$，这是肯定是要超时的，这和是直接 $O(\sqrt{a})$ 枚举 $a$ 的约数一个个判断没什么两样，但是好像由于数据比较水，结果能 $AC$。

我的做法是：首先前面预处理素数是一样的，但是是$O(\frac{\sqrt{a}}{log\sqrt{a}})$把 $a$ 的质因子，质因子出现的次数保存下来。然后由倍数法的推论 $1$ ~ $d$ 中平均每个数的约数个数大约只有 $logd。$即一个数 $x$ 的约数个数远远达不到 $\sqrt{x}$ 这个上界，经实际测算 $int$ 以内的数约数个数最多的只有 $1600$ 个，所以我们可以推测 $10^{12}$ 以内的数约数个数最多也是很少的。因此可以 $dfs$ 搜索组成 $a$ 的所有约数。这部分的时间复杂度是不会很高的，因此总体时间复杂度可以看作 $O(t\ast \frac{\sqrt{a}}{log\sqrt{a}})$

代码

#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define pii pair<int,int>
#define pll pair<ll,ll>
#define pli pair<ll,int>
#define Min(a,b,c)  min(a,min(b,c))
#define Max(a,b,c)  max(a,max(b,c))

typedef long long ll;
typedef unsigned long long ull;

const double pi = 3.141592653589793;
const double eps = 1e-8;
const int INF = 0x3f3f3f3f;

const int N = 1000010;

ll a, b;
int prime[N], v[N], m;
ll p[100], c[100], cnt;
ll ans = 0;

void primes(int n)
{
    for (int i = 2; i <= n; i++)
    {
        if (v[i] == 0)
            v[i] = i, prime[++m] = i;
        for (int j = 1; j <= m; j++)
        {
            if (prime[j] > v[i] || prime[j] * i > n) break;
            v[prime[j] * i] = prime[j];
        }
    }
}

void dfs(int x, ll val)
{
    
    if (x > cnt)
    {
        if (val >= b && a / val >= b && val != a / val) {ans++; }
        return;
    }
    for (int i = 0; i <= c[x]; i++)
    {
        dfs(x + 1, val);
        val *= p[x];
    }
}

int main()
{
    primes(1000000);
    int t, tt = 0;
    scanf("%d", &t);
    while (t--)
    {
        scanf("%lld%lld", &a, &b);
        cnt = 0;
        ans = 0;
        ll temp = a;
        for (int i = 1; i <= m && (ll)prime[i] * prime[i] <= a; i++)
        {
            if (a % prime[i] == 0)
            {
                p[++cnt] = prime[i], c[cnt] = 0;
                while (a % prime[i] == 0)
                {
                    a /= prime[i];
                    c[cnt]++;
                }
            }
        }
        if (a > 1)
            p[++cnt] = a, c[cnt] = 1;
        a = temp;
        dfs(1, 1);
        printf("Case %d: %lld\n", ++tt, ans / 2);
    }
    return 0;
}