分析:找到每一个点的左边离他最近的不互质数,记录下标(L数组),右边一样如此(R数组),预处理

这个过程需要分解质因数O(n*sqrt(n))

然后离线,按照区间右端点排序

然后扫一遍,对于当前拍好顺序的第i个询问,将所有小于r的点加入更新

更新的过程是这样的

(1)对于刚加入点x,树状数组L[x]位置+1   把这个定义为左更新

(2)对于所有R[i]=x的点,树状数组L[i]位置-1,i位置+1   把这个定义为右更新

(3)查询是询问区间 l->r的和

时间复杂度分析,因为每个点左更新,右更新各一次,所以单组测试用例是O(nlogn)的

下面来解释为啥这样更新

查看当前询问 l , r,对于所有小于 l 的点 i,它的所有更新(左更新和右更新)不会影响到树状数组区间 l 到 r 的 和

对于在l,r区间的点 i,如果 R[i]<=r,那么对于区间 l ->r 有 1 的贡献

如果 R[i]>r,对于这个点,只有左更新L[i]+1;

如果 L[i]>=l ,对于查询区间有 1 的贡献

如果 L[i]<l,对于查询区间没有贡献

不难发现,这样对于查询区间有贡献的点,都是会和别人打架的点

所以该查询的答案就是  查询区间长度 - 树状数组区间和

然后以下看代码

#include<cstdio>
#include<cstring>
#include<queue>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;
typedef long long LL;
const int N=2e5;
const LL mod=1e9+;
int a[N+],c[N+];
int f[N+],h[N+],k[N+];
vector<int>g[N+];
vector<int>b[N+];
struct Que
{
int l,r,id;
bool operator<(const Que &e)const
{
return r<e.r;
}
}q[N+];
int n,m;
bool vis[N+];
int prime[N+],cnt,res[N+];
void add(int x,int t)
{
for(;x<=n+;x+=(x&(-x)))
c[x]+=t;
}
int get(int x)
{
int ans=;
for(;x>;x-=(x&(-x)))
ans+=c[x];
return ans;
}
int main()
{
for(int i=; i*i<=N; ++i)
{
if(vis[i])continue;
for(int j=i*i; j<=N; j+=i)
vis[j]=;
}
for(int i=; i<=N; ++i)
if(!vis[i])prime[cnt++]=i;
for(int i=; i<=N; ++i)
{
int t=i;
for(int j=; j<cnt&&prime[j]*prime[j]<=i; ++j)
{
if(t%prime[j])continue;
g[i].push_back(prime[j]);
while(t%prime[j]==)t/=prime[j];
if(t==)break;
}
if(!vis[t]&&t!=)
g[i].push_back(t);
}
while(~scanf("%d%d",&n,&m),n)
{
memset(f,,sizeof(f));
memset(c,,sizeof(c));
for(int i=; i<=n+; ++i)
{
scanf("%d",&a[i]);
h[i]=;
k[i]=n+;
}
for(int i=; i<=n+; ++i)
{
for(int j=; j<g[a[i]].size(); ++j)
{
int x=g[a[i]][j];
h[i]=max(h[i],f[x]);
f[x]=i;
}
}
for(int i=; i<=N; ++i)
f[i]=n+;
for(int i=n+; i>; --i)
{
for(int j=; j<g[a[i]].size(); ++j)
{
int x=g[a[i]][j];
k[i]=min(k[i],f[x]);
f[x]=i;
}
}
for(int i=;i<=m;++i)
scanf("%d%d",&q[i].l,&q[i].r),q[i].id=i;
sort(q+,q++m);
for(int i=;i<=n;++i)
b[i].clear();
for(int i=;i<=n+;++i)b[k[i]].push_back(i);
int x=;
for(int i=;i<=m;++i)
{
while(x<=n+&&x<=q[i].r+)
{
add(h[x],);
for(int j=;j<b[x].size();++j)
{
int y=b[x][j];
add(h[y],-);
add(y,);
}
++x;
}
int t=get(q[i].r+)-get(q[i].l);
res[q[i].id]=q[i].r-q[i].l+-t;
}
for(int i=;i<=m;++i)
printf("%d\n",res[i]); }
return ;
}

HDU 4777 Rabbit Kingdom 树状数组的更多相关文章

  1. HDU 4777 Rabbit Kingdom(树状数组)

    HDU 4777 Rabbit Kingdom 题目链接 题意:给定一些序列.每次询问一个区间,求出这个区间和其它数字都互质的数的个数 #include <cstdio> #include ...

  2. HDU 4777 Rabbit Kingdom --容斥原理+树状数组

    题意: 给一个数的序列,询问一些区间,问区间内与区间其他所有的数都互质的数有多少个. 解法: 直接搞有点难, 所谓正难则反,我们求区间内与其他随便某个数不互质的数有多少个,然后区间长度减去它就是答案了 ...

  3. HDU 3333 | Codeforces 703D 树状数组、离散化

    HDU 3333:http://acm.hdu.edu.cn/showproblem.php?pid=3333 这两个题是类似的,都是离线处理查询,对每次查询的区间的右端点进行排序.这里我们需要离散化 ...

  4. HDU 3333 - Turing Tree (树状数组+离线处理+哈希+贪心)

    题意:给一个数组,每次查询输出区间内不重复数字的和. 这是3xian教主的题. 用前缀和的思想可以轻易求得区间的和,但是对于重复数字这点很难处理.在线很难下手,考虑离线处理. 将所有查询区间从右端点由 ...

  5. HDU 3333 Turing Tree (树状数组)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3333 题意就是询问区间不同数字的和. 比较经典的树状数组应用. //#pragma comment(l ...

  6. HDU 4325 Flowers(树状数组+离散化)

    http://acm.hdu.edu.cn/showproblem.php?pid=4325 题意:给出n个区间和m个询问,每个询问为一个x,问有多少个区间包含了x. 思路: 因为数据量比较多,所以需 ...

  7. hdu 5775 Bubble Sort 树状数组

    Bubble Sort 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=5775 Description P is a permutation of t ...

  8. HDU - 1541 Stars 【树状数组】

    题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=1541 题意 求每个等级的星星有多少个 当前这个星星的左下角 有多少个 星星 它的等级就是多少 和它同一 ...

  9. HDU 3854 Glorious Array(树状数组)

    题意:给一些结点,每个结点是黑色或白色,并有一个权值.定义两个结点之间的距离为两个结点之间结点的最小权值当两个结点异色时,否则距离为无穷大.给出两种操作,一种是将某个结点改变颜色,另一个操作是询问当前 ...

随机推荐

  1. 分类-Category

    1.基本用途 如何在不改变原来类模型的前提下,给类扩充一些方法?有2种方式 继承 分类(Category) 2.格式 分类的声明 @interface 类名 (分类名称) // 方法声明 @end 分 ...

  2. 【BZOJ 1070】[SCOI2007]修车

    Description 同一时刻有N位车主带着他们的爱车来到了汽车维修中心.维修中心共有M位技术人员,不同的技术人员对不同的车进行维修所用的时间是不同的.现在需要安排这M位技术人员所维修的车及顺序,使 ...

  3. spring dataSourceRouter自动切换数据源

    spring多数据源的切换,主要用到的是AbstractRoutingDataSource这个路由类,当我们的自定义的一个路由分发类继承AbstractRoutingDataSource类后,重写de ...

  4. JavaScript与DOM的关系

    JavaScript与浏览器的工作 1.浏览器获取并加载你的页面,从上至下解析它的内容. 遇到JavaScript时,浏览器会解析代码,检查它的正确性,然后执行代码. 浏览器还会建立一个HTML页面的 ...

  5. 相似元素存在的意义---HTML&CSS

    1.<q> 效果: 告诉浏览器这是一段短引用,让浏览器以合适的方法来显示 注: 不能直接以双引号直接代替<q>,因为有些浏览器<q>的效果不是双引号. 不要忘了移动 ...

  6. c++ string char* const char*

    #include <iostream> #include <string> #include <cstring> using namespace std; int ...

  7. 谈谈怎么实现Oracle数据库分区表

    谈谈怎么实现Oracle数据库分区表 数据库的读写分离 SQLSERVER性能监控级别步骤 Oracle索引问题诊断与优化(1)

  8. linux下c++實現簡單的生產者消費者隊列模式

    引言 生產者消費者是一個經典的模式 利用生產者,消費者和緩衝區降低了生產者和消費者之間的的耦合度 便於對生產者和消費者的修改 下面記錄的是一個經典的單一生產者多消費者的模式 設計思路 以隊列做為緩衝區 ...

  9. MetadataType的使用

    MetadataType的使用,MVC的Model层数据验证指定要与数据模型类关联的元数据类 using System.ComponentModel.DataAnnotations; //指定要与数据 ...

  10. Java中的try、catch、finally块简单的解析

    package com.wangzhu; import java.util.HashMap; import java.util.Map; /** * 在try.catch.finally块中,若try ...