区间最小公倍数

题目链接：Codeforces Hello 2015 (Div.2) - LCM Query

求出某个序列所有给定长度的连续子序列的最小公倍数中的最小值。考虑到有大量询问、无需进行修改，故采用 稀疏表 + 离线处理 的方法。

题目描述

出题人喜欢最小公倍数 (lcm)({\rm lcm})，他提供了一个数列 a1,a2,…,ana_1,a_2,\ldots,a_n，并给出了一些 xx，希望你对每个 xx 都求出以下式子的值：

min⁡i=1n−x+1lcm(ai,ai+1,…,ai+x−1)\min_{i=1}^{n-x+1}{\rm lcm}(a_i,a_{i+1},\ldots,a_{i+x-1})

简单来说求出这个数列所有长度为 xx 的连续子序列的最小公倍数中的最小值（简单来说？

数列 ai{a_i} 的长度为 nn，一个 xx 就是一次询问，共有 mm 次询问，并有以下限制：

1≤n≤2×1041\leq n\leq2\times10^4
1≤m≤1061\leq m\leq10^6
1≤ai≤601\leq a_i\leq60
1≤x≤n1\leq x\leq n

题解

对于一个数 xx，可以将其分解质因数：x=p1b1×p2b2×…×pkbkx=p_1^{b_1}\times p_2^{b_2}\times\ldots\times p_k^{b_k}，其中 pip_i 指第 ii 个素数。于是 xx 就可以表示为一个数列 bb 。用这种表示方法，可以很方便地表示最小公倍数。假设有两个数 x=p1b1×p2b2×…×pkbkx=p_1^{b_1}\times p_2^{b_2}\times\ldots\times p_k^{b_k} 和 y=p1c1×p2c2×…×pkcky=p_1^{c_1}\times p_2^{c_2}\times\ldots\times p_k^{c_k}，那么就有：

lcm(x,y)=p1max⁡(b1,c1)×p2max⁡(b2,c2)×…×pkmax⁡(bk,ck){\rm lcm}(x,y)=p_1^{\max(b_1,c_1)}\times p_2^{\max(b_2,c_2)}\times\ldots\times p_k^{\max(b_k,c_k)}

由于题目要求让结果对 109+710^9 + 7 取模，算出的结果不能直接比较大小。所以利用对数的性质，取对数后再比较对数的大小：

log⁡(x)=b1log⁡(p1)+b2log⁡(p2)+…+bklog⁡(pk)log⁡(y)=c1log⁡(p1)+c2log⁡(p2)+…+cklog⁡(pk)log⁡(x)<log⁡(y)  ⟹  x<y\begin{aligned}\log(x)&=b_1\log(p_1)+b_2\log(p_2)+\ldots+b_k\log(p_k)\\\log(y)&=c_1\log(p_1)+c_2\log(p_2)+\ldots+c_k\log(p_k)\\\log(x)&\lt\log(y)\implies x\lt y\end{aligned}

最后是查询，对于每一个 xx 都重新遍历一遍数列是肯定会 TLE 的，所以要先将每一个区间长度的结果处理出来。大概过程是这样的：创建两个 “指针” ll 和 rr 。一开始，ll 和 rr 都指向第一个元素，之后让 rr 不断向后移动，直至指向第一个 lcm{\rm lcm} 值lcm(al,al+1,…,ar) {\rm lcm}(a_l,a_{l+1},\ldots,a_r) 与之前的 lcm{\rm lcm} 值不同的元素，并用之前的 lcm{\rm lcm} 值更新 ansr−lans_{r-l} 的值。再向右移动 rr 并如此更新，直至数列尾后元素。之后 ll 和 rr 同时指向第二个元素，重复上述操作，再同时指向第三个元素…… 直至同时指向数列的尾部元素。

以数列 12,3,6,25,5,712,3,6,25,5,7 举例：

  1. l=0,r=3l=0,r=3，用 lcm(12,3,6)=12{\rm lcm}(12,3,6)=12 更新 ans3ans_3
  2. l=0,r=5l=0,r=5，用 lcm(12,3,6,25,5)=300{\rm lcm}(12,3,6,25,5)=300 更新 ans5ans_5
  3. l=0,r=6l=0,r=6，用 lcm(12,3,6,25,5,7)=2100{\rm lcm}(12,3,6,25,5,7)=2100 更新 ans6ans_6
  4. l=1,r=3l=1,r=3，用 lcm(3,6)=6{\rm lcm}(3,6)=6 更新 ans2ans_2
  5. ……

这部分比较难搞懂，代码中已经用 6464 道反斜杠划出了重点。

参考代码

lcm-query.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
#include <bits/stdc++.h>
using namespace std;
typedef vector<int> vi;
typedef long long ll;
const int maxn = 2e4 + 1;
const int lg2maxn = 15;
const ll mod = 1e9 + 7;

int lg2[maxn], pw2[lg2maxn], ansi[maxn];
vi rmq[maxn][lg2maxn], ans[maxn], p;
vector<double> lnp;

void init()
{
    for (int i = 2; i != maxn; ++i)
        lg2[i] = lg2[i >> 1] + 1;
    for (int i = 0; i != lg2maxn; ++i)
        pw2[i] = 1 << i;
    for (int i = 2; i <= 60; ++i) {
        bool prime = true;
        for (int j = 2; prime && j * j <= i; ++j)
            if (!(i % j))
                prime = false;
        if (prime)
            p.push_back(i), lnp.push_back(log((double)i));
    }
    fill(ans, ans + maxn, vi(p.size(), 7));
}

vi lcm(const vi& a, const vi& b)
{
    vi res(p.size());
    for (int i = 0, i_ = p.size(); i != i_; ++i)
        res[i] = max(a[i], b[i]);
    return res;
}

vi intvlcm(int l, int r)
{
    int s = lg2[r - l];
    return lcm(rmq[l][s], rmq[r - pw2[s]][s]);
}

vi itov(int x)
{
    vi res(p.size(), 0);
    for (int i = 0, i_ = p.size(); i != i_; ++i)
        while (!(x % p[i]))
            ++res[i], x /= p[i];
    return res;
}

int vtoi(const vi& x)
{
    ll res = 1ll;
    for (int i = 0, i_ = p.size(); i != i_; ++i)
        for (int j = x[i]; j--;)
            res = (res * p[i]) % mod;
    return res;
}

bool cmp(const vi& a, const vi& b)
{
    double x = 0.0, y = 0.0;
    for (int i = 0, i_ = p.size(); i != i_; ++i)
        x += lnp[i] * a[i], y += lnp[i] * b[i];
    return x < y;
}

int main()
{
    init();
    int n, q;
    scanf("%d%d", &n, &q);
    for (int i = 0, t; i != n; ++i)
        scanf("%d", &t), rmq[i][0] = itov(t);
    for (int j = 1; j <= lg2[n]; ++j)
        for (int i = n - pw2[j]; i != -1; --i)
            rmq[i][j] = lcm(rmq[i][j - 1], rmq[i + pw2[j - 1]][j - 1]);
    ////////////////////////////////
    for (int l = 0; l != n; ++l) {
        vi cur = rmq[l][0];
        for (int r = l; r != n;) {
            cur = lcm(cur, rmq[r][0]);
            for (int k = lg2[n - r]; k != -1; --k)
                if (r + pw2[k] <= n && lcm(cur, rmq[r][k]) == cur)
                    r += pw2[k];
            ans[r - l] = min(ans[r - l], cur, cmp);
        }
    }
    ////////////////////////////////
    for (int x; q--;) {
        scanf("%d", &x);
        printf("%d\n", ansi[x] ? ansi[x] : (ansi[x] = vtoi(ans[x])));
    }
}
#​#
Miller–Rabin 素性测试
稀疏表及其应用

Comments

Powered By Valine
v1.5.2