NKOJ 9106

最新推荐文章于 2024-09-26 00:03:11 发布

zl12345Q

最新推荐文章于 2024-09-26 00:03:11 发布

阅读量51

点赞数

分类专栏：数论 dfs 文章标签：深度优先算法 Powered by 金山文档

本文链接：https://blog.csdn.net/zl12345Q/article/details/129220494

版权

数论同时被 2 个专栏收录

3 篇文章 0 订阅

订阅专栏

dfs

1 篇文章 0 订阅

订阅专栏

文章介绍了如何使用深度优先搜索（DFS）来枚举质因子个数，以解决最大约数个数的问题。通过优化策略，如从小到大考虑质因子，避免乘法溢出等，来提高算法效率。文中提供了两种版本的DFS实现，一种是基础版，另一种是考虑了上界的优化版。尽管优化版考虑了log常数，但在某些情况下仍可能较慢导致超时。

摘要由CSDN通过智能技术生成

A - 反素数加强版 | CQNK (nks.edu.cn)

dfs枚举质因子个数

根据d(n)(约数个数) = (r1+1)(r2+1)(r3+1)...(rn+1)，求最大约数个数就行了

性质及优化：

一定是先考虑小质因子，再考虑大的

p1 < p2 < p3 < p4...<pn && r1 >= r2 >= r3 >= r4... >= rn, 把上一层的指数存下来, 这一层从1到上层指数遍历

判断是否大于n时，想想乘会不会溢出，把x * y >n改成x > n / y，或者提前算好最多可以乘多少，当前质因子最大指数为logpi(n/now)

#define int long long
int n;
int ans;
int anscnt;
int prime[20] = {0, 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 47, 53, 59};
inline void dfs(int num, int now, int nxt, int cnt) {
    if(prime[now] == 0) return ;
    if(cnt > anscnt) {
        ans = num;
        anscnt = cnt;
    }
    if(cnt == anscnt && num < ans) {
        ans = num;
    }
    for(re int i = 1; i <= nxt; ++ i) {
        if(num > n / prime[now]) return ;
        num *= prime[now];
        dfs(num, now + 1, i, cnt * (i + 1));
    }
}

判上界版

#define int long long
int prime[20] = {0, 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 47, 53, 59};
inline void dfs(int num, int now, int nxt, int cnt) {
    if(prime[now] == 0) return ;
    if(cnt > anscnt) {
        ans = num;
        anscnt = cnt;
    }
    if(cnt == anscnt && num < ans) {
        ans = num;
    }
    int op = min(nxt, (int)(log(n / now) / log(prime[now])));
    for(re int i = 1; i <= op; ++ i) {
        if(num > n / prime[now]) return ;
        num *= prime[now];
        dfs(num, now + 1, i, cnt * (i + 1));
    }
}

log常数较大，整体其实还跑的慢些（tle90）