[CSP-S模拟测试]:matrix(DP)
题目描述
求出满足以下条件的$n\times m$的$01$矩阵个数:
(1)第$i$行第$1~l_i$列恰好有$1$个$1$。
(2)第$i$行第$r_i~m$列恰好有$1$个$1$。
(3)每列至多有$1$个$1$。
输入格式
第一行两个整数$n,m$。接下来$n$行每行$2$个整数$l_i,r_i$。
输出格式
一行一个整数表示答案。对998244353取模。
样例
样例输入
2 6
2 4
5 6
样例输出
12
数据范围与提示
对于$20\%$的数据,$n,m\leqslant 12$。
对于$40\%$的数据,$n,m\leqslant 50$。
对于$70\%$的数据,$n,m\leqslant 300$。
对于$100\%$的数据,$n,m\leqslant 3000$,$1\leqslant l_i<r_i\leqslant m$。
题解
看到这道题,首先应该想到组合数,然后……
我一开始想的是容斥,正好能模过样例,但是忽然发现是一个多步容斥,无语……
$20\%$算法:
使劲搜就好了,别搜错了,记得别用clock(),递归函数里clock()返回值玄学(考场上被这东西干没了40分……)。
时间复杂度:$\Theta(m^n)~\Theta(m^{\frac{n^2}{4}})$。
期望得分:$20$分。
$100\%$算法:
考虑$DP$,定义有点意思,定义$dp[i][j]$表示当前到了第$i$列,已经有$j$行在右侧区间放$1$的方案数。
下面来解释一下,注意右侧区间不是第$i$列的右侧,而是$r_i$,理解这东西我用了半个小时……
下面在来看一下如何转移:
首先,要预处理一些东西,用一个$fl$数组表示在第$i$列以前结束的左区间的个数,$fr$数组表示在地$i$列以前开始的右区间的个数,注意这里的左区间和右区间也指的是$l_i$和$r_i$。
然后就来明确两个值:$i-j-l[i-1]$就是第$i$列左侧$1$的个数,$l[i]-l[i-1]$就是在当前列结束的左区间的个数。
那么,将有$4$种转移方式:
$\alpha.$如果$(i-j-l[i-1])<l[i]-l[i-1]$那么当前的$i$列$j$行以及往后更大的$j$都不可能贡献方案数,所以直接$break$即可。
$\beta.$$dp[i][j]=dp[i][j]\times A_{i-j-l[i-1]}^{l[i]-l[i-1]}$,需要用到这种转移是因为我们在下面两种转移的时候不方便执行这种操作。
$\gamma.$$dp[i-1][j]+=dp[i][j]$,这种就是不在这列放$1$,进行直接转移。
$\delta.$$dp[i][j]=dp[i-1][j-1]\times (r[i]-j+1)$这种就是放$1$,方案数为$r[i]-j+1$。
排列数$A$可以选择直接用$\Theta(n\times m)$进行打表即可。
时间复杂度:$\Theta(n\times m)$。
期望得分:$100$分。
代码时刻
$20\%$算法:
#include<bits/stdc++.h>
#define mod 998244353
using namespace std;
int n,m;
struct rec{int l,r;}q[3001];
bool vis[3001];
int ans;
int lmax,rmin=10000;
int qpow(int x,int y)
{
int res=1;
while(y)
{
if(y&1)res=res*x%mod;
x=x*x%mod;
y>>=1;
}
return res;
}
void dfs(int x)
{
if(x>n)
{
ans=(ans+1)%mod;
return;
}
for(int i=1;i<=q[x].l;i++)
for(int j=q[x].r;j<=m;j++)
if(!vis[i]&&!vis[j])
{
vis[i]=vis[j]=1;
dfs(x+1);
vis[i]=vis[j]=0;
}
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%d%d",&q[i].l,&q[i].r);
lmax=max(lmax,q[i].l);
rmin=min(rmin,q[i].r);
}
dfs(1);
printf("%d",ans);
return 0;
}
$100\%$算法:
#include<bits/stdc++.h>
using namespace std;
int n,m;
int l[5000],r[5000],A[5000][5000],dp[5000][5000];
void pre_work()
{
A[0][0]=1;
for(int i=1;i<=m;i++)
{
A[i][0]=1;
for(int j=1;j<=i;j++)
A[i][j]=(1LL*A[i][j-1]*(i-j+1))%998244353;
}
}
int main()
{
scanf("%d%d",&n,&m);
pre_work();
for(int i=1;i<=n;i++)
{
int x,y;
scanf("%d%d",&x,&y);
l[x]++;r[y]++;
}
for(int i=1;i<=m;i++)
{
l[i]+=l[i-1];
r[i]+=r[i-1];
}
dp[1][0]=1;
for(int i=1;i<=m;i++)
{
for(int j=0;j<=r[i];j++)
{
if(i-j<l[i])break;
dp[i][j]=1LL*dp[i][j]*A[i-j-l[i-1]][l[i]-l[i-1]]%998244353;
dp[i+1][j]=(dp[i+1][j]+dp[i][j])%998244353;
dp[i+1][j+1]=(dp[i+1][j+1]+1LL*dp[i][j]*(r[i+1]-j)%998244353)%998244353;
}
}
cout<<dp[m][n];
return 0;
}
rp++
[CSP-S模拟测试]:matrix(DP)的更多相关文章
- noi2019模拟测试赛(四十七)
noi2019模拟测试赛(四十七) T1与运算(and) 题意: 给你一个序列\(a_i\),定义\(f_i=a_1\&a_2\&\cdots\&a_i\),求这个序列的所 ...
- [考试反思]1109csp-s模拟测试106:撞词
(撞哈希了用了模拟测试28的词,所以这次就叫撞词吧) 蓝色的0... 蓝色的0... 都该联赛了还能CE呢... 考试结束前15分钟左右,期望得分300 然后对拍发现T2伪了写了一个能拿90分的垃圾随 ...
- [考试反思]1003csp-s模拟测试58:沉淀
稳住阵脚. 还可以. 至少想拿到的分都拿到了,最后一题的确因为不会按秩合并和线段树分治而想不出来. 对拍了,暴力都拍了.挺稳的. 但是其实也有波折,险些被卡内存. 如果内存使用不连续或申请的内存全部使 ...
- [考试反思]0814NOIP模拟测试21
前两名是外校的240.220.kx和skyh拿到了190的[暴力打满]的好成绩. 我第5是170分,然而160分就是第19了. 在前一晚上刚刚爆炸完毕后,心态格外平稳. 想想前一天晚上的挣扎: 啊啊啊 ...
- csp-s模拟测试98
csp-s模拟测试98 $T1$??不是我吹我轻松手玩20*20.$T2$装鸭好像挺可做?$T3$性质数据挺多提示很明显? $One$ $Hour$ $Later$ 这$T1$什么傻逼题真$jb$难调 ...
- csp-s模拟测试97
csp-s模拟测试97 猿型毕露.水题一眼秒,火题切不动,还是太菜了. $T1$看了一会儿感觉$woc$期望题$T1??$假的吧??. $T2$秒. $T3$什么玩意儿. 40 01:24:46 00 ...
- csp-s模拟测试95
csp-s模拟测试95 去世场祭. $T1$:这不裸的除法分块吗. $T2$:这不裸的数据结构优化$Dp$吗. $T3$:这不裸的我什么都不会搜索骗$30$分吗. 几分钟后. 这除法分块太劲了..(你 ...
- csp-s模拟测试93
csp-s模拟测试93 自闭场. $T1$想到$CDQ$,因为复杂度少看见一个$0$打了半年还用了$sort$直接废掉,$T2$,$T3$直接自闭暴力分都没有.考场太慌了,心态不好. 02:07:34 ...
- csp-s模拟测试92
csp-s模拟测试92 关于$T1$:最短路这一定建边最短路. 关于$T2$:傻逼$Dp$这一定线段树优化$Dp$. 关于$T3$:最小生成树+树P+换跟一定是这样. 深入(?)思考$T1$:我是傻逼 ...
随机推荐
- 初识HTML标签
web概念概述 JavaWeb: 使用Java语言开发基于互联网的项目 软件架构: C/S: Client/Server 客户端/服务器端 在用户本地有一个客户端程序,在远程有一个服务器端程序 如:Q ...
- 数据分析 - 缺失值、异常值、一致性分析方法及Python实现
1.数据质量分析 数据质量分析主要任务:检查原始数据是否存在脏数据. 脏数据: 缺失值 异常值 不一致的值 重复数据及含有特殊符号(如:#.¥.*)的数据 1.1 缺失值分析 数据缺失主要包括:记录的 ...
- 解决IDEA中自动生成返回值带final修饰的问题
修改配置文件: Editor--Code Style--Java--Code Generation--将Make generated local variables final勾选上
- C++中的函数重载分析(一)
1,重载是 C 语言到 C++ 语言的一个飞跃,C 语言中没有重载的概念,所有的函数 名是不允许有重复的,在 C++ 中因为引进了重载,所以函数名可以重复: 2,自然语言中的上下文: 1,你知道上面词 ...
- python3—廖雪峰之练习(三)
列表生成式练习 请修改列表生成式,通过添加if语句保证列表生成式能正确执行: L1 = ['Hello', 'World', 18, 'Apple', None] L2 = [] for x in L ...
- dfs(最佳路径)
http://acm.hdu.edu.cn/showproblem.php?pid=1242 Rescue Time Limit: 2000/1000 MS (Java/Others) Memo ...
- 【LGR-063】洛谷11月月赛 I & MtOI2019 Ex Div.2
[MtOI2019]黑蚊子多: 送分向水题,直接模拟即可. #include<iostream> #include<cstdio> #define N 1505 using n ...
- 批处理遍历文件夹执行git pull
echo off & color 0A for /d %%f in (D:\www\*) do ( D: cd %%f chdir git pull ) pause 遍历D:\www\这个文件 ...
- TMS320F28335——SCI串口
一.IO配置 以SCIA为例:使用的是GPIO35--SCITXDA GPIOA36--SCIRXDA 使用寄存器: GPBPUD :设置上拉 GPIO32-GPIO63 对应位0 使能上 ...
- ECMAScript 6 学习笔记(一)
ECMAScript 6简介 ECMAScript 6.0(以下简称ES6)是JavaScript语言的下一代标准,已经在2015年6月正式发布了.它的目标,是使得JavaScript语言可以用来编写 ...