【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 ...
随机推荐
- vue 模版组件用法
第一种 //首先,别忘了引入vue.js <div id="user_name_01"></div> <script src="../nod ...
- LeetCode387First Unique Character in a String字符串中第一个唯一字符
给定一个字符串,找到它的第一个不重复的字符,并返回它的索引.如果不存在,则返回 -1. 案例: s = "leetcode" 返回 0. s = "loveleetcod ...
- ubuntu 安装samba共享文件夹
安装samba sudo apt-get install samba smbclient 配置samba sudo cp /etc/samba/smb.conf /etc/samba/smb.conf ...
- Python之MySQL语法(增删改查)
-- ID: 新闻的唯一标示 -- title: 新闻的标题 -- content: 新闻的内容 -- created_at: 新闻添加的时间 -- types: 新闻的类型 -- image: 新的 ...
- Vue实例属性/方法/生命周期
一 实例属性 var app = new Vue(); app.$el ===>返回当前挂在DOM对象. app.$data ===>返回当前实例data里面的数据 app.$option ...
- ArcGIS Data Interoperability 的使用(1)
今天在用OneMap的时候,发现OneMap中注册过后的WFS服务无法在skyline中加载,于是想知道OneMap注册后的WFS服务与server中的原生态WFS服务有啥区别.首先想到是否能在Arc ...
- Python技巧—list与字符串互相转换
Python技巧-list与字符串互相转换 在Python的编程中,经常会涉及到字符串与list之间的转换问题,下面就将两者之间的转换做一个梳理. 1.list转换成字符串 命令:list() 例子: ...
- 挑逗 Java 程序员的那些 Scala 绝技
有个问题一直困扰着 Scala 社区,为什么一些 Java 开发者将 Scala 捧到了天上,认为它是来自上帝之吻的完美语言:而另外一些 Java 开发者却对它望而却步,认为它过于复杂而难以理解.同样 ...
- python 当文件目录不存在时,如何自动创建
import os if not os.path.exists('foldername'): os.mkdir('foldername')
- 在银行业中,BP是指什么?
基点 Basis Point(BP)债券和票据利率改变量的度量单位.一个基点等于0.01个百分点,即0.01%,因此,100个基点等于1%.[例]一浮动利率债券的利率可能比LIBOR高10个基点,10 ...