题目大概说给一个包含n个数的序列,多次询问有多少个区间GCD值等于某个区间的gcd值。
任何一个区间不同的GCD个数是log级别的,因为随着右端点向右延伸GCD是单调不增的,而每次递减GCD至少除以2。
考虑固定左端点,最多就nlogn种GCD,可以直接把所有区间GCD值预处理出来,用map存储各种GCD值的个数,查询时直接输出。
具体是这样处理的:枚举左端点,进行若干次二分查找,看当前GCD值最多能延伸到哪儿,进而统计当前GCD值的数量。
而求区间GCD,用ST表,预处理一下,就能在O(1)时间复杂度求出任意区间的gcd了。
注意第二个结果个数有可能是long long, 还有输出别写错了,二分的时候应该是找最右边的那个, 别加错了
#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 10;
int n; int a[maxn];
int dp[maxn][50];
void st(){
for(int j = 1; (1 << j) <= n; j ++){
for(int i = 1; i + (1 << j) - 1 <= n ; i ++){
dp[i][j] = __gcd(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1]);
}
}
}
int query(int l , int r){
int len = r - l + 1;
int k = 0 ;
while((1 << (k + 1)) <= len) k ++;
return __gcd(dp[l][k], dp[r - (1 << k) + 1][k]);
}
map<int, long long> mp;
int binary_search(int pos, int l, int r, int num){
int s = r + 1;
while(l <= r){
int mid = (l + r) / 2;
if(query(pos, mid) >= num){
l = mid + 1; s = mid;
}
else r = mid - 1;
}
return s;
}
void solve(){
for(int i = 1 ; i <= n ; i ++){
int num = a[i];
int l = i; int r = n; int pre = i;
while(l <= n){
l = binary_search(i, l ,r ,num);
// cout << "i = " << i << " l =" << l << endl;
// system("pause");
mp[num] += (1ll * (l - pre + 1));
l ++; num = __gcd(num, a[l]);
pre = l;
}
}
}
int main(){
// freopen("data.out", "r", stdin);
// freopen("data1.out", "w", stdout);
int t ; cin >> t; int cas = 1;
while(t --){
scanf("%d", &n); mp.clear();
for(int i = 1; i <= n ; i++){
scanf("%d", &a[i]); dp[i][0] = a[i];
}
st();
solve();
int q; scanf("%d", &q);printf("Case #%d:\n", cas ++);
while(q --){
int l, r; scanf("%d%d", &l, &r);
int gcd = query(l, r);
printf("%d %lld\n", gcd, mp[gcd]);
}
}
return 0;
}