hdu 4777 Rabbit Kingdom 离线+树状数组+数论

http://acm.hdu.edu.cn/showproblem.php?pid=4777

这题又一次地证明了我的SB，模拟赛的时候就想nsqrt(n)分块离线，结果交上去果断T了，呵呵呵呵

可以把询问按右端点排序，然后维护结果区间。网上有其他人的题解~我也不详细说了。

但是这个离线方法明明见过，但却没想到！这已经不能忍了！太弱了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn = 200001;
int num[maxn],L[maxn],R[maxn];
int gcd(int x,int y){return y==0?x:gcd(y,x%y);}
struct query{
	int l,r,id;
}q[maxn];
bool cmp(const query &a,const query &b){
	return a.r<b.r;
}
vector<int>rec[maxn];
int ans[maxn];
int res[maxn];
int n,m;
void add(int x,int v){
	if(!x)return;
	for(int i=x;i<=n;i+=i&(-i)){
		res[i]+=v;
	}
}
int get(int x){
	int ans=0;
	for(int i=x;i>0;i-=i&(-i)){
		ans+=res[i];
	}
	return ans;
}
bool flag[maxn];
vector<int>fac[maxn];
int prime[20000],top;
void get_prime(){
	for(int i=2;i<maxn;i++){
		if(!flag[i]){
			prime[top++]=i;
			for(int j=i;j<maxn;j+=i){
				fac[j].push_back(top-1);
				flag[j]=true;
			}
		}
	}
}
int pos[20000];
int main(){
	get_prime();
	while(~scanf("%d%d",&n,&m)){
		if(n==0&&m==0)break;
		memset(pos,-1,sizeof(pos));
		for(int i=1;i<=n;i++){
			scanf("%d",&num[i]);
			res[i]=0;
			rec[i].clear();
			L[i]=1;
			R[i]=n;
		}
		for(int i=1;i<=n;i++){
			int x=num[i];
			for(int j=0;j<fac[x].size();j++){
				if(pos[fac[x][j]]!=-1){
					L[i]=max(pos[fac[x][j]]+1,L[i]);
				}
				pos[fac[x][j]]=i;
			}
		}
		memset(pos,-1,sizeof(pos));
		for(int i=n;i>=1;i--){
			int x=num[i];
			for(int j=0;j<fac[x].size();j++){
				if(pos[fac[x][j]]!=-1){
					R[i]=min(pos[fac[x][j]]-1,R[i]);
				}
				pos[fac[x][j]]=i;
			}
		}
		for(int i=1;i<=n;i++){
			rec[R[i]+1].push_back(i);
		}
		for(int i=0;i<m;i++){
			scanf("%d%d",&q[i].l,&q[i].r);
			q[i].id=i;
		}
		sort(q,q+m,cmp);
		int i=0;
		for(int r=1;r<=n;r++){
			add(r,1);
			add(L[r]-1,-1);
			for(int j=0;j<rec[r].size();j++){
				add(rec[r][j],-1);
				add(L[rec[r][j]]-1,1);
			}
			while(i<m&&q[i].r==r){
				ans[q[i].id]=get(r)-get(q[i].l-1);
				i++;
			}
		}
		for(int i=0;i<m;i++){
			printf("%d\n",ans[i]);
		}
	}
}