【BZOJ】【3612】【HEOI 2014】平衡
DP
唉我还是too naive
这是个整数划分题……
我想的DP方式是f[i][j][k]表示前 i 个数拼出 j 用了 k 个数的方案数……
转移当然是比较直观……
但是只能得30分QAQ
正确的DP姿势:http://blog.csdn.net/Vmurder/article/details/42551603
分析:
数据范围不大,我们可以写整数划分。
f[i][j]表示将i划分成j个互不相同的正整数,且最大不超过n 的划分方案数。
这里说一下这道题的整数划分。
我们不妨先来反向思考一下。---------------------------------------------------------------------
首先考虑f[i][j](下图每一列都代表一个数,高度就是数值)
它可以在底下添加一行,进行转移
分为两种情况:
Ⅰ. 转移过后最小数不为1
Ⅱ. 转移过后最小数为1
然后f[i][j]可以向这两个方向转移。
而我们还需要保证最大数不能大于n,那么如下图
在最后加入一层,使得当前所有整数都+1时,发现出现了一个101,而不妨当作n是100
那么显然我们可以很方便地清掉这个数带来的情况。
即当前是f[i][j],那么最后一列是大于n的情况显然只能是有一个整数n+1
不看这个n+1,情况数是f[i-(n+1)][j-1] ,我们把这个情况集删掉就好了。
正向考虑:-----------------------------------------------------------------------------------
首先不妨把刚才的图片按顺序记作图1、2、3、4。
我们把f[i][j](图1)这么多方案分成两种情况:
Ⅰ. 最小的数不为1:
好说。 直接由f[i-j][j]在底下加一行得到。就是图2。
此时原来划分出来的整数不同,新的这些整数显然依然不同。
Ⅱ. 最小的数为1:
那么显然它可以由f[i-j][j-1]转移得到,
即在f[i-j][j-1]代表图形下面整体+1,最后加上一个整数1,即图3。
注意此时f[i-j][j-1]代表的所有图形整数都不同(性质/定义),那么新加1后所有整数依然不同,且均>=2
这个时候再来个整数1,依然满足所有整数不同。
而这两种情况显然互补,即这两种情况的转移包含了f[i][j]的所有情况(两个命题“最小数是1”,“最小数不是1”,显然包含全部情况),也就是说转移完成。
但是我们注意到还需要让最大数不能超过n,
所以有了图4。
也就是我们要减去最大数超过n的情况,方法前文图下有说明。
这道题难点解决了。
现在说一下其它细节:
f[i][j]算出来后直接暴力枚举两边的权值,及用点个数(不要忘了中心支点)
然后check。end。
/**************************************************************
Problem: 3612
User: Tunix
Language: C++
Result: Accepted
Time:3260 ms
Memory:15336 kb
****************************************************************/ //Huce #2 B
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
using namespace std; int getint(){
int v=,sign=; char ch=getchar();
while(ch<''||ch>'') {if (ch=='-') sign=-; ch=getchar();}
while(ch>=''&&ch<='') {v=v*+ch-''; ch=getchar();}
return v*sign;
}
typedef long long LL;
const int N=,INF=~0u>>;
/*******************tamplate********************/
LL f[N][];
int n,k,P,sum[],mx[N][]; int main(){
#ifndef ONLINE_JUDGE
freopen("B.in","r",stdin);
freopen("B.out","w",stdout);
#endif
int T=getint();
f[][]=;
while(T--){
n=getint(); k=getint(); P=getint();
if (k==){puts("");continue;}
int w=n*(k-);
F(i,,w)
F(j,,k-){
f[i][j]=i>=j ? (f[i-j][j]+f[i-j][j-])%P : ;
f[i][j]=i>=n+ ? (f[i][j]-f[i-(n+)][j-]+P)%P : f[i][j];
}
LL ans=;
F(i,,w) F(j,,k-)
(ans+=f[i][j]*f[i][k-j]+f[i][j]*f[i][k-j-])%=P;
printf("%lld\n",ans);
}
return ;
}
3612: [Heoi2014]平衡
Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 110 Solved: 72
[Submit][Status][Discuss]
Description
Input
第一行,一个正整数,表示数据组数 T(萱萱向你询问的次数)。
Output
共 T 行,每行一个正整数,代表你得出的对应问题的答案。
Sample Input
6 5 10000
4 1 10000
9 6 10000
4 6 10000
5 1 10000
8318 10 9973
9862 9 9973
8234 9 9973
9424 9 9973
9324 9 9973
Sample Output
1
920
8
1
4421
2565
0
446
2549
HINT
T <= 20,1 <= n <= 10000,1 <= k <= 10,2 <= p <= 10000,且 k <= 2n+1。
Source
【BZOJ】【3612】【HEOI 2014】平衡的更多相关文章
- 洛谷 P4106 / bzoj 3614 [ HEOI 2014 ] 逻辑翻译 —— 思路+递归
题目:https://www.luogu.org/problemnew/show/P4106 https://www.lydsy.com/JudgeOnline/problem.php?id=3614 ...
- [BZOJ 3167][HEOI 2013]SAO
[BZOJ 3167][HEOI 2013]SAO 题意 对一个长度为 \(n\) 的排列作出 \(n-1\) 种限制, 每种限制形如 "\(x\) 在 \(y\) 之前" 或 & ...
- BZOJ 3612: [Heoi2014]平衡( dp )
枚举Fl, 就变成一个整数划分的问题了...f(i,j) = f(i-j,j-1)+f(i-j,j)-f(i-N-1,j-1)递推.f(i,j)表示数i由j个不同的数组成,且最大不超过N的方案数 -- ...
- BZOJ 3612: [Heoi2014]平衡
3612: [Heoi2014]平衡 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 283 Solved: 219[Submit][Status][ ...
- bzoj 3612 [Heoi2014]平衡——整数划分(dp)
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3612 因为力矩的缘故,变成了整数划分. 学习到了整数划分.就是那个图一样的套路.https: ...
- bzoj 3612: [Heoi2014]平衡【整数划分dp】
其实就是-n~n中求选k个不同的数,和为0的方案数 学到了新姿势叫整数划分,具体实现是dp 详见:https://blog.csdn.net/Vmurder/article/details/42551 ...
- BZOJ 3672 NOI 2014 购票
题面 Description 今年夏天,NOI在SZ市迎来了她30周岁的生日.来自全国 n 个城市的OIer们都会从各地出发,到SZ市参加这次盛会. 全国的城市构成了一棵以SZ市为根的有根树,每个城市 ...
- [BZOJ 3530][Sdoi 2014]数数
阿拉~好像最近总是做到 AC 自动机的题目呢喵~ 题目的算法似乎马上就能猜到的样子…… AC 自动机 + 数位 dp 先暴力转移出 f[i][j] :表示从 AC 自动机上第 j 号节点走 i 步且不 ...
- [BZOJ 3503][Cqoi 2014]和谐矩阵
我觉得这一题的样例输出一点都不和谐,大家千万别像我一样被坑了…… 题目不算难,果然是进错省系列555,不过搞出 O(n*m*2m) 的还是不要挣扎的比较好 我们暴力地推出第 n 行 第 m 列中每个数 ...
随机推荐
- Loadrunner乱码问题
在LoadRunner中录制脚本时,出现乱码的问题解决 我在录制一个Web的脚本时,出现中文乱码. 原因为Web中采用的是UTF-8编码,而录制脚本的选项默认没有把支持UTF8选中. 方法:1. To ...
- NetworkX 使用(二)
官方教程 博客:NetworkX %pylab inline import networkx as nx Populating the interactive namespace from numpy ...
- Java中面向对象的理解
按照惯例,先做一个简单的介绍,现在开始学习 Thinging in Java 4 ,一边看,一边记录,我都不想给自己设定时间安排了,毕竟很少实现过.所以就这样吧!不定期的更新,我都会放到博客中的. 所 ...
- php 安全模式限制函数
表 42-2. 安全模式限制函数 函数名 限制 dbmopen() 检查被操作的文件或目录是否与正在执行的脚本有相同的 UID(所有者). dbase_open() 检查被操作的文件或目录是否与正在执 ...
- PHP渗透中的奇淫技巧--检查相等时的漏洞
PHP是现在网站中最为常用的后端语言之一,是一种类型系统 动态.弱类型的面向对象式编程语言.可以嵌入HTML文本中,是目前最流行的web后端语言之一,并且可以和Web Server 如apache和n ...
- hdu 4403 爆搜
题意:给一串数字,在其间加入若干加号和一个等号,问使等式成立的方案总数 if the digits serial is "1212", you can get 2 equation ...
- Codeforces Round #349 (Div. 1) A. Reberland Linguistics 动态规划
A. Reberland Linguistics 题目连接: http://www.codeforces.com/contest/666/problem/A Description First-rat ...
- Codeforces Round #298 (Div. 2) E. Berland Local Positioning System 构造
E. Berland Local Positioning System Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://codeforces.c ...
- BZOJ 1196: [HNOI2006]公路修建问题 Kruskal/二分
1196: [HNOI2006]公路修建问题 Time Limit: 1 Sec Memory Limit: 162 MB 题目连接 http://www.lydsy.com/JudgeOnline ...
- 采用LinkedList来模拟栈数据结构的集合--先进后出
三.用LinkedList来模拟栈数据结构的集合 /* * 自定义一个数据结构为LinkedList的集合类*/public class MyCollection_LinkedList { publi ...