题解 matrix
无比毒瘤的dp题,而且伪装地好像很可做的样子
考场上我给它氪了差不多一个小时最后还是只能扔了个20pts状压走人
以下思路基本均来源于题解:
对于此题,题面中三个限制条件:
(1)第 i 行第 1~li 列恰好有 1 个 1。 (li+1到ri-1不能放1)
(2)第 i 行第 ri~m 列恰好有 1 个 1。
(3)每列至多有 1 个 1。
根本注意不到注意到条件(3)相对比较好转移,所以此题考虑dp列而不是dp行
那么考虑从中选择一列该如何转移
显然每列只有放1或不放1两种情况,但条件(1)(2)在干扰我们转移
所以尝试消除干扰: 仅考虑单一区间(以右区间为例)
显然只有左端点到达这一列的区间需要被考虑
那么令\(dp[i][j]\)为从左到右第i列,所有跨越第i列的右区间中有j个已放过1
对于不放1的情况,\(dp[i][j] += dp[i-1][j]\)
对于放1,先定义数组\(cntl[i]\)为第i列及其左侧左区间的右端点的个数,\(cntr[i]\)同理
则第i列可以放1的左端点有\(cntr[i]-(j-1)\)个,所以\(dp[i][j] += dp[i-1][j-1]\times (cntr[i]-(j-1))\)
上面仅考虑了跨越第i列的右区间的方案数,那么下面处理左区间
首先定义\(f[i][j]\)时j指的是「所有跨越第i列的右区间中有j个已放过1」
那就不必考虑第i列的1应该给左区间还是右区间了
对于这种类似左右两边对抗的方案数dp
(就是说类似只能 在某条线以左/以右/以此线为分界左右同时 进行某种操作)
(为什么我一想到这就想起alpha-beta对抗搜索啊,好像将答案区间划分的思路差不多?)
考虑拆分区间,拆分出对左侧产生影响的右区间
则其对左区间影响已知,可计算出左区间
对于此题,i列左侧总共能放\(i-(j-1)\)个1,需要放\(cntl[i]-cntl[i-1]\)个
这个\(cntl[i]-cntl[i-1]\)其实是有多少个左区间在i位置结束,即「新增加的必须放1的区间个数」
是在满足i列以前必须放1的区间都满足的条件下转移
\(cntl[i]-cntl[i-1]\)个1放到\(i-(j-1)\)个行中,考虑不同方案,应该是排列数
所以\(dp[i][j] *= A^{cntl[i]-cntl[i-1]}_{i-(j-1)}\)
就可以转移了
一大坑点: 这类n, m混杂的题一定要分清n, m!我预处理阶乘逆元的时候习惯性打了i<=n直接调了一晚上最后还是战神帮忙指出来的谢谢战神小可爱啦大雾逃
Code:
#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define N 3010
#define ll long long
#define ld long double
#define usd unsigned
#define ull unsigned long long
//#define int long long
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf, 1, 1<<21, stdin)), p1==p2?EOF:*p1++)
char buf[1<<21], *p1=buf, *p2=buf;
inline int read() {
int ans=0, f=1; char c=getchar();
while (!isdigit(c)) {if (c=='-') f=-f; c=getchar();}
while (isdigit(c)) {ans=(ans<<3)+(ans<<1)+(c^48); c=getchar();}
return ans*f;
}
int n, m;
int l[N], r[N], lcnt[N], rcnt[N];
ll dp[N][N], fac[N], inv[N];
const ll mod=998244353;
inline ll A(ll n, ll k) {return fac[n]*inv[n-k]%mod;}
inline ll md(ll a) {return a>=mod?a-mod:a;}
signed main()
{
#ifdef DEBUG
freopen("1.in", "r", stdin);
#endif
n=read(); m=read();
if ((n<<1)>m) {puts("0"); return 0;}
for (int i=1; i<=n; ++i) l[i]=read(), r[i]=read();
sort(l+1, l+n+1); sort(r+1, r+n+1);
fac[0]=fac[1]=1; inv[0]=inv[1]=1;
for (int i=2; i<=m; ++i) fac[i]=fac[i-1]*i%mod;
for (int i=2; i<=m; ++i) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
for (int i=2; i<=m; ++i) inv[i]=inv[i-1]*inv[i]%mod;
for (int i=1,p=1; i<=m; ++i) {rcnt[i]=rcnt[i-1]; while (p<=n&&r[p]==i) ++rcnt[i],++p;}
for (int i=1,p=1; i<=m; ++i) {lcnt[i]=lcnt[i-1]; while (p<=n&&l[p]==i) ++lcnt[i],++p;}
dp[0][0]=1;
for (int i=1; i<=m; ++i) {
dp[i][0] = dp[i-1][0]*A(i-lcnt[i-1], lcnt[i]-lcnt[i-1])%mod;
for (int j=1; j<=min(i, n); ++j)
dp[i][j] = md(dp[i-1][j]+dp[i-1][j-1]*max(rcnt[i]-(j-1), 0)%mod)*A(i-j-lcnt[i-1], lcnt[i]-lcnt[i-1])%mod; //, cout<<"A: "<<dp[i][j]<<endl;
}
printf("%lld\n", dp[m][n]);
return 0;
}
题解 matrix的更多相关文章
- 算法与数据结构基础 - 数组(Array)
数组基础 数组是最基础的数据结构,特点是O(1)时间读取任意下标元素,经常应用于排序(Sort).双指针(Two Pointers).二分查找(Binary Search).动态规划(DP)等算法.顺 ...
- [LeetCode 题解] Spiral Matrix
前言 [LeetCode 题解]系列传送门: http://www.cnblogs.com/double-win/category/573499.html 题目链接 54. Spiral Matrix ...
- 【题解】Sonya and Matrix Beauty [Codeforces1080E]
[题解]Sonya and Matrix Beauty [Codeforces1080E] 传送门:\(Sonya\) \(and\) \(Matrix\) \(Beauty\) \([CF1080E ...
- 【题解】#6622. 「THUPC 2019」找树 / findtree(Matrix Tree+FWT)
[题解]#6622. 「THUPC 2019」找树 / findtree(Matrix Tree+FWT) 之前做这道题不理解,有一点走火入魔了,甚至想要一本近世代数来看,然后通过人类智慧思考后发现, ...
- [题解]UVa 11082 Matrix Decompressing
开始眨眼一看怎么也不像是网络流的一道题,再怎么看也觉得像是搜索.不过虽然这道题数据范围很小,但也不至于搜索也是可以随随便便就可以过的.(不过这道题应该是special judge,因为一题可以多解而且 ...
- [LeetCode]题解(python):074-Search a 2D Matrix
题目来源 https://leetcode.com/problems/search-a-2d-matrix/ Write an efficient algorithm that searches fo ...
- [LeetCode]题解(python):059-Spiral Matrix II
题目来源 https://leetcode.com/problems/spiral-matrix-ii/ Given an integer n, generate a square matrix fi ...
- [LeetCode]题解(python):054-Spiral Matrix
题目来源 https://leetcode.com/problems/spiral-matrix/ Given a matrix of m x n elements (m rows, n column ...
- [LeetCode]题解(python):073-Set Matrix Zeroes
题目来源: https://leetcode.com/problems/set-matrix-zeroes/ 题意分析: 输入一个m×n矩阵,如果出现有0,那么将对应的行和列都变成0. 题目思路: 简 ...
随机推荐
- centos7 U盘安装及Raid划分的完整流程
目录 一.Centos7的新特性: 二.安装方法与准备工作(U盘镜像) 1. 安装方法介绍 2. Centos iso 常用镜像下载地址: 3. UltraISO制作U盘系统镜像 3.1 准备工作: ...
- Selenium的Css Selector使用方法
什么是Css Selector? Css Selector定位实际就是HTML的Css选择器的标签定位 工具 Css Selector的练习建议大家安装火狐浏览器后,下载插件,FireFinder 或 ...
- RabbitMQ 安装与界面管理
RabbitMQ 安装与界面管理 RabbitMQ概述 官网:https://www.rabbitmq.com/ RabbitMQ是部署最广泛的开源消息代理. RabbitMQ拥有成千上万的用户,是最 ...
- MapReduce处理简单数据
首先要说明的是,关于老师给的实验要求,我在网上看到了原文,原文地址:https://blog.csdn.net/qq_41035588/article/details/90514824,有兴趣的同学可 ...
- [刘阳Java]_Spring IOC程序代码如何编写_第3讲
第2讲我们介绍了Spring IOC的基本原理,这篇文章告诉大家Spring IOC程序代码是如何编写的,从而可以更好的理解IOC和DI的概念(所有的Java类的初始化工作扔给Spring框架,一个J ...
- Mybatis order by 动态传参出现的一个小bug
大家好,我是老三,一个平平无奇的CRUD仔. 今天,我正在愉快地CRUD,突然发现出现一个Bug,我们来看看是怎么回事吧! 问题由来 一个简单的需求,要求把和当前用户相关的数据置顶展示. 这里,我用了 ...
- 开源百宝箱《HelloGitHub》第 64 期
兴趣是最好的老师,HelloGitHub 让你对编程感兴趣! 简介 分享 GitHub 上有趣.入门级的开源项目. HelloGitHub 有实战.教程.黑科技.开源书籍.企业级开源项目,涵盖多种编程 ...
- 解决org.hibernate.LazyInitializationException的正确姿势
项目运行过程中,一个报错信息,报错信息如下: org.hibernate.LazyInitializationException: could not initialize proxy [xxx.do ...
- 构建前端第3篇之---使用scss
张艳涛 写于2021-1-20 主要解决俩个问题: 在单个vue文件中 <style rel="stylesheet/scss" lang="scss"& ...
- 微信小程序 -- 英语词典 (小程序插件)
英语词典小程序 基于英语词典小程序插件 - 提供开源地址 项目地址 英语词典小程序插件: 微信小程序 词典 真题基础服务插件(gitee.com) 功能特色 [x] 全面详实的经典词库,详细释义覆盖约 ...