Codeforces Round #236 (Div. 2)__Upgrading Array

题目链接

• 思路：
  求出所有前缀区间的gcd，然后倒着处理。对于当前位置的前缀gcd，求一下这个数的f函数值，如果是负的，那么肯定去掉是对的。然后对之前的位置，gcd值要除以它之后的区间已经除过的数（不用依次去除，只用记录一下之前除的数是多少，然后计算到一个区间的时候，直接gcd除以这个数即可），然后依次向前处理直到结束即可
• 分析：
  1.为什么要倒着处理呢，因为题目要求每次都是前缀操作，那么对于最后一个区间（最长的区间），没有人会对他产生影响。
  2.方法的正确性：由于是前缀区间的gcd，那么最后一个区间的gcd肯定包括在其他所有区间中，而且题目的操作单位是gcd而不是质因子。那么对于前边的区间，可以分成两个部分，一个是最后一个区间的gcd值，另一个是剩下的部分。
  考虑两种情况：（1）最后区间的gcd为正：不妨设在i位置时候需要除去，那么在最后一个区间不除肯定比较优，因为少除去了一部分正值；（2）如果为负，那么同理会多除去一部分负值，也最优。所以这个策略是没错的

int na, nb;int ipt[MAXN], gcd[MAXN], bp[MAXN];void solve(int& n, int v){    if (binary_search(bp, bp + nb, v))        n--;    else        n++;}int fun(int n){    int ret = 0;    for (int i = 2; i * i <= n; i++)//这里使用 i <= n / i的话会超时...无奈    {        while (n % i == 0)        {            solve(ret, i);            n /= i;        }    }    if (n > 1)        solve(ret, n);    return ret;}int main(){//    freopen("in.txt", "r", stdin);    while (~RII(na, nb))    {        int ans = 0, used = 1, t;        gcd[0] = 0;        FE(i, 1, na)        {            RI(ipt[i]);            gcd[i] = __gcd(gcd[i - 1], ipt[i]);        }        REP(i, nb) RI(bp[i]);        sort(bp, bp + nb);        FE(i, 1, na)            ans += fun(ipt[i]);        FED(i, na, 1)        {            if ((t = fun(gcd[i] / used)) < 0)            {                ans -= t * i;                used = gcd[i];            }        }        cout << ans << endl;    }    return 0;}