最快速题解 线性筛+素数性质| #kotori和素因子#

kotori和素因子

https://www.nowcoder.com/practice/7b1c858a3e7a41ed8364178979eaae67

/*
首先用线性筛求出1000以内的质数(这里可以打表做到更快,但实在不必要)
然后求出每个数的所有因子,这一步有很多可以优化的地方。1.先看本身是不是素数,2.将已知素数从小向大比较,3.如果是因子,在计入的同时,还要让这个数本身除掉这个因子(如果是多重还可以除多次)。
然后是深搜过程,这里用到一个特性:素数作为因子是唯一的,当我们要用到一个素数时,在总乘积中乘以这个素数,再想用这个素数就会发现乘积中已经存在,则不可使用。换句话说,就是用乘法,省去了 !use.count(prim) 的比较过程。
在深搜中还有一步剪枝,如果现在的和已经大于了最小值,就把这个枝减掉,免去无意义的计算。
这个方法应该把每一步都优化到极限了,如果还有更优解欢迎回复。
*/

#include <vector>
#include <iostream>

#define N 1000

using namespace std;

vector<int> prime;
bool mark[N+5];
vector<int> factor[12];
int n, MinSum=2147483647;

void getPrime() { //线性筛法求素数
    for(int i=2; i<=N; i++) {
        if(!mark[i]) prime.push_back(i);
        for(int num : prime) {
            if(num * i > N) break;
            mark[num*i] = 1;
            if(i % num == 0) break;
        }
    }
}

void dfs(int i, long long mul, int sum) {
    if(sum > MinSum) return; //剪枝
    if(i>=n) {
        if(sum < MinSum) MinSum = sum;
        return;
    }

    for(int num: factor[i]) {
        if(mul%num != 0) dfs(i+1, mul*num, sum+num); //余数为0说明没被使用过,然后就在乘积中乘以
    }
}

int main() {
    getPrime();
    int num;
    cin >> n;
    for(int i=0; i<n ;i++) {
        scanf("%d", &num);
        if(!mark[num]) { //如果本身是素数则直接结束
            factor[i].push_back(num); 
            continue;
        }
        for(int p: prime) {
            if(p > num) break;
            if(num % p == 0) {
                factor[i].push_back(p);
                while(num % p == 0) num /= p; //快速得到因子
            }
        }
    }
    dfs(0, 1, 0);
    if(MinSum < 2147483647) 
        cout << MinSum << endl;
    else
        puts("-1");
    return 0;
}

首先用线性筛求出1000以内的质数(这里可以打表做到更快,但实在不必要)

然后求出每个数的所有因子,这一步有很多可以优化的地方。1.先看本身是不是素数,2.将已知素数从小向大比较,3.如果是因子,在计入的同时,还要让这个数本身除掉这个因子(如果是多重还可以除多次)。

然后是深搜过程,这里用到一个特性:素数作为因子是唯一的,当我们要用到一个素数时,在总乘积中乘以这个素数,再想用这个素数就会发现乘积中已经存在,则不可使用。换句话说,就是用乘法,省去了 !use.count(prim) 的比较过程。

在深搜中还有一步剪枝,如果现在的和已经大于了最小值,就把这个枝减掉,免去无意义的计算。

这个方法应该把每一步都优化到极限了,如果还有更优解欢迎回复。

全部评论

相关推荐

10-05 23:02
东北大学 Java
我说句实话啊:那时候看三个月培训班视频,随便做个项目背点八股,都能说3 40w是侮辱价
点赞 评论 收藏
分享
10-15 09:13
已编辑
天津大学 soc前端设计
点赞 评论 收藏
分享
评论
点赞
收藏
分享
牛客网
牛客企业服务