洛咕

题意:原题面见链接,简单来说就是给出一个\(n*m\)的矩阵,每一行代表同一种烹饪方法,每一列代表同一种食材,\(a_{i,j}\)表示使用第i种烹饪方法第j种食材能做出多少种菜,要求至少做一道菜,而且每种烹饪方法至多用一次,并且每种食材不超过总菜数的一半,问有多少种做菜的方案?

分析:直观来讲就是在矩阵上选数,每一行要么不选数要么选一个数,每一列选的数不超过所有选数的一半。

如果只考虑前两个限制条件,设\(g[i][j]\)表示前i行选了j个数的方案数,则\(g[i][j]=g[i-1][j]+g[i-1][j-1]*sum[i]\),其中\(sum[i]\)表示第i行所有数(菜的种类)的和,则总方案数为\(\sum_{j=1}^{n}g[n][j]\)。

再考虑第三个限制条件,拿上面所求的总方案数减掉其中不合法的方案数(不满足第三个限制条件的),对于第\(k\)列,设\(f[i][j]\)表示前i行选的数中,在第k列的总数减去不在第k列的总数为j的情况下 的方案数,则\(f[i][j]=f[i-1][j]+f[i-1][j-1]*a[i][k]+f[i-1][j+1]*(sum[i]-a[i][k])\),注意\(j\)的范围应该是\([-n,n]\),所以可以统一右移\(n\)个单位,则不满足第三个限制条件的方案数为\(\sum_{k=1}^{m}\sum_{j=n+1}^{2n}f[n][j]\)。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
inline int read(){
int x=0,o=1;char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
if(ch=='-')o=-1,ch=getchar();
while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
return x*o;
}
const int N=105;
const int M=2005;
const int mod=998244353;
int a[N][M],sum[N],g[N][N],f[N][N*2];
int main(){
int n=read(),m=read();
for(int i=1;i<=n;++i){
for(int j=1;j<=m;++j){
a[i][j]=read();
sum[i]=(sum[i]+a[i][j])%mod;
}
}
g[0][0]=1;
for(int i=1;i<=n;++i){
g[i][0]=1;
for(int j=1;j<=i;++j){
g[i][j]=g[i-1][j]+(1ll*g[i-1][j-1]*sum[i])%mod;
g[i][j]%=mod;
}
}
ll tot1=0,tot2=0;
for(int j=1;j<=n;++j)tot1=(tot1+g[n][j])%mod;
for(int k=1;k<=m;++k){
memset(f,0,sizeof(f));
f[0][n]=1;//f[0][0]=1
for(int i=1;i<=n;++i)
for(int j=n-i;j<=n+i;++j){
f[i][j]=f[i-1][j];
f[i][j]=(f[i][j]+(1ll*f[i-1][j-1]*a[i][k])%mod)%mod;
f[i][j]=(f[i][j]+(1ll*f[i-1][j+1]*(sum[i]-a[i][k]))%mod)%mod;
}
for(int j=n+1;j<=n*2;++j)tot2=(tot2+f[n][j])%mod;
}
cout<<((tot1-tot2)%mod+mod)%mod<<endl;
return 0;
}

[CSP-S2019] Emiya 家今天的饭的更多相关文章

  1. 洛谷P5664 Emiya 家今天的饭 问题分析

    首先来看一道我编的题: 安娜写宋词 题目背景 洛谷P5664 Emiya 家今天的饭[民间数据] 的简化版本. 题目描述 安娜准备去参加宋词大赛,她一共掌握 \(n\) 个 词牌名 ,并且她的宋词总共 ...

  2. 洛谷P5664 Emiya 家今天的饭 题解 动态规划

    首先来看一道题题: 安娜写宋词 题目背景 洛谷P5664 Emiya 家今天的饭[民间数据] 的简化版本. 题目描述 安娜准备去参加宋词大赛,她一共掌握 \(n\) 个 词牌名 ,并且她的宋词总共有 ...

  3. Emiya家今天的饭 NOIP2019 (CSP?) 类DP好题 luoguP5664

    luogu题目传送门! 首先,硬求可行方案数并不现实,因为不好求(去年考场就这么挂的,虽然那时候比现在更蒟). 在硬搞可行方案数不行之后,对题目要求的目标进行转换: 可行方案数 = 总方案数 - 不合 ...

  4. 【CSP-S 2019】【洛谷P5664】Emiya 家今天的饭【dp】

    题目 题目链接:https://www.luogu.org/problem/P5664 Emiya 是个擅长做菜的高中生,他共掌握 \(n\) 种烹饪方法,且会使用 \(m\) 种主要食材做菜.为了方 ...

  5. 【CSP-S 2019】D2T1 Emiya 家今天的饭

    Description 传送门 Solution 算法1 32pts 爆搜,复杂度\(O((m+1)^n)\) 算法2 84pts 裸的dp,复杂度\(O(n^3m)\) 首先有一个显然的性质要知道: ...

  6. CSP2019 Emiya 家今天的饭 题解

    这题在考场上只会O(n^3 m),拿了84分.. 先讲84分,考虑容斥,用总方案减去不合法方案,也就是枚举每一种食材,求用它做超过\(\lfloor \frac{k}{2} \rfloor\) 道菜的 ...

  7. Emiya 家今天的饭

    \(dp_{i,j,k}\)表示前\(i\)种烹饪方法,假设最多的是食材\(j\),食材\(j\)比其他食材多\(k\)次出现 其中\(i \in [1,n],j \in [1,m],k \in [- ...

  8. 【NOIP/CSP2019】D2T1 Emiya 家今天的饭

    这个D2T1有点难度啊 原题: 花了我一下午的时间,作为D2T1的确反常 条件很奇怪,感觉不太直观,于是看数据范围先写了个暴力 写暴力的时候我就注意到了之前没有仔细想过的点,烹饪方式必须不同 虽然a很 ...

  9. csp2019 Emiya家今天的饭题解

    qwq 由于窝太菜了,实在是不会,所以在题解的帮助下过掉了这道题. 写此博客来整理一下思路 正文 传送 简化一下题意:现在有\(n\)行\(m\)列数,选\(k\)个数的合法方案需满足: 1.一行最多 ...

  10. 【JZOJ6433】【luoguP5664】【CSP-S2019】Emiya 家今天的饭

    description analysis 首先可以知道不符合要求的食材仅有一个,于是可以容斥拿总方案数减去选不合法食材的不合法方案数 枚举选取哪一个不合法食材,设\(f[i][j]\)表示到第\(i\ ...

随机推荐

  1. 2.2 在resources目录下,新建applicationContext.xml

    <?xml version="1.0" encoding="UTF-8"?> <beans> <!-- 这个bean标签的作用是, ...

  2. VMware-SSH协议的认证方式

    SSH1协议支持非对称密钥认证方式.口令认证,无法保证连接的完整性. SSH2协议支持SSH1协议支持的所有认证方式,增加数据保密性. 基于主机的认证方式[!不安全!] 当本地计算机收到执行远程命令的 ...

  3. UISelector

    1.UiSelector的基本方法 UiSelector对象可以理解为一种条件对象,描述的是一种条件,可以配合UiObject使用得到某个符合条件的控件对象. 所有的方法都是public的,且都返回U ...

  4. codeforce C. Maximal Intersection

    http://codeforces.com/contest/1029/problem/C 从第一天吃晚饭做到第二天吃完饭--你无法想象我的代码曾经150行 o( ̄┰ ̄*)ゞ 找到所有线段最远的左边和最 ...

  5. mybatis_pagehelper_selectOne的SQL语句被莫名的增加分页相关设置,暂定解决办法

    在使用mybatis.以及其分页插件pagehelper时,原本的一个selectOne的sql语句被莫名(原因未知)的加上了分页相关,引起如下异常: exception is org.apache. ...

  6. 自动删除几天前的备份集文件脚本 for windows

    将以下内容根据自己想要备份的库修改后保存为del.bat,加入到计划任务中 @echo 删除指定路径下指定N天前的文件 @echo off@echo.------------------------- ...

  7. v-if和v-show最重要一点

    vif是新建与销毁标签 show是显示与隐藏标签 如果涉及到权限操作的话 必须使用vif 因为vshow中把标签改成display:none的话可以手动改成显示

  8. lombok.config

    # 声明该配置文件是一个根配置文件,从该配置文件所在的目录开始扫描 config.stopBubbling=true # 全局配置 equalsAndHashCode 的 callSuper 属性为t ...

  9. JS form表单数据校验及失效情况下的解决方案

    如下图,当执行提交操作之前,我们需要对序号,要求完成时间,责任人,措施内容四项进行非空,字符长度及输入内容的类型进行校验. 直接贴样式代码 <div class="wrapper an ...

  10. Unity动态修改材质球RenderingMode属性

    Material四个模式 动态修改代码 using System.Collections; using System.Collections.Generic; using UnityEngine; p ...