洛谷传送门
BZOJ传送门
题目描述
传说很久以前,大地上居住着一种神秘的生物:地精。
地精喜欢住在连绵不绝的山脉中。具体地说,一座长度为 N N 的山脉可分为从左到右的 N N 段,每段有一个独一无二的高度,其中 Hi H i 是 1 1 到之间的正整数。
如果一段山脉比所有与它相邻的山脉都高,则这段山脉是一个山峰。位于边缘的山脉只有一段相邻的山脉,其他都有两段(即左边和右边)。
类似地,如果一段山脉比所有它相邻的山脉都低,则这段山脉是一个山谷。
地精们有一个共同的爱好——饮酒,酒馆可以设立在山谷之中。地精的酒馆不论白天黑夜总是人声鼎沸,地精美酒的香味可以飘到方圆数里的地方。
地精还是一种非常警觉的生物,他们在每座山峰上都可以设立瞭望台,并轮流担当瞭望工作,以确保在第一时间得知外敌的入侵。
地精们希望这N段山脉每段都可以修建瞭望台或酒馆的其中之一,只有满足这个条件的整座山脉才可能有地精居住。
现在你希望知道,长度为 N N 的可能有地精居住的山脉有多少种。两座山脉和 B B 不同当且仅当存在一个,使得 Ai≠Bi A i ≠ B i 。由于这个数目可能很大,你只对它除以 P P 的余数感兴趣。
输入输出格式
输入格式:
输入文件goblin.in仅含一行,两个正整数。
输出格式:
输出文件goblin.out仅含一行,一个非负整数,表示你所求的答案对 P P 取余之后的结果。
输入输出样例
输入样例#1:
4 7
输出样例#1:
3
说明
【数据规模和约定】
对于20%的数据,满足;
对于40%的数据,满足 N≤18 N ≤ 18 ;
对于70%的数据,满足 N≤550 N ≤ 550 ;
对于100%的数据,满足 3≤N≤4200,P≤1e9 3 ≤ N ≤ 4200 , P ≤ 1 e 9 。
解题分析
又是一道 dp d p 神题… 不会做只好膜了一波其他大佬的博客。
我们设 dp[i][j] d p [ i ] [ j ] 表示前 i i 个数的排列, 满足为第一个且为山峰的方案数。
这里有个很显然的结论:如果 i,i+1 i , i + 1 不相邻, 那么交换这两个数仍然是一个合法的排列(因为相对于它们其他数的大小是没有变化的), 所以 dp[i][j]+=dp[i][j−1] d p [ i ] [ j ] + = d p [ i ] [ j − 1 ] 。
考虑 j,j−1 j , j − 1 相邻的情况, 那么 j−1 j − 1 应该在山谷的位置。 由于一个我们用 1+n 1 + n 去减一个排列的所有数也能得到一个合法的排列, 且山峰山谷的位置相反, 而剩下的 i−1 i − 1 个数中 j−1 j − 1 还是排第 j−1 j − 1 位, 所以 j−1 j − 1 排第一个, 且为山谷的方案数相当于 (i−1)−(j−1)+1 ( i − 1 ) − ( j − 1 ) + 1 排第一个, 且为山峰的方案数。
综上, dp[i][j]=dp[i][j−1]+dp[i−1][i−j+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);
}