Solution -「APIO 2016」「洛谷 P3643」划艇
\(\mathcal{Description}\)
给定 \(n\) 个区间 \([a_i,b_i)\)(注意原题是闭区间,这里只为方便后文描述),求 \(\{c_n\}\) 的个数,使得:
- \(\forall i~~~~c_i=0\lor c_i\in[a_i,b_i)\)。
- \(\forall i<j~~~~c_i\not=0\land c_j\not=0\Rightarrow c_i<c_j\)。
对 \(10^9+7\) 取模。
\(n\le500\),\(1\le a_i\le b_i\le10^9\)。
\(\mathcal{Solution}\)
一个很 naive 的 DP 想法,\(f(i,j)\) 表示考虑前 \(i\) 个位置,\(c_i=j~(j\not=0)\) 时的方案数。问题在于第二维开销过大,考虑离散化所有端点坐标。
先来一个引理,取值在 \([a,b)\),长度为 \(n\) 的上升整数序列的个数为 \(\binom{b-a}{n}\),显然选 \(n\) 个数就可以了。
再来一个引理,取值在 \([a,b)\cup\{0\}\),长度为 \(n\),非零的位置是上升整数序列的序列个数为 \(\binom{b-a+n}{n}\),证明也很显然,有几个 \(0\) 可以选,虽然不同的 \(0\) 可以任意排列,但看上去都是一样的。所以钦定 \(0\) 的大小关系后就等价于令区间为 \([a-n,b]\),长度为 \(n\) 时的上一引理。
接着刚才的思路,离散化时,排过序的端点们把坐标轴分为若左闭右开的区间,从左开始第 \(t\) 个区间称作第 \(t\) 段。令 \(f(i,j)\) 表示考虑前 \(i\) 个位置,\(c_i\) 属于前 \(j\) 段时的方案数。设 \([j,j+1)\) 实际映射 \([a,b)\),枚举 \(k<i\),转移:
\]
组合数运用了引理二。注意钦定 \(c_i\) 不为 \(0\),所以上面 \(-1\)。
扫 \(f\) 就结束了,交换枚举顺序,第二维还可以滚掉。复杂度 \(\mathcal O(n^3)\)。
\(\mathcal{Code}\)
/* Clearink */
#include <cstdio>
#include <algorithm>
const int MAXN = 500, MOD = 1e9 + 7;
int n, a[MAXN + 5], b[MAXN + 5], tmp[MAXN * 2 + 5], inv[MAXN + 5];
int f[MAXN + 5], comb[MAXN + 5];
inline int mul ( long long a, const int b ) { return a * b % MOD; }
inline int& addeq ( int& a, const int b ) { return ( a += b ) < MOD ? a : a -= MOD; }
int main () {
scanf ( "%d", &n ), inv[1] = 1;
for ( int i = 1; i <= n; ++ i ) {
if ( i > 1 ) inv[i] = mul ( MOD - MOD / i, inv[MOD % i] );
scanf ( "%d %d", &a[i], &b[i] ), ++ b[i];
tmp[2 * i - 1] = a[i], tmp[i << 1] = b[i];
}
std::sort ( tmp + 1, tmp + ( n << 1 | 1 ) );
int m = std::unique ( tmp + 1, tmp + ( n << 1 | 1 ) ) - tmp - 1;
for ( int i = 1; i <= n; ++ i ) {
a[i] = std::lower_bound ( tmp + 1, tmp + m + 1, a[i] ) - tmp;
b[i] = std::lower_bound ( tmp + 1, tmp + m + 1, b[i] ) - tmp;
}
f[0] = 1;
for ( int j = 1, len; j < m; ++ j ) {
len = tmp[j + 1] - tmp[j], comb[0] = 1;
for ( int i = 1; i <= n; ++ i ) {
comb[i] = mul ( mul ( comb[i - 1], len + i - 1 ), inv[i] );
}
for ( int i = n; i; -- i ) {
if ( j < a[i] || b[i] <= j ) continue;
for ( int k = i - 1, c = 1; ~k; -- k ) {
addeq ( f[i], mul ( comb[c], f[k] ) );
if ( a[k] <= j && j < b[k] ) ++ c;
}
}
}
int ans = 0;
for ( int i = 1; i <= n; ++ i ) addeq ( ans, f[i] );
printf ( "%d\n", ans );
return 0;
}
\(\mathcal{Details}\)
目前洛谷最优解,兔的代码吸口氧真的快到飞起 www。
Solution -「APIO 2016」「洛谷 P3643」划艇的更多相关文章
- [洛谷P3643] [APIO2016]划艇
洛谷题目链接:[APIO2016]划艇 题目描述 在首尔城中,汉江横贯东西.在汉江的北岸,从西向东星星点点地分布着 \(N\) 个划艇学校,编号依次为 \(1\) 到 \(N\).每个学校都拥有若干艘 ...
- 洛谷P3643 [APIO2016]划艇(组合数学)
题面 传送门 题解 首先区间个数很少,我们考虑把所有区间离散化(这里要把所有的右端点变为\(B_i+1\)代表的开区间) 设\(f_{i,j}\)表示考虑到第\(i\)个学校且第\(i\)个学校必选, ...
- 洛谷 P3643 - [APIO2016]划艇(dp)
题面传送门 一道难度中等的 \(dp\)(虽然我没有想出来/kk). 首先一眼 \(dp_{i,j}\) 表示考虑到第 \(i\) 个学校,第 \(i\) 个学校派出了 \(j\) 个划艇的方案数,转 ...
- 「区间DP」「洛谷P1043」数字游戏
「洛谷P1043」数字游戏 日后再写 代码 /*#!/bin/sh dir=$GEDIT_CURRENT_DOCUMENT_DIR name=$GEDIT_CURRENT_DOCUMENT_NAME ...
- Solution -「JSOI 2019」「洛谷 P5334」节日庆典
\(\mathscr{Description}\) Link. 给定字符串 \(S\),求 \(S\) 的每个前缀的最小表示法起始下标(若有多个,取最小的). \(|S|\le3\time ...
- Solution -「洛谷 P4372」Out of Sorts P
\(\mathcal{Description}\) OurOJ & 洛谷 P4372(几乎一致) 设计一个排序算法,设现在对 \(\{a_n\}\) 中 \([l,r]\) 内的元素排 ...
- Solution -「POI 2010」「洛谷 P3511」MOS-Bridges
\(\mathcal{Description}\) Link.(洛谷上这翻译真的一言难尽呐. 给定一个 \(n\) 个点 \(m\) 条边的无向图,一条边 \((u,v,a,b)\) 表示从 ...
- 「洛谷4197」「BZOJ3545」peak【线段树合并】
题目链接 [洛谷] [BZOJ]没有权限号嘤嘤嘤.题号:3545 题解 窝不会克鲁斯卡尔重构树怎么办??? 可以离线乱搞. 我们将所有的操作全都存下来. 为了解决小于等于\(x\)的操作,那么我们按照 ...
- 「洛谷3338」「ZJOI2014」力【FFT】
题目链接 [BZOJ] [洛谷] 题解 首先我们需要对这个式子进行化简,否则对着这么大一坨东西只能暴力... \[F_i=\sum_{j<i} \frac{q_iq_j}{(i-j)^2}-\s ...
随机推荐
- Google Java 风格指南(Google Java Style Guide)
官方地址 google.github.io 本文档作为 Google 的 Java 编程语言源代码编码标准的完整定义.当且仅当它遵守此处的规则时,Java 源文件才被描述为 Google 风格. 前言 ...
- VictoriaMerics学习笔记(2):核心组件
核心组件 1. 单机版 victoria-metrics-prod 单一二进制文件 读写都在一个节点上 作者推荐单机版 特性 merge方式配置 通过HTTP协议提供服务 内存限制(防止OOM) 使用 ...
- 【解决了一个小问题】golang的go.mod中出现版本错误
代码中的这一句使用prometheus2.28.0版本的代码: import "github.com/prometheus/prometheus/prompb" 我把require ...
- 关于3G移动通信网络中用户ip的配置过程的研究(中国电信cdma2000)
在RP口对ppp过程进行研究 PPP协商过程,如下图所示: 在建立ppp过程中pdsn需要与FAAA.HAAA交互.同时在分组数据业务进行过程中这种交互更加频繁,介绍如下,分为两种情况,简单ip,移动 ...
- 新手如何入门linux,linux原来还可以这么学
前言 在这个只有cangls和小白两人的小房间中,展开了一次关于学习方法的讨论. 小白:cangls啊,我想请教一个问题,您是如何记住那么多linux命令的. cangls:我啊,别人都看我的小电影, ...
- node.js和vue cli脚手架下载安装配置方法
一.node.js安装以及环境配置 1.下载vue.js 下载地址: https://nodejs.org/en/ 2.安装node.js 下载完成后,双击安装包开始安装.安装地址最好换成自己指定的地 ...
- python开发之函数
转:https://www.tuicool.com/wx/vEVrqeR 06 python开发之函数 博客园精华区12-12 20:56 06 python开发之函数 目录 6.2 调用函数与函数返 ...
- Java中的方法内联
Java中的方法内联 1. 什么是方法内联 例如有下面的原始代码: static class B { int value; final int get() { return value; } } pu ...
- 学习JAVAWEB第四天
# 今日内容 1. JDBC基本概念 2. 快速入门 3. 对JDBC中各个接口和类详解 ## JDBC: 1. 概念:Java DataBase Connectivity Java 数据库连接, J ...
- l线程池抓取lianjia
1. 线程池 的应用 from multiprocessing.dummy import Pool import requests from lxml import etree url="h ...