题目
https://atcoder.jp/contests/nomura2020/tasks/nomura2020_f?lang=en
思路
考虑给你一个序列如何判断合不合法。
一个很显然的结论(1):如果存在 i < j i<j i<j且 a i > a j a_i>a_j ai>aj且 a i a_i ai 和 a j a_j aj 不相同的位数>1,那么这个序列就不合法。
另外一个结论(2)是:捣乱的那个人只会从高位到低位删。
证明:
首先我们只需要存在一对逆序对不满足结论(1)这个序列就不合法,所以我们只用考虑一对。
对于一对,我们可以把相同的缩成一段,这些段删了是没有意义的(删不删都行)。考虑最高位不同的地方:如果是1、0,那这一对已经是逆序对了,下面的就没必要再删;如果是0、1那这一行必须删,直到找到一个1、0.
考虑如何求解:
用DP统计,设
d
p
i
,
j
dp_{i,j}
dpi,j表示
n
=
i
,
m
=
j
n=i,m=j
n=i,m=j时的答案,第一种情况就枚举0,1分界点的位置,第二种情况就枚举一个长度为
k
<
j
k<j
k<j的序列,选择其中的一个数复制
j
−
k
j-k
j−k次,这
j
−
k
+
1
j-k+1
j−k+1个数满足第
i
i
i位,它们中最左边的数为1,它们的左边的数为0,它们中最右边的数为0,它们右边的数为1,其他的
j
−
k
−
1
j-k-1
j−k−1个数任意
d
p
i
,
j
=
d
p
i
−
1
,
j
(
j
+
1
)
+
∑
k
<
j
d
p
i
−
1
,
k
k
⋅
2
j
−
k
−
1
dp_{i,j}=dp_{i-1,j}(j+1)+\sum_{k<j} dp_{i-1,k}k \cdot 2^{j-k-1}
dpi,j=dpi−1,j(j+1)+k<j∑dpi−1,kk⋅2j−k−1
注意到可以用前缀和优化,时间复杂度 O ( n ) O(n) O(n)
代码
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=5077,mod=1e9+7;
int n,m;
ll f[N][N];
int main()
{
scanf("%d%d",&n,&m);
for(int i=1; i<=m; i++) f[0][i]=1;
for(int i=1; i<=n; i++)
{
ll sum=0;
for(int j=1; j<=m; j++)
{
f[i][j]=(sum+(j+1)*f[i-1][j]%mod)%mod;
sum=(sum*2%mod+f[i-1][j]*j%mod)%mod;
}
}
printf("%lld",f[n][m]);
}