CF402D 【Upgrading Array】

题目链接：

CF402D

题目分析：

首先考虑一下怎么求每个数的分数。把每个数分解到最后会发现它的坏质因子对它分数的贡献是\(-1\)，好质因子对它分数的贡献是\(1\)，那么最后的分数就是好质因数-坏质因数

然后想一想怎么操作。我们的最优答案是把所有能除掉的负数分数的\(gcd\)全部除掉，一个很显然的贪心是从后往前操作，因为前面操作了之后后面的\(gcd\)就恒为\(1\)，操作不下去了

另一个\(point\)是操作是显然正确的，因为后面的操作会除掉前面操作原本的\(gcd\)的一部分，而剩下的一部分应该被除掉的会在之后的贪心中被除去

举个例子
假设区间\([1, 5]\)需要被除掉的\(gcd\)是\(12\)，区间\([1, 9]\)需要被除掉的\(gcd\)是\(4\)，那么在\(9\)号位上除去\(4\)后再在\(5\)号位上除去\(3\)即可

处理一下前缀\(gcd\)，从后向前贪心，并记录当前已经除去了多少，每次进行操作之前要先除去这个量。（另：素数只需筛到\(sqrt(1e9)\)即可，后面的可以暴力判断），对于每一位，质因数分解一下当前位置上的前缀\(gcd\)并计算它的分数，如果\(<0\)就除掉它

用\(bitset\)记录坏质数

代码：

#include<bits/stdc++.h>
#define N (100000 + 5)
using namespace std;
int n, m, prime[N], tot, ans, g[N];
bool vis[N];
bitset <1000000005> S;inline int read() {
    int cnt = 0, f = 1; char c = getchar();
    while (!isdigit(c)) {if (c == '-') f = -f; c = getchar();}
    while (isdigit(c)) {cnt = (cnt << 3) + (cnt << 1) + (c ^ 48); c = getchar();}
    return cnt * f;
}
int a[N], b[N];
int gcd(int a, int b) {return b ? gcd(b, a % b) : a;}
void init() {
    for (register int i = 2; i <= N - 5; ++i) {
        if (!vis[i]) vis[i] = 1, prime[++tot] = i;
        for (register int j = 1; j <= tot && prime[j] * i <= N - 5; ++j) {
            vis[i * prime[j]] = 1;
            if (i % prime[j] == 0) break;
        }
    }
}

int divide(int x) {
    int ans = 0;
    for (register int i = 1; i <= tot && prime[i] * prime[i] <= x; ++i) {
        while (x % prime[i] == 0) ans += S[prime[i]] ? -1 : 1, x /= prime[i];
    }
    if (x > 1) ans += S[x] ? -1 : 1;
    return ans;
}
int main() {
    n = read(), m = read();
    for (register int i = 1; i <= n; ++i) a[i] = read();
    for (register int i = 1; i <= m; ++i) S[read()] = 1;
    init();
    for (register int i = 1; i <= n; ++i) ans += divide(a[i]);
    for (register int i = 1; i <= n; ++i) g[i] = gcd(g[i - 1], a[i]);
    int div = 1;
    for (register int i = n; i >= 1; --i) {
        g[i] /= div;
        int x = divide(g[i]);
        if (x < 0) ans += i * (-x), div *= g[i];
    }
    printf("%d\n", ans);
    return 0;
}