题目
题目背景
在
o
e
i
s
\tt oeis
oeis 中已经收录了这个数列……
题目描述
世人不患寡而患不均。而
H
a
n
d
I
n
D
e
v
i
l
\sf HandInDevil
HandInDevil 只愿天下大乱!
请你找到一个长度为 n n n 的排列,使得相邻两个数的差值的绝对值都不小于 2 2 2 。
数据范围与提示
n
≤
1
0
3
n\le 10^3
n≤103,输出取模给定的数
m
m
m 。(估计是为了防止打表。)
思路
限制只发生在相邻的两个数之间?容易想到连续段 d p \tt dp dp 。用 f ( i , j , 0 / 1 / 2 , 0 / 1 ) f(i,j,0/1/2,0/1) f(i,j,0/1/2,0/1) 表示,已经填了前 i i i 个数字、形成了 j j j 个连续段、 i i i 是几个端点( 0 0 0 在段中间, 1 1 1 是段端点, 2 2 2 是单独成段)、 i i i 所在的区间是否是最靠左或最靠右的,存储方案数。这样做明显比较麻烦嘛……
另一个想法是,模仿这道题,直接依次插入。其实是完全可以的,只要记录形如 ⟨ i , i + 1 ⟩ \langle i,i+1\rangle ⟨i,i+1⟩ 的 “危险缝隙” 的数量,然后用数字插入来隔开它们就行。复杂度 O ( n 2 ) \mathcal O(n^2) O(n2),具体看代码吧。
代码
#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;
}
后记
o
e
i
s
\tt oeis
oeis 如是说:初值
F
(
0
)
=
F
(
1
)
=
1
,
F
(
2
)
=
F
(
3
)
=
0
F(0)=F(1)=1,\;F(2)=F(3)=0
F(0)=F(1)=1,F(2)=F(3)=0,当
n
≥
4
n\ge 4
n≥4 时有
F
(
n
)
=
(
n
+
1
)
F
(
n
−
1
)
−
(
n
−
2
)
F
(
n
−
2
)
−
(
n
−
5
)
F
(
n
−
3
)
+
(
n
−
3
)
F
(
n
−
4
)
F(n)=\\ (n+1)F(n-1)\\-(n-2)F(n-2)\\-(n-5)F(n-3)\\+(n-3)F(n-4)
F(n)=(n+1)F(n−1)−(n−2)F(n−2)−(n−5)F(n−3)+(n−3)F(n−4)
原因未知。但是
S
t
a
t
e
:
A
p
p
r
o
v
e
d
.
\rm State:Approved.
State:Approved.