【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 列中每个数 ...
随机推荐
- caffe中根据 *_train_test.prototxt文件生成 *_deploy.prototxt文件 (转载)
见博客:http://blog.csdn.net/u010417185/article/details/52137825
- 基于Laravel开发博客应用系列 —— 构建博客后台管理系统
一个完整的博客应用不能没有后台管理系统.所以在本节中我们将继续完善博客应用 —— 开发后台管理系统. 1.创建路由 在上一节十分钟创建博客项目中,已经设置过了 app/Http/routes.php, ...
- PHP 博客收集
https://lvwenhan.com/ www.chrisyue.com https://silex.symfony.com/ https://www.chrisyue.com/translati ...
- QString::arg()//用字符串变量参数依次替代字符串中最小数值
QString i = "iTest"; // current file's number QString total = "totalTest&qu ...
- @react-native-community/async-storage在Android上的手动link问题
PS C:\Users\linjin\Desktop\RN_APP> react-native link @react-native-community/async-storage error ...
- POJ - 2456 Aggressive cows 二分 最大化最小值
Aggressive cows Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 18099 Accepted: 8619 ...
- 层级目录结构的Makefile递归编译方法
层级目录结构的Makefile编写方法. 层级目录结构的Makefile编写方法. 0.前言 1.如何编译整个工程 2.过滤每层不需要编译的目录 3将所有输出文件定向输出. 0.前言 假如现在有这样一 ...
- CF582A GCD Table
A. GCD Table time limit per test 2 seconds memory limit per test 256 megabytes input standard input ...
- 选择 React Native 的理由
转载:选择 React Native 的理由 从开始知道 React Native 到现在已经过了5个月,真实的试用也经历了三个月的时间.阅读文档开始,了解是什么,到简单的理解为什么,都是在聆听不同的 ...
- JS 常用库汇总收集
本文不定期更新, 用于汇总记录一些看着 ok 的 JS 库. 库名 简介 项目地址 macy.js 仅 4 kb的 原生 流布局插件 http://macyjs.com/ Driver.js 仅 4 ...