20250216 总结
20250216 总结
A - Codeforces - 1547F
简化题意：给定一个长度为 \(n\) 的数组 \(a\)，每次操作会将 \(a_i(1\le i<n)\) 变为 \(\gcd (a_i,a_{i+1})\)，将 \(a_n\) 变为 \(\gcd (a_n,a_1)\)，求最少的使得 \(a\) 中所有数相等的操作次数。
思路：设开始时 \(b\) 数组等于 \(a\) 数组 观察操作过程：
第 \(1\) 次操作：
第 \(2\) 次操作
以此类推，将 \(b\) 数组倍长，发现一个规律：
第 \(t\) 次操作后，\(a_i=\gcd (b_i,\cdots,b_{i+t})\) 。
这就好办了，ST 表维护区间 gcd，二分操作次数并 check \(a_i\) 即 \(\gcd (b_i,\cdots,b_{i+t})\) 是否相等即可，时间复杂度 \(\Omicron (n\log_2 n)\)。
细节及难点：如果没有倍长 \(b\) 数组，那么二分时 check 函数需要分讨，较为麻烦。所以题目没有刻意卡空间限制时能倍长则倍长。
B - Codeforces - 1548B
简化题意：给定一个长度为 \(n\) 的数组 \(a\)，求 \(a\) 数组中的一个子数组 \(a_l,a_{l+1},\cdots,a_{r}\)，使得存在一个 \(m\ge 2\) 满足 \(a_l\bmod m = a_{l+1}\bmod m =\cdots = a_r\bmod m\)，求其最大长度。
思路：思考一个 \(m\) 满足什么条件才能使 \(x\bmod m = y\bmod m\)。满足 \(x\bmod m = y\bmod m\) 即满足 \(x\equiv y (\bmod\ m)\)，则 \(|x-y|\equiv 0 (\bmod\ m)\)，则 \(m\) 为 \(|x-y|\) 的因数。
对于一个子数组 \(a_l,a_{l+1},\cdots,a_{r}\)，满足 \(a_l\bmod m = a_{l+1}\bmod m =\cdots = a_r\bmod m\) 即满足 \(m\) 为 \(|a_{l+1}-a_l|,|a_{l+2}-a_{l+1}|,\cdots,|a_r-a_{r-1}|\) 的因数，即 \(m\) 为 \(\gcd (a_{l+1}-a_l,a_{l+2}-a_{l+1},\cdots,a_r-a_{r-1})\) 的因数，则当 \(\gcd (a_{l+1}-a_l,a_{l+2}-a_{l+1},\cdots,a_r-a_{r-1})\ge 2\) 时 \(m\) 满足所有条件。
这也好办，ST 表维护 \(a\) 的差分数组的区间 gcd，枚举区间起始点，二分区间长度并 check 是否满足上述条件即可，时间复杂度 \(\Omicron (n\log_2 n)\)。
细节及难点：本题需分析同余式的本质与性质，其余实现不难。
C - Codeforces - 689D
简化题意：给定长度为 \(n\) 的两数组 \(a,b\)，求满足 \(\max\limits_{i=l}^{r} a_i =\min\limits_{i=l}^{r} b_i\) 的整数对 \((l,r)\) 的数量。
思路：我们对于每个左端点 \(l\)，求出第一个和最后一个满足条件的右端点 \(r_0,r_1\)，则当前左端点对答案的贡献为 \(r_1-r_0+1\)，暴力求的复杂度为 \(\Omicron (n^2)\)，不可取。
发现 \(\max\limits_{i=l}^{r} a_i\) 的值随 \(i\) 变大越来越大，\(\min\limits_{i=l}^{r} b_i\) 的值随 \(i\) 变大越来越小，则 \(\max\limits_{i=l}^{r} a_i -\min\limits_{i=l}^{r} b_i\) 单调不降，ST 表维护 \(a\) 的区间最大值和 \(b\) 的区间最小值，二分出 \(r_0,r_1\) 即可 \(\Omicron (n\log_2 n)\) 解决。
细节及难点：无细节及难点，发现单调性后维护就比较简单了。
D - Codeforces - 863E
简化题意：给定 \(n\) 个整数区间 \([l_i,r_i]\)，求满足删除该区间前后被覆盖的整点集合相同的任意一个区间。
思路：考虑一个区间满足什么条件才可删，容易发现如果把删前的 \(n\) 个区间对于每个整点的覆盖次数记录进数组 \(t\)，则满足 \(\min\limits_{i=l}^{r} t_i\ge 2\) 的区间可删。对于 \(1\le l_i\le r_i\le 10^9\) 的数据，显然暴力统计覆盖次数以及判断区间合法性都不现实。对于覆盖次数的统计，因为没有修改，可以离散化后差分；对于判断区间合法性，可以 ST 表维护 \(\min t_i\) 并 \(\Omicron (1)\) 求解 \(\min\limits_{i=l}^{r} t_i\)，总时间复杂度 \(\Omicron (n\log_2 n)\)。
细节及难点：本题细节及难点在于离散化的细节。因为差分时需对 \(l_i\) 和 \(r_i+1\) 进行修改，而判断区间合法性时需对 \([l_i,r_i]\) 这个区间进行判断，所以需要离散化 \(l_i,r_i,r_{i+1}\) 三个值，开始时 WA 错在没有离散化 \(r_{i+1}\) 导致 \(r_i+1\) 对应了错误的值。
E - Codeforces - 1175E
简化题意：给定 \(n\) 个线段 \((l,r)\) 和 \(q\) 组询问 \((x,y)\)，每次询问求区间 \((x,y)\) 需要至少几条线段才能完全覆盖。
思路：朴素的贪心是按右端点将线段排序，并选择能与当前已覆盖区间重合且右端点最远的线段，复杂度 \(\Omicron (nq)\)，需优化。
考虑使用 ST 表以倍增优化。设 \(f_{i,j}\) 表示从 \(1\sim i\) 中的节点开始选 \(2^j\) 条线段能完全覆盖到的最远右端点，则转移方程为 \(f_{i,j}\gets f_{f_{i,j-1},j-1}\)，询问时则 \(\Omicron (\log_2 n)\) 的倍增跳。开始时 \(f_{i,0}\gets\max\limits_{i=1}^{n} r_i [l_i\le i]\)，这一步可以先 \(f_{i,0}\gets\max\limits_{i=1}^{n} r_i [l_i =i]\) 然后 \(f_{i,0}\gets\max (f_{i,0},f_{i-1,0})\) 来 \(\Omicron (n)\) 预处理。总时间复杂度 \(\Omicron (n\log_2 n)\)。
细节及难点：细节及难点在于想到倍增优化和 ST 表的状态定义，实现不难。
My Code
A
#include <bits/stdc++.h>
using namespace std;
namespace rab {
int t,n,a[200010],lg[200010];
class ST {
public:
int f[200010][20];
void init () {
for (int i=1;i<=n;i++) f[i][0]=a[i];
for (int j=1;j<=lg[n];j++)
for (int i=1;i+(1<<j)-1<=n;i++)
f[i][j]=__gcd (f[i][j-1],f[i+(1<<j-1)][j-1]);
}
int query (int l,int r) {
if (l>r) return 0;
int k=lg[r-l+1];
return __gcd (f[l][k],f[r-(1<<k)+1][k]);
}
} st;
bool check (int x) {
int r1=st.query (1,x);
for (int i=2;i<=n;i++) {
int m=min (n,i+x-1)-i+1;
int p1=st.query (i,min (n,i+x-1));
int p2=st.query (1,x-m);
if (__gcd (p1,p2)!=r1) return 0;
}
return 1;
}
void solve () {
cin>> n;
for (int i=1;i<=n;i++) cin>> a[i];
st.init ();
int l=1,r=n,ans=0;
while (l<=r) {
int mid=l+r>>1;
if (check (mid)) r=(ans=mid-1);
else l=mid+1;
}
cout<< ans<< "\n";
}
int main () {
cin>> t;
for (int i=2;i<=200000;i++) lg[i]=lg[i>>1]+1;
while (t--) solve ();
return 0;
}
}
int main () {
ios::sync_with_stdio (0);
cin.tie (0);cout.tie (0);
return rab::main ();
}
B
#include <bits/stdc++.h>
using namespace std;
#define ll long long
namespace rab {
int t,n,lg[200010];
ll a[200010];
class ST {
public:
ll f[200010][20];
void init () {
for (int i=1;i<n;i++) f[i][0]=abs (a[i+1]-a[i]);
for (int j=1;j<=lg[n-1];j++)
for (int i=1;i+(1<<j)-1<n;i++)
f[i][j]=__gcd (f[i][j-1],f[i+(1<<j-1)][j-1]);
}
ll query (int l,int r) {
if (l>r) return 0;
int k=lg[r-l+1];
return __gcd (f[l][k],f[r-(1<<k)+1][k]);
}
} st;
bool check (int l,int x) {return st.query (l,l+x-1)>1;}
void solve () {
cin>> n;
for (int i=1;i<=n;i++) cin>> a[i];
st.init ();
int ans=1;
for (int i=1;i<n;i++) {
int l=1,r=n-i,rs=1;
while (l<=r) {
int mid=l+r>>1;
if (check (i,mid)) l=(rs=mid+1);
else r=mid-1;
}
ans=max (ans,rs);
}
cout<< ans<< "\n";
}
int main () {
cin>> t;
for (int i=2;i<=200000;i++) lg[i]=lg[i>>1]+1;
while (t--) solve ();
return 0;
}
}
#undef ll
int main () {
ios::sync_with_stdio (0);
cin.tie (0);cout.tie (0);
return rab::main ();
}
C
#include <bits/stdc++.h>
using namespace std;
#define ll long long
namespace rab {
int n,a[200010],b[200010],lg[200010];
ll ans=0;
class ST {
public:
int f[200010][20],g[200010][20];
void init () {
for (int i=1;i<=n;i++) f[i][0]=a[i];
for (int i=1;i<=n;i++) g[i][0]=b[i];
for (int j=1;j<=lg[n];j++) {
for (int i=1;i+(1<<j)-1<=n;i++) {
f[i][j]=max (f[i][j-1],f[i+(1<<j-1)][j-1]);
g[i][j]=min (g[i][j-1],g[i+(1<<j-1)][j-1]);
}
}
}
int query (int l,int r,bool fl) {
int k=lg[r-l+1];
if (fl) return max (f[l][k],f[r-(1<<k)+1][k]);
else return min (g[l][k],g[r-(1<<k)+1][k]);
}
} st;
int main () {
cin>> n;
for (int i=2;i<=n;i++) lg[i]=lg[i>>1]+1;
for (int i=1;i<=n;i++) cin>> a[i];
for (int i=1;i<=n;i++) cin>> b[i];
st.init ();
for (int i=1;i<=n;i++) {
if (a[i]>b[i]) continue;
int l1=i,r1=n,k1=-1;
while (l1<=r1) {
int mid=l1+r1>>1;
int q1=st.query (i,mid,1);
int q0=st.query (i,mid,0);
if (q1>=q0) {
r1=mid-1;
if (q1==q0) k1=mid;
}
else l1=mid+1;
}
if (!~k1) continue;
int l2=k1,r2=n,k2=-1;
while (l2<=r2) {
int mid=l2+r2>>1;
int q1=st.query (i,mid,1);
int q0=st.query (i,mid,0);
if (q1>q0) r2=mid-1;
else l2=(k2=mid)+1;
}
ans+=(k2-k1+1);
}
cout<< ans;
return 0;
}
}
#undef ll
int main () {
ios::sync_with_stdio (0);
cin.tie (0);cout.tie (0);
return rab::main ();
}
D
#include <bits/stdc++.h>
using namespace std;
namespace rab {
int n,k,fr,h[600010],lg[600010];
class L {
public:
int l,r;
} l[200010];
class ST {
public:
int f[600010][20];
void init () {
for (int i=1;i<=fr;i++)
f[i][0]+=f[i-1][0];
for (int j=1;j<=lg[fr];j++)
for (int i=1;i+(1<<j)-1<=fr;i++)
f[i][j]=min (f[i][j-1],f[i+(1<<j-1)][j-1]);
}
int query (int l,int r) {
int k=lg[r-l+1];
return min (f[l][k],f[r-(1<<k)+1][k]);
}
} st;
int main () {
cin>> n;
for (int i=1,x,y;i<=n;i++) {
cin>> x>> y;
l[i]={h[++k]=x,h[++k]=y};
h[++k]=y+1;
}
sort (h+1,h+k+1);
k=unique (h+1,h+k+1)-h-1;
for (int i=1;i<=n;i++) {
int r1=lower_bound (h+1,h+k+1,l[i].r+1)-h;
l[i].l=lower_bound (h+1,h+k+1,l[i].l)-h;
l[i].r=lower_bound (h+1,h+k+1,l[i].r)-h;
fr=max (fr,l[i].r);
st.f[l[i].l][0]++;
st.f[r1][0]--;
}
for (int i=2;i<=fr;i++) lg[i]=lg[i>>1]+1;
st.init ();
for (int i=1;i<=n;i++) {
if (st.query (l[i].l,l[i].r)>1) {
cout<< i;
return 0;
}
}
cout<< -1;
return 0;
}
}
int main () {
ios::sync_with_stdio (0);
cin.tie (0);cout.tie (0);
return rab::main ();
}
E
#include <bits/stdc++.h>
using namespace std;
namespace rab {
int n,m,fr,lg[1000010];
class ST {
public:
int f[1000010][30];
void init () {
for (int i=1;i<=fr;i++)
f[i][0]=max (f[i][0],f[i-1][0]);
for (int j=1;j<=lg[fr];j++)
for (int i=0;i<=fr;i++)
f[i][j]=f[f[i][j-1]][j-1];
}
int query (int x,int y) {
int p=x,ans=0;
for (int i=lg[fr];~i;i--) {
if (f[p][i]<y) {
p=f[p][i];
ans+=(1<<i);
}
}
if (f[p][0]>=y) return ans+1;
else return -1;
}
} st;
int main () {
cin>> n>> m;
for (int i=1,l,r;i<=n;i++) {
cin>> l>> r;
fr=max (fr,r);
st.f[l][0]=max (st.f[l][0],r);
}
for (int i=2;i<=fr;i++) lg[i]=lg[i>>1]+1;
st.init ();
while (m--) {
int x,y;
cin>> x>> y;
cout<< st.query (x,y)<< "\n";
}
return 0;
}
}
int main () {
ios::sync_with_stdio (0);
cin.tie (0);cout.tie (0);
return rab::main ();
}
浙公网安备 33010602011771号