BZOJ 1925 地精部落 DP

最新推荐文章于 2021-01-25 14:29:59 发布

原创最新推荐文章于 2021-01-25 14:29:59 发布 · 195 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#DP #BZOJ

总结同时被 2 个专栏收录

37 篇文章

订阅专栏

8 篇文章

订阅专栏

本文探讨了波动数列的特性及其在动态规划中的应用。通过三个关键性质，揭示了如何利用DP解决这类问题，并提供了详细的转移方程及代码实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

By Bartholomew

~~其实不要小看一道地精部落，有比较大的思维量在里面！~~

小结论：

我们首先知道 3 个性质：

如果有想看证明的，请自动转到代码片下面，因为考虑有些人不想看证明

$First$ :
对于每一个数字 i 和 i+1 ，如果这两个数不是相邻的，那么交换两个数字的对应的方案数是一样的！
比如有波动序列 $3 2 4 1 5$ 和 $1 2 4 3 5$
$Second$ :
我们由 1~n 的波动数列的任意一种，将每一个 $a_{i}$ 都可以用（n+1）减去，那么得到的新的序列其实还是合法的,而且相对的山谷和山峰会改变！
比如有波动序列 $32415$ 和 $34251$
$Third$ :
波动数列具有对称性……

DP转移：

其实就是变成了 $dp[i][j] =$ $dp[i-1][j] +$ $dp[i-1][i+1-j]$
$dp[i][j]$ 表示可以用当前离散化之后的1~i 号的以 $j$ 为开头的山峰的方案个数
我们可以思考一下：
一.
我们假设有 $j$ 与 $j+1$ 并不相邻—>那么我们就是方案数为 $dp[i-1][j]$
为什么？

因为 $dp[i][j-1]$ 表示的是将j放入第一个的方案数,而且 $j-1$ 就是山峰,那么我们就知道了 $j$ 肯定不是第二个,那么这样 $j-1$ 就不是山峰了!

二.
如果 $j$ 与 $j-1$ 是相邻的,那么此时的方案数就是 dp[i+1-j]
因为我们可以考虑就是 $j$ , $j-1$ , $????...$
那么就是相当于转换成了 $1$ ~ $i-1$ 的数字 $j-1$ 为山谷的情况了! 因为 $j>j-1$ 那么j-1 的右边只能是比 $j-1$ 要大的数字,也就是它是山谷的时候了!
那么我们根据性质二可以知道就是等同于给你 $[1,j-1] U [j+1,i]$ 的数,让你求以 $j-1$ 为开头的山谷的方案数.那么其实又可以发现 $[j+1,i]$ 的数其实也就是相当于 $[j,i-1]$ ,可以自行思考一发!
那么就是求 1~i-1 以 $j-1$ 为开头的山谷的方案数. 也就是同时用 $i-1+1$ 减去之后的方案数,即 dp[i-1][i-1+1-(j-1)] = dp[i-1][i+1-j]的方案数.!

代码:

用了滚动数组!

#pragma GCC optimize(2)
#include <iostream>
#include <algorithm>
#include <queue>
#include <vector>
#include <map>
#include <math.h>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#define MAXN 4201
#define N 1000050
#define INF 0x3f3f3f3f
#define P pair<int,int>
#define ll long long
using namespace std;
int n,MOD,dp[2][MAXN];
inline ll read()
{
    ll x=0;
    char c=getchar();
    bool flag=0;
    while(c<'0'||c>'9'){if(c=='-')flag=1;   c=getchar();}
    while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+c-'0';c=getchar();}
    return flag?-x:x;
}
int main(int argc,char const* argv[])
{
    n = read(); MOD = read();
    dp[0][2] = 1;
    for(register int i=3;i<=n;++i) 
        for(register int j=1;j<=i;++j)
            dp[i&1][j] = (dp[i&1][j - 1] + dp[(i-1)&1][i-j+1]) % MOD;// dp[i-1][i-j+1] = dp[i-1][(i-1+1) - (j-1)];
    int ans = 0;
    for(int i=2;i<=n;i++) ans = (ans + dp[n&1][i])%MOD;
    printf("%d",(ans<<1) % MOD);
    return 0;
}