Written with StackEdit.

Description

传说很久以前,大地上居住着一种神秘的生物:地精。 地精喜欢住在连绵不绝的山脉中。具体地说,一座长度为 \(N\) 的山脉 \(H\)可分 为从左到右的\(N\) 段,每段有一个独一无二的高度 \(H_i\),其中\(H_i\)是\(1\)到\(N\) 之间的正 整数。 如果一段山脉比所有与它相邻的山脉都高,则这段山脉是一个山峰。位于边 缘的山脉只有一段相邻的山脉,其他都有两段(即左边和右边)。 类似地,如果一段山脉比所有它相邻的山脉都低,则这段山脉是一个山谷。 地精们有一个共同的爱好——饮酒,酒馆可以设立在山谷之中。地精的酒馆 不论白天黑夜总是人声鼎沸,地精美酒的香味可以飘到方圆数里的地方。 地精还是一种非常警觉的生物,他们在每座山峰上都可以设立瞭望台,并轮 流担当瞭望工作,以确保在第一时间得知外敌的入侵。 地精们希望这\(N\) 段山脉每段都可以修建瞭望台或酒馆的其中之一,只有满足 这个条件的整座山脉才可能有地精居住。 现在你希望知道,长度为\(N\) 的可能有地精居住的山脉有多少种。两座山脉\(A\)和\(B\)不同当且仅当存在一个 \(i\),使得 \(A_i≠B_i\)。由于这个数目可能很大,你只对它除以\(P\)的余数感兴趣。

Input

仅含一行,两个正整数 \(N\), \(P\)。

Output

仅含一行,一个非负整数,表示你所求的答案对\(P\)取余之后的结果。

Sample Input

4 7

Sample Output

3

HINT



对于 \(20\%\)的数据,满足 \(N≤10\);

对于 \(40\%\)的数据,满足 \(N≤18\);

对于 \(70\%\)的数据,满足 \(N≤550\);

对于 \(100\%\)的数据,满足 \(3≤N≤4200,P≤10^9.\)

Solution

  • 题面废话很多,就是求长度为\(N\)的波动序列的个数.
  • 注意到对于任意一个波动序列\(a\),将其中的每个\(a_i\)变为\(N+1-a_i\),仍是一个波动序列,且由先升后降变为先降后升,或由先降后升变为先升降.
  • eg.\(1,3,2,4\)->\(4,2,3,1\).
  • 那么我们只需要求出先降后升的序列个数,再乘\(2\)即为答案.
  • 考虑令\(f[i]\)表示\(i\)个不同的数组成的先降后升的序列个数.(不同即可,与具体大小无关)
  • 计算\(f[i]\)时,考虑枚举最大的数所在的位置\(j\),显然是一个奇数.
  • 那么它前面还有\(j-1\)个数,后面还有\(i-j\)个数.
  • 我们从剩余的\(i-1\)个数中选出\(j-1\)个数放在前面,剩余的放在后面.
  • 有转移方程\(f[i]=\sum_{j=2k+1,k\in N}^{i} C_{i-1}^{j-1}*f[j-1]*f[i-j]\).
  • 时间复杂度为\(O(n^2)\).
  • 由于\(P\)是读入的,组合数的计算使用递推.但直接开\(n^2\)的数组会\(MLE\),需要滚动数组计算(第一次见到组合数还要滚动数组的).
#include<bits/stdc++.h>
using namespace std;
typedef long long LoveLive;
inline int read()
{
int out=0,fh=1;
char jp=getchar();
while ((jp>'9'||jp<'0')&&jp!='-')
jp=getchar();
if (jp=='-')
{
fh=-1;
jp=getchar();
}
while (jp>='0'&&jp<='9')
{
out=out*10+jp-'0';
jp=getchar();
}
return out*fh;
}
int n,P;
inline int add(int a,int b)
{
return (a + b) % P;
}
inline int mul(int a,int b)
{
return 1LL * a * b % P;
}
int fpow(int a,int b)
{
int res=1;
while(b)
{
if(b&1)
res=mul(res,a);
b>>=1;
a=mul(a,a);
}
return res;
}
const int MAXN=4396;
int f[MAXN];
int C[2][MAXN];
int main()
{
n=read(),P=read();
f[0]=1;
int id=1;
for(int i=1;i<=n;++i)
{
C[id^1][0]=1;
C[id][0]=1;
for(int j=1;j<=i;++j)
{
C[id][j]=add(C[id^1][j],C[id^1][j-1]);
if(j&1)
f[i]=add(f[i],mul(C[id^1][j-1],mul(f[j-1],f[i-j])));
}
id^=1;
}
int ans=mul(f[n],2);
printf("%d\n",ans);
return 0;
}

bzoj 1925 地精部落的更多相关文章

  1. [sdoi 2010][bzoj 1925]地精部落(神仙dp)

    Description 传说很久以前,大地上居住着一种神秘的生物:地精. 地精喜欢住在连绵不绝的山脉中.具体地说,一座长度为 N 的山脉 H可分 为从左到右的 N 段,每段有一个独一无二的高度 Hi, ...

  2. [SDOI2010][BZOJ 1925]地精部落

    Description 传说很久以前,大地上居住着一种神秘的生物:地精. 地精喜欢住在连绵不绝的山脉中.具体地说,一座长度为 N 的山脉 H可分 为从左到右的 N 段,每段有一个独一无二的高度 Hi, ...

  3. BZOJ 1925 地精部落(DP)

    一道很经典的DP题. 题意:求n排列中波动排列的种数. 不妨考虑DP,令dp1[i][j],表示1-j的排列中,第一项为i之后递增的波动排列种数.dp2[i][j]表示1-j的排列中,第一项为i之后递 ...

  4. BZOJ 1925地精部落题解

    题目链接 一道神仙题,有很多思考的方式,这里选择最好理解的一种来讲 我们将序列分为两种,一种开头递增,一种开头递减,显然这两种序列的数目是一样的 现在我们只用考虑开头递增的情况 f[i][j]表示前i ...

  5. BZOJ 1925: [Sdoi2010]地精部落( dp )

    dp(i,j)表示1~i的排列中, 以1~j为开头且开头是下降的合法方案数 这种数列具有对称性, 即对于一个满足题意且开头是上升的n的排列{an}, 令bn = n-an+1, 那么{bn}就是一个满 ...

  6. 1925: [Sdoi2010]地精部落

    1925: [Sdoi2010]地精部落 Time Limit: 10 Sec Memory Limit: 64 MB Submit: 1929 Solved: 1227 [Submit][Statu ...

  7. BZOJ-1925 地精部落 烧脑DP+滚动数组

    1925: [Sdoi2010]地精部落 Time Limit: 10 Sec Memory Limit: 64 MB Submit: 1053 Solved: 633 [Submit][Status ...

  8. 【BZOJ1925】[SDOI2010]地精部落(动态规划)

    [BZOJ1925][SDOI2010]地精部落(动态规划) 题面 BZOJ 洛谷 题解 一道性质\(dp\)题.(所以当然是照搬学长PPT了啊 先来罗列性质,我们称题目所求的序列为抖动序列: 一个抖 ...

  9. 【bzoj1925】地精部落[SDOI2010](dp)

    题目传送门:1925: [Sdoi2010]地精部落 这道题,,,首先可以一眼看出他是要我们求由1~n的排列组成,并且抖来抖去的序列的方案数.然后再看一眼数据范围,,,似乎是O(n^2)的dp?然后各 ...

随机推荐

  1. 【Java】仿真qq尝试:聊天界面 && 响应用户输入

    需求分析: 逐步完善一个“qq仿真”程序. 参考: 1.文本框与文本区:http://www.weixueyuan.net/view/6062.html 2.java布局:http://www.cnb ...

  2. Forward Proxy vs Reverse Proxy

    Overview We've talked about reverse proxy servers and how they can really be good at protecting the ...

  3. shell 中的 eval 及 crontab 命令

    eval eval会对后面的命令进行两遍扫描,如果第一遍扫描后,命令是个普通命令,则执行此命令:如果命令中含有变量的间接引用,则保证间接引用的语义.也就是说,eval命令将会首先扫描命令行进行所有的置 ...

  4. Ctrl+Z 暂停程序及重启程序【转】

    本文转自:https://blog.csdn.net/duyiwuer2009/article/details/43191799 Ctrl+Z - 暂停进程并放入后台 jobs - 显示当前暂停的进程 ...

  5. Monkey测试命令【学习笔记】

    monkey --ignore-crashes --ignore-timeouts --ignore-security-exceptions --ignore-native-crashes -v 30 ...

  6. 一个线程知识点, 一个MongoDB的知识点

    //WINForm窗体中切换前后台线程执行任务: protected void RunOnUI(Action action) { Invoke(action); } protected void Ru ...

  7. Mysql导出数据结构 or 数据

    如果我们单单只想导出mysql数据表结构,通过navcat工具还不行,这时我们可以用mysqldump工具 在mysql server的安装目录:C:\Program Files\MySQL\MySQ ...

  8. shell script语法高亮和自动缩进的配置

    编辑/etc/profile文件,在文件末尾加一下内容: export TERM=xterm-color 接着让其变为全局变量 source /etc/profile 编辑/etc/vimrc文件,在 ...

  9. ThinkPHP关于模板的一些嵌套、IF判断使用

    > ##### 前言,现在有一组数据(涉及到3个数据表,order订单表,order_process办理流程表,process_status流程描述表),根据当前订单,展示相应信息 1.办理流程 ...

  10. NO.4 Android开发中常用框架及工具

    android-pulltorefresh 一个强大的拉动刷新开源项目,支持各种控件下拉刷新ListView.ViewPager.WevView.ExpandableListView.GridView ...