小学生放假了

总时限 26s 内存限制 256MB
出题人 zsyzzsoft 提交情况 16/150
初始分值 1500 锁定情况

背景

我们能见到的最可怕的事情,莫过于小学生放假了!

描述

小学生要放假了!MT学校一共有N个小学生,学校旁边的ET小卖部希望在小学生放假之前做好坑蒙小学生的准备!ET小卖部一共有M个不同的商品,每个商品的价格可以定位任意非负整数,每个商品的数量是无限的。每个小学生有Ci RMB,每人只能购买一个商品,他们希望他们购买的商品尽量贵。小卖部应该如何设定每个商品的价格,使得他们坑蒙小学生的收入尽可能多呢?请输出最多的收入。

输入格式

第一行两个用空格隔开的整数N,M。

紧接着N行,第i+1行一个整数,表示Ci(见题目描述)

输出格式

一个整数,表示最多的收入。

样例输入

5 3
1
3
5
7
9

样例输出

22

样例解释

三个商品的价格分别设置为3RMB,7RMB和9RMB。

第一个小学生由于没有足够的RMB,不购买任何商品;

第二个小学生和第三个小学生只能购买3RMB的商品;

第四个小学生可以购买7RMB的商品;

第五个小学生可以购买9RMB的商品。

3 + 3 + 7 + 9 = 22,所以这种方案获得了22RMB的收入。

可以证明,没有更优的方案。

数据范围与约定

对于100%的数据,1 <= Ci <= 109,1 <= N <= 10000,1 <= M <= 2000。

单点时间限制2s。

来源

原创

斜率优化易错点。。。

1.不等式*(-1),2边都要变号。

2.注意叉积乘爆,切忌bool强转 (bool) t >0

优先级高

3.初值P(0,a[1],0)

4.eps乱用

5.a.y-b.y

a.x=b.x的情况  *inf 注意+,-

否则long double 可能出现真的正负无穷   然后T得惨惨的,还查不出错。。。

能出正解的其实,上面。。。↑

写了一天各种错。。无语了。。后天考试。。考文化。。。坐等爆0。。。。

顺便说一下clock()的用法

clock()/CLOCKS_PER_SEC 返回程序运行开始到执行的时间(单位:s)

卡时间专用。。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define MAXN (10000+10)
#define MAXM (2000+10)
#define eps 1e-13
#define Read(x) { \
while (!isdigit(c=getchar())); \
x=c-48; \
while (isdigit(c=getchar())) x=x*10+c-48; \
}
long long mul(long long a,long long b){return (a*b)%F;}
long long add(long long a,long long b){return (a+b)%F;}
long long sub(long long a,long long b){return (a-b+(a-b)/F*F+F)%F;}
typedef long long ll;
int n,m;
char c;
ll a[MAXN];
ll f[MAXM][MAXN]={0},T[MAXN]={0},h[MAXN];
struct P
{
int i;
long double x,y;
P(int _i,ll _x,ll _y):i(_i),x(_x),y(_y){}
P(ll _x,ll _y):x(_x),y(_y){}
P(){}
friend long double kk(P a,P b){if (abs(a.x-b.x)<eps) return (b.y-a.y)*INF;return (b.y-a.y)/(b.x-a.x); }
}st[MAXN];
struct V
{
long double x,y;
V(ll _x,ll _y):x(_x),y(_y){}
V(){}
V(P a,P b):x(b.x-a.x),y(b.y-a.y){}
//friend bool operator*(V a,V b){return ((long double)a.x/a.y-(long double)b.x/b.y)>eps;}
friend long double operator*(V a,V b){/*cout<<a.x*b.y-a.y*b.x<<endl;*/return a.x*b.y-a.y*b.x; }
};
int main()
{
// freopen("input","r",stdin);
Read(n) Read(m)
For(i,n)
Read(a[i])
sort(a+1,a+1+n);
For(i,n-1) h[i]=(ll)i*a[i+1];
// For(j,n) f[0][j]=0;
For(i,m)
{
For(j,n-1) T[j]=f[i-1][j]-h[j];
int head=1,tail=1;
st[1]=P(0,a[1],0);
For(j,n)
{
// cout<<head<<' '<<tail<<endl;
while (head^tail&&kk(st[head],st[head+1])>=-j) head++;
// cout<<kk(st[head],st[head+1])<<' '<<-j<<endl;
int k=st[head].i;
f[i][j]=f[i-1][k]+(j-k)*a[k+1];
// cout<<i<<' '<<j<<':'<<k<<' '<<f[i][j]<<' '<<f[i-1][k]<<' '<<(j-k)*a[k+1]<<endl;
// if (j<n)
// {
P A=P(j,a[j+1],T[j]);
while (head^tail&&V(st[tail-1],st[tail])*V(st[tail],A)>=0) tail--;
st[++tail]=A;
// }
}
}
//ll ans=0;
/*
For(i,m)
{
For(j,n) ans=max(ans,f[i][j]);//,cout<<f[i][j]<<' ';cout<<endl;
}*/
// cout<<ans<<endl;
printf("%lld\n",f[m][n]);
//For (j,n) cout<<f[m][j]<<endl;
// cout<<clock()/CLOCKS_PER_SEC<<endl;
return 0;
}

CH BR8(小学生放假了-clock()/CLOCKS_PER_SEC-斜率优化常错集锦)的更多相关文章

  1. CH BR8(小学生在上课-逆元和互质数一一对应关系)

    小学生在上课 总时限 11s 内存限制 256MB 出题人 jzc 提交情况 66/277 初始分值 600 锁定情况 背景 小学生在学校上数学课…… 描述 数学课上,小学生刚学会了乘除法.老师问了他 ...

  2. CF1179D Fedor Runs for President [DP,斜率优化]

    Codeforces 思路 考虑把连的那两个点中间的链提出来,那么就会变成一条链,链上的每个点挂着一棵子树的形式. 设那些子树的大小为\(S_1,S2,\cdots\),那么新加的简单路径个数就是 \ ...

  3. 【BZOJ-1096】仓库建设 斜率优化DP

    1096: [ZJOI2007]仓库建设 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 3719  Solved: 1633[Submit][Stat ...

  4. 【BZOJ-4518】征途 DP + 斜率优化

    4518: [Sdoi2016]征途 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 230  Solved: 156[Submit][Status][ ...

  5. 【BZOJ-3675】序列分割 DP + 斜率优化

    3675: [Apio2014]序列分割 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 1420  Solved: 583[Submit][Statu ...

  6. 【BZOJ-3437】小P的牧场 DP + 斜率优化

    3437: 小P的牧场 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 705  Solved: 404[Submit][Status][Discuss ...

  7. 【BZOJ-3156】防御准备 DP + 斜率优化

    3156: 防御准备 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 951  Solved: 446[Submit][Status][Discuss] ...

  8. 【BZOJ-1010】玩具装箱toy DP + 斜率优化

    1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 8432  Solved: 3338[Submit][St ...

  9. 【BZOJ-1597】土地购买 DP + 斜率优化

    1597: [Usaco2008 Mar]土地购买 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 2931  Solved: 1091[Submit] ...

随机推荐

  1. Eclipse 导入项目乱码问题(中文乱码)

    1.编码不对 a.对某文件或某工程更改编码:   鼠标移到工程名或文件名,右键->Properties->Resource->Text file enCoding ->更改编码 ...

  2. VM下Linux网卡丢失(pcnet32 device eth0 does not seem to be ...)解决方案

    系统启动日志:Bringing up interface eth0: pcnet32 device eth0 does not seepresent, delaying initialization. ...

  3. Oracle中包的创建

    包是过程和函数的集合体,包包括创建包和创建包体,创建包的时候在可以定义过程和函数,包体中则具体实现过程和函数. eg: --创建包 create  or replace package mypac1 ...

  4. WebBasic-表单

    用来提交数据 <form></form> 属性:action:提交的url  method:表单数据提交的方式  enctype:表单数据的编码方式 表单控件 --input元 ...

  5. 武汉科技大学ACM:1005: 华科版C语言程序设计教程(第二版)例题5.8

    Problem Description 老师给小豪出了一道题目:给你两个整数x和n(-10<=x<=10,1<=n<=10),让你求出x^1+x^2+x^3+……+x^n的结果 ...

  6. STM8S和STM8L调试串口中断的注意点

    1. STM8L串口中断注意点 在调试PM2.5传感器GP2Y1051的时候,发现在仿真的时候开始能够进行数据的接受,但是如果暂停之后就不能接受数据,其实只是接收了一次完整的数据. 问题程序 解决方法 ...

  7. jquery获取浏览器高度、宽度和滚动条高度(来自网络)

    Jquery代码: alert($(window).height()); //浏览器时下窗口可视区域高度 alert($(document).height()); //浏览器时下窗口文档的高度 ale ...

  8. PHPCMS v9 在windows2008系统 IIS7 下设置伪静态的方法

    安装环境:windows2008+IIS7.0+PHP5+MYSQL5 一.安装phpcms v9程序,设置伪静态.如图: 二.安装IIS7官方Url重写模块 1.先到IIS官方下载模块 下载地址:h ...

  9. 1109HTML学习

    <div><!--face里面用逗号隔开表示 字体优先选择.size是字体1到7 --><font color="red" face="微软 ...

  10. arm linux上的第一个应用程序 BOA移植

    1. 首先, linux在开发板上能跑起来了. 包括网络驱动也有了, ifconfig之后, 能看到在rcS里面设置的IP, 也能ping通windows主机了, 当然, 也要window关掉防火墙才 ...