题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4777
题目大意:
给n个数,求LR区间内有多少个数和当前区间所有其他数互质。
解题思路:
训练赛完全想不到,首先我们需要预处理出 对于一个数 i 它最左和最右到哪不互质,记为 pl pr,
那么当我们查询区间 L R时,实际上就是查询区间中有多少 pl<L pr>R的数。
首先基本没有涉及修改操作的题目,几乎都是需要的离线的,这题也不例外,做法也比较多,我们先按照R排序,
当我们遇到i的时候,在pl处+1,在pr处-1,遇到pr的时候,在i处+1,pl处-1。
这其中个人感觉很难理解,最好自己推一下样例,应该会容易理解很多。
Ac代码:
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+5;
bool vis[maxn];
vector<int> pri[maxn],vec[maxn]; //储存每个数的素因子
int n,m,a[maxn],pos[maxn],ans[maxn];
struct node
{
int l,r;
int idx;
bool operator<(const node &p) const
{
return r<p.r;
}
}pa[maxn],q[maxn];
void init(int si) //得到每个数的素因子
{
for(int i=0;i<maxn;i++) pri[i].clear();
for(int i=2;i<maxn;i++)
if(pri[i].size()==0)
for(int j=i;j<maxn;j+=i)
pri[j].push_back(i);
}
void vs() //找出每个数左边界和右边界
{
for(int i=1;i<=n;i++) pa[i].l=0,pa[i].r=n+1;
for(int i=0;i<maxn;i++) pos[i]=0;
for(int i=1;i<=n;i++)
{
int now=a[i];
for(int j=0;j<pri[now].size();j++)
{
int t=pri[now][j];
pa[i].l=max(pa[i].l,pos[t]);
pa[pos[t]].r=min(pa[pos[t]].r,i);
pos[t]=i;
}
}
}
int bit[maxn];
int lowbit(int x) { return x&(-x); }
void add(int x,int val)
{
if(x==0) return ;
while(x<=n)
{
bit[x]+=val;
x+=lowbit(x);
}
}
int sum(int x)
{
int s=0;
while(x)
{
s+=bit[x];
x-=lowbit(x);
}
return s;
}
int main()
{
init(maxn);
while(scanf("%d%d",&n,&m)!=EOF)
{
if(n==0&&m==0) break;
memset(bit,0,sizeof bit);
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].idx=i;
vs();
sort(q+1,q+1+m);
for(int i=0;i<n+1;i++) vec[i].clear();
for(int i=1;i<=n;i++) vec[pa[i].r].push_back(i);
int id=1;
for(int i=1;i<=m;i++) //区间处理
{
while(id<=n&&id<=q[i].r)
{
add(pa[id].l,1);
int sz=vec[id].size();
for(int j=0;j<sz;j++)
{
int v=vec[id][j];
add(pa[v].l,-1);
add(v,1);
}
id++;
}
ans[q[i].idx]=sum(q[i].r)-sum(q[i].l-1);
ans[q[i].idx]=q[i].r-q[i].l+1-ans[q[i].idx];
}
for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
}
}