hdu 5381 The sum of gcd 莫队+预处理
The sum of gcd
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)
Let f(l,r)=∑ri=l∑rj=igcd(ai,ai+1....aj)
First line has one integers n
Second line has n integers Ai
Third line has one integers Q,the number of questions
Next there are Q lines,each line has two integers l,r
1≤T≤3
1≤n,Q≤104
1≤ai≤109
1≤l<r≤n
5
1 2 3 4 5
3
1 3
2 3
1 4
4
4 2 6 9
3
1 3
2 4
2 3
6
16
18
23
10
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<iostream>
#include<cstdio>
#include<cmath>
#include<string>
#include<queue>
#include<algorithm>
#include<stack>
#include<cstring>
#include<vector>
#include<list>
#include<set>
#include<map>
using namespace std;
#define ll long long
#define pi (4*atan(1.0))
#define eps 1e-14
#define bug(x) cout<<"bug"<<x<<endl;
const int N=1e4+,M=4e6+,inf=;
const ll INF=1e18+,mod=1e9+; /// 数组大小 int a[N];
int n,pos[N],k;
vector<pair<int,int> >l[N],r[N];
struct is
{
int l,r,now;
bool operator <(const is &b)const
{
if(pos[l]!=pos[b.l])
return pos[l]<pos[b.l];
return r<b.r;
}
}p[N];
ll out[N],ans;
void prex(int L,int R,int flag)
{
ll sum=;
int p=L;
for(int i=;i<l[L].size();i++)
{
int k=min(R,l[L][i].first);
if(k>=p)
sum+=1LL*(k-p+)*l[L][i].second;
if(k>=R)break;
p=k+;
}
if(flag==)ans+=sum;
else ans-=sum;
} void nexx(int L,int R,int flag)
{
ll sum=;
int p=R;
for(int i=;i<r[R].size();i++)
{
int k=max(r[R][i].first,L);
if(p>=k)
sum+=1LL*(p-k+)*r[R][i].second;
if(k<=L)break;
p=k-;
}
if(flag==)ans+=sum;
else ans-=sum;
}
void init()
{
for(int i=;i<=n;i++)
l[i].clear(),r[i].clear();
ans=;
// 预处理l
l[n].push_back(make_pair(n,a[n]));
for(int i=n-;i>=;i--)
{
int g=a[i];
int p=i;
for(int j=;j<l[i+].size();j++)
{
int k=__gcd(g,l[i+][j].second);
if(g!=k)l[i].push_back(make_pair(p,g));
g=k;
p=l[i+][j].first;
}
l[i].push_back(make_pair(p,g));
}
//预处理r
r[].push_back(make_pair(,a[]));
for(int i=;i<=n;i++)
{
int g=a[i];
int p=i;
for(int j=;j<r[i-].size();j++)
{
int k=__gcd(g,r[i-][j].second);
if(g!=k)r[i].push_back(make_pair(p,g));
g=k;
p=r[i-][j].first;
}
r[i].push_back(make_pair(p,g));
}
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
k=sqrt(n);
for(int i=; i<=n; i++)
scanf("%d",&a[i]),pos[i]=(i-)/k+;
init();
int q;
scanf("%d",&q);
for(int i=; i<=q; i++)
scanf("%d%d",&p[i].l,&p[i].r),p[i].now=i;
sort(p+,p++q);
int L=,R=;
for(int i=; i<=q; i++)
{
while(L<p[i].l)
{
prex(L,R,);
L++;
}
while(L>p[i].l)
{
L--;
prex(L,R,);
}
while(R>p[i].r)
{
nexx(L,R,);
R--;
}
while(R<p[i].r)
{
R++;
nexx(L,R,);
}
out[p[i].now]=ans;
}
for(int i=; i<=q; i++)
printf("%lld\n",out[i]);
}
return ;
}
hdu 5381 The sum of gcd 莫队+预处理的更多相关文章
- HDOJ 5381 The sum of gcd 莫队算法
大神题解: http://blog.csdn.net/u014800748/article/details/47680899 The sum of gcd Time Limit: 2000/1000 ...
- hdu 5381 The sum of gcd(线段树+gcd)
题目链接:hdu 5381 The sum of gcd 将查询离线处理,依照r排序,然后从左向右处理每一个A[i],碰到查询时处理.用线段树维护.每一个节点表示从[l,i]中以l为起始的区间gcd总 ...
- HDU 5381 The sum of gcd (技巧,莫队算法)
题意:有一个含n个元素的序列,接下来有q个询问区间,对每个询问区间输出其 f(L,R) 值. 思路: 天真单纯地以为是道超级水题,不管多少个询问,计算量顶多就是O(n2) ,就是暴力穷举每个区间,再直 ...
- hdu 5381 The sum of gcd 2015多校联合训练赛#8莫队算法
The sum of gcd Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) T ...
- hdu 4676 Sum Of Gcd 莫队+phi反演
Sum Of Gcd 题目连接: http://acm.hdu.edu.cn/showproblem.php?pid=4676 Description Given you a sequence of ...
- hdu 4676 Sum Of Gcd 莫队+数论
题目链接 给n个数, m个询问, 每个询问给出[l, r], 问你对于任意i, j.gcd(a[i], a[j]) L <= i < j <= R的和. 假设两个数的公约数有b1, ...
- hdu5381 The sum of gcd]莫队算法
题意:http://acm.hdu.edu.cn/showproblem.php?pid=5381 思路:这个题属于没有修改的区间查询问题,可以用莫队算法来做.首先预处理出每个点以它为起点向左和向右连 ...
- HDU-4676 Sum Of Gcd 莫队+欧拉函数
题意:给定一个11~nn的全排列AA,若干个询问,每次询问给出一个区间[l,r][l,r],要求得出∑l≤i<j≤r gcd(Ai,Aj)的值. 解法:这题似乎做的人不是很多,蒟蒻当然不会做只 ...
- hdu 5381 The sum of gcd
知道对于一个数列,如果以x为左(右)端点,往右走,则最多会有log(a[x])个不同的gcd,并且有递减性 所以会分成log段,每一段的gcd相同 那我们可以预处理出对于每一个位置,以这个位置为左端点 ...
随机推荐
- [LeetCode] 152. Maximum Product Subarray_Medium tag: Dynamic Programming
Given an integer array nums, find the contiguous subarray within an array (containing at least one n ...
- JS参差不齐的数组
<html><head> <title>参差不齐的数组</title> <meta charset="utf-8"> & ...
- Java缓存学习之六:Spring单独使用EhCache
通过在Application Context中配置EhCacheManagerFactoryBean和EhCacheFactoryBean,我们就可以把对应的EhCache的CacheManage ...
- 生产者消费者模型——wait/notify/notifyAll使用
告警系统架构如下 1. 数据处理系统处理完原始数据并入库后,发送消息到kafka系统: 2. 告警生产者从kafka系统查询消息存入告警消息队列: 3. 告警消费者从告警消息队列查询消息进行处理. 这 ...
- Linux基础命令---more
more 将内容较长的文本文件内容分屏显示,支持定位关键字.此命令的适用范围:RedHat.RHEL.Ubuntu.CentOS.SUSE.openSUSE.Fedora. 1.语法 ...
- 干货:Java并发编程系列之synchronized(一)
1. 使用方法 synchronized 是 java 中最常用的保证线程安全的方式,synchronized 的作用主要有三方面: 确保线程互斥的访问代码块,同一时刻只有一个方法可以进入到临界区 保 ...
- VM三种网络连接方式
bridge:这种方式最简单,直接将虚拟网卡桥接到一个物理网卡上面,和Linux下一个网卡 绑定两个不同地址类似,实际上是将网卡设置为混杂模式,从而达到侦听多个IP的能力. 在此种模式下,虚拟机内部的 ...
- elasticsearch分词器ik
1. 下载和es配套的版本 git clone https://github.com/medcl/elasticsearch-analysis-ik 2. 编译 cd elasticsearch-an ...
- Class<T>
首先,什么是类类型? 可见: https://www.cnblogs.com/yanze/p/9717658.html Class<T>即T的类类型 如何获取Class<T>? ...
- 20145104张家明 《Java程序设计》第三次实验设计
合作伙伴是20145103 下面是我们的git成果 首先下载他托管上去的代码 然后运行下载的代码 之后对下载的代码进行修改 然后推送上去 下载修改后的代码并运行 •软件工程是把系统的.有序的.可量化的 ...