bzoj4306: 玩具厂
Description
Input
Output
Sample Input
50
25
25
50
1
2
3
2
1
Sample Output
0.000
0.000
0.090
0.821
HINT
#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#define maxn 100005
#define eps 1E-9
#define inf 1E100
using namespace std;
int n,flag,cnt,siz,pos[maxn];
double w[maxn],d[maxn],dist[maxn],tot,tmp,ans,t1,t2,t3,val[maxn];
double a,b,l,r,f[maxn],p[maxn];
struct Line{
double k,b;
int id;
}line[maxn],tl[maxn];
bool cmp(Line x,Line y){
if (abs(x.k-y.k)>eps) return x.k>y.k;
return x.b<y.b;
}
double calc(Line l1,Line l2){return -(l1.b-l2.b)/(l1.k-l2.k);}
struct DATA{
double x;
int id;
}list[maxn];
void solve(){
memcpy(tl,line,sizeof(line));
sort(tl+,tl+n+,cmp);
for (int i=,j,t;i<=n;i=j){
line[++cnt]=tl[i],t=;
for (j=i;tl[j].k==tl[i].k;j++) if (tl[j].b==tl[i].b) pos[tl[j].id]=cnt,t++;
for (j=i;tl[j].k==tl[i].k;j++) if (tl[j].b==tl[i].b) p[tl[j].id]=1.0/t;
}
list[siz=]=(DATA){-inf,};
for (int i=;i<=cnt;i++){
double t=calc(line[i],line[list[siz].id]);
while (t<list[siz].x) siz--,t=calc(line[i],line[list[siz].id]);
list[++siz]=(DATA){t,i};
}
list[++siz]=(DATA){inf,cnt+};
l=a,r=b,flag=;
for (int i=;i<=siz&&flag!=;i++){
if (list[i].x>l){
if (!flag) flag=,f[list[i-].id]=(list[i].x-l)/(r-l);
else if (flag==&&list[i].x<r) f[list[i-].id]=(list[i].x-list[i-].x)/(r-l);
else if (list[i].x>=r) f[list[i-].id]=(r-list[i-].x)/(r-l),flag=;
}
}
for (int i=;i<=n;i++) p[i]=f[pos[i]]*p[i];
for (int i=;i<=n;i++) printf("%.3lf\n",p[i]);
}
int main(){
scanf("%d%lf%lf",&n,&a,&b);
for (int i=;i<n;i++) scanf("%lf",&w[i]);
for (int i=;i<=n;i++) scanf("%lf",&d[i]);
for (int i=;i<=n;i++) dist[i]=d[i];
for (int i=;i<=n;i++) dist[i]+=dist[i-];
tot=dist[n];
if (a==b){
w[n]=a;
t1=w[]*tot,t3=w[];
for (int i=;i<=n;i++){
tmp+=w[i]*(dist[i]-dist[])*(tot-dist[i]+dist[]);
t1+=w[i]*(dist[i]-dist[]);
t2+=w[i]*(tot-dist[i]+dist[]);
t3+=w[i];
}
val[]=ans=tmp;
for (int i=;i<=n;i++){
tmp+=t1*d[i]-t2*d[i]-t3*d[i]*d[i];
ans=min(ans,tmp);
val[i]=tmp;
t1+=w[i]*tot,t2-=w[i]*tot;
t1-=t3*d[i],t2+=t3*d[i];
}
cnt=;
for (int i=;i<=n;i++) if (abs(val[i]-ans)<=eps) cnt++;
for (int i=;i<=n;i++) if (abs(val[i]-ans)<=eps) printf("%.3lf\n",1.0/cnt);
else puts("0.000");
return ;
}
l=,r=tot;
for (int i=;i<=n;i++) l+=d[i],r-=d[i],line[i].k=l*r;
t1=w[]*tot,t3=w[];
for (int i=;i<=n;i++){
tmp+=w[i]*(dist[i]-dist[])*(tot-dist[i]+dist[]);
t1+=w[i]*(dist[i]-dist[]);
t2+=w[i]*(tot-dist[i]+dist[]);
t3+=w[i];
}
line[].b=tmp;
for (int i=;i<=n;i++){
tmp+=t1*d[i]-t2*d[i]-t3*d[i]*d[i];
line[i].b=tmp;
t1+=w[i]*tot,t2-=w[i]*tot;
t1-=t3*d[i],t2+=t3*d[i];
}
for (int i=;i<=n;i++) line[i].id=i;
solve();
return ;
}
bzoj4306: 玩具厂的更多相关文章
- C#设计模式之装饰者
IronMan之装饰者 前言 上一篇的文章我们讲到要给"IronMan"配备"武器",并且还使用了"武器",效果还是不错的,对于多种环境.多 ...
- C#设计模式之工厂
IronMan之工厂 前言 实用为主,学一些用得到的技术更能在如今的社会里保命. 虽然在日常的工作中设计模式不是经常的用到,但是呢,学习它只有好处没有坏处. 设计模式像是一种“标签”,它是代码编写者思 ...
- 2017.07.11【NOIP提高组】模拟赛B组
Summary 今天的比赛打得还不错,第一题被同桌灌输的贪心,纯模拟洗脑了,然后steal的看了一下,发现怎么也对不了,一直在检查.最后10分钟才找出反例,推出动态规划方程,没有想到怎么转移,比赛就结 ...
- python 全栈开发,Day132(玩具管理页面,控制玩具通讯录,基于请求的好友关系建立)
先下载github代码,下面的操作,都是基于这个版本来的! https://github.com/987334176/Intelligent_toy/archive/v1.5.zip 注意:由于涉及到 ...
- python 全栈开发,Day129(玩具开机提示语,为多个玩具发送点播,聊天界面,app录音,app与服务器端文件传输,简单的对话)
一.玩具开机提示语 先下载github代码,下面的操作,都是基于这个版本来的! https://github.com/987334176/Intelligent_toy/archive/v1.2.zi ...
- 45.更新一下scrapy爬取工商信息爬虫代码
这里是完整的工商信息采集代码,不过此程序需要配合代理ip软件使用.问题:1.网站对ip之前没做限制,但是采集了一段时间就被检测到设置了反爬,每个ip只能访问十多次左右就被限制访问.2.网站对请求头的检 ...
- VR与AR的发展趋势分析
概要 你是否想象过与神秘的深海生物近距离接触?你是否梦想过穿戴钢铁侠那样的超先进科技装备成为超级英雄?你又是否幻想过与梦中的女神面对面的交流?这些可能在以前都只能是存在于脑海中的幻想,可是在如今有一项 ...
- 《深入浅出数据分析》-利用Excel的Solver求橡皮玩具的最大利润
背景:一玩具厂可以生产两种浴盆玩具,分别是橡皮鸭和橡皮鱼,并且原材料和生产时间都有所限制,求如何才分配生产两种玩具的数量才可以让厂商达到最大利润. 假设条件如下: 产品名称 数量 duck a ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
随机推荐
- [github] 创建个人网页
创建 github.io 的个人子域名网页,其实非常简单,只需要两步: 1)创建名为 username.github.io 的代码仓库 2)把网页代码上传到代码库中 在 username.github ...
- 简单约瑟夫环的循环单链表实现(C++)
刚刚接触C++以及数据结构,今天做了第一次尝试用C++和数据结构解决问题,问题是基于约瑟夫环问题的简单版. 先来看看约瑟夫环问题的介绍: 约瑟夫环是一个数学的应用问题:已知n个人(以编号1,2,3.. ...
- Java Annotation 必须掌握的特性
什么是Annotation? Annotation翻译为中文即为注解,意思就是提供除了程序本身逻辑外的额外的数据信息.Annotation对于标注的代码没有直接的影响,它不可以直接与标注的代码产生交互 ...
- Vs 2008 对 OpenMP 的 支持 以及 OpenMP的环境变量及库函数
Visual C++® 2008对OpenMP的支持 VC++2008根据项目属性配置的指示进行 /openmp编译器切换,当配置了OpenMP支持后,编译器会提供_OPENMP定义,可以使用#ifd ...
- Lucene中string docvalues使用utf-16的优化
原来的string docvalues使用utf-8编码,载入时转码花费大量时间,我们把转码实现从new String(bytes, "UTF-8")改用lucene的bytesR ...
- QT unit test code coverage
准备环境: qt-creator5.2.1 , gcov(gcc 默认安装),lcov(gcov 的图形化显示界面),qt_testlib 各环境介绍: 1.gcov gcov 是一个可用于C/C ...
- 使用CSS、JS的几种方法
引用CSS到html网页里的4种方法 1.行内式 在标记的style属性中设定CSS样式,这种方式本质上没有体现出CSS的优势,因此不推荐使用 <div style="widt ...
- 转载:C#中事件的由来
原文地址 http://www.tracefact.net/CSharp-Programming/Delegates-and-Events-in-CSharp.aspx 感谢博主分享! 我们继续思考转 ...
- 5事件DOM零级事件跟DOM二级事件
事件的行为传播,行为本身跟事件绑定没有关系:1.全新认识事件(某一个具体的行为)->行为本身:浏览器天生自带的一些行为操作->click,mouseover(mouseenter),mou ...
- 不同浏览器创建XMLHttpRequest对象
function getXHR() { if (XMLHttpRequest) { return new XMLHttpRequest(); } else { return new ActiveXOb ...