[HDU5726]GCD

题目描述
给出一个有N（N≤100,000）个整数的序列A1,…,An(0<Ai≤1000,000,000)。然后回答 Q(Q≤100,000)个询问，每个询问一个区间[L, R]的GCD值是多少，即GCD(Al, Al+1,…, Ar)，并回答在整个序列上有多少个区间[L’, R’](1<=L’ < R’ <= N)满足GCD(Al’, …Ar’) = GCD(Al, Al+1,…, Ar)

输入
第1行：1个整数T，表示测试数据的组数。每组数据的格式为：

第1行1个整数N，表示序列中元素的个数

第2行：N个整数Ai(0<Ai≤1000,000,000)

第3行：1个整数Q

接下来Q行，每行2个整数L，R，表示询问的区间的左、右端点

输出
每组测试数据，先在一行上输入 “Case #:t”，t是第几组测试数据，从1开始编号

接下来输出每组测试数据的答案，2个整数，第1行是询问的区间[L, R]的GCD值，第2个整数表示有多少个区间的GCD值与它相等。

这题主要是区间GCD的套路，固定一个点为左（右）区间的端点，共有logn种不同的区间GCD，于是共有O(nlogn)种GCD。
本来可以利用上一个点的答案，计算GCD可以做到O(Nlogn),但是本人很懒，直接二分套RMQ O（nlog^3n）玄学卡常。

#include<bits/stdc++.h>
#define maxn 100005
#define lim 17
#define LL long long
using namespace std;

int n,Q,a[maxn];
int st[lim][maxn],lg[maxn];
inline int gcd(int a,int b){ return !b ? a : gcd(b , a % b); }
map<int,LL>M;
int Query(int u,int v)
{
	int tmp = lg[v-u+1];
	return gcd(st[tmp][u], st[tmp][v-(1<<tmp)+1]);
}

int main()
{
	int T;scanf("%d",&T);
	for(int t=1;t<=T;t++)
	{
		M.clear();
		scanf("%d",&n);
		for(int i=1;i<=n;i++) scanf("%d",&a[i]),st[0][i] = a[i];
		for(int j=1;j<lim;j++)
			for(int i=1;i<=n;i++)
				if((i+(1<<j-1)) <= n)
					st[j][i] = gcd(st[j-1][i] , st[j-1][i+(1<<j-1)]);
		for(int i=1,sum=2,LG=0;i<=n;lg[i]=LG,i++) if(i==sum) sum<<=1,LG++;
		for(int i=1,j;i<=n;i++)
		{
			j=i;
			for(int nowgcd=a[i];j<=n;)
			{
				int L = j, R = n , Mid;
				for(;L<R;)
				{
					Mid = (L+R+1) >> 1;
					if(Query(i,Mid) < nowgcd) R = Mid - 1;
					else L = Mid;
				}
				M[nowgcd] += (L-j+1);
				j = L + 1;
				if(j<=n) nowgcd = Query(i,j);
			}
		}
		
		int Q;
		printf("Case #%d:\n",t);
		scanf("%d",&Q);
		for(;Q--;)
		{
			int x,y;
			scanf("%d%d",&x,&y);
			int gcd = Query(x,y);
			printf("%d %lld\n",gcd,M[gcd]);
		}
	}
}