Codeforces Round #717 (Div. 2) D. Cut

最新推荐文章于 2022-04-02 22:53:43 发布

TURNINING

最新推荐文章于 2022-04-02 22:53:43 发布

阅读量168

点赞数 1

CC 4.0 BY-SA版权

分类专栏：数据结构

本文链接：https://blog.csdn.net/TURNINING/article/details/118723277

数据结构专栏收录该内容

8 篇文章

订阅专栏

本文介绍了一种解决数组区间划分问题的方法，通过贪心策略和双指针技巧，计算在保持区间内元素互质的前提下，划分次数的最小值。涉及的算法包括线段树求公共素数最大值和动态规划求最长区间。适合于对算法竞赛、数据结构和计算机数学感兴趣的读者。

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

传送门

题意：数组有n个数，有q次询问，每次询问给你一个区间l，r，问[l, r]最少的划分次数是多少。划分的规则为：每个划分后的区间每个数都互质。

思路：对于每次询问，要想划分次数最少显然我们得贪心的令子区间的长度最长。求子区间的最长长度思路很简单，从i开始一直往后遍历一直到gcd不为1为止(互质)。现在我们就可以得到以i开头，最长的划分区间了。对于每次询问我们就可以直接算次数了。但我们能直接暴力的算次数吗？显然不行，当每个数都为1时可以把我们卡上天。那怎么办？倍增

现在我们思路都已经很清楚了，那么怎么求每个数的最长区间呢？(n^2暴力求肯定不行的吧？) 对于这种区间有限制的贪心问题，考虑双指针，这样我们就可以尽可能的减少移动次数。设双指针左右区间为l, r当加上arr[r]后若公共素数最大值为1（求公共素数最大值我套了个线段树）那么我们就移动右指针，大于1就说明l最长右点就是r-1(为什么是r-1，因为我们每次都会把r+1)，到此问题就全部解决了。这个双指针我没想到，我好菜。。。。

#include<bits/stdc++.h>
using namespace std;

#define lsn (u << 1)
#define rsn (u << 1 | 1)
#define mid (l + r >> 1)

typedef long long ll;
typedef unsigned long long ull;

typedef pair<int, int> P;
typedef pair<double, double> PD;

const int MAXN = 1e5 + 10;
const int MAX_LEN = 100000 + 10;
const int MAX_LOG_V = 22;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const double eps = 1e-7;
const ull B = 100000007;

int n, q; 
int arr[MAXN];
int nxt[MAXN];
int f[MAXN][21];

int tr[4*MAXN];

void pushup(int u) { tr[u] = max(tr[lsn], tr[rsn]); }

void add(int u, int l, int r, int p, int v) {
    if(l == r && l == p) tr[u] += v;
    else {
        if(p <= mid) add(lsn, l, mid, p, v);
        else add(rsn, mid+1, r, p, v);
        pushup(u);
    }
}

void add(int x) {
    int last = -1;
    while(x > 1) {
        if(nxt[x] != last) {
            add(1, 1, 1e5, nxt[x], 1);
            last = nxt[x];
        }
        x /= nxt[x];
    }
}

void del(int x) {
    int last = -1;
    while(x > 1) {
        if(nxt[x] != last) {
            add(1, 1, 1e5, nxt[x], -1);
            last = nxt[x];
        }
        x /= nxt[x];
    }
}

void solve() {
    scanf("%d %d", &n, &q);
    for(int i = 1; i <= n; i++) scanf("%d", &arr[i]);
    int r = 1;
    for(int i = 1; i <= n; i++) {
        while(r <= n && tr[1] < 2) {
            add(arr[r++]);
        }
        f[i][0] = r - 1;
        if(tr[1] < 2 && r == n + 1) f[i][0] = n + 1;
        del(arr[i]); 
    }
    for(int k = 0; k < 20; k++) {
        for(int i = 1; i <= n; i++) {
            if(f[i][k] == n+1) f[i][k+1] = n+1;
            else f[i][k+1] = f[f[i][k]][k]; 
        }
    }
    for(int i = 1; i <= q; i++) {
        int l, r; scanf("%d %d", &l, &r);
        int ans = 0;
        for(int k = 20; k >= 0; k--) {
            if(f[l][k] <= r) {
                l = f[l][k];
                ans += (1 << k);
            }
        }
        printf("%d\n", ans + 1);
    }
}   


int main() {
    //ios::sync_with_stdio(false);
    for(int i = 2; i < MAXN; i++) {
        if(!nxt[i]) {
            for(int j = i; j < MAXN; j += i) {
                nxt[j] = i;
            }
        }
    }
    int t = 1; //scanf("%d", &t);
    while(t--) {
        solve();
    }
    return 0;
}