[DTOI 2023] C. 不见故人（贪心）

[DTOI 2023] C. 不见故人

题目背景

虽然 luanmenglei 已经是成熟的高中生了，但每次提起 luanmenglei 八年级的女朋友时，luanmenglei 都会沉浸在美好的回忆中，不可自拔。

题目描述

给定 $n,k$ 和序列 { a n } \{a_n\} {an​}，你同时有一个临时变量 $x$，你可以进行以下操作若干次（也可以是 $0$ 次），一次操作的流程是：

1. 选定一个区间 $[l,r]$，$\forall i\in [l,r]$，$x\leftarrow \gcd (a_l,a_{l+1},\cdots ,a_r)$。
2. $\forall i\in [l,r]$，$a_i\leftarrow x$。

简而言之，你每次可以选定一个区间并将其中每个数变成这个区间的 $\gcd$。

一次操作的代价是 $r-l+1+k$，现在你希望把这个序列的每个数都变成相等的，求最小代价和。

---

如果您不了解 $\gcd$ 或者多元 $\gcd$ 的含义，可以参照如下定义：

• $\gcd (a_1,a_2,\dots ,a_k)$ 表示 $a_1,a_2,\dots ,a_k$ 的最大公约数，即最大的能同时整除 $a_1,a_2,\dots ,a_k$ 的正整数。

输入格式

第一行两个非负整数 $n,k$。

第二行 $n$ 个数，表示 { a n } \{a_n\} {an​}。

输出格式

一行一个数，表示答案。

样例 #1

样例输入 #1

10 3
2 2 2 1 2 2 2 1 2 2

样例输出 #1

13

样例 #2

样例输入 #2

10 0
2 2 2 1 2 2 2 1 2 3

样例输出 #2

9

样例 #3

样例输入 #3

11 0
2 2 2 1 2 2 2 1 1 3 3

样例输出 #3

10

提示

【样例 1 解释】

操作一次，选择区间 $[1,10]$。

【样例 4】

见附加文件中的 old/old4.in 与 old/old4.out。

该样例满足测试点 $9\sim 12$ 的限制。

【样例 5】

见附加文件中的 old/old5.in 与 old/old5.out。

该样例满足测试点 $13\sim 16$ 的限制。

【数据范围与提示】

对于所有数据，保证 $1\le n\le 4\times 10^6$，$0\le k\le 10^9$，$1\le a_i\le 10^9$。

每个测试点的具体限制见下表：

测试点编号	$n\le$	$k,a_i\le$	特殊性质
$1$	$10^6$	$10^9$	所有数都相等
$2\sim 4$	$4$	$10^9$	无
$5\sim 8$	$100$	$10^9$	无
$9\sim 12$	$1000$	$10^9$	无
$13\sim 16$	$10^6$	$10^9$	无
$17\sim 20$	$4\times 10^6$	$10^9$	无

本题的读入量较大，请选择较快的读入方式，下面提供一种读入策略：

请在代码的开头加入此行：std::ios::sync_with_stdio(false);std::cin.tie(0);。（我就不加）

请注意，加入本行后 cin/cout 的效率将大幅提高，保证其能在 250 ms 内读入所有数据，但使用后你仅能使用 cin/cout 流读入数据。

思路

• 首先，我们会发现，这道题要让所有数最终相等，那么就得让它等于 $t=gcd(a_1,a_2...a_n)$ 。
• 其次，我们算出了 $t$ 是吧，然后我们把数列中等于 $t$ 的作为区间的分界点，比如：

10 3
2 2 2 1 2 2 2 1 2 2 

我们就划分成：

(2 2 2) 1 (2 2 2) 1 (2 2) 

然后我们一顿区间合并（就是看如果合并的话能然答案更小，我们就要了）

• 细节1：本人发现 vector 的 size() 函数，它是 unsigned 类型的，因此如果它等于 0 再减去 1 ，会出问题，因此，这道题我们最好用数组来做。
• 细节2：开long long。
• 细节3：每个组我们要维护 $l,r,w$ ，分别表示左端点，右端点，代价值。
• 细节4：因为有的区间的 $gcd$ 不等于 $t$ 的话，我们此时就在此区间的 $w+1$ 就行了。（为什么呢？因为此时我们就得扩展一下区间到那个边界，此时代价值就得+1）
• 细节5：我们合并完区间后，要把原区间的 $w$ 清零。

代码

//这道题就枚举分界点呗，看看下一步是否要区间合并

#include<iostream>
#include<algorithm>
#include<cstring>

#define int long long

using namespace std;

const int N = 5e6+10;

int w[N],a[N];
int n,m;
int cnt1,cnt2;
int ans;
int t;//因为最后的结果肯定是全为一个数
struct E{
    int l,r,w;
}b[N];

int gcd(int a,int b){
    return b?gcd(b,a%b):a;
}

signed main(){
    scanf("%lld%lld",&n,&m);
    
    for(int i=1;i<=n;i++){
        scanf("%lld",&w[i]);
        t=gcd(t,w[i]);        
    }
    
    for(int i=1;i<=n;i++){
        if(t==w[i]){
            a[++cnt1]=i;//记录分界点
        }
    }
    
    //(222)1(222)1(23)
    //边界特判
    if(a[1]>=2)b[++cnt2]={1,a[1]-1,a[1]-1+m};
    
    for(int i=1;i<cnt1;i++){
        if(a[i+1]-a[i]==1)continue;
        b[++cnt2]={a[i]+1,a[i+1]-1,a[i+1]-a[i]-1+m};
    }
    
    if(a[cnt1]<n)b[++cnt2]={a[cnt1]+1,n,n-a[cnt1]+m};
    
    //下面记录每一段内的gcd是不是答案，如果不是，还要加上左边（或者右边）的全局gcd的一个贡献
    
    for(int i=1;i<=cnt2;i++){
        int g=w[b[i].l];
        for(int j=b[i].l;j<=b[i].r;j++){
            g=gcd(g,w[j]);
        }
        
        if(g!=t){
            b[i].w++;
        }
    }
    
    for(int i=1;i<cnt2;i++){
        if(b[i+1].r-b[i].l+1+m<b[i+1].w+b[i].w){//如果合并更优，我们就合并
            b[i].w=0;
            b[i+1]={b[i].l,b[i+1].r,b[i+1].r-b[i].l+1+m};
        }
    }
    
    
    for(int i=1;i<=cnt2;i++){
        ans+=b[i].w;
    }
    
    printf("%lld",ans);
    
    return 0;
    
}