[ACNOI2021]OEIS yyds

最新推荐文章于 2023-09-24 19:53:51 发布

OneInDark

最新推荐文章于 2023-09-24 19:53:51 发布

阅读量236

点赞数

分类专栏： C++ 动态规划

本文链接：https://blog.csdn.net/qq_42101694/article/details/115488489

版权

OEIS序列动态规划排列问题差值绝对值数学算法

关键词由CSDN通过智能技术生成

C++ 同时被 2 个专栏收录

547 篇文章 3 订阅

订阅专栏

动态规划

220 篇文章 2 订阅

订阅专栏

题目

题目背景
在 $\tt oeis$ 中已经收录了这个数列……

题目描述
世人不患寡而患不均。而 $\sf HandInDevil$ 只愿天下大乱！

请你找到一个长度为 $n$ 的排列，使得相邻两个数的差值的绝对值都不小于 $2$ 。

数据范围与提示
$n\le 10^3$ ，输出取模给定的数 $m$ 。（估计是为了防止打表。）

思路

限制只发生在相邻的两个数之间？容易想到连续段 $\tt dp$ 。用 $f (i, j, 0 / 1 / 2, 0 / 1)$ 表示，已经填了前 $i$ 个数字、形成了 $j$ 个连续段、 $i$ 是几个端点（ $0$ 在段中间， $1$ 是段端点， $2$ 是单独成段）、 $i$ 所在的区间是否是最靠左或最靠右的，存储方案数。这样做明显比较麻烦嘛……

另一个想法是，模仿这道题，直接依次插入。其实是完全可以的，只要记录形如 $\langle i,i+1\rangle$ 的 “危险缝隙” 的数量，然后用数字插入来隔开它们就行。复杂度 $\mathcal O(n^2)$ ，具体看代码吧。

代码

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}

const int MaxN = 1005;
int dp[2][MaxN][2], M;
inline void add(int&x,const int&y){
	((x += y) >= M) ? (x -= M) : 0;
}

int main(){
	for(int n; ~scanf("%d",&n); ){
		M = readint(); // mod
		memset(dp[1],0,MaxN*2<<2);
		dp[1][0][0] = 1; // only 1 itself
		for(int i=2,fr=1; i<=n; ++i,fr^=1){
			memset(dp[i&1],0,MaxN*2<<2);
			rep(j,0,i-1) rep(d,0,1){
				/* case 1: put in a dangerous gap without creating new */ ;
				if(j > d) add(dp[i&1][j-1][0],dp[fr][j][d]*(j-d)%M);
				/* case 2: put in a dangerous gap meanwhile creating new */ ;
				if(j && d) add(dp[i&1][j][1],dp[fr][j][d]); // *d
				/* case 3: put in a safe gap without creating new */ ;
				/* i gap in total, j is dangerous, 2-d will create */ ;
				add(dp[i&1][j][0],dp[fr][j][d]*(i-j-2+d)%M);
				/* case 4: put in a safe gap meanwhile creating new */ ;
				add(dp[i&1][j+1][1],dp[fr][j][d]*(2-d)%M);
			}
		}
		printf("%d\n",dp[n&1][0][0]); // no dangerous gap
	}
	return 0;
}