HDU 5288 OO’s Sequence 区间个数

本文介绍了一种解决特定区间内非倍数计数问题的高效算法。通过预处理约数并采用二分查找,实现了O(n*log(n))的时间复杂度。文章详细解释了算法思路与实现步骤。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题意:对一个给定的数列,定义一个函数f(l, r),表示在[l, r]这个区间内,满足不是区间内任何其他数的倍数的数的个数。比如一个区间内的数为2,3,4,那么函数值为2。因为2和3都不是其他任何数的倍数,但4是2的倍数。现在要求所有的区间的函数值的和。






看似比较复杂,不过我们可以换一种思路。与其求区间的函数值之和,不如计算每个数被算了多少次。打个比方,第一个样例是1,2,3,4,5,在[2,3]、[2,4]、[3,4]这个区间内,3都是符合要求的数,于是3就被计算了3次(当然并不是只有这3个区间,只是举个例子)。所以我们可以计算数列里每个数在多少个区间内会被计算。然后再加起来,这样外层的复杂度就确定为O(n)了。

接下来就是怎么求每个数被计算多少次的问题。直观上看,如果我们要求某个数所能被计算的区间,首先这个区间不能含有它的约数。进一步想,如果我们能求出这个数左边离它最近的且是它的约数的数的位置L,以及右边离它最近的且是它的约数的数的位置R,这样在[L + 1, R - 1]这个区间内,任意一个包含该数的子区间均可以让这个数被计算一次。设这个数的位置为pos,于是这个数被计算了(pos - L)*(R - pos)次。现在问题转化为了怎么求出每个数的上述L和R值。

考虑到每个数最大也只有10000,不由得想到去预处理。首先预处理出1到10000每个数的约数,然后,对于数列里某个数,枚举它的所有约数(最多64个),用vector记录每个约数出现在哪些位置上(即下标),再二分找到左边最近的和右边最近的位置,更新离它最近的约数的位置即可。这样只有log(n)复杂度。注意左边更新是取最大值,右边更新是取最小值。初始化为左边为0,右边为n + 1。最终复杂度为O(n*log(n))。





#include <cstdio>
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>
#include <stack>
#include <vector>
#include <map>
#include <set>
using namespace std;

const int MAX = 1e5 + 5;
const int MOD = 1e9 + 7;
vector <int> yueshu[MAX], pos[MAX];

int n, a[MAX];

void initial() //预处理约数
{
    for(int i = 0; i < MAX; i++)
        yueshu[i].clear();
    for(int i = 1; i <= 10000; i++)
    {
        for(int j = 1; j*j <= i; j++)
        {
            if(i%j == 0)
            {
                yueshu[i].push_back(j);
                if(j*j != i)
                    yueshu[i].push_back(i/j);
            }
        }
    }
}

void input()
{
    for(int i = 1; i <= n; i++)
        scanf("%d", &a[i]);
}

void solve()
{
    for(int i = 0; i < MAX; i++)
        pos[i].clear();
    for(int i = 1; i <= n; i++) //记录每个数出现的所有下标
        pos[a[i]].push_back(i);
    long long ans = 0;
    for(int i = 1; i <= n; i++)
    {
        int maxl = 0, minr = n + 1;
        for(int j = 0; j < yueshu[a[i]].size(); j++) //枚举约数
        {
            int temp = yueshu[a[i]][j];
            if(pos[temp].size() == 0)
                continue;
            if(pos[temp][0] > i)
            {
                minr = min(minr, pos[temp][0]);
                continue;
            }
            if(pos[temp][pos[temp].size() - 1] < i)
            {
                maxl = max(maxl, pos[temp][pos[temp].size() - 1]);
                continue;
            }
            int l = 0, r = pos[temp].size() - 1, mid;
            while(l < r - 1) //二分找离a[i]左边和右边最近的当前约数
            {
                mid = (l + r)/2;
                if(pos[temp][mid] < i)
                    l = mid;
                else
                    r = mid;
            }
            if(r + 1 < pos[temp].size()) //更新两边最近约数的位置
                minr = min(minr, pos[temp][r + 1]);
            if(pos[temp][r] > i)
                minr = min(minr, pos[temp][r]);
            if(pos[temp][l] < i)
                maxl = max(maxl, pos[temp][l]);
        }
        ans = (ans + (long long)(i - maxl)*(minr - i))%MOD;
    }
    printf("%I64d\n", ans);
}

int main()
{
    initial();
    while(scanf("%d", &n) != EOF)
    {
        input();
        solve();
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值