http://uoj.ac/problem/204

肯定要离散化的,先离散化出\(O(n)\)个取值区间。

设\(f(i,j)\)表示第\(i\)所学校派出的划艇数量在\(j\)区间中。

\(f(i,j)=\sum\limits_{k=0}^{i-1}\left(\sum\limits_{t=1}^{j-1}f(k,t)\right)\times Cal(k+1,i,j)\)

\(Cal(l,r,j)\)表示\([l,r)\)中的每所学校要不然不派出划艇,要不然派出数量在\(j\)区间中的划艇,第\(r\)所学校一定要派出数量在\(j\)区间中的划艇,且满足划艇数递增的方案个数。

假设\([l,r)\)中只有\(m\)所学校能满足派出数量在\(j\)区间中的划艇,设\(j\)区间的大小为\(l\),那么\(Cal(l,r,j)=\sum\limits_{i=1}^{m+1}{l\choose i}\times{m\choose i-1}={l+m\choose m+1}\)。

利用这个组合数,再记录一下dp的前缀和,时间复杂度\(O(n^3)\)。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll; const int N = 503;
const int p = 1000000007; int l[N], a[N], b[N], H[N << 1], cnt = 0, n, f[N][N << 1], ni[N]; int main() {
scanf("%d", &n);
for (int i = 1; i <= n; ++i) {
scanf("%d%d", a + i, b + i);
H[++cnt] = a[i];
H[++cnt] = ++b[i];
} stable_sort(H + 1, H + cnt + 1);
cnt = unique(H + 1, H + cnt + 1) - H; for (int i = 1; i <= n; ++i) {
a[i] = lower_bound(H + 1, H + cnt, a[i]) - H;
b[i] = lower_bound(H + 1, H + cnt, b[i]) - H;
} cnt -= 2;
for (int i = 1; i <= cnt; ++i)
l[i] = H[i + 1] - H[i]; ni[1] = 1;
for (int i = 2; i <= n; ++i)
ni[i] = 1ll * (p - p / i) * ni[p % i] % p; for (int i = 0; i <= cnt; ++i) f[0][i] = 1;
for (int i = 1; i <= n; ++i) {
for (int j = a[i], up = b[i]; j < up; ++j) {
int C = l[j], r = l[j], c = 1;
for (int k = i - 1; k >= 0; --k) {
(f[i][j] += 1ll * f[k][j - 1] * C % p) %= p;
if (a[k] <= j && j < b[k]) {
++r; ++c;
C = 1ll * C * r % p * ni[c] % p;
}
}
} for (int j = 2; j <= cnt; ++j)
(f[i][j] += f[i][j - 1]) %= p;
} int ans = 0;
for (int i = 1; i <= n; ++i) (ans += f[i][cnt]) %= p;
printf("%d\n", ans);
return 0;
}

【UOJ #204】【APIO 2016】Boat的更多相关文章

  1. 【UOJ #205】【APIO 2016】Fireworks

    http://uoj.ac/problem/205 好神的题啊. dp[i][j]表示以i为根的子树调整成长度j需要的最小代价. 首先要观察到dp值是一个下凸壳. 因为从儿子合并到父亲时要把所有儿子的 ...

  2. 【UOJ #206】【APIO 2016】Gap

    http://uoj.ac/problem/206 对于T=1,直接从两端往中间跳可以遍历所有的点. 对于T=2,先求出最小值a和最大值b,由鸽巢原理,答案一定不小于\(\frac{b-a}{N-1} ...

  3. UOJ 275. 【清华集训2016】组合数问题

    UOJ 275. [清华集训2016]组合数问题 组合数 $C_n^m $表示的是从 \(n\) 个物品中选出 \(m\) 个物品的方案数.举个例子,从$ (1,2,3)(1,2,3)$ 三个物品中选 ...

  4. UOJ #269. 【清华集训2016】如何优雅地求和

    UOJ #269. [清华集训2016]如何优雅地求和 题目链接 给定一个\(m\)次多项式\(f(x)\)的\(m+1\)个点值:\(f(0)\)到\(f(m)\). 然后求: \[ Q(f,n,x ...

  5. [UOJ#276]【清华集训2016】汽水

    [UOJ#276][清华集训2016]汽水 试题描述 牛牛来到了一个盛产汽水的国度旅行. 这个国度的地图上有 \(n\) 个城市,这些城市之间用 \(n−1\) 条道路连接,任意两个城市之间,都存在一 ...

  6. UOJ #274. 【清华集训2016】温暖会指引我们前行 [lct]

    #274. [清华集训2016]温暖会指引我们前行 题意比较巧妙 裸lct维护最大生成树 #include <iostream> #include <cstdio> #incl ...

  7. 【模拟】【数学】CSU 1803 2016 (2016湖南省第十二届大学生计算机程序设计竞赛)

    题目链接: http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1803 题目大意: 给定n,m(n,m<=109)1<=i<=n,1& ...

  8. 【UOJ】67 新年的毒瘤 &【BZOJ】1123 BLO

    [UOJ 67] 题目链接: 传送门 题解: 第一眼很懵逼……这什么鬼. 思考什么点复合条件……(o(>﹏<)o 1.树,也就是说还剩n-2条边,等价于要删去一个度数为m-n+2的点. 2 ...

  9. 【UOJ#236】[IOI2016]railroad(欧拉回路,最小生成树)

    [UOJ#236][IOI2016]railroad(欧拉回路,最小生成树) 题面 UOJ 题解 把速度看成点,给定的路段看成边,那么现在就有了若干边,然后现在要补上若干边,以及一条\([inf,\) ...

随机推荐

  1. ASP.NET读取RSS

    从网上找的一段读取RSS的代码,经测能用: /// <summary> /// 加载RSS /// </summary> /// <param name="Rs ...

  2. 51nod1110 距离之和最小 V3

    基准时间限制:1 秒 空间限制:131072 KB 分值: 40  X轴上有N个点,每个点除了包括一个位置数据X[i],还包括一个权值W[i].该点到其他点的带权距离 = 实际距离 * 权值.求X轴上 ...

  3. 【洛谷 P1772】 [ZJOI2006]物流运输(Spfa,dp)

    题目链接 \(g[i][j]\)表示不走在\(i\text{~}j\)时间段中会关闭的港口(哪怕只关\(1\)天)从\(1\)到\(m\)的最短路. \(f[i]\)表示前\(i\)天的最小花费.于是 ...

  4. 【leetcode 简单】第三十六题 最小栈

    设计一个支持 push,pop,top 操作,并能在常数时间内检索到最小元素的栈. push(x) -- 将元素 x 推入栈中. pop() -- 删除栈顶的元素. top() -- 获取栈顶元素. ...

  5. python设计模式之单例模式(一)

    前言 单例模式是创建模式中比较常见和常用的模式,在程序执行的整个生命周期只存在一个实例对象. 系列文章 python设计模式之单例模式(一) python设计模式之常用创建模式总结(二) python ...

  6. linux下的僵尸进程处理SIGCHLD信号【转】

    转自:http://www.cnblogs.com/wuchanming/p/4020463.html 什么是僵尸进程? 首先内核会释放终止进程(调用了exit系统调用)所使用的所有存储区,关闭所有打 ...

  7. puppet practice

    目标 试验环境有两台主机(VM)构成,一台是master,一台是agent,完成以下工作: 新建用户newuser; 安装 ubuntu-cloud-keyring package,更改文件/etc/ ...

  8. Spring Cloud Feign 输出日志

    还需要在application 文件中配置: #feign调用日志输出logging.level.cn.XXX=DEBUG Logger.Level下面有几种级别. BASIC : 只输出 请求URL ...

  9. springboot + swagger2 生成api文档

    直接贴代码: <dependency> <groupId>io.springfox</groupId> <artifactId>springfox-sw ...

  10. 虚拟机NAT网络设置

    1. 虚拟机设置 2. 本地网络设置 3. 本地虚拟网卡设置 4. 安装虚拟机,设置网络为NAT方式即可访问外网.