hdu5726 GCD 多校1

题目链接：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;
}