CF 757 B Bash's Big Day gcd 质因子分解

题意: 给出一个数组, 问最多能取出多少个数使得其gcd不为1.

思路:

若干个数的gcd!=1等价于这些数有一个公共质因子.

这样, 问题就变为: 一个质数被数组中多少个数用到, 找到被用到最多的那个质数, 输出它被用了多少次.

只要对数组中的每个数分解质数因子, 统计每个质数被用了的次数即可.

质数分解:

    int x = a[i];
    int pos = 0;
    while (x != 1 && !is_prime[x]) {
        if (x % prime[pos] == 0)son[prime[pos]]++;
        while (x % prime[pos] == 0)x /= prime[pos];
        pos++;
    }
    if (x != 1)
        son[x]++;

可以证明, 求1个数的所有质因子的时间复杂度为O(s), s为n以下质数的个数. (先筛质数)

本题代码:

#include<bits/stdc++.h>
using namespace std;
const int M = 1e5 + 5;

bool is_prime[M + 5];
int prime[M + 5], cnt = 0;// 素数个数

void jp() {
    fill(is_prime, is_prime + M, 1);
    is_prime[0] = is_prime[1] = 0;
    for (int i = 2; i <= M; ++i) {
        if (is_prime[i]) {
            prime[cnt++] = i;
        }
        for (int j = 0; j < cnt && prime[j] * i <= M; ++j) {
            is_prime[prime[j] * i] = 0;
            if (i % prime[j] == 0)
                break;
        }
    }
}

int gcd(int a, int b) { return b == 0 ? a : gcd(b, a % b); }

int n;
int a[M];

map<int, int> son;

void init() {
    jp();
    scanf("%d", &n);
    for (int i = 0; i < n; i++) {
        scanf("%d", &a[i]);
    }
    for (int i = 0; i < n; i++) {
        if (is_prime[a[i]])son[a[i]]++;
        else {
            int x = a[i];
            int pos = 0;
            while (x != 1 && !is_prime[x]) {
                if (x % prime[pos] == 0)son[prime[pos]]++;
                while (x % prime[pos] == 0)x /= prime[pos];
                pos++;
            }
            if (x != 1)
                son[x]++;
        }
    }
    int ans = 1;
    for (auto x:son) {
        ans = max(ans, x.second);
    }
    printf("%d\n", ans);
}


int main() {
    init();
    return 0;
}