HDU - 1299 Diophantus of Alexandria（因数分解求方程整数解的个数）

传送门

---

题目大意

给定$n$，求出$\frac{1}{x}+\frac{1}{y}=\frac{1}{n}$的正整数解的个数

解题思路

一开始我以为$x,y$其中有一个一定是$n$的倍数，但是后来才发现如果这样的话答案应该就是$n$的因子个数，但是第二个并不是这样的。又从其他地方下手找找，然鹅也没有找到正确解法

正解是：不难发现$x,y$均大于$n$且一定有一个较大者。设$x\le y$，令$y=n+k$，代入得到$x=\frac{n^2}{k}+n$，这样问题就转化为了求$n^2$的所有因子。又因为$x\le y$，那么可以得出只需要求出$n^2$的大于$n$的因子个数

下面代码我省事写了$set$暴力存，正解应该是先对$n$质因数分解，然后对每个质因数的指数乘上$2$即为$n^2$的质因数分解结果，最后再搜索求其所有的因子数

//
// Created by Happig on 2020/10/18
//
#include <bits/stdc++.h>
#include <unordered_map>
#include <unordered_set>

using namespace std;
#define fi first
#define se second
#define pb push_back
#define ins insert
#define Vector Point
#define ENDL "\n"
#define lowbit(x) (x&(-x))
#define mkp(x, y) make_pair(x,y)
#define mem(a, x) memset(a,x,sizeof a);
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef pair<double, double> pdd;
const double eps = 1e-8;
const double pi = acos(-1.0);
const int inf = 0x3f3f3f3f;
const double dinf = 1e300;
const ll INF = 1e18;
const int Mod = 1e9 + 7;
const int maxn = 2e5 + 10;


set<ll> s;
ll fac[maxn];
int tot;

void getFac(ll n) {
    tot = 0;
    for (ll i = 1; i * i <= n; i++) {
        if (n % i == 0) {
            fac[++tot] = i;
            if (n / i != i) fac[++tot] = n / i;
        }
    }
}

int main() {
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0);
    int T, kase = 0;
    ll x, y, n;
    cin >> T;
    while (T--) {
        cin >> n;
        getFac(n);
        s.clear();
        for (int i = 1; i <= tot; i++) {
            for (int j = 1; j <= tot; j++) {
                s.insert(fac[i] * fac[j]);
            }
        }
        auto pos = s.lower_bound(n);
        ll cnt = 0;
        for (auto i = s.begin(); i != pos; i++) {
            cnt++;
        }
        cout << "Scenario #" << ++kase << ":" << ENDL;
        cout << s.size() - cnt << ENDL;
        cout << ENDL;
    }
    return 0;
}