2016多校第四题 GCD(hdu 5726)(*)

题意:给一个数组a,大小为n,接下来有m个询问,每次询问给出l、r,定义f[l,r]=gcd(al,al+1,...,ar),问f[l,r]的值 和 有多少对(l',r')使得f[l',r']=f[l,r]。n<=10万,m<=10万,1<=l<=r<=n,1<=l'<=r'。

题解:官方题解很6啊,对于每个点维护以这个点为右端点的所有区间,将这些区间按照GCD一段段合并在一块,然后递推维护下一个点的信息,将结果保存在map中,时间复杂度就是O(nlog(n))。

#include <set>
#include <map>
#include <queue>
#include <deque>
#include <cmath>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cassert>
#include <iostream>
#include <algorithm>
using namespace std;
#define L(i) i<<1
#define R(i) i<<1|1
#define INF  0x3f3f3f3f
#define pi acos(-1.0)
#define eps 1e-9
#define maxn 100010
#define MOD 1000000007
#define zero(x) (((x)>0?(x):-(x))<eps)

vector <pair <int, int> > gcd[maxn];
map <int, long long> mp;

int n;
int num[maxn];

int main()
{
    int t,C = 1;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i = 1; i <= n; i++)
            scanf("%d",&num[i]);
        mp.clear();
        for(int k = 1; k <= n; k++)
        {
            int last = -1;
            for(int i = 0; i < gcd[k-1].size(); i++)
            {
                int u = gcd[k-1][i].first;
                int v = gcd[k-1][i].second;
                int g = __gcd(u,num[k]);
                if(g == last)
                    continue;
                last = g;
                gcd[k].push_back(make_pair(g,v));
            }
            if(num[k] != last)
                gcd[k].push_back(make_pair(num[k],k));
            for(int i = 0; i < gcd[k].size(); i++)
                mp[gcd[k][i].first] += (i + 1 == gcd[k].size()?k+1:gcd[k][i+1].second) - gcd[k][i].second;
        }
        int q;
        scanf("%d",&q);
        printf("Case #%d:\n",C++);
        while(q--)
        {
            int s,t,k;
            scanf("%d%d",&s,&t);
            for(k = 0; k < gcd[t].size(); k++)
                if(gcd[t][k].second > s)
                    break;
            printf("%d %I64d\n",gcd[t][k-1].first,mp[gcd[t][k-1].first]);
        }
        for (int i = 0; i <= n; ++i)
			gcd[i].clear();
    }
    return 0;
}

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值