hdu4777 Rabbit Kingdom 离线树状数组 求询问区间内的区间数

2024-08-26 13:58

本文主要是介绍hdu4777 Rabbit Kingdom 离线树状数组 求询问区间内的区间数,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题意:询问区间内有多少个数与区间中其他的数都互质

分析:易得,一个区间内的数的个数减去,与其他数不互质的数即可——即离当前数i左边最近的不互质的数的位置(设为L[i])和右边最近的不互质的数的位置(设为R[i])有一个在区间[L,R]内。那么问题就变成统计:1.区间[L,R]中有多少个数的L[i]或R[i]在区间[L,R]内。2.多少个数的L[i]且R[i]在区间[L,R]内。对于每个询问,答案就是区间内的数减去1的结果,再加上2的结果。

2的结果其实就是询问有多少个区间[L[i],R[i]]完全在给定区间[L,R]内。其实1也可以转化为相同的问题,即区间[L[i],i]或[R[i],i],是否在给定区间内。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<ostream>
#include<istream>
#include<algorithm>
#include<queue>
#include<string>
#include<cmath>
#include<set>
#include<map>
#include<stack>
#include<vector>
#define fi first
#define se second
#define ll long long
#define pii pair<int,int>
#define inf (1<<30)
#define eps 1e-8
#define pb push_back
using namespace std;
const int maxn=200010;
int n,m;
int d[maxn];
int L[maxn];
int R[maxn];
int Div[maxn];
int prime[maxn];
int num;
bool vis[maxn];
void getPrime()
{int N=200000;int m=sqrt(N+0.5);memset(vis,0,sizeof(vis));for(int i=2;i<=m;i++) {if(vis[i]==0) {for(int j=i*i;j<=N;j+=i)vis[j]=1;}}num=0;for(int i=2;i<=N;i++) {if(!vis[i])prime[num++]=i;}
}
struct Node
{int l,r;int id;bool operator<(const Node& u)const{return r<u.r;}
}q[maxn],s[3][maxn],ans[3][maxn];
int sum[maxn];
int lowbit(int x)
{return x&-x;
}
int Sum(int x)
{int tot=0;while(x) {tot+=sum[x];x-=lowbit(x);}return tot;
}
void update(int x)
{while(x<=n) {sum[x]++;x+=lowbit(x);}
}
void pre()
{memset(Div,0,sizeof(Div));memset(L,0,sizeof(L));for(int i=1;i<=n;i++) {int g=d[i];for(int j=0;j<num;j++) {if((ll)prime[j]*prime[j]>g)break;if(g%prime[j]==0) {L[i]=max(L[i],Div[prime[j]]);Div[prime[j]]=i;while(g%prime[j]==0)g/=prime[j];}}if(g>1) {L[i]=max(L[i],Div[g]);Div[g]=i;}// cout<<Div[2]<<"==="<<endl;
//        int st=sqrt(d[i]+0.5);
//        for(int j=2;j<=st;j++) {
//            if(d[i]%j==0) {
//                L[i]=max(L[i],Div[j]);
//                if(d[i]/j!=j) {
//                    L[i]=max(L[i],Div[d[i]/j]);
//                }
//                Div[j]=i;
//                Div[d[i]/j]=i;
//            }
//        }
//        if(d[i]>1) {
//            L[i]=max(L[i],Div[d[i]]);
//            Div[d[i]]=i;
//        }}int N=200000;for(int i=1;i<=N;i++)Div[i]=n+1;for(int i=1;i<=n;i++)R[i]=n+1;for(int i=n;i>=1;i--) {int g=d[i];for(int j=0;j<num;j++) {if((ll)prime[j]*prime[j]>g)break;if(g%prime[j]==0) {R[i]=min(R[i],Div[prime[j]]);Div[prime[j]]=i;while(g%prime[j]==0)g/=prime[j];}}if(g>1) {R[i]=min(R[i],Div[g]);Div[g]=i;}
//        int st=sqrt(d[i]+0.5);
//        for(int j=2;j<=st;j++) {
//            if(d[i]%j==0) {
//                R[i]=min(R[i],Div[j]);
//                if(d[i]/j!=j) {
//                    R[i]=min(R[i],Div[d[i]/j]);
//                }
//                Div[j]=i;
//                Div[d[i]/j]=i;
//            }
//        }
//        if(d[i]>1) {
//            R[i]=min(R[i],Div[d[i]]);
//            Div[d[i]]=i;
//        }}for(int i=1;i<=n;i++) {s[0][i-1].l=L[i];s[0][i-1].r=R[i];s[1][i-1].l=L[i];s[1][i-1].r=i;s[2][i-1].l=i;s[2][i-1].r=R[i];}
}
void solve(int c)
{memset(sum,0,sizeof(sum));int cnt=0;for(int i=0;i<m;i++) {int st=q[i].l;int ed=q[i].r;//cout<<st<<" "<<ed<<endl;while(cnt<n && s[c][cnt].r<=ed) {if(s[c][cnt].l>0)update(s[c][cnt].l);cnt++;}//cout<<cnt<<endl;ans[c][i].l=Sum(ed)-Sum(st-1);ans[c][i].id=q[i].id;}
}
bool cmp(const Node& a,const Node& b)
{return a.id<b.id;
}
int main()
{getPrime();while(~scanf("%d%d",&n,&m)) {if(n==0 && m==0) break;for(int i=1;i<=n;i++) {scanf("%d",&d[i]);}pre();for(int i=0;i<m;i++) {scanf("%d%d",&q[i].l,&q[i].r);q[i].id=i;}// cout<<1<<endl;
//        cout<<s[0][0].l<<" "<<s[0][0].r<<endl;
//        cout<<s[0][1].l<<" "<<s[0][1].r<<endl;
//        cout<<s[0][2].l<<" "<<s[0][2].r<<endl;sort(q,q+m);sort(s[0],s[0]+n);sort(s[1],s[1]+n);sort(s[2],s[2]+n);for(int i=0;i<3;i++)solve(i);// cout<<2<<endl;sort(q,q+m,cmp);sort(ans[0],ans[0]+m,cmp);sort(ans[1],ans[1]+m,cmp);sort(ans[2],ans[2]+m,cmp);for(int i=0;i<m;i++) {printf("%d\n",(q[i].r-q[i].l+1)-ans[1][i].l-ans[2][i].l+ans[0][i].l);}}return 0;
}


这篇关于hdu4777 Rabbit Kingdom 离线树状数组 求询问区间内的区间数的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1108761

相关文章

Java数组初始化的五种方式

《Java数组初始化的五种方式》数组是Java中最基础且常用的数据结构之一,其初始化方式多样且各具特点,本文详细讲解Java数组初始化的五种方式,分析其适用场景、优劣势对比及注意事项,帮助避免常见陷阱... 目录1. 静态初始化:简洁但固定代码示例核心特点适用场景注意事项2. 动态初始化:灵活但需手动管理代

C++中初始化二维数组的几种常见方法

《C++中初始化二维数组的几种常见方法》本文详细介绍了在C++中初始化二维数组的不同方式,包括静态初始化、循环、全部为零、部分初始化、std::array和std::vector,以及std::vec... 目录1. 静态初始化2. 使用循环初始化3. 全部初始化为零4. 部分初始化5. 使用 std::a

shell编程之函数与数组的使用详解

《shell编程之函数与数组的使用详解》:本文主要介绍shell编程之函数与数组的使用,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录shell函数函数的用法俩个数求和系统资源监控并报警函数函数变量的作用范围函数的参数递归函数shell数组获取数组的长度读取某下的

C++原地删除有序数组重复项的N种方法

《C++原地删除有序数组重复项的N种方法》给定一个排序数组,你需要在原地删除重复出现的元素,使得每个元素只出现一次,返回移除后数组的新长度,不要使用额外的数组空间,你必须在原地修改输入数组并在使用O(... 目录一、问题二、问题分析三、算法实现四、问题变体:最多保留两次五、分析和代码实现5.1、问题分析5.

Java中数组转换为列表的两种实现方式(超简单)

《Java中数组转换为列表的两种实现方式(超简单)》本文介绍了在Java中将数组转换为列表的两种常见方法使用Arrays.asList和Java8的StreamAPI,Arrays.asList方法简... 目录1. 使用Java Collections框架(Arrays.asList)1.1 示例代码1.

Python依赖库的几种离线安装方法总结

《Python依赖库的几种离线安装方法总结》:本文主要介绍如何在Python中使用pip工具进行依赖库的安装和管理,包括如何导出和导入依赖包列表、如何下载和安装单个或多个库包及其依赖,以及如何指定... 目录前言一、如何copy一个python环境二、如何下载一个包及其依赖并安装三、如何导出requirem

C++一个数组赋值给另一个数组方式

《C++一个数组赋值给另一个数组方式》文章介绍了三种在C++中将一个数组赋值给另一个数组的方法:使用循环逐个元素赋值、使用标准库函数std::copy或std::memcpy以及使用标准库容器,每种方... 目录C++一个数组赋值给另一个数组循环遍历赋值使用标准库中的函数 std::copy 或 std::

C++初始化数组的几种常见方法(简单易懂)

《C++初始化数组的几种常见方法(简单易懂)》本文介绍了C++中数组的初始化方法,包括一维数组和二维数组的初始化,以及用new动态初始化数组,在C++11及以上版本中,还提供了使用std::array... 目录1、初始化一维数组1.1、使用列表初始化(推荐方式)1.2、初始化部分列表1.3、使用std::

C++ Primer 多维数组的使用

《C++Primer多维数组的使用》本文主要介绍了多维数组在C++语言中的定义、初始化、下标引用以及使用范围for语句处理多维数组的方法,具有一定的参考价值,感兴趣的可以了解一下... 目录多维数组多维数组的初始化多维数组的下标引用使用范围for语句处理多维数组指针和多维数组多维数组严格来说,C++语言没

Java 字符数组转字符串的常用方法

《Java字符数组转字符串的常用方法》文章总结了在Java中将字符数组转换为字符串的几种常用方法,包括使用String构造函数、String.valueOf()方法、StringBuilder以及A... 目录1. 使用String构造函数1.1 基本转换方法1.2 注意事项2. 使用String.valu