图论 + 数论 ---- CF1325E E. Ehab‘s REAL Number Theory Problem (约数个数 + 枚举 + bfs找最小环)[从图结构优化搜索]

题目链接

---

题目大意：

给你$n$个数，每个数保证约数个数不超过$7$个。问你从中选出若
干个数的乘积是完全平方数，最少要选择多少个数？，不存在输出-1.

---

题目思路：

1.每个数约数不超过7个$\to$每个数最多只有两个不同的质因子
$$n=p_1^{a1}{p}_2^{a2}...p_n^{an}$$
$$\text{Number of divisors}=(a_1+1)\ast (a_2+1)..\ast (an+1)\le 7$$
那么可以解除$n\le 2$

2. 我知道平方数里面的质因子的指数都是偶数，那么我们先把偶数的质因子给去掉那么每个质因子要么没有，要么只剩下一个那么$n$可以表示为$n=p_1\times p_2$我们对$p_1和p_2$直接连接无向边，那么要使得是是完全平方数，就是p1,p2,p3要出现偶数个$\to$每个点的度数是偶数
3. 我们要求最小那么每个点的度数最好为2,那么本质上就是找最小环，每个点的度数都是$2$
4. 无向图最小环$\to FloyedO(n^3)$直接T飞了！！
5. 我们观察一下图结构$p_1\le \sqrt{n},p_2>\sqrt{n}$就是小点和大点交错连边
6. 那么我们可以枚举起点，用$bfs$去找最小环 ，$bfs$找最小环是保证起点一定在环上面，但是不能枚举所有的点，我们发现$小点连大点，就是一条边的两端肯定是一个\ge \sqrt{n},一个小于\sqrt{n}$
7. 那么我们只要枚举小于$\sqrt{n}$就可以了$也就最多也就1000$
   那么就可以A了

---

AC code

#include <bits/stdc++.h>
#define mid ((l + r) >> 1)
#define Lson rt << 1, l , mid
#define Rson rt << 1|1, mid + 1, r
#define ms(a,al) memset(a,al,sizeof(a))
#define log2(a) log(a)/log(2)
#define lowbit(x) ((-x) & x)
#define IOS std::ios::sync_with_stdio(0); cin.tie(0); cout.tie(0)
#define INF 0x3f3f3f3f
#define LLF 0x3f3f3f3f3f3f3f3f
#define f first
#define s second
#define endl '\n'
using namespace std;
const int N = 2e6 + 10, mod = 1e9 + 9;
const int maxn = N;
const long double eps = 1e-5;
const int EPS = 500 * 500;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef pair<ll,ll> PLL;
typedef pair<double,double> PDD;
template<typename T> void read(T &x) {
   x = 0;char ch = getchar();ll f = 1;
   while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
   while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
template<typename T, typename... Args> void read(T &first, Args& ... args) {
   read(first);
   read(args...);
}
inline void out(int x) {
    if (x > 9) out(x / 10);
    putchar(x % 10 + '0');
}
vector<int> edge[maxn];
int dist[maxn], fa[maxn];
int ans = INF;
void bfs(int u) {
    memset(dist, INF, sizeof(dist));
    dist[u] = 0;
    queue<int> q;
    q.push(u);
    while (!q.empty()) {
        int u = q.front();
        q.pop();
        for (int i = 0; i < edge[u].size(); i++) {
            int v = edge[u][i];
            if (v != fa[u]) {
                if (dist[v] != INF) ans = min(ans, (dist[u] + dist[v] + 1));
                else fa[v] = u, dist[v] = dist[u] + 1, q.push(v);;
            }
        }
    }
}
int main() {
    int n, x;
    read(n);
    for (int i = 0; i < n; i++) {
        read(x);
        for (int i = 2; i <= sqrt(x); i++) {
            while (x % (i * i) == 0) x /= i * i;
        }
        if (x == 1) {
            out(1);
            return 0;
        }
        int u = x, v = 1;
        edge[u].push_back(v);
        edge[v].push_back(u);
        for (int i = 2; i <= sqrt(x); i++) {
            if (x % i == 0) {
                u = x / i, v = i;
                edge[u].push_back(v);
                edge[v].push_back(u);
            }
        }
    }
    for (int i = 1; i <= sqrt(maxn); i++) bfs(i);
    if (ans != INF) out(ans);
    else printf("-1");
    return 0;
}