HDU 4630 No Pain No Game 离线+线段树

最新推荐文章于 2019-09-21 03:13:32 发布

little_胖

最新推荐文章于 2019-09-21 03:13:32 发布

阅读量497

点赞数

分类专栏：线段树

本文链接：https://blog.csdn.net/wangjie_wang/article/details/10468845

版权

线段树专栏收录该内容

5 篇文章 0 订阅

订阅专栏

题意：给出n个数，以及一些区间，求区间中两两数字间的最大公约数。

思路：用pre[j]表示约数j上一次出现的位置，则当约数 j 出现时，它上一次出现的位置所能得到的最大公约数也就是max(j,原来的值)。所以求区间的最大公约数也就成了求区间的最大值，先将所要访问的区间按r由小到大排序，用一颗线段树来维护，求出这个区间的最大值即可。(求约数也可在logn的时间内求出来)。

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <map>
#include <set>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#include <cstdlib>
#include <algorithm>
using namespace std;
typedef __int64 int64;
typedef long long ll;
#define M 50005
#define N 1000005
#define max_inf 0x7f7f7f7f
#define min_inf 0x80808080
#define mod 1000000007
#define lc rt<<1
#define rc rt<<1|1

struct Node
{
	int l , r , pos;
}q[M];
struct node
{
	int l , r , val;
}tree[M<<2];
int arr[M] , n , pre[M] , m , ans[M];

void Build(int rt , int l , int r)
{
	tree[rt].l = l;
	tree[rt].r = r;
	tree[rt].val = 0;
	if (l == r)return;
	int mid = (l+r)>>1;
	Build(lc , l , mid);
	Build(rc , mid+1 , r);
}

bool cmp(Node a , Node b)
{
	return a.r < b.r;
}

void Updata(int rt , int k , int val)
{
	if (tree[rt].l == k && tree[rt].r == k)
	{
		tree[rt].val = max(tree[rt].val,val);
		return;
	}
	int mid = (tree[rt].l+tree[rt].r)>>1;
	if (k <= mid)Updata(lc , k , val);
	else Updata(rc , k , val);
	tree[rt].val = max(tree[lc].val,tree[rc].val);
}

int Query(int rt , int l , int r)
{
	if (tree[rt].l == l && tree[rt].r == r)return tree[rt].val;
	int mid = (tree[rt].l+tree[rt].r)>>1;
	if (r <= mid)return Query(lc,l,r);
	if (l > mid)return Query(rc,l,r);
	return max(Query(lc,l,mid),Query(rc,mid+1,r));
}

int main()
{
	int i , j , t;
	scanf("%d",&t);
	while (t--)
	{
		scanf("%d",&n);
		Build(1,1,n);
		for (i = 1 ; i <= n ; i++)scanf("%d",arr+i);
		scanf("%d",&m);
		for (i = 0 ; i < m ; i++)
		{
			q[i].pos = i;
			scanf("%d%d",&q[i].l,&q[i].r);
		}
		sort(q , q+m , cmp);
		memset(pre , -1 , sizeof pre);
		int cnt = 0;
		for (i = 1 ; i <= n ; i++)
		{
			for (j = 1 ; j*j <= arr[i] ; j++)
			{
				if (arr[i]%j == 0)
				{
					if (pre[j] != -1)Updata(1,pre[j],j);
					pre[j] = i;

					int temp = arr[i]/j;
					if (temp != j && pre[temp] != -1)Updata(1,pre[temp],temp);
					pre[temp] = i;
				}
			}

			while (cnt < m && i == q[cnt].r)
			{
				ans[q[cnt].pos] = Query(1,q[cnt].l,q[cnt].r);
				cnt++;
			}
		}
		for (i = 0 ; i < m ; i++)printf("%d\n",ans[i]);
	}
	return 0;
}