#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int MAXN=100009;
int st[MAXN][20];
map<int,LL>ans;
void RMQ(int n){ //RMQ一般就是指ST表
int Len=log2(n); //得到了要维护的深度
for(int i=1;i<=Len;i++){ //i深度,在输入的时候st[j][0]已经存好,这里i直接从1开始
for(int j=1;j<=n;j++){ //j位置,注意输入时存的下标是从1~N
if(j-1+(1<<i)>n)break; //退出条件:从j位起长为(1<<i)区间最后一位置j+(1<<i)-1,故>n跳出
st[j][i]=__gcd(st[j][i-1],st[j+(1<<(i-1))][i-1]);//st[j][i]指第j位起长为2^i长度区间GCD
}
}
}
int gcd(int l,int r){ //传入查找区间
int Len=log2(r-l+1); //在第几深度查找
return __gcd(st[l][Len],st[r-(1<<Len)+1][Len]);//左右区间GCD的GCD即为所求(区间可部分覆盖)
}
//gcd从左端点开始要么是一长段不变 , 要么是直接减少一半以上
void init(int n){ //就是用来更新ANS
for(int i=1;i<=n;i++){ //枚举左界,二分求出每次变化的点
int now=i; //在以i为左界下,now是当前判到的变化点,即由i开始往右扫区间GCD有变化的点,左界i初始化为第一个
while(now<=n){ //变化点未去到尽头就要继续,因为要把以I为左界的右边的区间全部扫完
int l=now,r=n,x=now; //解释x:在以i为左界下,以当前变化点now为右界的GCD记为G1
//则x是在以i为左界下,能使GCD也是G1的最大右界
while(l<=r){ //二分未结束
int mid=l+r>>1; //二分区间中点
if(gcd(i,mid)==gcd(i,now)) //如果i,mid区间与i,now区间GCD下相同
x=mid,l=mid+1; //表明下一变化点在右半区间,则最大右界移到中点,然后左界右移
else r=mid-1; //否则下一变化点在左半区间,则右界左移,最大右界不动
}
ans[gcd(i,now)]+=x-now+1; //一整段gcd都相同,最大右界减当前变化点的长度加进对应答案里计数
now=x+1; //下一次变化点就是当前变化点下GCD不变的最大右界
}
}
}
int main(){
int t,ca=0;cin>>t;
while(t--){
printf("Case #%d:\n",++ca);ans.clear();
int n;cin>>n;for(int i=1;i<=n;i++)cin>>st[i][0];//输入N个数,注意是从1开始存的
RMQ(n); //处理ST表
init(n); //处理GCD映射MAP(GCD值,GCD数量)
int q;cin>>q; //询问数
while(q--){
int l,r;cin>>l>>r;int g=gcd(l,r); //找出区间GCD
printf("%d %lld\n",g,ans[g]); //输出区间找到的GCD及一样GCD的区间数
}
}
}
RMQ预处理gcd(hdu 5726 GCD)
题意:n的数组 , q次询问 , 每次询问(l,r) , 求区间内所有数的gcd , 以及有多少个区间的gcd和这个区间相同(包括自己)
Sample Input
1
5
1 2 4 6 7
4
1 5
2 4
3 4
4 4
Sample Output
Case #1:
1 8
2 4
2 4
6 1