洛谷 P5664 Emiya 家今天的饭 题解

最新推荐文章于 2022-06-06 18:16:57 发布
最新推荐文章于 2022-06-06 18:16:57 发布
阅读量339 收藏 1
点赞数
文章标签： 动态规划
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xsaploser/article/details/103652358
版权
  • 纯暴力（32 pts）
    • 开始的第一思路，可以过测试点1-8
    • 直接枚举出每种烹饪方法所用的主要食材
    • 用搜索枚举，传入两个参数：现在枚举的烹饪方法编号以及在此之前的方案数
    • 代码：
#include <iostream>
#include <cstdio>

const long long mod = 998244353 ;

using namespace std ;

int n , m ;
int a[105][2005] ;
int cnt[2005] ;
long long ans ;

void input() 
{
	scanf("%d%d" , &n , &m) ;
	for ( int i = 1 ; i <= n ; i++ )
	{
		for ( int j = 1 ; j <= m ; j++ )
		{
			scanf("%d" , &a[i][j]) ;
		}
	}
}

void dfs(int x , long long sum)
{
	if ( x > n )
	{
		long long num = 0 ;
		for ( int i = 1 ; i <= m ; i++ )
		{
			num += cnt[i] ;
		}
		num /= 2 ;
		for ( int i = 1 ; i <= m ; i++ )
		{
			if ( cnt[i] > num )
				return ;
		}
		ans += sum ;
		ans = ans % mod ;
		return ;
	}
	dfs(x + 1 , sum) ;
	for ( int i = 1 ; i <= m ; i++ )
	{
		if ( a[x][i] == 0 )
			continue ;
		cnt[i]++ ;
		dfs(x + 1 , sum * a[x][i] % mod) ;
		cnt[i]-- ;
	}
}

int main()
{
	input() ;
	dfs(1 , 1) ;
	printf("%lld\n" , ans - 1) ;
	return 0 ;
}
  • 用dp过9-12（48pts）
    • 在写完32pts后，想先考虑m = 2（显然，m = 2会比其他简单），发现只有两种主要食材，所以必须选第一种食材的次数等于选第二种食材的次数
    • 想到包裹分拣，可以开三位分别表示现在决策的烹饪方法编号、现在选的主要食材一的个数以及现在选的主要食材二的个数
    • 之后动态规划一下就可以了
    • 代码：
#include <iostream>
#include <cstdio>

const long long mod = 998244353 ;

using namespace std ;

int n , m ;
int a[105][2005] ;
int cnt[2005] ;
long long ans ;
long long dp[105][55][55] ;

void input() 
{
	scanf("%d%d" , &n , &m) ;
	for ( int i = 1 ; i <= n ; i++ )
	{
		for ( int j = 1 ; j <= m ; j++ )
		{
			scanf("%d" , &a[i][j]) ;
		}
	}
}

void dfs(int x , long long sum)
{
	if ( x > n )
	{
		long long num = 0 ;
		for ( int i = 1 ; i <= m ; i++ )
		{
			num += cnt[i] ;
		}
		num /= 2 ;
		for ( int i = 1 ; i <= m ; i++ )
		{
			if ( cnt[i] > num )
				return ;
		}
		ans += sum ;
		ans = ans % mod ;
		return ;
	}
	dfs(x + 1 , sum) ;
	for ( int i = 1 ; i <= m ; i++ )
	{
		if ( a[x][i] == 0 )
			continue ;
		cnt[i]++ ;
		dfs(x + 1 , sum * a[x][i] % mod) ;
		cnt[i]-- ;
	}
}

int main()
{
	input() ;
	if ( n <= 10 )
	{
		dfs(1 , 1) ;
		printf("%lld\n" , ans - 1) ;
		return 0 ;
	}
	if ( m == 2 ) 
	{
		ans = 0 ;
		dp[0][0][0] = 1 ;
		for ( int i = 1 ; i <= n ; i++ )
		{
			for ( int j = 0 ; j <= n / 2 ; j++ )
			{
				for ( int k = 0 ; k <= n / 2 ; k++ )
				{
					dp[i][j][k] = dp[i - 1][j][k] ;
					if ( j > 0 )
					{
						dp[i][j][k] = (dp[i][j][k] + (dp[i - 1][j - 1][k] * a[i][1]) % mod) % mod ; 
					}
					if ( k > 0 )
					{
						dp[i][j][k] = (dp[i][j][k] + (dp[i - 1][j][k - 1] * a[i][2]) % mod) % mod ; 
					}
				}
			}
		}
		for ( int i = 1 ; i <= n / 2 ; i++ )
		{
			ans = ( ans + dp[n][i][i] ) % mod ;
		}
		printf("%lld\n" , ans) ;
		return 0 ;
	}
	return 0 ;
}
  • m = 3（64pts）
    • 由m = 2的dp想到m = 3也可用dp，开四维处理即可
    • 直接dp一下即可
    • 代码几乎与m = 2相同
#include <iostream>
#include <cstdio>

const long long mod = 998244353 ;

using namespace std ;

int n , m ;
int a[105][2005] ;
int cnt[2005] ;
long long ans ;
long long dp[105][55][55] ;
long long dp2[105][55][55][55] ;

void input() 
{
	scanf("%d%d" , &n , &m) ;
	for ( int i = 1 ; i <= n ; i++ )
	{
		for ( int j = 1 ; j <= m ; j++ )
		{
			scanf("%d" , &a[i][j]) ;
		}
	}
}

void dfs(int x , long long sum)
{
	if ( x > n )
	{
		long long num = 0 ;
		for ( int i = 1 ; i <= m ; i++ )
		{
			num += cnt[i] ;
		}
		num /= 2 ;
		for ( int i = 1 ; i <= m ; i++ )
		{
			if ( cnt[i] > num )
				return ;
		}
		ans += sum ;
		ans = ans % mod ;
		return ;
	}
	dfs(x + 1 , sum) ;
	for ( int i = 1 ; i <= m ; i++ )
	{
		if ( a[x][i] == 0 )
			continue ;
		cnt[i]++ ;
		dfs(x + 1 , sum * a[x][i] % mod) ;
		cnt[i]-- ;
	}
}

int main()
{
	input() ;
	if ( n <= 10 )
	{
		dfs(1 , 1) ;
		printf("%lld\n" , ans - 1) ;
		return 0 ;
	}
	if ( m == 2 ) 
	{
		ans = 0 ;
		dp[0][0][0] = 1 ;
		for ( int i = 1 ; i <= n ; i++ )
		{
			for ( int j = 0 ; j <= n / 2 ; j++ )
			{
				for ( int k = 0 ; k <= n / 2 ; k++ )
				{
					dp[i][j][k] = dp[i - 1][j][k] ;
					if ( j > 0 )
					{
						dp[i][j][k] = (dp[i][j][k] + (dp[i - 1][j - 1][k] * a[i][1]) % mod) % mod ; 
					}
					if ( k > 0 )
					{
						dp[i][j][k] = (dp[i][j][k] + (dp[i - 1][j][k - 1] * a[i][2]) % mod) % mod ; 
					}
				}
			}
		}
		for ( int i = 1 ; i <= n / 2 ; i++ )
		{
			ans = ( ans + dp[n][i][i] ) % mod ;
		}
		printf("%lld\n" , ans) ;
		return 0 ;
	}
	else if ( m == 3 )
    {
    	dp2[0][0][0][0] = 1 ;
        for ( int i = 1 ; i <= n ; i++ )
            for ( int j = 0; j <= n; j++ )
                for ( int k = 0 ; k <= i - j ; k++ )
                	for ( int l = 0 ; l <= i - j - k ; l++ )
	                {
	                    dp2[ i ][ j ][ k ][ l ] = 
							dp2[ i - 1 ][ j ][ k ][ l ] ;
	                    if ( j >= 1 ) 
							dp2[ i ][ j ][ k ][ l ] += 
							dp2[ i - 1 ][ j - 1 ][ k ][ l ] * 
							a[ i ][ 1 ] % mod ;
	                    if ( k >= 1 ) 
							dp2[ i ][ j ][ k ][ l ] += 
							dp2[ i - 1 ][ j ][ k - 1 ][ l ] * 
								a[ i ][ 2 ] % mod ;
	                    if ( l >= 1 ) 
							dp2[ i ][ j ][ k ][ l ] += 
							dp2[ i - 1 ][ j ][ k ][ l - 1 ] * 
								a[ i ][ 3 ] % mod ;
	                }
        for ( int i = 0 ; i <= n / 2 ; i++ ) 
        	for ( int j = 0 ; j <= n / 2 ; j++ ) 
        		for ( int k = 0 ; k <= n / 2 ; k++ ) 
        		{
        			if ( i <= j + k && j <= i + k && k <= i + j ) 
        				ans += dp2[ n ][ i ][ j ][ k ] ;
        			ans = ans % mod ;
        		}
        printf( "%lld\n" , ans - 1 ) ;
    }
	return 0 ;
}
  • n = 40 && m = 500 (84pts)
    • 发现Yazid的条件是相对较难做到，从Yazid的条件下手讨论
    • 做完64pts后，感觉不太容易用顺向思维p做（m = 500肯定会炸空间）
    • 于是想到逆推（以）考虑，这是先算出总方案（定义sum[i]表示用第i种烹饪方法所能做出的总方案数，则总方案数为（sum[1] + 1）×（sum[2] + 1）× …… ×（sum[n] + 1））
    • 如果不满足Yazid的条件，则必然有一种菜品作为主要食材用了至少n / 2 + 1次，枚举这种菜品，之后将菜品分为主菜（即枚举出的菜品）与非主菜，显然这其中每类没有交集，所以可直接求出不满足Yazid的条件的总和,相减再减一即可
    • 其中枚举过程中只要用m = 2的解法即可
    • 代码：
#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std ;

const int mod = 998244353 ; 
const int N = 2100 ;
int n , m ;
int x[105][2005] ;
long long sum[105] ;
long long a[2005] , b[2005] ;
long long dp[105][105][105] ;

void input()
{
	scanf("%d%d" , &n , &m) ;
	for ( int i = 1 ; i <= n ; i++ )
	{
		for ( int j = 1 ; j <= m ; j++ )
		{
			scanf("%d" , &x[i][j]) ;
		}
	} 
	return ;
}

void doit()
{
	long long tot = 1 ;
	for ( int i = 1 ; i <= n ; i++ )
	{
		sum[i] = 0 ;
		for ( int j = 1 ; j <= m ; j++ )
		{
			sum[i] += x[i][j] ;
			sum[i] %= mod ;
		}
		tot = ( tot * ( sum[i] + 1 ) ) % mod ;
	}
	long long wrong = 0 ;
	for ( int food = 1 ; food <= m ; food++ )
	{
		memset(dp , 0 , sizeof(dp)) ;
		for ( int i = 1 ; i <= n ; i++ )
			a[i] = x[i][food] , b[i] = (sum[i] - a[i] + mod) % mod ;
		dp[0][0][0] = 1 ;
		for ( int i = 1 ; i <= n ; i++ )
			for ( int j = 0 ; j <= i ; j++ )
				for ( int k = 0 ; k <= i - j ; k++ )
				{
					dp[i][j][k] = dp[i - 1][j][k] ;
					if ( j > 0 )
					{
						dp[i][j][k] = (dp[i][j][k] + (dp[i - 1][j - 1][k] * a[i]) % mod) % mod ; 
					}
					if ( k > 0 )
					{
						dp[i][j][k] = (dp[i][j][k] + (dp[i - 1][j][k - 1] * b[i]) % mod) % mod ; 
					}
				}
		for ( int i = 1 ; i <= n ; i++ )
			for ( int j = 0 ; j <= n - i ; j++ )
				if ( i > j )
					wrong += dp[n][i][j] ;
		wrong %= mod ;
	}
	printf("%lld\n" , (tot - 1 - wrong + mod) % mod) ;
} 

int main()
{
	input() ;
	doit() ;
	return 0 ;
}
  • 100pts
    • 思路与84pts的几乎一样，只是发现dp其实关心的只有差，所以可以用差值dp优化为O（n²），总复杂度O（n³）
    • 注意差值dp时必须平移，因为可能差为负数
    • 代码:

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std ;

const int mod = 998244353 ; 
const int N = 2100 ;
int n , m ;
int x[105][2005] ;
long long sum[105] ;
long long a[2005] , b[2005] ;
long long dp[105][2105] ;

void input()
{
	scanf("%d%d" , &n , &m) ;
	for ( int i = 1 ; i <= n ; i++ )
	{
		for ( int j = 1 ; j <= m ; j++ )
		{
			scanf("%d" , &x[i][j]) ;
		}
	} 
	return ;
}

void doit()
{
	long long tot = 1 ;
	for ( int i = 1 ; i <= n ; i++ )
	{
		sum[i] = 0 ;
		for ( int j = 1 ; j <= m ; j++ )
		{
			sum[i] += x[i][j] ;
			sum[i] %= mod ;
		}
		tot = ( tot * ( sum[i] + 1 ) ) % mod ;
	}
	long long wrong = 0 ;
	for ( int food = 1 ; food <= m ; food++ )
	{
		memset(dp , 0 , sizeof(dp)) ;
		for ( int i = 1 ; i <= n ; i++ )
			a[i] = x[i][food] , b[i] = (sum[i] - a[i] + mod) % mod ;
		dp[0][N] = 1 ;
		for ( int i = 1 ; i <= n ; i++ )
		{
			for ( int j = -1 * i ; j <= i ; j++ )
			{
				dp[i][N + j] = dp[i - 1][N + j] ;
				if ( N + j - 1 >= 0 )
					dp[i][N + j] += ( dp[i - 1][N + j - 1] * a[i] % mod ) ;
				dp[i][N + j] += ( dp[i - 1][N + j + 1] * b[i] % mod ) ;
				dp[i][N + j] %= mod ; 
			}
		}
		for ( int i = 1 ; i <= n ; i++ )
			wrong += dp[n][N + i] ;
		wrong %= mod ;
	}
	printf("%lld\n" , (tot - 1 - wrong + mod) % mod) ;
} 

int main()
{
	input() ;
	doit() ;
	return 0 ;
}
xsaploser
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
1989：【19CSPS提高组】Emiya家今天的饭
11-10
c++noip完全AC代码（信奥赛昵称Kirito榜四可查），这题不算特别难，但也不简单，大概花了我两个小时左右吧。。。学习忙，没去别的平台上测试，请见谅
P5664 [CSP-S2019] Emiya 家今天的饭
GDW07的博客
06-22 779
P5664 [CSP-S2019] Emiya 家今天的饭
Emiya 家今天的饭
C202044zxy的博客
11-19 647
一、题目 点此看题 二、解法 0x01 暴力 可以考虑搜索，时间复杂度O(mn)O(m^n)O(mn)，可以拿到323232分的高分。 考虑m≤3m\leq 3m≤3的部分分，直接四维dpdpdp，定义dp[i][j][k][l]dp[i][j][k][l]dp[i][j][k][l]为前iii行分别选j,k,lj,k,lj,k,l个，转移也不难，时间复杂度O(n4m)O(n^4m)O(n4m)，...
NOIP2019 Emiya家今天的饭
SUNSET
11-17 2017
NOIP2019 Emiya家今天的饭 ACM退役选手远程口胡 考虑使用容斥方法: 1 采用动态规划,先求出在无限制情况下,安排kkk种烹饪方法总的方案数. 记dp2[i][j]dp2[i][j]dp2[i][j]表示已经考虑完前iii种烹饪方法,共做了jjj个菜的方案数. 那么显然,决策分2种情况,用或不用第iii种烹饪方法,用的话就只能选一种主要食材. dp2[i][j]=dp2[i−1][j...
CSP-S2 Emiya 家今天的饭
fengqiyuka的博客
11-30 329
题目描述 题目分析 这一道题本身算法并不难，但思考难度也不低。 “烹饪方法互不相同”这个条件用DP可以很好解决，关键是每一个食材出现次数不超过k2\frac{k}{2}2k​，接下来我们所说的不合法均指不满足该条件。 本题最关键的一点便是——容斥。考虑先求出总的方案数，再减去不合法的方案数。 总数求法很简单，那么如何求不合法的方案数呢？我们考虑先枚举超出限制的食材。因为超限的食材会出现超过k...
[csp2019]Emiya家今天的饭
NephrenRuqInsania
11-26 813
作为提高组d2t1d2t1d2t1，比去年难 所以这道题我打的特别的差 32pts 这道题我们很显然可以看到可以打一个暴力 复杂度o(n∗n!)o(n*n!)o(n∗n!) 我考场上就达到了这里——我太菜了 void dfs(int u,ll plus){ if(u==n+1){ if(!tot)return; Rep(i,1,m)if(cnt[i]>tot/2)return; ...
博客系统solo的第三方主题emiya
08-11
博客系统solo的第三方主题 emiya
语音合成代码matlab-peass-software:来自http://bass-db.gforge.inria.fr/peass/PEAS
05-27
Emiya和Emmanuel Vincent，法国INRIA。 这是Valentin Emiya（法国INRIA），Emmanuel Vincent（法国INRIA），Niklas Harlander（德国奥尔登堡大学），Volker Hohmann（德国奥尔登堡大学）的联合工作。 减少Audionamix...
MATLAB小波软硬阈值去噪代码-HPSS:高压钠灯
05-20
另一个目标是使用PEASS（Emiya等，2011）和APE（De Man和Reiss，2014）工具包评估一些结果。 但是，后者没有得到维护，因此与MATLAB的较新版本不兼容。 因此，使用了另一个名为PEAQ的感知评估工具包（Kabal 2002）。...
clnaga.github.io
03-30
Emiya的摘录室1.优化算法于2021.3.29更新
luoguP5664 Emiya 家今天的饭
到处遛鸟刑警
03-14 122
考虑维护一个差值 每次枚举是那个食材爆炸。。。然后就ok了 嘛，建议看题解这个题。。。 #include<bits/stdc++.h> using namespace std; long long n,m,mod = 998244353; long long a[105][2005],s[105],f[205][205],sum; void init(){ memset(f,0,...
[CSP-S2019] Emiya 家今天的饭
qq_43649491的博客
06-06 117
Day 6+7 DP（那一年的day2题1，完全被支配了┭┮﹏┭┮）代码如下（示例）：
P5664 [CSP-S2019] Emiya 家今天的饭（dp+容斥）
m0_51841071的博客
09-14 338
[题目传送门]
「CSP-S 2019」 Emiya 家今天的饭 题解
Clever_Hard的博客
03-09 326
题面 样例 2 3 1 0 1 0 1 1 3 分析 考虑正着限制最大的数不超过一半不好做，那我们可以反着来。 令 dp[i][j][k]dp[i][j][k]dp[i][j][k] 为第 iii 行，此列指定数 pointpointpoint 数量为 jjj 个，选的非指定数的个数为 kkk 的方案数。 则一共有三种情况 Case 1:Case\ 1:Case 1: dp[i][j][k]=dp[i−1][j][k]dp[i][j][k] = dp[i - 1][j][k]d
【CSP-S 2019】【洛谷P5664】Emiya 家今天的饭【dp】
stoorz
11-23 682
菜到爆炸qwq，84pts的出血大放送都没拿到
Emiya 家今天的饭's 题解
Belief_yfly的博客
11-30 258
又是一道令人无语的D2T1。。。 考场搞了半天还是只能敲暴力。。。 题目链接 题目太繁琐。。。简单点说就是找满足条件的组合数。 直接进暴力 深搜暴力。。。（因为我正解的转移推不出来） 首先，第一维表示目前枚举到了第几种做法，第二维表示选了几种做法，第三维表示能组成的组合数。 到当前层，你可以枚举每个不为0的来统计，如果发现可以满足条件，累乘上去，搜索下一层。最后再不要这层枚举，再有一个简单的剪枝，...
【luoguP5664】 Emiya 家今天的饭 动态规划
EM LGH
01-02 176
现在看来，这道题的 84 pts 随便拿啊. 100 pts 也不是很难，但是考场上太过于紧张加上心理素质不好吧... 84pts 暴力的话设状态 $f[i][j][k]$ 表示决策到第 $i$ 行，枚举的拿了 $j$ 个，其余拿了 $k$ 个. 但是我们把那个不等式拆开后发现其实表示的就是钦定的比其余的多拿就行. 所以我们只关注那个差值就好了. 令 $f[i][j]...
洛谷P5664 Emiya 家今天的饭
lzx_xzl_______的博客
12-25 218
题目描述 Emiya 是个擅长做菜的高中生，他共掌握nn种烹饪方法，且会使用mm种主要食材做菜。为了方便叙述，我们对烹饪方法从1 \sim n1∼n编号，对主要食材从1 \sim m1∼m编号。 Emiya 做的每道菜都将使用恰好一种烹饪方法与恰好一种主要食材。更具体地，Emiya 会做a_{i,j}ai,j​道不同的使用烹饪方法ii和主要食材jj的菜（1 \leq i \leq n, 1 \leq j \leq m1≤i≤n,1≤j≤m），这也意味着 Emiya 总共会做...
【题解】LuoGu5664：Emiya 家今天的饭
ModestCoder_的博客
01-20 340
原题传送门 每行最多取一个告诉我们可以枚举行 所以这道题目总体复杂度里肯定有行的复杂度O(n)O(n)O(n) 考场上写的是m=2/3m=2/3m=2/3的暴力，直接把每一列分别取了几个写到状态里面去 正解需要考虑正难则反 我们写直接将每行最多一个当作大前提 要求的是没有一列取了超过一半的东西 那么这个答案其实 =所有方案-有一列取了超过一半的方案 这里就有一个重要性质：最多只有一列会超过一半（虽...
用c语言设计一个程序，用于维吉尼亚密码破译，以输入密文：M mu rhi nwle sr uw saazb，输入密钥：emiya，输出明文：I am the bone of my sword为例
最新发布
05-18
char key[] = "emiya"; // 密钥 char plaintext[100]; // 存储解密后的明文 int len = strlen(ciphertext); // 密文的长度 int i, j, k; int key_len = strlen(key); // 密钥的长度 // 对每一组密文进行解密...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值