杆子的排列

题目背景:

10.22 NOIP模拟作业T2

分析:DP

定义状态dp[i][j][k]表示,目前枚举到第i大的数(即n - i + 1)那么显然如果这一个数放在左边,可以在左边被看到,放在右边可以在右边看到,放在中间任何位置都不能被看见,(也就是说当前枚举的数是目前最小的一个)那么转移就很显然了

dp[i][j][k] = (dp[i - 1][j - 1][k] + dp[i - 1][j][k - 1] + dp[i -1][j][k] * (i - 2))

dp[i - 1][j - 1][k] à 这一个放在最左边

dp[i - 1][j][k - 1] à 这一个放在最右边

dp[i - 1][j][k] à 这一个放在中间,(一共有i - 1个,中间位置有i - 2个)

上面的DP是O(n3)的,那么总复杂度是O(n3+ T)

然而,刚才又知道了一种也很妙妙的方法,dp[i][j]表示,当前枚举到第i大的数,然后左边的上升序列长度为j的方案数,显然,右边的上升序列长度为j的方案数也是dp[i][j],那么对于每一个n,我们枚举最高的杆子n所在的位置x,然后在剩下的n - 1个书中选择x - 1个放在左边,并且上升序列为l - 1,然后右边的n - x个数满足上升序列为r - 1那么显然方案数位dp[x - 1][l - 1] *dp[n - x][r - 1] * c(n - 1, x - 1),(c为组合数),那么我们就可以在O(n2+ Tn)的时间里面解决了(是不是可以该数据范围啊,T = 1000000, n =200或者T = 1000, n = 1000呵呵哒)

O(n3 + T)

Source:

/*created by scarlyw
*/
#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <cctype>
#include <set>
#include <map>
#include <vector>
#include <queue>const int MAXN = 200 + 10;
const int mod = 998244353;int n, l, r, t;
int dp[MAXN][MAXN][MAXN];inline void add(int &x, int t) {x += t, (x >= mod) ? (x -= mod) : (0);
}inline void pre_work() {dp[1][1][1] = 1;for (int i = 2; i < MAXN; ++i)for (int j = 1; j < MAXN; ++j)for (int k = 1; k < MAXN; ++k) {if (j + k > i + 1) break ;if (j + k < 3) continue ;add(dp[i][j][k], dp[i - 1][j - 1][k]);add(dp[i][j][k], dp[i - 1][j][k - 1]);add(dp[i][j][k], (long long)dp[i - 1][j][k] * (i - 2) % mod);}
//  for (int i = 1; i <= 5; ++i)
//      for (int j = 1; j <= 5; ++j)
//          for (int k = 1; k <= 5; ++k)
//              std::cout << i << " " << j << " " << k << " "
//                  << dp[i][j][k] << '\n';
}int main() {scanf("%d", &t), pre_work();while (t--) scanf("%d%d%d", &n, &l, &r), printf("%d\n", dp[n][l][r]);return 0;
}

O(n2 + Tn)

Source:

/*created by scarlyw
*/
#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <cctype>
#include <set>
#include <map>
#include <vector>
#include <queue>const int MAXN = 200 + 10;
const int mod = 998244353;int n, l, r, t, ans;
int dp[MAXN][MAXN], c[MAXN][MAXN];inline void add(int &x, int t) {x += t, (x >= mod) ? (x -= mod) : (0);
}inline void pre_work() {dp[0][0] = 1;for (int i = 1; i < MAXN; ++i)for (int j = 1; j <= i; ++j) {add(dp[i][j], dp[i - 1][j - 1]);if (i > j) add(dp[i][j], (long long)dp[i - 1][j] * (i - 1) % mod);}
//  for (int i = 1; i <= 4; ++i)
//      for (int j = 1; j <= 4; ++j)
//          std::cout << i << " " << j << " " << dp[i][j] << '\n';for (int i = 0; i < MAXN; ++i) c[i][i] = 1, c[i][0] = 1;for (int i = 2; i < MAXN; ++i)for (int j = 1; j < i; ++j)c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;
//  for (int i = 1; i <= 5; ++i)
//      for (int j = 1; j <= 5; ++j)
//          for (int k = 1; k <= 5; ++k)
//              std::cout << i << " " << j << " " << k << " "
//                  << dp[i][j][k] << '\n';
}int main() {scanf("%d", &t), pre_work();while (t--) {scanf("%d%d%d", &n, &l, &r), ans = 0;for (int i = 0; i <= n - 1; ++i)add(ans, (long long)dp[i][l - 1] * (long long)dp[n - 1 - i][r - 1] % mod * (long long)c[n - 1][i] % mod);printf("%d\n", ans);}return 0;
}

NOIP模拟(10.22)T2 杆子的排列相关推荐

  1. JZOJ 3885. 【长郡NOIP2014模拟10.22】搞笑的代码

    Description 在OI界存在着一位传奇选手--QQ,他总是以风格迥异的搞笑代码受世人围观 某次某道题目的输入是一个排列,他使用了以下伪代码来生成数据 while 序列长度< n do { ...

  2. JZOJ 3886. 【长郡NOIP2014模拟10.22】道路维护

    Description 最近徆多人投诉说C国的道路破损程度太大,以至亍无法通行 C国的政府徆重视这件事,但是最近财政有点紧,丌可能将所有的道路都进行维护,所以他们决定按照下述方案进行维护 将C国抽象成 ...

  3. 【NOIP2015模拟10.22】最小代价

    Description 给出一幅由n个点m条边构成的无向带权图. 其中有些点是黑点,其他点是白点. 现在每个白点都要与他距离最近的黑点通过最短路连接(如果有很多个黑点,可以选取其中任意一个),我们想要 ...

  4. 【NOIP2015模拟10.22】最大子矩阵

    Description 我们将矩阵A中位于第i行第j列的元素记作A[i,j].一个矩阵A是酷的仅当它满足下面的条件: A[1,1]+A[r,s]<=A[1,s]+A[r,1] (r,s>1 ...

  5. 【NOIP2015模拟10.22】矩形

    Description 给定一个由数字(0-9)构成的字符串s.我们可以由此定义出size(s) * size(s) 大 小的矩阵b,其中b[i][j] = s[i] * s[j]:请问在这个矩阵b中 ...

  6. 【JZOJ3885】【长郡NOIP2014模拟10.22】搞笑的代码

    ok 在OI界存在着一位传奇选手--QQ,他总是以风格迥异的搞笑代码受世人围观 某次某道题目的输入是一个排列,他使用了以下伪代码来生成数据 while 序列长度<n do { 随机生成一个整数属 ...

  7. {题解}[jzoj3885]【长郡NOIP2014模拟10.22】搞笑的代码

    传送门 Analysis 其实吧 出题人已经够友善啦 出题人已经把方程给了 我们尝试写出来 设f(i)为生成i排列的期望次数 设f(i)为生成i排列的期望次数 f(i)=[(f(i−1)+1)∗n−i ...

  8. JZOJ5920. 【NOIP2018模拟10.22】风筝

    题意: 当一阵风吹来,风筝飞上天空,为了你,而祈祷,而祝福,而感动-- oyiya 在 AK 了 IOI 之后来到了乡下,在田野中玩耍,放松身心. 他发现前面有一排小朋友在放风筝,每一个风筝有一个高度 ...

  9. NOIP模拟测试22「位运算」

    范围n-----$100000$   m $30$ 输出方案 这是一个很好的$dp$题 首先我们应该看出来一条性质只要你最后有方案达到$n$个$1$,那么你可以达到任何`一种$n$个$1$的情况 例如 ...

最新文章

  1. 项目实战-解决AES(java.security.InvalidKeyException: Illegal key size)
  2. 打印tensorflow恢复模型中所有变量与操作节点
  3. 对java面向对象的三大特征的理解_Java面向对象的三大特征是什么?
  4. SQL Server的镜像是基于物理块变化的复制 镜像Failover之后数据的预热问题
  5. 正则表达式里转义字符_五分钟搞定正则表达式,如果没搞定,再加两分钟
  6. UNP Chapter 27 - 客户-服务器程序其他设计方法
  7. 为什么物理不能理解爱情?
  8. excel文件存入mysql_解析excel文件并将数据导入到数据库中
  9. 简述python的安装过程_python3+ selenium3开发环境搭建-手把手教你安装python(详细)...
  10. airpods pro连接安卓声音小_airpods pro 麦克风声音小
  11. CSS + DIV 让页脚始终保持在页面底部
  12. 【整理】牛客网编程题前端篇(中等难度)
  13. yoga ins android.ios,to Yoga Lin:
  14. 最近做的一些文件破解的结果
  15. 【UOJ 710】魔塔 OL(贪心)(四毛子分块)
  16. 设置此计算机的用户已选择屏蔽此网站,已阻止此网站显示有安全证书错误的内容,详细教您已阻止此网站显示有安全证书错误的内容怎么解决...
  17. Oracle首席工程师四火:技术面试中,怎样的问题才是好问题?
  18. 前端白屏问题_小程序白屏问题和内存解决方法
  19. 软件不能用Surface Pro触控笔 微软这有解决办法
  20. const指针的引用

热门文章

  1. Unity3d静动态批处理Batching优化帧率、Batches(DrawCall) 初使用
  2. solidworks钣金折弯扣除表_钣金折弯人员必备知识
  3. ansible剧本功能编写
  4. ZUI 动态加载树形菜单+动态打开Tab
  5. 城市智慧泊车整体解决方案
  6. [听侯捷演讲]关于IT人士个人发展的一点思索
  7. 浏览器的同源策略和跨域请求_学习版
  8. TSC 342E Pro标签打印机设置方法
  9. Tampermonkey油猴
  10. jsplumb 点击节点或者连线高亮显示