【JZOJ4710】【NOIP2016提高A组模拟8.17】Value
题目描述
输入
输出
样例输入
5
8 2
10 7
5 1
11 8
13 3
样例输出
27
数据范围
解法
选定一些物品a[1],a[2],a[3]…a[num],尝试交换a[i],a[j],那么对答案的贡献是:
-(v[i]-xsum)-(v[j]-w[i]-xsum)+(v[j]-xsum)+(v[i]-w[j]-xsum) (xsum为sigma(w[1..i-1]));
整理后可得+w[i]-w[j],所以说如果w[i]>w[j],那么我们交换a[i],a[j]是对答案有积极贡献的。
所以我们以w为关键字从小到大排序,然后就变成了经典的背包问题。
但是直接套上背包会超时,只可以拿40%的分。
接下来会有两个选择支:
A.优化背包,去除冗余状态
排序后,经典背包是这样的f[j+w[i]]=max{f[j]+v[i]-j};
如果想让f[j]对f[j+w[i]]做出贡献,那么就必须保证v[i]-j>0,所以所以枚举j时不能超过v[i]。
B.更改动态规划方法(标准解法)
排序后,设f[i][j]表示倒序动态规划到了第i个物品,已经选了j个物品。
那么f[i][j]=max{f[i+1][j],f[i+1][j-1]+v[i]-j*w[i]}
这样就利用了选择一个物品只会对这个物品后的其他物品造成影响的性质。
时间复杂度为O(n^2)。
代码
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define ll long long
#define sqr(x) ((x)*(x))
#define ln(x,y) int(log(x)/log(y))
using namespace std;
const char* fin="aP2.in";
const char* fout="aP2.out";
const int inf=0x7fffffff;
const int maxn=5007,maxm=100007;
int n,i,j,k,maxa,ans;
struct node{
int v,w;
}a[maxn];
int f[maxn][maxn];
bool cmp(node a,node b){
return a.w<b.w;
}
int main(){
scanf("%d",&n);
for (i=1;i<=n;i++){
scanf("%d%d",&a[i].v,&a[i].w);
maxa=max(maxa,a[i].v);
}
sort(a+1,a+n+1,cmp);
for (i=n;i>0;i--){
for (j=1;j<=n;j++){
f[i][j]=max(f[i+1][j],f[i+1][j-1]-a[i].w*(j-1)+a[i].v);
ans=max(ans,f[i][j]);
}
}
printf("%d",ans);
return 0;
}
启发
1.尝试优化动态规划的冗余状态,剔除不必要的转移。
2.挖掘题目性质,顺着性质设置状态。
【JZOJ4710】【NOIP2016提高A组模拟8.17】Value的更多相关文章
- 【JZOJ4787】【NOIP2016提高A组模拟9.17】数格子
题目描述 输入 输出 样例输入 1 10000 3 10000 5 10000 0 0 样例输出 1 11 95 数据范围 每个测试点数据组数不超过10组 解法 状态压缩动态规划. 设f[i][j]表 ...
- 【NOIP2016提高A组模拟9.17】序列
题目 分析 首先用\(a_i\)表示达到目标的步数\(B_i-A_i(mod 4)\) 根据粉刷栅栏,先不管mod 4的情况,答案就是\(\sum\max(a_i-a_{i+1},0)\) 那我们刚才 ...
- NOIP2016提高A组模拟9.17总结
第一题,典型的隔板问题, 但是我忘记隔板问题怎么打,一开始在花了1小时,还是没想出来,果断弃疗, 最后的40分钟,我打完了第二题,接着又用了20分钟推敲出一种极其猥琐的式子来代替,可惜预处理的阶乘忘记 ...
- 【NOIP2016提高A组模拟9.17】数格子
题目 分析 设表示每一行的状态,用一个4位的二进制来表示,当前这一行中的每一个位数对下一位有没有影响. 设\(f_{i,s}\)表示,做完了的i行,其状态为s,的方案数. 两个状态之间是否可以转移就留 ...
- 【NOIP2016提高A组模拟9.17】小a的强迫症
题目 分析 题目要求第i种颜色的最后一个珠子要在第i+1种颜色的最后一个珠子之前, 那么我们从小到大枚举做到第i种,把第i种的最后一颗珠子取出,将剩下的\(num(i)-1\)个珠子插入已排好的前i- ...
- 【NOIP2016提高A组模拟8.17】(雅礼联考day1)总结
考的还ok,暴力分很多,但有点意外的错误. 第一题找规律的题目,推了好久.100分 第二题dp,没想到. 第三题树状数组.比赛上打了个分段,准备拿60分,因为时间不够,没有对拍,其中有分段的20分莫名 ...
- 【NOIP2016提高A组模拟8.17】(雅礼联考day1)Binary
题目 分析 首先每个数对\(2^i\)取模.也就是把每个数的第i位以后删去. 把它们放进树状数组里面. 那么当查询操作, 答案就位于区间\([2^i-x,2^{i-1}-1-x]\)中,直接查询就可以 ...
- 【NOIP2016提高A组模拟8.17】(雅礼联考day1)Value
题目 分析 易证,最优的答案一定是按\(w_i\)从小到大放. 我们考虑dp, 先将w从小到大排个序,再设\(f_{i,j}\)表示当前做到第i个物品,已选择了j个物品的最大值.转移就是\[f_{i, ...
- 【NOIP2016提高A组模拟8.17】(雅礼联考day1)Matrix
题目 分析 假设,我们从\(F_{i,2}\)出发,那么对\(F_{n,n}\)的贡献就是\(某个系数乘以a^{n-i}b^{n-1}r_i\): 同理,如果从\(F_{2,i}\)出发,那么对\(F ...
随机推荐
- typeof 、Object.prototype.toString和 instanceof
数据类型 js 基本类型包括:Undefined symbol null string boolean number js 引用类型包括:object array Date RegExp typeo ...
- np一些基本操作2
import numpy as nparr1 = np.arange(32).reshape(8,4)print(arr1)arr1 = arr1.reshape(-1);print(arr1)arr ...
- iPadOS 更新日志 - 持续更新中
本文只是为了简单记录一下每个正式版本发布时间和更新内容,只有这个初衷,从2019年9月25日开始,将会持续更新. iPadOS 13.1 - 2019年9月25日 经全新命名的 iPadOS 是一款强 ...
- HBase+Redis
- mysqlbinlog: unknown variable 'default-character-set=utf8'
[xxx@dbhost log]$ mysqlbinlog mysql-bin. mysqlbinlog: unknown variable 'default-character-set=utf8' ...
- Python学习之enumerate
enumerate还可以接收第二个参数,用于指定索引起始值 2. 注意open返回文件对象,可迭代,而os.open返回的是文件指针,int类型, <wiz_tmp_tag ...
- git commit规范工具
npm install -g commitizen commitizen init cz-conventional-changelog --save --save-exact 以后,凡是用到git c ...
- 左神算法书籍《程序员代码面试指南》——1_08构造数组的MaxTree
[题目] 将一个没有重复数字的数组中的数据构造一个二叉树 每个节点都是该子树的最大值 [要求] 时间复杂度为O(N)[题解] 使用单调栈,栈的顺序是维持从大到小排序 通过使用单调栈,将数组中中所有数的 ...
- Python eval()函数的用法
Python eval()函数的用法 eval(str)函数很强大,官方解释为:将字符串str当成有效的表达式来求值并返回计算结果.所以,结合math当成一个计算器很好用. eval()函数常见作用有 ...
- Codeforces Round #466 (Div. 2) 题解940A 940B 940C 940D 940E 940F
Codeforces Round #466 (Div. 2) 题解 A.Points on the line 题目大意: 给你一个数列,定义数列的权值为最大值减去最小值,问最少删除几个数,使得数列的权 ...