HDU 5726 GCD(DP+线段树)

给定一个序列,对于每次查询,有多少个连续的子序列gcd等于x到y的序列的gcd


an记录不同gcd的数目,cot为滚动map用来记录第i-1位为终止位置的子区间的gcd的种类和数目,然后处理出所有lr范围内的gcd的数目

然后对于每次查询的lr的gcd,用线段树查询一下就好了


#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <map>
#include <stack>
#include <queue>
#include <string.h>
#include <string>
#include <iomanip>
#include <cstring>
#include <vector>
#include <functional>
#include <cmath>
#include <sstream>
#include <set>

using namespace std;
#define PI acos(-1.0)
#define sp system("pause")
typedef long long ll;
typedef unsigned long long ull;
const ll mod = 1000000007;
#define PI acos(-1.0)

int tree[100000 * 4 + 20];
vector<int>v;
int gcd(int a, int b)
{
	return b == 0 ? a : gcd(b, a%b);
}

void build_tree(int l, int r, int now)
{
	if (l == r)tree[now] = v[l-1];
	else
	{
		int mid = (l + r) >> 1;
		build_tree(l, mid, now << 1);
		build_tree(mid + 1, r, now << 1 | 1);
		tree[now] = gcd(tree[now << 1], tree[now << 1 | 1]);
	}
}

int query(int l, int r,int now,int ql,int qr)
{
	if (l == ql&&r == qr)return tree[now];
	else
	{
		int mid = (l + r) >> 1;
		if (qr <= mid)return query(l, mid, now << 1, ql, qr);
		else if (ql > mid)return query(mid + 1, r, now << 1 | 1, ql, qr);
		else return gcd(query(l, mid, now << 1, ql, mid), query(mid + 1, r, now << 1 | 1, mid + 1, qr));
	}
}


int main()
{
	int T;
	cin >> T;
	int cot = 1;
	while (T--)
	{
		printf("Case #%d:\n", cot++);
		int n;
		map<int, ll>cot[2];
		map<int, ll>an;
		v.clear();
		cin >> n;
		int pos = 0;
		for (int i = 0; i < n; i++)
		{
			int x;
			scanf("%d", &x);
			v.push_back(x);
			map<int, ll>::iterator po;
			for (po = cot[pos].begin(); po != cot[pos].end(); po++)
			{
				int cc = gcd(po->first, x);
				cot[pos^1][cc] += cot[pos][po->first];
				an[cc] += cot[pos][po->first];
			}
			cot[pos].clear();
			cot[pos^1][x]++;
			an[x]++;
			pos ^= 1;
		}
		build_tree(1, n, 1);
		int m;
		cin >> m;
		while (m--)
		{
			int x, y;
			scanf("%d%d", &x, &y);
			int ans = query(1, n, 1, x, y);
			printf("%d %lld\n", ans, an[ans]);
		}
	}
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值