題目鏈接: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]);
}
}