题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5726
题目大意:给你一个数列,问有多少个子区间的GCD等于给定区间的GCD
思路:区间GCD用线段树,查询多少个子集可以暴力预处理出所有区间的GCD,处理时,利用滚动的map记录到达当前位置的gcd和数量,代码里面用f表示
#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <cmath>
#include <stack>
#include <queue>
#include <algorithm>
#include <vector>
#include <map>
#include <set>
#include <stdlib.h>
#include <iomanip>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
#define maxn 100005
#define MOD 1000000007
#define mem(a , b) memset(a , b , sizeof(a))
#define LL long long
#define ULL unsigned long long
#define FOR(i , n) for(int i = 1 ; i<= n ; i ++)
typedef pair<int , int> pii;
int t , n , q , a[maxn];
map<int , LL>M;
map<int , LL>f[maxn];
map<int , LL>::iterator it;
struct node
{
int l , r , g;
}Tr[maxn*3];
int gcd(int a , int b)
{
if(b == 0) return a;
return gcd(b , a % b);
}
void build(int ro , int l , int r)
{
Tr[ro].l = l , Tr[ro].r = r;
if(l == r)
{
Tr[ro].g = a[l];
return ;
}
int mid = (l + r) >> 1;
build(ro << 1 , l , mid);
build(ro << 1 | 1 , mid + 1 , r);
Tr[ro].g = gcd(Tr[ro<<1].g , Tr[ro<<1|1].g);
}
int query(int ro , int l , int r)
{
if(Tr[ro].r < l || Tr[ro].l > r) return 0;
int mid = (Tr[ro].l + Tr[ro].r) >> 1;
if(Tr[ro].l == l && Tr[ro].r == r)
{
return Tr[ro].g;
}
else if(l > mid)
return query(ro << 1 | 1 , l , r);
else if(r <= mid)
return query(ro << 1 , l , r);
else
return gcd(query(ro << 1 , l , mid) , query(ro << 1 | 1 , mid + 1 , r));
}
void solve()
{
for(int i = 0 ; i <= n ; i ++) f[i].clear();
M[a[1]] = 1 , f[1][a[1]] = 1;
for(int i = 2 ; i <= n ; i ++)
{
M[a[i]] += 1;
f[i][a[i]] += 1;
for(it = f[i-1].begin() ; it != f[i-1].end() ; it++)
{
int tg = gcd(a[i] , it -> first);
f[i][tg] += it -> second;
M[tg] += it -> second;
}
}
}
int main()
{
scanf("%d" , &t);
int k = 1;
while(t --)
{
scanf("%d" , &n);
mem(Tr , 0);M.clear();
for(int i = 1 ; i <= n ; i ++) scanf("%d" , &a[i]);
build(1 ,1 , n);
solve();
// cout << query(1 , 4 , 4) << endl;
/*
for(int i = 1 ; i <= n ; i ++)
{
int j = i , g = a[i];
while(j <= n)
{
int l = j , r = n;
while(l < r)
{
int mid = (l + r) >> 1;
int tg = query(1 , i , mid);
if(tg == g) l = mid + 1;
else r = mid;
}
if(query(1 , i , l) != g) l--;
M[g] += l - j + 1;
g = query(1 , i , l + 1);
j = l + 1;
}
}*/
//cout << "YESS" << endl;
scanf("%d" , &q);
printf("Case #%d:\n" , k ++);
for(int i = 1 ; i <= q ; i ++)
{
int l , r;
scanf("%d %d" , &l , &r);
int tg = query(1 , l , r);
printf("%d %lld\n" , tg , M[tg]);
}
}
return 0;
}