博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ3874:[AHOI2014&JSOI2014]宅男计划(爬山法)
阅读量:6564 次
发布时间:2019-06-24

本文共 2422 字,大约阅读时间需要 8 分钟。

Description

 【故事背景】
自从迷上了拼图,JYY就变成了个彻底的宅男。为了解决温饱问题,JYY
不得不依靠叫外卖来维持生计。
【问题描述】
外卖店一共有N种食物,分别有1到N编号。第i种食物有固定的价钱Pi和保质期Si。第i种食物会在Si天后过期。JYY是不会吃过期食物的。
比如JYY如果今天点了一份保质期为1天的食物,那么JYY必须在今天或
者明天把这个食物吃掉,否则这个食物就再也不能吃了。保质期可以为0天,这
样这份食物就必须在购买当天吃掉。
JYY现在有M块钱,每一次叫外卖需要额外付给送外卖小哥外送费F元。
送外卖的小哥身强力壮,可以瞬间给JYY带来任意多份食物。JYY想知道,在
满足每天都能吃到至少一顿没过期的外卖的情况下,他可以最多宅多少天呢?

Input

第一行包含三个整数M,F和N。
接下来N行,第i行包含两个整数Pi和Si。

Output

输出仅包含一行一个整数表示JYY可以宅的最多的天数。

Sample Input

32 5 2
5 0
10 2

Sample Output

3

HINT

 【样例说明】

JYY的最佳策略是:
第一天买一份食物1和一份食物2并且吃一份食物1;
第二天吃一份食物2;
第三天买一份食物1并且吃掉。
【数据规模与约定】
对于100%的数据满足0<=Si<=10^18,1<=F,Pi,M<=10^18,1<=N<=200

Solution

调了一晚上真的难受……
首先我们可以先将食物sort一下然后去重,使之成为一个价格递增且保质期递增的数列
被我们从食物中舍弃的肯定都是原序列中有比它更优的
然后我们可以发现,对于确定的点外卖次数,我们可以通过贪心法来求更多生存天数
具体贪心方法是对于我们要订的t次外卖,如果食物i(1<=i<=n)能订的话那显然比i+1更优,
不超保质期的情况下优先用i把t次订单填满。
把订外卖次数当成x轴,然后就可以用模拟退火乱搞了
然而搞来搞去只有九十,我们发现答案始终比正确答案小一点怎么办呢?
由于购买次数少的时候买的东西贵,购买次数多的时候配送费贵,
我们可以大胆猜测结果单峰
然后把源程序注释掉一点就能过了

Code

1 #include
2 #include
3 #include
4 #include
5 #include
6 #include
7 #include
8 using namespace std; 9 10 struct Node{ long long cost,life;}Food[205];11 bool cmp(Node a,Node b){ return a.cost
b.life;}12 13 long long M,F,N,cnt,ans;14 15 double Rand(){ return rand()%10001/10000.0;}16 17 long long Calc(long long t)18 {19 long long Left=M-t*F,sum=0;20 for (int i=1; i<=N; ++i)21 {22 if (Left*1.0/Food[i].cost/t>=(Food[i].life-Food[i-1].life+(i==1)))23 {24 Left-=Food[i].cost*(Food[i].life-Food[i-1].life+(i==1))*t;25 sum+=(Food[i].life-Food[i-1].life+(i==1))*t;26 continue;27 }28 else29 {30 sum+=Left/Food[i].cost;31 break;32 }33 }34 if (sum>ans) ans=sum;35 return sum;36 }37 38 void Simulate_Anneal()39 {40 double T=M/F+1000;41 long long now=1;42 while (T>=1)43 {44 long long nxt=now+T*(2*Rand()-1);45 while (nxt<=0) nxt=T*Rand();46 nxt%=(M/F+1);47 48 long long delta=Calc(nxt)-Calc(now);49 if (delta>0/*|| exp(delta/T)>Rand()*/)50 now=nxt;51 T*=0.93;52 }53 }54 55 int main()56 {57 scanf("%lld%lld%lld",&M,&F,&N);58 for (int i=1; i<=N; ++i)59 scanf("%lld%lld",&Food[i].cost,&Food[i].life);60 sort(Food+1,Food+N+1,cmp);61 for (int i=2; i<=N; ++i)62 if (Food[i].life<=Food[i-1].life)63 {64 for (int j=i; j<=N; ++j)Food[j]=Food[j+1];65 N--; i--;66 }67 for (int i=1; i<=20; ++i)68 Simulate_Anneal();69 printf("%lld",ans);70 }

转载于:https://www.cnblogs.com/refun/p/9286962.html

你可能感兴趣的文章
android 实现ListView中添加RaidoButton单选
查看>>
WS-Security 中文问题&Stax(Streaming API for XML) (二)
查看>>
dos 分页显示及查看应用程序占用端口
查看>>
Oracle数据库:启动操作
查看>>
限制VLAN之间互访实例
查看>>
关于selenium中断言判断url获取错误解决
查看>>
Ubuntu12下挂载硬盘(9TB)
查看>>
好用的PHP分页类
查看>>
linux下的防火墙
查看>>
简练软考知识点整理-创建工作分解结构过程
查看>>
NVisionXR_iOS教程一 —— NVisionXR从零搭建一个AR项目
查看>>
oracle 12c ins-30131 执行安装程序验证所需的初始设置失败
查看>>
SNAT与DNAT
查看>>
BGP十三条规则
查看>>
Linux 修改密码“ Authentication token manipulation err”
查看>>
openstack
查看>>
【顶】(与同事合作的快乐)技术人员也需要先学会做人,再学会做事,再是能成事,最后是成名得利...
查看>>
Lync Server 2013 安装体验(一)
查看>>
Hadoop2.6.0学习笔记(五)自定义InputFormat和RecordReader
查看>>
EBB-24、DNS2
查看>>