HDU5726-GCD 区间GCD+二分

GCD

Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 4611    Accepted Submission(s): 1655

Problem Description

Give you a sequence of  N(N≤100,000) integers :  a1,...,an(0<ai≤1000,000,000). There are  Q(Q≤100,000) queries. For each query  l,r you have to calculate  gcd(al,,al+1,...,ar) and count the number of pairs (l′,r′)(1≤l<r≤N)such that  gcd(al′,al′+1,...,ar′) equal  gcd(al,al+1,...,ar).

 

Input

The first line of input contains a number  T, which stands for the number of test cases you need to solve.

The first line of each case contains a number  N, denoting the number of integers.

The second line contains  N integers,  a1,...,an(0<ai≤1000,000,000).

The third line contains a number  Q, denoting the number of queries.

For the next  Q lines, i-th line contains two number , stand for the  li,ri, stand for the i-th queries.

 

Output

For each case, you need to output “Case #:t” at the beginning.(with quotes,  t means the number of the test case, begin from 1).

For each query, you need to output the two numbers in a line. The first number stands for  gcd(al,al+1,...,ar) and the second number stands for the number of pairs (l′,r′) such that  gcd(al′,al′+1,...,ar′) equal  gcd(al,al+1,...,ar).

链接

一、题意

        给定一个长度为N的数组和Q个查询，对于每个查询，要求输出数组中区间[L,R]的区间GCD和数组中有多少子区间的区间GCD等于该值（包括[L,R]自身）。

二、思路

        对于区间GCD，可以使用RMQ来实现计算。记dp[i][j]为i开始，2的j次幂个数的GCD，则有递推关系dp[i][j] = GCD(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1])。查询时利用RMQ查询。

        对于第二个问题，因为当左界不变时，随着区间的增大，区间GCD不可能增加，所以可以利用二分进行打表处理。枚举左界i，右界j从i开始，然后通过二分找到最大的右界j'，使得[i,j']的区间GCD等于[i,j]的区间GCD，易知这之间的GCD都相等，记gcd为区间[i,j]的区间GCD，map[gcd]为区间GCD等于gcd的区间数，则有map[gcd] += j' - j + 1。

三、代码

#include <iostream>
#include <sstream>
#include <iomanip>
#include <string>
#include <numeric>
#include <vector>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <utility>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>

using namespace std;
typedef long long LL;
const int MAXN = 100100;
const int MOD7 = 1e9 + 7;
const int MOD9 = 1e9 + 9;
const int INF = 2e9;
const double EPS = 1e-6;
const double PI = 3.14159265358979;
const int dir_4r[] = { -1, 1, 0, 0 };
const int dir_4c[] = { 0, 0, -1, 1 };
const int dir_8r[] = { -1, -1, -1, 0, 0, 1, 1, 1 };
const int dir_8c[] = { -1, 0, 1, -1, 1, -1, 0, 1 };

int input[MAXN];
int dp[MAXN][20];
map<int, LL> mp;

int GCD(int a, int b) {
	int r;
	while (b) {
		r = a % b;
		a = b;
		b = r;
	}
	return a;
}

//RMQ初始化
void init(int n) {
	for (int i = 0; i < n; ++i)
		dp[i][0] = input[i];
	for (int j = 1; (1 << j) <= n; ++j)
		for (int i = 0; i + (1 << j) - 1 < n; ++i)
			dp[i][j] = GCD(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1]);
}

int getGcd(int l, int r) {
	int k = 0;
	while ((1 << (k + 1)) <= r - l + 1)
		k++;
	return GCD(dp[l][k], dp[r - (1 << k) + 1][k]);
}

int main() {
	int t, n, q, l, r;
	scanf("%d", &t);
	for (int kase = 1; kase <= t; ++kase) {
		scanf("%d", &n);
		for (int i = 0; i < n; ++i)
			scanf("%d", input + i);
		init(n);

		//二分
		mp.clear();
		for (int i = 0; i < n; ++i) {//左界i
			int j = i;
			while (j < n) {//右界j
				int gcd = getGcd(i, j);
				l = j;
				r = n - 1;
				while (l < r) {
					int mid = (l + r + 1) >> 1;
					if (gcd > getGcd(i, mid))
						r = mid - 1;
					else
						l = mid;
				}
				mp[gcd] += l - j + 1;
				j = l + 1;
			}
		}

		scanf("%d", &q);
		printf("Case #%d:\n", kase);
		for (int i = 0; i < q; ++i) {
			scanf("%d%d", &l, &r);
			l--;
			r--;
			int gcd = getGcd(l, r);
			printf("%d %lld\n", gcd, mp[gcd]);
		}
	}

	//system("pause");
	return 0;
}