[Luogu P2467] [BZOJ 1925] [SDOI2010]地精部落

洛谷传送门
BZOJ传送门

题目描述

传说很久以前,大地上居住着一种神秘的生物:地精。

地精喜欢住在连绵不绝的山脉中。具体地说,一座长度为 N N 的山脉H可分为从左到右的 N N 段,每段有一个独一无二的高度Hi,其中 Hi H i 1 1 N之间的正整数。

如果一段山脉比所有与它相邻的山脉都高,则这段山脉是一个山峰。位于边缘的山脉只有一段相邻的山脉,其他都有两段(即左边和右边)。

类似地,如果一段山脉比所有它相邻的山脉都低,则这段山脉是一个山谷。

地精们有一个共同的爱好——饮酒,酒馆可以设立在山谷之中。地精的酒馆不论白天黑夜总是人声鼎沸,地精美酒的香味可以飘到方圆数里的地方。

地精还是一种非常警觉的生物,他们在每座山峰上都可以设立瞭望台,并轮流担当瞭望工作,以确保在第一时间得知外敌的入侵。

地精们希望这N段山脉每段都可以修建瞭望台或酒馆的其中之一,只有满足这个条件的整座山脉才可能有地精居住。

现在你希望知道,长度为 N N 的可能有地精居住的山脉有多少种。两座山脉A B B 不同当且仅当存在一个i,使得 AiBi A i ≠ B i 。由于这个数目可能很大,你只对它除以 P P 的余数感兴趣。

输入输出格式

输入格式:

输入文件goblin.in仅含一行,两个正整数N,P

输出格式:

输出文件goblin.out仅含一行,一个非负整数,表示你所求的答案对 P P 取余之后的结果。

输入输出样例

输入样例#1:

4 7

输出样例#1:

3

说明

img

【数据规模和约定】

对于20%的数据,满足N10

对于40%的数据,满足 N18 N ≤ 18

对于70%的数据,满足 N550 N ≤ 550

对于100%的数据,满足 3N4200P1e9 3 ≤ N ≤ 4200 , P ≤ 1 e 9

解题分析

又是一道 dp d p 神题… 不会做只好膜了一波其他大佬的博客。

我们设 dp[i][j] d p [ i ] [ j ] 表示前 i i 个数的排列, 满足j为第一个且为山峰的方案数。

这里有个很显然的结论:如果 i,i+1 i , i + 1 不相邻, 那么交换这两个数仍然是一个合法的排列(因为相对于它们其他数的大小是没有变化的), 所以 dp[i][j]+=dp[i][j1] d p [ i ] [ j ] + = d p [ i ] [ j − 1 ]

考虑 j,j1 j , j − 1 相邻的情况, 那么 j1 j − 1 应该在山谷的位置。 由于一个我们用 1+n 1 + n 去减一个排列的所有数也能得到一个合法的排列, 且山峰山谷的位置相反, 而剩下的 i1 i − 1 个数中 j1 j − 1 还是排第 j1 j − 1 位, 所以 j1 j − 1 排第一个, 且为山谷的方案数相当于 (i1)(j1)+1 ( i − 1 ) − ( j − 1 ) + 1 排第一个, 且为山峰的方案数。

综上, dp[i][j]=dp[i][j1]+dp[i1][ij+1] d p [ i ] [ j ] = d p [ i ] [ j − 1 ] + d p [ i − 1 ] [ i − j + 1 ]

注意特判 n=1 n = 1 的情况。

代码如下:

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <algorithm>
#define R register
#define IN inline
#define W while
#define gc getchar()
#define ll long long
#define MX 5005
template <class T>
IN void in(T &x)
{
    x = 0; R char c = gc;
    for (; !isdigit(c); c = gc);
    for (;  isdigit(c); c = gc)
    x = (x << 1) + (x << 3) + c - 48;
}
ll dp[2][MX];
ll mod, ans;
int n, cur;
int main(void)
{
    in(n), in(mod);
    if(n == 1) printf("%lld", 1ll * n % mod), exit(0);
    dp[cur][2] = 1;
    for (R int i = 3; i <= n; ++i)
    {
        cur ^= 1;
        for (R int j = 2; j <= i; ++j)
        dp[cur][j] = (dp[cur][j - 1] + dp[cur ^ 1][i - j + 1]) % mod;
    }
    for (R int i = 2; i <= n; ++i) ans = (ans + dp[cur][i]) % mod;
    printf("%lld", (ans << 1) % mod);
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值