13杭州区域赛现场赛Rabbit Kingdom（树状数组+离线）

题意：给你一个长度数列，再给你m个询问（一个区间），问你在这个区间里面有多少个数与其他的数都互质。

解题思路：你看这种类型的题目都可以肯定这是 离线+树状数组（线段树）。主要就是他的更新信息。这里我的处理是先把1-200000（每个数的范围）数里面所有的质因子求出来。然后从后往前遍历数组。会出现以下几种情况

          1.a[k]的质因子在后面出现过而【更新标记】 和【被更新标记】 都为假

          2.a[k]的质因子在后面出现过的那个位置 I   【更新标记】为 真 。

          3.a[k]的质因子在后面出现过且那个位置 I  【被更新标记】  为真 。

    如果只出现 1 这种情况，那么这个位置   +1，【更新标记】 变为真；

    如果只出现 2 这种情况， 那么对这个位置不更新。

    如果出现 3 这种情况， 那么  I  的被更新位置 -1  ，I位置 +1  ，且 I的被更新标记变为假。

    在这个过程中，a[k]的所有质因子的 最新位置 都变为K ，再找到到离 a【K】 最近的和它 不互质的数 a【X】 把K的 【被更新标记】变为  X，再在X位置 +1

    （上面所有的操作都是为了使得记数不 错误）

      最后的答案就是   区间内数的个数 - 区间求的和。

解题代码：

// File Name: c.cpp
// Author: darkdream
// Created Time: 2014年05月07日 星期三 17时07分23秒

#include<vector>
#include<list>
#include<map>
#include<set>
#include<deque>
#include<stack>
#include<bitset>
#include<algorithm>
#include<functional>
#include<numeric>
#include<utility>
#include<sstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<ctime>

using namespace std;
#define MAXN  200005


int a[MAXN];
struct lr{
    int l,r,num;
}q[MAXN];
int cmp(lr a, lr b)
{
    return a.l > b.l;
}
int tree[MAXN];
int lowbit(int x)
{
    return x&(-x);
}
int getsum(int l,int k)
{
    int sum = 0 ; 
    while(k >= l )
    {
        sum += tree[k];
        k -= lowbit(k);
    }
    return sum ; 
}
int n;
int hs[MAXN];
void update(int k  ,int v )
{
    //    printf("%d
",k);
    while(k <= n)
    {
        tree[k] += v; 
        k += lowbit(k);
    }

}
int ans[MAXN];
int prime[MAXN][20];
int num[MAXN];
int last[MAXN];
void solve(){
    memset(num,0,sizeof(num));
    memset(hs,0,sizeof(hs));
    for(int i = 2 ;i <= 200000;i ++)
    {
        if(hs[i] == 0 )
        {
            int k = i;
            while(k <= 200000)
            {
                num[k] ++ ; 
                prime[k][num[k]] = i ; 
                hs[k] = 1; 
                k += i ;
            }
        }
    }
}
int isupdate[MAXN];
int cmp1(int a, int b)
{
    return a < b ;
}
int main(){
    //freopen("intpu.in","r",stdin);
    int m ;
    solve();
    //printf("***
");
    while(scanf("%d %d", &n,&m) != EOF)
    {
        if(n == 0 && m == 0 )
            break;
        memset(tree,0,sizeof(tree));
        memset(q,0,sizeof(q));
        
        
        for(int i = 1;i <= n;i ++)
            scanf("%d",&a[i]);
        for(int i = 1;i <= m;i ++)
        {
            scanf("%d %d",&q[i].l,&q[i].r);
            q[i].num = i; 
        }
        sort(q+1,q+m+1,cmp);
        
        int p = 1;
        memset(hs,0,sizeof(hs));
        memset(last,0,sizeof(last));
        memset(isupdate,0,sizeof(isupdate));
        //memset(q,0,sizeof(q));
        for(int i = n;i >= 1 ;--i)
        {
            int ta[20];
            int tanum = 0 ; 
            memset(ta,0,sizeof(ta));
            for(int j = 1; j <= num[a[i]]; ++j )
            {
                int temp = prime[a[i]][j];
                if(hs[temp])
                {
                    tanum ++ ; 
                    ta[tanum] = hs[temp];
                }
                hs[temp] = i ; 
            }
            if(tanum)
            {
                sort(ta+1,ta+tanum+1,cmp1);
                int min = ta[1];
                for(int j = 1;j <= tanum;j ++)
                {
                    if(ta[j]!= ta[j-1])
                    {
                        if(last[ta[j]])
                        {
                            update(last[ta[j]],-1);
                            update(ta[j],1);
                            last[ta[j]] = 0 ; 
                           }else {
                              if(!isupdate[ta[j]])
                              {
                                 update(ta[j],1);
                                 isupdate[ta[j]] = 1 ; 
                              }
                        }
                       if(ta[j] < min)
                          min = ta[j];
                    }
                }
                update(min,1);
                isupdate[min] = 1; 
                last[i] = min;
                isupdate[i] = 1;

            }
            while(q[p].l == i )
            {
                ans[q[p].num] = (q[p].r - q[p].l +1) - getsum(q[p].l,q[p].r);
                p++;
            }
        }
        for(int i = 1;i <= m;++i)
            printf("%d
",ans[i]);
    }
    return 0;
}