[ AT_agc009_c]Division into Two

题目传送门

大家都知道集合吧,原谅我喜欢说废话

解法

先钦定 A > B A>B A>B,

先把数排好序得到数组 a a a,考虑先解决集合 X X X的问题,

设计状态:

明显只有一维( n < = 1 e 5 n<=1e5 n<=1e5),所以有:
f i : X 集合最后放的是第 i 个数,对于前 i 个数来说 X , Y 都合法的方案数 f_i:X集合最后放的是第i个数,对于前i个数来说X,Y都合法的方案数 fi:X集合最后放的是第i个数,对于前i个数来说X,Y都合法的方案数

发现后 ( n − i ) (n-i) (n−i)个放在 Y Y Y的数我们未保证合法,但是我们之后特判一下就好了

状态转移

设 f i f_i fi由 f j f_j fj转移而来

考虑有哪些 j j j可以贡献:

  1. i > j i>j i>j(废话)

  2. a i − a j ≥ A a_i-a_j\ge A ai−aj≥A(差值大于 A A A)

  3. max ⁡ { a k − a k − 1 } ≥ B , k ∈ ( j , i ) \max\{a_k-a_{k-1}\}\ge B,k\in(j,i) max{ak−ak−1}≥B,k∈(j,i)(满足 Y Y Y集合的要求)

则: f i = ∑ f j f_i=\sum f_j fi=∑fj

code:

复制代码
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int mod = 1e9 + 7,N=1e5+7;
int n,l,r;
int f[N],sum[N];
ll A,B;
ll a[N];
int main(){
	scanf("%d%lld%lld",&n,&A,&B);
	if(A<B) swap(A,B);
	for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
	for (int i=1;i+2<=n;++i)
        if (a[i+2]-a[i]<B) return puts("0"),0;
	a[0]=0;
	f[0]=sum[0]=1;
	for(int i=1;i<=n;i++) {
		while(a[i]-a[r+1]>=A) r++;
		if(l<=r) {
			f[i]=(1ll*f[i]+1ll*sum[r])%mod;
			if(l) f[i]=(1ll*f[i]-1ll*sum[l-1]+mod)%mod;
		}
		sum[i]=(1ll*sum[i-1]+1ll*f[i])%mod;
		if(a[i]-a[i-1]<B) l=i-1;
	}
	int ans=0;
	for(int i=n;~i;i--){	
		ans=(1ll *ans+1ll*f[i])	%mod;
		if(i<n && a[i+1]-a[i] <B) break;
	}
	printf("%d\n",ans);
}

TXL

相关推荐
小O的算法实验室几秒前
2024年IEEE TMC SCI1区TOP,面向无人机辅助 MEC 系统的轨迹规划与任务卸载的双蚁群算法,深度解析+性能实测
算法·无人机·论文复现·智能算法·智能算法改进
无才顽石1 小时前
什么是数学
算法·数理象
AlexMercer10121 小时前
【操作系统】操作系统期末考试 简答题 焚决
c语言·经验分享·笔记·操作系统
CoderCodingNo1 小时前
【GESP】C++五级真题(数论, 贪心思想考点) luogu-B4070 [GESP202412 五级] 奇妙数字
开发语言·c++·算法
百***58841 小时前
MATLAB高效算法实战技术文章大纲1
人工智能·算法·matlab
学嵌入式的六子1 小时前
如何使用VScode开发STM32【喂饭级教程】-全过程讲解
c语言·ide·vscode·stm32·单片机·嵌入式硬件
墨辰JC1 小时前
C语言可变参数讲解:stdarg.h应用
c语言·开发语言·蓝桥杯·内存·蓝桥杯嵌入式
C语言小火车1 小时前
C++右值引用与转移语义详解
c语言·开发语言
hans汉斯1 小时前
【人工智能与机器人研究】自动移液设备多轴运动控制系统设计
算法·机器学习·3d·自然语言处理·机器人·硬件架构·汉斯出版社
guygg882 小时前
经典信道估计MATLAB实现(含LSMMSE算法)
深度学习·算法·matlab