CF 402D Upgrading Array

暴力搞过。。。处理出从1 到 i的公约数。从右往左开始判断，如果公约数的值小于零直接约掉。然后计算总和即可。注意处理素数大于sqrt(1e9)的情况。当给出素数中存在大于素数表最大值的素数时,加入素数表。v记录是否为坏素数。当除以约数结束时，判断剩下的值是否存在于素数表中，存在则根据v数组的记录判断加或减，否则加1.

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <cmath>
using namespace std;
typedef long long LL;
#define INF 100007
#define N 5005
int n, m;
int a[N];
int b[N];
int p[40000];
int r[100000];
int v[40000];
int gd[N];
int ptr;
int gcd(int a, int b){
    int r = a % b;
    while(r){
        a = b;
        b = r;
        r = a % b;
    }
    return b;   
}
void prime(){
    memset(r, 0, sizeof(r));
    ptr = 0;
    for(int i = 2; i < 32000; i++){
        if(r[i] == 0)p[ptr ++] = i;
        for(int j = 0; j < ptr; j++){
            if(i * p[j] < 32000)r[i * p[j]] = 1;
        }
    }
}
int jud(int x){
    int ul = sqrt(x);
    int sum = 0;
    int i;
    for(i = 0; p[i] <= ul; i++){
        while(x % p[i] == 0){
            x /= p[i];
            if(v[i])sum--;
            else sum++;
        }
        if(x == 1)break;
    }
    if(x > 1){
        while(p[i] != x && i < ptr)i++;
        if(i == ptr){//不存在于素数表中
            sum++;
        }
        else {
            if(v[i])sum--;
            else sum++;
        }
    }
    return sum;
}
int main(){
    prime();
    scanf("%d%d", &n, &m);
    for(int i = 0; i < n; i++){
        scanf("%d", &a[i]);
    }
    for(int i = 0; i < m; i++){
        scanf("%d", &b[i]);
        if(b[i] > p[ptr - 1]){
            p[ptr ++] = b[i];
            v[ptr] = 1;
        }
        else
        for(int j = 0; j < 32000; j++){
            if(b[i] == p[j]){
                v[j] = 1;break;
            }
        }
    }
    gd[0] = a[0];
    for(int i = 1; i < n; i++){
        gd[i] = gcd(a[i], gd[i - 1]);
    }
    int rec = 1;
    for(int i = n - 1; i >= 0; i--){
        gd[i] /= rec;
        if(jud(gd[i]) < 0){
            rec *= gd[i];
        }
        a[i] /= rec;
    }
    int ans = 0;
    for(int i = 0; i < n; i++){
        int ul = sqrt(a[i]);
        int j;
        for(j = 0; p[j] <= ul; j++){
            if(a[i] == 1)break;
            while(a[i] % p[j] == 0){
                a[i] /=  p[j];
                if(v[j])ans --;
                else ans ++;
            }
        }
        if(a[i] > 1){
            while(p[j] != a[i] && j < ptr) j++;
            if(j == ptr){//不存在于素数表中
                ans ++;
            }
            else {
                if(v[j])ans --;
                else  ans ++;
            }

        }
    }
    printf("%d\n", ans);
    return 0;
}