题目传送门

盒子和小球之二:N个有差别的盒子(1<=N<=20)。你有A个红球和B个蓝球。0 <= A <= 15, 0 <= B <= 15。球除了颜色没有任何区别。你可以将球放进盒子。一个盒子可以同时放进两种球,也可以只放一种,也可以空着。球不必全部放入盒子中。编程计算有多少种放置球的方法。

考虑动态规划。一个盒子可以同时放进两种球,但是每个盒子我们可以对他进行很多操作,所以也就可以近似看成每个盒子可以放无数球,但是这并不重要。设f[i][j][k]表示我们当前放到第i个盒子,已经放了j个红球,k个篮球的状态方案数。(因为我们要求在状态设计的时候要求顾及到题目中的重要量,而这个状态恰好达到了这点)另外,由于本题数据范围较小,所以我们思考的不用太复杂。

则我们可以理所当然地想到转移:f[i][j][k]=sigmaf[i-1][j-x][k-x],其中x是我们枚举的在当前盒子中放几个。

时间复杂度为O(n*A^2*B^2),可以轻松通过。

目标答案显然我们可以累加f[n][i][j],枚举ij。

 #include<cstdio>
#include<algorithm> using namespace std;
typedef long long ll;
//complexity n*A^2*B^2
int n,A,B;
ll ans,f[][][]; int main()
{
scanf("%d%d%d",&n,&A,&B);
f[][][]=;//初值,很重要
for(int i=;i<=n;i++)
for(int j=;j<=A;j++)//注意放几个从0开始枚举
for(int k=;k<=B;k++)
for(int a=;a<=j;a++)
for(int b=;b<=k;b++)
f[i][j][k]+=f[i-][j-a][k-b];
for(int i=;i<=A;i++)
for(int j=;j<=B;j++)
ans+=f[n][i][j];//球不必全部放入盒子中
printf("%lld",ans);
return ;
}

盒子与小球之三:有N个相同的球,M个不同的盒子,每个盒子最多放K个球 ,请计算将这N个球全部放入盒子中的方案数模1000007后的结果 ,N<=5000,M<=5000。

我们可以继续考虑动态规划,相似的状态与转移。设f[i][j]为当前放到第i个盒子,已经放了j个球的方案数。显然有转移:

f[i][j]=sigmaf[i-1][j-k],其中k枚举在当前盒子里放了几个。但由于本题数据比上一题大很多,这样的n^3转移会超时。状态设计的没有什么可以优化的了,考虑在转移优化。我们仔细观察转移方程,发现我们要加的其实是连续的一段,这就启示我们可以用前缀和优化,在转移前提前算出在i-1状态下的j前缀和。

 #include<cstdio>
#include<algorithm>
#define maxn 5090 using namespace std;
typedef long long ll; int n,m,k;
ll moder=,f[maxn][maxn],sum[maxn]; int main()
{
scanf("%d%d%d",&n,&m,&k);
for(int i=;i<=m;i++) f[i][]=;
for(int i=;i<=m;i++)
{
for(int j=;j<=n;j++)
(sum[j]=sum[j-]+f[i-][j])%=moder;
for(int j=;j<=n;j++)
{
//if(i==1&&j==1) printf("~~~~%d %d %d\n",sum[j],sum[max(j-k-1,0)],max(j-k-1,0));
if(j-k-<) (f[i][j]=sum[j]+moder)%=moder;
else (f[i][j]=sum[j]-sum[j-k-]+moder)%=moder;
//printf("%d %d:=%d\n",i,j,f[i][j]);
}
}
printf("%lld",f[m][n]);
return ;
}

盒子与小球之四:给定N个各不相同的小球,和M个不同的BOX,有多少种不同的放球方法,使得每个BOX里的小球个数不小于K。N,M,K均小于15。

我们还可以继续考虑动态规划。设f[i][j]表示当前放到第i个盒子,已经放了j个球的方案数。小球个数不小于k,只需要在转移时从k开始枚举。由于本题小球各不相同,所以我们需要组合数。本题数据在15,非常小可以递推预处理出组合数。于是有转移如下:

f[i][j]=sigmaf[i-1][j-k]*C(n-(j-k),k).

 #include<cstdio>
#include<algorithm>
#include<cstring> using namespace std;
typedef long long ll; int n,m,jue;
ll f[][],C[][]; ll Combines()
{
for(int i=;i<=;i++) C[i][]=,C[i][i]=;
for(int i=;i<=;i++)
for(int j=;j<=i;j++)
C[i][j]=C[i-][j]+C[i-][j-];
} int main()
{
Combines();
while(scanf("%d%d%d",&n,&m,&jue)!=EOF&&n!=)
{
memset(f,,sizeof f);
for(int i=jue;i<=n;i++) f[][i]=C[n][i];
for(int i=;i<=m;i++)
for(int j=jue;j<=n;j++)
for(int k=jue;k<=j;k++)
f[i][j]+=C[n-(j-k)][k]*f[i-][j-k];
printf("%lld\n",f[m][n]);
}
return ;
}

NOI题库--盒子和小球系列 By cellur925的更多相关文章

  1. NOI题库刷题日志 (贪心篇题解)

    这段时间在NOI题库上刷了刷题,来写点心得和题解 一.寻找平面上的极大点 2704:寻找平面上的极大点 总时间限制:  1000ms  内存限制:  65536kB 描述 在一个平面上,如果有两个点( ...

  2. NOI题库 1768最大子矩阵 题解

    NOI题库 1768最大子矩阵  题解     总时间限制: 1000ms 内存限制: 65536kB   描述   已知矩阵的大小定义为矩阵中所有元素的和.给定一个矩阵,你的任务是找到最大的非空(大 ...

  3. NOI题库 09:图像旋转翻转变换

    NOI题库开始的题,也是略水,当然也是大水,所以彼此彼此 09:图像旋转翻转变换 总时间限制: 1000ms 内存限制: 65536kB 描述 给定m行n列的图像各像素点灰度值,对其依次进行一系列操作 ...

  4. NOI题库-小学奥赛QwQ

    今天Loli教育我们让我们来看看NOI题库的奥赛部分,不过,为何是小学的( ⊙ o ⊙ )啊!感觉智商被各种侮辱. 余数相同问题: 描述 已知三个正整数 a,b,c. 现有一个大于1的整数x,将其作为 ...

  5. noi题库(noi.openjudge.cn) 1.7编程基础之字符串T31——T35

    T31 字符串P型编码 描述 给定一个完全由数字字符('0','1','2',-,'9')构成的字符串str,请写出str的p型编码串.例如:字符串122344111可被描述为"1个1.2个 ...

  6. NOI题库192 生日蛋糕

    192:生日蛋糕 总时间限制: 5000ms 内存限制: 65536kB 描述 7月17日是Mr.W的生日,ACM-THU为此要制作一个体积为Nπ的M层生日蛋糕,每层都是一个圆柱体. 设从下往上数第i ...

  7. NOI 题库 9272 题解

    9272   偶数个数字3 描述 在所有的N位数中,有多少个数中有偶数个数字3? 输入 一行给出数字N,N<=1000 输出 如题 样例输入 2 样例输出 73 Solution : 令f ( ...

  8. noi题库(noi.openjudge.cn) 1.5编程基础之循环控制T36——T45

    T36 计算多项式的值 描述 假定多项式的形式为xn+xn-1+-+x2+x+1,请计算给定单精度浮点数x和正整数n值的情况下这个多项式的值. 输入 输入仅一行,包括x和n,用单个空格隔开.x在flo ...

  9. noi题库(noi.openjudge.cn) 1.7编程基础之字符串T21——T30

    T21:单词替换 描述 输入一个字符串,以回车结束(字符串长度<=100).该字符串由若干个单词组成,单词之间用一个空格隔开,所有单词区分大小写.现需要将其中的某个单词替换成另一个单词,并输出替 ...

随机推荐

  1. 2017-10-04-morning

    改题面只有1改为0 .. #include <cstring> #include <cstdio> inline void read(int &x) { x=; reg ...

  2. 洛谷—— P2812 校园网络

    P2812 校园网络 题目背景 浙江省的几所OI强校的神犇发明了一种人工智能,可以AC任何题目,所以他们决定建立一个网络来共享这个软件.但是由于他们脑力劳动过多导致全身无力身体被♂掏♂空,他们来找你帮 ...

  3. Date日期模式

    package cn.zmh.Date; import java.text.SimpleDateFormat; import java.util.Date; public class DateDemo ...

  4. jquery显示和隐藏元素

    1.$('#id').show()/$('#id').hide()/$('#id').toggle() 2.$('#id').css('display','none')/$('#id').css('d ...

  5. jenkins的代理设置,方便下载插件

    jenkins在下载插件的时候,总是网络不通,需要设置代理跨越长城 java.net.SocketTimeoutException: connect timed out Caused: java.ne ...

  6. 2003 -Can't connection to mysql server on | navicat for mysql Access denied for user 'root'@''ip'(using password :yes)

    用本机windows上的Navicat for mysql链接虚拟机Linux的mysql数据库时,第一次连接的时候报的错误是 2003 -Can't connection to mysql serv ...

  7. POJ2752 Seek the Name, Seek the Fame 【KMP】

    Seek the Name, Seek the Fame Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 11602   Ac ...

  8. [Tools] Create a Chrome Extension

    Creating a Chrome extension requires a manifest.json file which defines how your extension will beha ...

  9. oracle统计信息

    手工刷ORACLE统计信息  select count(1) from LOG_TRX_DETAIL;  select * from user_tab_statistics where table_n ...

  10. CxImage的编译及简单使用举例

    1.  从http://sourceforge.net/projects/cximage/下载最新的CxImage 702源代码. 2.  解压缩后,以管理员身份打开CxImageFull_vc10. ...